cp's OEIS Frontend

This is a front-end for the Online Encyclopedia of Integer Sequences, made by Christian Perfect. The idea is to provide OEIS entries in non-ancient HTML, and then to think about how they're presented visually. The source code is on GitHub.

Showing 1-4 of 4 results.

A052331 Inverse of A052330; A binary encoding of Fermi-Dirac factorization of n, shown in decimal.

Original entry on oeis.org

0, 1, 2, 4, 8, 3, 16, 5, 32, 9, 64, 6, 128, 17, 10, 256, 512, 33, 1024, 12, 18, 65, 2048, 7, 4096, 129, 34, 20, 8192, 11, 16384, 257, 66, 513, 24, 36, 32768, 1025, 130, 13, 65536, 19, 131072, 68, 40, 2049, 262144, 258, 524288, 4097, 514, 132, 1048576, 35
Offset: 1

Views

Author

Christian G. Bower, Dec 15 1999

Keywords

Comments

Every number can be represented uniquely as a product of numbers of the form p^(2^k), sequence A050376. This sequence is a binary representation of this factorization, with a(p^(2^k)) = 2^(i-1), where i is the index (A302778) of p^(2^k) in A050376. Additive with a(p^e) = sum a(p^(2^e_k)) where e = sum(2^e_k) is the binary representation of e and a(p^(2^k)) is as described above. - Franklin T. Adams-Watters, Oct 25 2005 - Index offset corrected by Antti Karttunen, Apr 17 2018

Examples

			n = 84 has Fermi-Dirac factorization A050376(5) * A050376(3) * A050376(2) = 7*4*3. Thus a(84) = 2^(5-1) + 2^(3-1) + 2^(2-1) = 16 + 4 + 2 = 22 ("10110" in binary = A182979(84)). - _Antti Karttunen_, Apr 17 2018

Links

Crossrefs

Cf. A050376, A052330, A064547, A302024, A302029, A302776, A302778, A302785, A302786, A302787, A302790, A302784.
Cf. A182979 (same sequence shown in binary).
One less than A064358.
Cf. also A156552.

Programs

  • PARI
    A052331=a(n)={for(i=1,#n=factor(n)~,n[2,i]>1||next; m=binary(n[2,i]); n=concat(n,Mat(vector(#m-1,j,[n[1,i]^2^(#m-j),m[j]]~)));n[2,i]%=2); n||return(0); m=vecsort(n[1,]); forprime(p=1,m[#m],my(j=0);while(p^2^j>1} \\ M. F. Hasler, Apr 08 2015
  • PARI
    up_to_e = 8192;
v050376 = vector(up_to_e);
ispow2(n) = (n && !bitand(n,n-1));
i = 0; for(n=1,oo,if(ispow2(isprimepower(n)), i++; v050376[i] = n); if(i == up_to_e,break));
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); }; \\ Antti Karttunen, Apr 12 2018

Formula

a(1)=0; a(n*A050376(k)) = a(n) + 2^k for a(n) < 2^k, k=0, 1, ... - Thomas Ordowski, Mar 23 2005
From Antti Karttunen, Apr 13 2018: (Start)
a(1) = 0; for n > 1, a(n) = A000079(A302785(n)-1) + a(A302776(n)).
For n > 1, a(n) = A000079(A302786(n)-1) * A302787(n).
a(n) = A064358(n)-1.
A000120(a(n)) = A064547(n).
A069010(a(n)) = A302790(n).
(End)

A302786 Index of the smallest Fermi-Dirac factor of n, a(1) = 0 by convention: a(n) = A302778(A223490(n)).

Original entry on oeis.org

0, 1, 2, 3, 4, 1, 5, 1, 6, 1, 7, 2, 8, 1, 2, 9, 10, 1, 11, 3, 2, 1, 12, 1, 13, 1, 2, 3, 14, 1, 15, 1, 2, 1, 4, 3, 16, 1, 2, 1, 17, 1, 18, 3, 4, 1, 19, 2, 20, 1, 2, 3, 21, 1, 4, 1, 2, 1, 22, 2, 23, 1, 5, 3, 4, 1, 24, 3, 2, 1, 25, 1, 26, 1, 2, 3, 5, 1, 27, 4, 28, 1, 29, 2, 4, 1, 2, 1, 30, 1, 5, 3, 2, 1, 4, 1, 31, 1, 6, 3, 32, 1, 33, 1, 2
Offset: 1

Views

Author

Antti Karttunen, Apr 13 2018

Keywords

Links

Crossrefs

A left inverse of A050376.
Cf. A001511, A052331, A223490, A302778, A302785, A302787, A302788 (ordinal transform), A302789.
Cf. also A055396.

Programs

  • Mathematica
    nn = 105; t = {}; k = 1;
While[lim = nn^(1/k); lim > 2,
     t = Join[t, Prime[Range[PrimePi[lim]]]^k]; k = 2 k];
