A182979 -id:A182979 - cp's OEIS frontend

A050376 "Fermi-Dirac primes": numbers of the form p^(2^k) where p is prime and k >= 0.

2, 3, 4, 5, 7, 9, 11, 13, 16, 17, 19, 23, 25, 29, 31, 37, 41, 43, 47, 49, 53, 59, 61, 67, 71, 73, 79, 81, 83, 89, 97, 101, 103, 107, 109, 113, 121, 127, 131, 137, 139, 149, 151, 157, 163, 167, 169, 173, 179, 181, 191, 193, 197, 199, 211, 223, 227, 229, 233, 239, 241

Offset: 1

Christian G. Bower, Nov 15 1999

Every number n is a product of a unique subset of these numbers. This is sometimes called the Fermi-Dirac factorization of n (see A182979). Proof: In the prime factorization n = Product_{j>=1} p(j)^e(j) expand every exponent e(j) as binary number and pick the terms of this sequence corresponding to the positions of the ones in binary (it is clear that both n and n^2 have the same number of factors in this sequence, and that each factor appears with exponent 1 or 0).
Or, a(1) = 2; for n>1, a(n) = smallest number which cannot be obtained as the product of previous terms. This is evident from the unique factorization theorem and the fact that every number can be expressed as the sum of powers of 2. - Amarnath Murthy, Jan 09 2002
Except for the first term, same as A084400. - David Wasserman, Dec 22 2004
The least number having 2^n divisors (=A037992(n)) is the product of the first n terms of this sequence according to Ramanujan.
According to the Bose-Einstein distribution of particles, an unlimited number of particles may occupy the same state. On the other hand, according to the Fermi-Dirac distribution, no two particles can occupy the same state (by the Pauli exclusion principle). Unique factorizations of the positive integers by primes (A000040) and over terms of A050376 one can compare with two these distributions in physics of particles. In the correspondence with this, the factorizations over primes one can call "Bose-Einstein factorizations", while the factorizations over distinct terms of A050376 one can call "Fermi-Dirac factorizations". - Vladimir Shevelev, Apr 16 2010
The numbers of the form p^(2^k), where p is prime and k >= 0, might thus be called the "Fermi-Dirac primes", while the classic primes might be called the "Bose-Einstein primes". - Daniel Forgues, Feb 11 2011
In the theory of infinitary divisors, the most natural name of the terms is "infinitary primes" or "i-primes". Indeed, n is in the sequence, if and only if it has only two infinitary divisors. Since 1 and n are always infinitary divisors of n>1, an i-prime has no other infinitary divisors. - Vladimir Shevelev, Feb 28 2011
{a(n)} is the minimal set including all primes and closed with respect to squaring. In connection with this, note that n and n^2 have the same number of factors in their Fermi-Dirac representations. - Vladimir Shevelev, Mar 16 2012
In connection with this sequence, call an integer compact if the factors in its Fermi-Dirac factorization are pairwise coprime. The density of such integers equals (6/Pi^2)*Product_{prime p} (1+(Sum_{i>=1} p^(-(2^i-1))/(p+1))) = 0.872497... It is interesting that there exist only 7 compact factorials listed in A169661. - Vladimir Shevelev, Mar 17 2012
The first k terms of the sequence solve the following optimization problem:
Let x_1, x_2,..., x_k be integers with the restrictions: 2<=x_1A064547(Product{i=1..k} x_i) >= k.  Let the goal function be Product_{i=1..k} x_i. Then the minimal value of the goal function is Product_{i=1..k} a(i). - Vladimir Shevelev, Apr 01 2012
From Joerg Arndt, Mar 11 2013: (Start)
Similarly to the first comment, for the sequence "Numbers of the form p^(3^k) or p^(2*3^k) where p is prime and k >= 0" one obtains a factorization into distinct factors by using the ternary expansion of the exponents (here n and n^3 have the same number of such factors).
The generalization to base r would use "Numbers of the form p^(r^k), p^(2*r^k), p^(3*r^k), ..., p^((r-1)*r^k) where p is prime and k >= 0" (here n and n^r have the same number of (distinct) factors).  (End)
The first appearance of this sequence as a multiplicative basis in number theory with some new notions, formulas and theorems may have been in my 1981 paper (see the Abramovich reference). - Vladimir Shevelev, Apr 27 2014
Numbers n for which A064547(n) = 1. - Antti Karttunen, Feb 10 2016
Lexicographically earliest sequence of distinct nonnegative integers such that no term is a product of 2 or more distinct terms. Removing the distinctness requirement, the sequence becomes A000040 (the prime numbers); and the equivalent sequence where the product is of 2 distinct terms is A026416 (without its initial term, 1). - Peter Munn, Mar 05 2019
The sequence was independently developed as a multiplicative number system in 1985-1986 (and first published in 1995, see the Uhlmann reference) using a proof method involving representations of positive integers as sums of powers of 2. This approach offers an arguably simpler and more flexible means for analyzing the sequence. - Jeffrey K. Uhlmann, Nov 09 2022

			Prime powers which are not terms of this sequence:
  8 = 2^3 = 2^(1+2), 27 = 3^3 = 3^(1+2), 32 = 2^5 = 2^(1+4),
  64 = 2^6 = 2^(2+4), 125 = 5^3 = 5^(1+2), 128 = 2^7 = 2^(1+2+4)
