|
%I #11 Apr 16 2018 22:11:22
%S 0,1,1,1,1,1,1,2,1,2,1,1,1,2,2,1,1,2,1,1,2,2,1,1,1,2,2,2,1,2,1,2,2,2,
%T 1,2,1,2,2,2,1,2,1,2,2,2,1,2,1,2,2,2,1,2,2,3,2,2,1,1,1,2,1,2,2,2,1,2,
%U 2,2,1,3,1,2,2,2,2,2,1,2,1,2,1,2,2,2,2,3,1,3,2,2,2,2,2,2,1,2,1,2,1,2,1,3,2
%N Number of runs of consecutive Fermi-Dirac factors of n (the runs are separated by gaps between indices of factors): a(n) = A069010(A052331(n)).
%H Antti Karttunen, <a href="/A302790/b302790.txt">Table of n, a(n) for n = 1..65537</a>
%F a(n) = A069010(A052331(n)).
%F a(n) = A069010(A302787(n)).
%F a(n) = A001221(A302024(n)).
%F For all n >= 1, a(n) <= A064547(n).
%e n = 84 has Fermi-Dirac factorization as A050376(2) * A050376(3) * A050376(5) = 3*4*7. Because there is a gap between A050376(3) and A050376(5), the factors occur in two separate runs (3*4 and 7), thus a(84) = 2.
%o (PARI)
%o up_to = 65537;
%o v050376 = vector(up_to);
%o ispow2(n) = (n && !bitand(n,n-1));
%o i = 0; for(n=1,oo,if(ispow2(isprimepower(n)), i++; v050376[i] = n); if(i == up_to,break));
%o A052331(n) = { my(s=0,e); while(n > 1, fordiv(n, d, if(((n/d)>1)&&ispow2(isprimepower(n/d)), e = vecsearch(v050376, n/d); if(!e, print("v050376 too short!"); return(1/0)); s += 2^(e-1); n = d; break))); (s); };
%o A069010(n) = ((1 + (hammingweight(bitxor(n, n>>1)))) >> 1); \\ From A069010
%o A302790(n) = A069010(A052331(n));
%Y Cf. A052331, A069010, A302787.
%K nonn
%O 1,8
%A _Antti Karttunen_, Apr 13 2018
|