A050376 = Union[t];
A223490[n_] := Table[{p, e} = pe; p^(2^IntegerExponent[e, 2]), {pe, FactorInteger[n]}] // Min;
a[n_] := If[n == 1, 0, FirstPosition[A050376, A223490[n]][[1]]];
Array[a, nn] (* Jean-François Alcover, Jan 08 2022, after T. D. Noe in A050376 *)
  • PARI
    up_to = 65537;
v050376 = vector(up_to);
ispow2(n) = (n && !bitand(n,n-1));
i = 0; for(n=1,oo,if(ispow2(isprimepower(n)), i++; v050376[i] = n); if(i == up_to,break));
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); };
A001511(n) = 1+valuation(n,2);
A302786(n) = if(1==n,0,A001511(A052331(n)));

Formula

a(n) = A302778(A223490(n)).
a(1) = 0; for n > 1, a(n) = A001511(A052331(n)).
For n >= 1, a(A050376(n)) = n.
For n > 1, A050376(a(n)) = A223490(n).

A302791 A filter sequence for Fermi-Dirac factorization: restricted growth sequence transform of A046523(A302024(n)).

Original entry on oeis.org

1, 2, 2, 2, 2, 3, 2, 4, 2, 4, 2, 3, 2, 4, 4, 2, 2, 4, 2, 3, 4, 4, 2, 5, 2, 4, 4, 4, 2, 6, 2, 4, 4, 4, 3, 4, 2, 4, 4, 6, 2, 6, 2, 4, 4, 4, 2, 4, 2, 4, 4, 4, 2, 6, 4, 7, 4, 4, 2, 5, 2, 4, 3, 4, 4, 6, 2, 4, 4, 6, 2, 7, 2, 4, 4, 4, 4, 6, 2, 4, 2, 4, 2, 6, 4, 4, 4, 7, 2, 7, 4, 4, 4, 4, 4, 6, 2, 4, 3, 4, 2, 6, 2, 7, 6
Offset: 1

Views

Author

Antti Karttunen, Apr 15 2018

Keywords

Comments

For all i, j: a(i) = a(j) => A064547(i) = A064547(j).
For all i, j: a(i) = a(j) => A302790(i) = A302790(j).
See also comments in A302024.

Links

Crossrefs

Cf. A037445, A046523, A050376 (gives the positions of 2's), A052331, A064547, A293442, A302024, A302787, A302790.

Programs

  • PARI
    allocatemem(2^30);
up_to = 65537;
v050376 = vector(up_to);
A050376(n) = v050376[n];
ispow2(n) = (n && !bitand(n,n-1));
i = 0; for(n=1,oo,if(ispow2(isprimepower(n)), i++; v050376[i] = n); if(i == up_to,break));
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); };
A005940(n) = { my(p=2, t=1); n--; until(!n\=2, if((n%2), (t*=p), p=nextprime(p+1))); t }; \\ Modified from code of M. F. Hasler
A302024(n) = A005940(1+A052331(n));
A046523(n) = { my(f=vecsort(factor(n)[, 2], , 4), p); prod(i=1, #f, (p=nextprime(p+1))^f[i]); };  \\ From A046523
Aux302791(n) = A046523(A302024(n));
rgs_transform(invec) = { my(om = Map(), outvec = vector(length(invec)), u=1); for(i=1, length(invec), if(mapisdefined(om,invec[i]), my(pp = mapget(om, invec[i])); outvec[i] = outvec[pp] , mapput(om,invec[i],i); outvec[i] = u; u++ )); outvec; };
write_to_bfile(start_offset,vec,bfilename) = { for(n=1, length(vec), write(bfilename, (n+start_offset)-1, " ", vec[n])); }
write_to_bfile(1,rgs_transform(vector(up_to,n,Aux302791(n))),"b302791.txt");

A302790 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)).

Original entry on oeis.org

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, 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, 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
Offset: 1

Views

Author

Antti Karttunen, Apr 13 2018

Keywords

Examples

			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.

Links

Crossrefs

Cf. A052331, A069010, A302787.

Programs

  • PARI
    up_to = 65537;
v050376 = vector(up_to);
ispow2(n) = (n && !bitand(n,n-1));
i = 0; for(n=1,oo,if(ispow2(isprimepower(n)), i++; v050376[i] = n); if(i == up_to,break));
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); };
A069010(n) = ((1 + (hammingweight(bitxor(n, n>>1)))) >> 1); \\ From A069010
A302790(n) = A069010(A052331(n));

Formula

a(n) = A069010(A052331(n)).
a(n) = A069010(A302787(n)).
a(n) = A001221(A302024(n)).
For all n >= 1, a(n) <= A064547(n).
Showing 1-4 of 4 results.

Started in 1964 by Neil J. A. Sloane | Maintained by The OEIS Foundation Inc.

Content is available under The OEIS End-User License Agreement.