"Fermi-Dirac factorizations":
  6 = 2*3, 8 = 2*4, 24 = 2*3*4, 27 = 3*9, 32 = 2*16, 64 = 4*16,
  108 = 3*4*9, 120 = 2*3*4*5, 121 = 121, 125 = 5*25, 128 = 2*4*16.

V. S. Abramovich, On an analog of the Euler function, Proceeding of the North-Caucasus Center of the Academy of Sciences of the USSR (Rostov na Donu) (1981) No. 2, 13-17 (Russian; MR0632989(83a:10003)).
S. Ramanujan, Highly Composite Numbers, Collected Papers of Srinivasa Ramanujan, p. 125, Ed. G. H. Hardy et al., AMS Chelsea 2000.
V. S. Shevelev, Multiplicative functions in the Fermi-Dirac arithmetic, Izvestia Vuzov of the North-Caucasus region, Nature sciences 4 (1996), 28-43 (in Russian; MR 2000f: 11097, pp. 3912-3913).
J. K. Uhlmann, Dynamic map building and localization: new theoretical foundations, Doctoral Dissertation, University of Oxford, Appendix 16, 1995.

T. D. Noe and Charles R Greathouse IV, Table of n, a(n) for n = 1..10000
Steven R. Finch, Unitarism and Infinitarism, February 25, 2004. [Cached copy, with permission of the author]
Simon Litsyn and Vladimir Shevelev, On factorization of integers with restrictions on the exponent, INTEGERS: Electronic Journal of Combinatorial Number Theory, 7 (2007), #A33, 1-36.
OEIS Wiki, "Fermi-Dirac representation" of n.
Vladimir Shevelev, Compact integers and factorials, Acta Arith. 126 (2007), no.3, 195-236.
J. K. Uhlmann, Appendix 16, Doctoral Dissertation, University of Oxford, page 243, 1995.

Cf. A000040 (primes, is a subsequence), A026416, A026477, A037992 (partial products), A050377-A050380, A052330, A064547, A066724, A084400, A176699, A182979.
Cf. A268388 (complement without 1).
Cf. A124010, subsequence of A000028, A000961, A213925, A223490.
Cf. A228520, A186945 (Fermi-Dirac analog of Ramanujan primes, A104272, and Labos primes, A080359).
Cf. also A268385, A268391, A268392.

a050376 n = a050376_list !! (n-1)
a050376_list = filter ((== 1) . a209229 . a100995) [1..]
-- Reinhard Zumkeller, Mar 19 2013

isA050376 := proc(n)
    local f,e;
    f := ifactors(n)[2] ;
    if nops(f) = 1 then
        e := op(2,op(1,f)) ;
        if isA000079(e) then
            true;
        else
            false;
        end if;
    else
        false;
    end if;
end proc:
A050376 := proc(n)
    option remember ;
    local a;
    if n = 1 then
        2 ;
    else
        for a from procname(n-1)+1 do
            if isA050376(a) then
                return a;
            end if;
        end do:
    end if;
end proc: # R. J. Mathar, May 26 2017

nn = 300; t = {}; k = 1; While[lim = nn^(1/k); lim > 2,  t = Join[t, Prime[Range[PrimePi[lim]]]^k]; k = 2 k]; t = Union[t] (* T. D. Noe, Apr 05 2012 *)

{a(n)= local(m, c, k, p); if(n<=1, 2*(n==1), n--; c=0; m=2; while( cMichael Somos, Apr 15 2005; edited by Michel Marcus, Aug 07 2021

lst(lim)=my(v=primes(primepi(lim)),t); forprime(p=2,sqrt(lim),t=p; while((t=t^2)<=lim,v=concat(v,t))); vecsort(v) \\ Charles R Greathouse IV, Apr 10 2012

is_A050376(n)=2^#binary(n=isprimepower(n))==n*2 \\ M. F. Hasler, Apr 08 2015

ispow2(n)=n && n>>valuation(n,2)==1
is(n)=ispow2(isprimepower(n)) \\ Charles R Greathouse IV, Sep 18 2015

isok(n)={my(e=isprimepower(n)); e && !bitand(e,e-1)} \\ Andrew Howroyd, Oct 16 2024

from sympy import isprime, perfect_power
def ok(n):
  if isprime(n): return True
  answer = perfect_power(n)
  if not answer: return False
  b, e = answer
  if not isprime(b): return False
  while e%2 == 0: e //= 2
  return e == 1
def aupto(limit):
  alst, m = [], 1
  for m in range(1, limit+1):
    if ok(m): alst.append(m)
  return alst
print(aupto(241)) # Michael S. Branicky, Feb 03 2021

from sympy import primepi, integer_nthroot
def A050376(n):
    def bisection(f,kmin=0,kmax=1):
        while f(kmax) > kmax: kmax <<= 1
        kmin = kmax >> 1
        while kmax-kmin > 1:
            kmid = kmax+kmin>>1
            if f(kmid) <= kmid:
                kmax = kmid
            else:
                kmin = kmid
        return kmax
    def f(x): return n+x-sum(primepi(integer_nthroot(x,1<Chai Wah Wu, Feb 18-19 2025

(define A050376 (MATCHING-POS 1 1 (lambda (n) (= 1 (A064547 n)))))
;; Requires also my IntSeq-library. - Antti Karttunen, Feb 09 2016

From Vladimir Shevelev, Mar 16 2012: (Start)
Product_{i>=1} a(i)^k_i = n!, where k_i = floor(n/a(i)) - floor(n/a(i)^2) + floor(n/a(i)^3) - floor(n/a(i)^4) + ...
Denote by A(x) the number of terms not exceeding x.
Then A(x) = pi(x) + pi(x^(1/2)) + pi(x^(1/4)) + pi(x^(1/8)) + ...
Conversely, pi(x) = A(x) - A(sqrt(x)). For example, pi(37) = A(37) - A(6) = 16-4 = 12. (End)
A209229(A100995(a(n))) = 1. - Reinhard Zumkeller, Mar 19 2013
From Vladimir Shevelev, Aug 31 2013: (Start)
A Fermi-Dirac analog of Euler product: Zeta(s) = Product_{k>=1} (1+a(k)^(-s)), for s > 1.
In particular, Product_{k>=1} (1+a(k)^(-2)) = Pi^2/6. (End)
a(n) = A268385(A268392(n)). [By their definitions.] - Antti Karttunen, Feb 10 2016
A000040 union A001248 union A030514 union A179645 union A030635 union .... - R. J. Mathar, May 26 2017

Edited by Charles R Greathouse IV, Mar 17 2010

More examples from Daniel Forgues, Feb 09 2011

A052330 Let S_k denote the first 2^k terms of this sequence and let b_k be the smallest positive integer that is not in S_k; then the numbers b_k*S_k are the next 2^k terms.

Original entry on oeis.org

1, 2, 3, 6, 4, 8, 12, 24, 5, 10, 15, 30, 20, 40, 60, 120, 7, 14, 21, 42, 28, 56, 84, 168, 35, 70, 105, 210, 140, 280, 420, 840, 9, 18, 27, 54, 36, 72, 108, 216, 45, 90, 135, 270, 180, 360, 540, 1080, 63, 126, 189, 378, 252, 504, 756, 1512, 315, 630, 945, 1890

Offset: 0

Christian G. Bower, Dec 15 1999

Inverse of sequence A064358 considered as a permutation of the positive integers. - Howard A. Landman, Sep 25 2001
This sequence is not exactly a permutation because it has offset 0 but doesn't contain 0. A052331 is its exact inverse, which has offset 1 and contains 0. See also A064358.
Are there any other values of n besides 4 and 36 with a(n) = n? - Thomas Ordowski, Apr 01 2005
4 = 100 = 4^1 * 3^0 * 2^0, 36 = 100100 = 9^1 * 7^0 * 5^0 * 4^1 * 3^0 * 2^0. - Thomas Ordowski, May 26 2005
Ordering of positive integers by increasing "Fermi-Dirac representation", which is a representation of the "Fermi-Dirac factorization", term implying that each prime power with a power of two as exponent may appear at most once in the "Fermi-Dirac factorization" of n. (Cf. comment in A050376; see also the OEIS Wiki page.) - Daniel Forgues, Feb 11 2011
The subsequence consisting of the squarefree terms is A019565. - Peter Munn, Mar 28 2018
Let f(n) = A050376(n) be the n-th Fermi-Dirac prime. The FDH-number of a strict integer partition (y_1,...,y_k) is f(y_1)*...*f(y_k). A binary index of n is any position of a 1 in its reversed binary expansion. The binary indices of n are row n of A048793. Then a(n) is the number whose binary indices are the parts of the strict integer partition with FDH-number n. - Gus Wiseman, Aug 19 2019
The set of indices of odd-valued terms has asymptotic density 0. In this sense (using the order they appear in this permutation) 100% of numbers are even. - Peter Munn, Aug 26 2019

			Terms following 5 are 10, 15, 30, 20, 40, 60, 120; this is followed by 7 as 6 has already occurred. - _Philippe Deléham_, Jun 03 2015
From _Antti Karttunen_, Apr 13 2018, after also _Philippe Deléham_'s Jun 03 2015 example: (Start)
This sequence can be regarded also as an irregular triangle with rows of lengths 1, 1, 2, 4, 8, 16, ..., that is, it can be represented as a binary tree, where each left hand child contains A300841(k), and each right hand child contains 2*A300841(k), when their parent contains k:
                                     1
                                     |
                  ...................2...................
                 3                                       6
       4......../ \........8                  12......../ \........24
      / \                 / \                 / \                 / \
     /   \               /   \               /   \               /   \
    /     \             /     \             /     \             /     \
   5       10         15       30         20       40         60      120
  7 14   21  42     28  56   84  168    35  70  105  210   140 280  420 840
  etc.
Compare also to trees like A005940 and A283477, and sequences A207901 and A302783.
(End)

Antti Karttunen, Table of n, a(n) for n = 0..8191 (Terms 0..1023 from T. D. Noe)
Michael De Vlieger, Fan style binary tree showing a(n), n = 0..16383, with a color function representing primes in red, perfect powers of primes in gold, squarefree composites in greens, and numbers neither squarefree nor prime powers in blue or purple. Purple represents powerful numbers that are not prime powers. This is a 15 level version of Karttunen's diagram in the example.
OEIS Wiki, Ordering of positive integers by increasing "Fermi-Dirac representation"
Index entries for sequences that are permutations of the natural numbers

Subsequences: A019565 (squarefree terms), A050376 (the left edge from 2 onward), A336882 (odd terms).
Cf. A050030, A052331 (inverse), A096111, A096113, A096114, A096115, A096116, A096118, A096119, A182979, A207901, A300841, A302023, A302783.
Cf. A000120, A029931, A048793, A064547, A070939, A213925, A299755, A299757, A327041.

a = {1}; Do[a = Join[a, a*Min[Complement[Range[Max[a] + 1], a]]], {n, 1, 6}]; a (* Ivan Neretin, May 09 2015 *)

up_to_e = 13; \\ Good for computing up to n = (2^13)-1
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));
A050376(n) = v050376[n];
A052330(n) = { my(p=1,i=1); while(n>0, if(n%2, p *= A050376(i)); i++; n >>= 1); (p); }; \\ Antti Karttunen, Apr 12 2018

a(0)=1; a(n+2^k)=a(n)*b(k) for n < 2^k, k = 0, 1, ... where b is A050376. - Thomas Ordowski, Mar 04 2005
The binary representation of n, n = Sum_{i=0..1+floor(log_2(n))} n_i * 2^i, n_i in {0,1}, is taken as the "Fermi-Dirac representation" (A182979) of a(n), a(n) = Product_{i=0..1+floor(log_2(n))} (b_i)^(n_i) where b_i is A050376(i), i.e., the i-th "Fermi-Dirac prime" (prime power with exponent being a power of 2). - Daniel Forgues, Feb 12 2011
From Antti Karttunen, Apr 12 & 17 2018: (Start)
a(0) = 1; a(2n) = A300841(a(n)), a(2n+1) = 2*A300841(a(n)).
a(n) = A207901(A006068(n)) = A302783(A003188(n)) = A302781(A302845(n)).
(End)

Entry revised Mar 17 2005 by N. J. A. Sloane, based on comments from several people, especially David Wasserman and Thomas Ordowski

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

Christian G. Bower, Dec 15 1999

nonn

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

			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

Antti Karttunen, Table of n, a(n) for n = 1..4096
Index entries for sequences that are permutations of the natural numbers

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.

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

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

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)

A279065 Fermi-Dirac primeth recurrence: a(0)=1; thereafter a(n+1) = a(n)-th number of the form p^(2^k) where p is prime and k>=0.

Original entry on oeis.org

1, 2, 3, 4, 5, 7, 11, 19, 47, 169, 907, 6829, 67931, 851891, 13034887, 237522877, 5057212439, 123890683831

Offset: 0

Gus Wiseman, Dec 10 2016

Daniel Forgues (see A182979) and Reinhard Zumkeller (see A213925) describe the increasing sequence of positive integers of the form p^(2^k) where p is prime and k>=0 (A050376 or A084400) as Fermi-Dirac primes, because any positive integer has a unique factorization into distinct terms.

OEIS Wiki, "Fermi-Dirac representation" of n.

Cf. A007097, A050376, A084400, A182979, A213925.

nn=10000;FDfactor[n_]:=If[n===1,{},Sort[Join@@Cases[FactorInteger[n],{p_,k_}:>Power[p,Cases[Position[IntegerDigits[k,2]//Reverse,1],{m_}->2^(m-1)]]]]];
FDprimeList=Array[FDfactor,nn,1,Union];
NestWhileList[Part[FDprimeList,#]&,1,#<=Length[FDprimeList]&]

lista(kmax) = {my(m = 1, c=0, isp); print1(1, ", "); for(k = 1, kmax, isp = isprimepower(k); if(isp && isp >> valuation(isp, 2) == 1, c++); if(c == m, print1(k,", "); m=k));} \\ Amiram Eldar, Oct 05 2023

a(15)-a(17) from Amiram Eldar, Oct 05 2023

A368781 The maximal exponent in the unique factorization of n in terms of distinct "Fermi-Dirac primes".

Original entry on oeis.org

0, 1, 1, 2, 1, 1, 1, 2, 2, 1, 1, 2, 1, 1, 1, 4, 1, 2, 1, 2, 1, 1, 1, 2, 2, 1, 2, 2, 1, 1, 1, 4, 1, 1, 1, 2, 1, 1, 1, 2, 1, 1, 1, 2, 2, 1, 1, 4, 2, 2, 1, 2, 1, 2, 1, 2, 1, 1, 1, 2, 1, 1, 2, 4, 1, 1, 1, 2, 1, 1, 1, 2, 1, 1, 2, 2, 1, 1, 1, 4, 4, 1, 1, 2, 1, 1, 1

Offset: 1

Amiram Eldar, Jan 05 2024

First differs from A335428 at n = 36. Differs from A050377, A344417 and A347437 at n = 1 and then at n = 36.

In the unique factorization of n in terms of distinct "Fermi-Dirac primes", n is represented as a product of prime powers (A246655) whose exponents are powers of 2 (A000079). a(n) is the maximal exponent of these prime powers (and not the maximal exponent of the exponents that are powers of 2). Thus, a(n) is a power of 2 for n >= 2.

			For n = 972 = 2^2 * 3^5, the unique factorization of 972 in terms of distinct "Fermi-Dirac primes" is 2^(2^1) * 3^(2^0) * 3^(2^2). Therefore, a(972) = 2^2 = 4.

Amiram Eldar, Table of n, a(n) for n = 1..10000
Wikipedia, Fermi-Dirac prime.

Cf. A000079, A050376, A051903, A053644, A182979, A246655, A299090, A353897.

Cf. A050377, A335428, A344417, A347437.

a[n_] := 2^Floor[Log2[Max[FactorInteger[n][[;; , 2]]]]]; a[1] = 0; Array[a, 100]

a(n) = if(n > 1, 2^exponent(vecmax(factor(n)[, 2])), 0);

from sympy import factorint
def A368781(n): return 1<1 else 0 # Chai Wah Wu, Apr 11 2025

a(n) = A053644(A051903(n)).
a(n) = 2^(A299090(n)-1) for n >= 2.
Asymptotic mean: Limit_{m->oo} (1/m) * Sum_{k=1..m} a(k) = 1 + Sum_{k>=1} (2^(k-1) * (1 - 1/zeta(2^k))) = 1.56056154773294953123... .
a(n) = A051903(A353897(n)). - Amiram Eldar, May 07 2024

cp's OEIS Frontend

A050376 "Fermi-Dirac primes": numbers of the form p^(2^k) where p is prime and k >= 0.

Views

Author

Keywords

Comments

Examples

References

Links

Crossrefs

Programs

Haskell

Maple

Mathematica

PARI

PARI

PARI

PARI

PARI

Python

Python

Scheme

Formula

Extensions

A052330 Let S_k denote the first 2^k terms of this sequence and let b_k be the smallest positive integer that is not in S_k; then the numbers b_k*S_k are the next 2^k terms.

Views

Author

Keywords

Comments

Examples

Links

Crossrefs

Programs

Mathematica

PARI

Formula

Extensions

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

Views

Author

Keywords

Comments

Examples

Links

Crossrefs

Programs

PARI

PARI

Formula

A279065 Fermi-Dirac primeth recurrence: a(0)=1; thereafter a(n+1) = a(n)-th number of the form p^(2^k) where p is prime and k>=0.

Views

Author

Keywords

Comments

Links

Crossrefs

Programs

Mathematica

PARI

Extensions

A368781 The maximal exponent in the unique factorization of n in terms of distinct "Fermi-Dirac primes".

Views

Author

Keywords

Comments

Examples

Links

Crossrefs

Programs

Mathematica

PARI

Python

Formula

A298480 Lexicographically earliest sequence of distinct positive terms such that the Fermi-Dirac factorizations of two consecutive terms differ by exactly one factor.

Views

Author

Keywords

Comments

Examples

Links