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-8 of 8 results.

A000720 pi(n), the number of primes <= n. Sometimes called PrimePi(n) to distinguish it from the number 3.14159...

Original entry on oeis.org

0, 1, 2, 2, 3, 3, 4, 4, 4, 4, 5, 5, 6, 6, 6, 6, 7, 7, 8, 8, 8, 8, 9, 9, 9, 9, 9, 9, 10, 10, 11, 11, 11, 11, 11, 11, 12, 12, 12, 12, 13, 13, 14, 14, 14, 14, 15, 15, 15, 15, 15, 15, 16, 16, 16, 16, 16, 16, 17, 17, 18, 18, 18, 18, 18, 18, 19, 19, 19, 19, 20, 20, 21, 21, 21, 21, 21, 21
Offset: 1

Views

Author

Keywords

Comments

Partial sums of A010051 (characteristic function of primes). - Jeremy Gardiner, Aug 13 2002
pi(n) and prime(n) are inverse functions: a(A000040(n)) = n and A000040(n) is the least number m such that A000040(a(m)) = A000040(n). A000040(a(n)) = n if (and only if) n is prime. - Jonathan Sondow, Dec 27 2004
See the additional references and links mentioned in A143227. - Jonathan Sondow, Aug 03 2008
A lower bound that gets better with larger N is that there are at least T prime numbers less than N, where the recursive function T is: T = N - N*Sum_{i=0..T(sqrt(N))} A005867(i)/A002110(i). - Ben Paul Thurston, Aug 23 2010
Number of partitions of 2n into exactly two parts with the smallest part prime. - Wesley Ivan Hurt, Jul 20 2013
Equivalent to the Riemann hypothesis: abs(a(n) - li(n)) < sqrt(n)*log(n)/(8*Pi), for n >= 2657, where li(n) is the logarithmic integral (Lowell Schoenfeld). - Ilya Gutkovskiy, Jul 05 2016
The second Hardy-Littlewood conjecture, that pi(x) + pi(y) >= pi(x + y) for integers x and y with min{x, y} >= 2, is known to hold for (x, y) sufficiently large (Udrescu 1975). - Peter Luschny, Jan 12 2021

Examples

			There are 3 primes <= 6, namely 2, 3 and 5, so pi(6) = 3.
		

References

  • M. Abramowitz and I. A. Stegun, eds., Handbook of Mathematical Functions, National Bureau of Standards Applied Math. Series 55, 1964 (and various reprintings), p. 870.
  • Tom M. Apostol, Introduction to Analytic Number Theory, Springer-Verlag, 1976, p. 8.
  • Raymond Ayoub, An Introduction to the Analytic Theory of Numbers, Amer. Math. Soc., 1963; p. 129.
  • Florian Cajori, A History of Mathematical Notations, Dover edition (2012), par. 409.
  • Richard Crandall and Carl Pomerance, Prime Numbers: A Computational Perspective, Springer, NY, 2001; see p. 5.
  • G. H. Hardy and E. M. Wright, An Introduction to the Theory of Numbers, 5th ed., Oxford Univ. Press, 1979, Theorems 6, 7, 420.
  • G. J. O. Jameson, The Prime Number Theorem, Camb. Univ. Press, 2003. [See also the review by D. M. Bressoud (link below).]
  • Władysław Narkiewicz, The Development of Prime Number Theory, Springer-Verlag, 2000.
  • Paulo Ribenboim, The Little Book of Bigger Primes, Springer-Verlag NY 2004. See pp. 132-133, 157-184.
  • József Sándor, Dragoslav S. Mitrinovic and Borislav Crstici, Handbook of Number Theory I, Springer Science & Business Media, 2005, Section VII.1. (For inequalities, etc.).
  • N. J. A. Sloane, A Handbook of Integer Sequences, Academic Press, 1973 (includes this sequence).
  • N. J. A. Sloane and Simon Plouffe, The Encyclopedia of Integer Sequences, Academic Press, 1995 (includes this sequence).
  • Gerald Tenenbaum and Michel Mendès France, Prime Numbers and Their Distribution, AMS Providence RI, 1999.
  • V. Udrescu, Some remarks concerning the conjecture pi(x + y) <= pi(x) + pi(y), Rev. Roumaine Math. Pures Appl. 20 (1975), 1201-1208.

Crossrefs

Closely related:
A099802: Number of primes <= 2n.
A060715: Number of primes between n and 2n (exclusive).
A035250: Number of primes between n and 2n (inclusive).
A038107: Number of primes < n^2.
A014085: Number of primes between n^2 and (n+1)^2.
A007053: Number of primes <= 2^n.
A036378: Number of primes p between powers of 2, 2^n < p <= 2^(n+1).
A006880: Number of primes < 10^n.
A006879: Number of primes with n digits.
A033270: Number of odd primes <= n.
A065855: Number of composites <= n.
For lists of large values of a(n) see, e.g., A005669(n) = a(A002386(n)), A214935(n) = a(A205827(n)).
Related sequences:
Primes (p) and composites (c): A000040, A002808, A065855.
Primes between p(n) and 2*p(n): A063124, A070046; between c(n) and 2*c(n): A376761; between n and 2*n: A035250, A060715, A077463, A108954.
Composites between p(n) and 2*p(n): A246514; between c(n) and 2*c(n): A376760; between n and 2*n: A075084, A307912, A307989, A376759.

Programs

  • Haskell
    a000720 n = a000720_list !! (n-1)
    a000720_list = scanl1 (+) a010051_list  -- Reinhard Zumkeller, Sep 15 2011
    
  • Magma
    [ #PrimesUpTo(n): n in [1..200] ];  // Bruno Berselli, Jul 06 2011
    
  • Maple
    with(numtheory); A000720 := pi; [ seq(A000720(i),i=1..50) ];
  • Mathematica
    A000720[n_] := PrimePi[n]; Table[ A000720[n], {n, 1, 100} ]
    Array[ PrimePi[ # ]&, 100 ]
    Accumulate[Table[Boole[PrimeQ[n]],{n,100}]] (* Harvey P. Dale, Jan 17 2015 *)
  • PARI
    A000720=vector(100,n,omega(n!)) \\ For illustration only; better use A000720=primepi
    
  • PARI
    vector(300,j,primepi(j)) \\ Joerg Arndt, May 09 2008
    
  • Python
    from sympy import primepi
    for n in range(1,100): print(primepi(n), end=', ') # Stefano Spezia, Nov 30 2018
  • Sage
    [prime_pi(n) for n in range(1, 79)]  # Zerinvary Lajos, Jun 06 2009
    

Formula

The prime number theorem gives the asymptotic expression a(n) ~ n/log(n).
For x > 1, pi(x) < (x / log x) * (1 + 3/(2 log x)). For x >= 59, pi(x) > (x / log x) * (1 + 1/(2 log x)). [Rosser and Schoenfeld]
For x >= 355991, pi(x) < (x / log(x)) * (1 + 1/log(x) + 2.51/(log(x))^2 ). For x >= 599, pi(x) > (x / log(x)) * (1 + 1/log(x)). [Dusart]
For x >= 55, x/(log(x) + 2) < pi(x) < x/(log(x) - 4). [Rosser]
For n > 1, A138194(n) <= a(n) <= A138195(n) (Tschebyscheff, 1850). - Reinhard Zumkeller, Mar 04 2008
For n >= 33, a(n) = 1 + Sum_{j=3..n} ((j-2)! - j*floor((j-2)!/j)) (Hardy and Wright); for n >= 1, a(n) = n - 1 + Sum_{j=2..n} (floor((2 - Sum_{i=1..j} (floor(j/i)-floor((j-1)/i)))/j)) (Ruiz and Sondow 2000). - Benoit Cloitre, Aug 31 2003
a(n) = A001221(A000142(n)). - Benoit Cloitre, Jun 03 2005
G.f.: Sum_{p prime} x^p/(1-x) = b(x)/(1-x), where b(x) is the g.f. for A010051. - Franklin T. Adams-Watters, Jun 15 2006
a(n) = A036234(n) - 1. - Jaroslav Krizek, Mar 23 2009
From Enrique Pérez Herrero, Jul 12 2010: (Start)
a(n) = Sum_{i=2..n} floor((i+1)/A000203(i)).
a(n) = Sum_{i=2..n} floor(A000010(n)/(i-1)).
a(n) = Sum_{i=2..n} floor(2/A000005(n)). (End)
Let pf(n) denote the set of prime factors of an integer n. Then a(n) = card(pf(n!/floor(n/2)!)). - Peter Luschny, Mar 13 2011
a(n) = -Sum_{p <= n} mu(p). - Wesley Ivan Hurt, Jan 04 2013
a(n) = (1/2)*Sum_{p <= n} (mu(p)*d(p)*sigma(p)*phi(p)) + sum_{p <= n} p^2. - Wesley Ivan Hurt, Jan 04 2013
a(1) = 0 and then, for all k >= 1, repeat k A001223(k) times. - Jean-Christophe Hervé, Oct 29 2013
a(n) = n/(log(n) - 1 - Sum_{k=1..m} A233824(k)/log(n)^k + O(1/log(n)^{m+1})) for m > 0. - Jonathan Sondow, Dec 19 2013
a(n) = A001221(A003418(n)). - Eric Desbiaux, May 01 2014
a(n) = Sum_{j=2..n} H(-sin^2 (Pi*(Gamma(j)+1)/j)) where H(x) is the Heaviside step function, taking H(0)=1. - Keshav Raghavan, Jun 18 2016
a(A014076(n)) = (1/2) * (A014076(n) + 1) - n + 1. - Christopher Heiling, Mar 03 2017
From Steven Foster Clark, Sep 25 2018: (Start)
a(n) = Sum_{m=1..n} A143519(m) * floor(n/m).
a(n) = Sum_{m=1..n} A001221(m) * A002321(floor(n/m)) where A002321() is the Mertens function.
a(n) = Sum_{m=1..n} |A143519(m)| * A002819(floor(n/m)) where A002819() is the Liouville Lambda summatory function and |x| is the absolute value of x.
a(n) = Sum_{m=1..n} A137851(m)/m * H(floor(n/m)) where H(n) = Sum_{m=1..n} 1/m is the harmonic number function.
a(n) = Sum_{m=1..log_2(n)} A008683(m) * A025528(floor(n^(1/m))) where A008683() is the Moebius mu function and A025528() is the prime-power counting function.
(End)
Sum_{k=2..n} 1/a(k) ~ (1/2) * log(n)^2 + O(log(n)) (de Koninck and Ivić, 1980). - Amiram Eldar, Mar 08 2021
a(n) ~ 1/(n^(1/n)-1). - Thomas Ordowski, Jan 30 2023
a(n) = Sum_{j=2..n} floor(((j - 1)! + 1)/j - floor((j - 1)!/j)) [Mináč, unpublished] (see Ribenboim, pp. 132-133). - Stefano Spezia, Apr 13 2025
a(n) = n - 1 - Sum_{k=2..floor(log_2(n))} pi_k(n), where pi_k(n) is the number of k-almost primes <= n. - Daniel Suteu, Aug 27 2025

Extensions

Additional links contributed by Lekraj Beedassy, Dec 23 2003
Edited by M. F. Hasler, Apr 27 2018 and (links recovered) Dec 21 2018

A002386 Primes (lower end) with record gaps to the next consecutive prime: primes p(k) where p(k+1) - p(k) exceeds p(j+1) - p(j) for all j < k.

Original entry on oeis.org

2, 3, 7, 23, 89, 113, 523, 887, 1129, 1327, 9551, 15683, 19609, 31397, 155921, 360653, 370261, 492113, 1349533, 1357201, 2010733, 4652353, 17051707, 20831323, 47326693, 122164747, 189695659, 191912783, 387096133, 436273009, 1294268491
Offset: 1

Views

Author

Keywords

Comments

See the links by Jens Kruse Andersen et al. for very large gaps.

References

  • B. C. Berndt, Ramanujan's Notebooks Part IV, Springer-Verlag, see p. 133.
  • D. E. Knuth, The Art of Computer Programming. Addison-Wesley, Reading, MA, Vol. 3, Sect 6.1, Table 1.
  • M. Kraitchik, Recherches sur la Théorie des Nombres. Gauthiers-Villars, Paris, Vol. 1, 1924, Vol. 2, 1929, see Vol. 1, p. 14.
  • N. J. A. Sloane, A Handbook of Integer Sequences, Academic Press, 1973 (includes this sequence).
  • N. J. A. Sloane and Simon Plouffe, The Encyclopedia of Integer Sequences, Academic Press, 1995 (includes this sequence).

Crossrefs

Cf. A000040, A001223, A000101 (upper ends), A005250 (record gaps), A000230, A111870, A111943.
See also A205827(n) = A000040(A214935(n)), A182514(n) = A000040(A241540(n)).

Programs

  • Mathematica
    s = {2}; gm = 1; Do[p = Prime[n]; g = Prime[n + 1] - p; If[g > gm, Print[p]; AppendTo[s, p]; gm = g], {n, 2, 1000000}]; s   (* Jean-François Alcover, Mar 31 2011 *)
    Module[{nn=10^7,pr,df},pr=Prime[Range[nn]];df=Differences[pr];DeleteDuplicates[ Thread[ {Most[ pr],df}],GreaterEqual[#1[[2]],#2[[2]]]&]][[All,1]] (* The program generates the first 26 terms of the sequence. *) (* Harvey P. Dale, Sep 24 2022 *)
  • PARI
    a(n)=local(p,g);if(n<2,2*(n>0),p=a(n-1);g=nextprime(p+1)-p;while(p=nextprime(p+1),if(nextprime(p+1)-p>g,break));p) /* Michael Somos, Feb 07 2004 */
    
  • PARI
    p=q=2;g=0;until( g<(q=nextprime(1+p=q))-p && print1(q-g=q-p,","),) \\ M. F. Hasler, Dec 13 2007

Formula

a(n) = A000101(n) - A005250(n) = A008950(n-1) - 1. - M. F. Hasler, Dec 13 2007
A000720(a(n)) = A005669(n).
a(n) = A000040(A005669(n)). - M. F. Hasler, Apr 26 2014

Extensions

Definition clarified by Harvey P. Dale, Sep 24 2022

A111870 Prime p with prime gap q - p of n-th record merit, where q is smallest prime larger than p and the merit of a prime gap is (q-p)/log(p).

Original entry on oeis.org

2, 3, 7, 113, 1129, 1327, 19609, 31397, 155921, 360653, 370261, 1357201, 2010733, 17051707, 20831323, 191912783, 436273009, 2300942549, 3842610773, 4302407359, 10726904659, 25056082087, 304599508537, 461690510011, 1346294310749, 1408695493609, 1968188556461, 2614941710599, 13829048559701, 19581334192423, 218209405436543, 1693182318746371
Offset: 1

Views

Author

N. J. A. Sloane, based on correspondence with Ed Pegg Jr, Nov 23 2005

Keywords

Comments

As I understand it, the sequence refers to "Smallest prime p whose following gap has bigger merit than the other primes smaller than p." If that is the case, then it has an error. The sequence starts: 2, 3, 7, 113, 1129, 1327, 19609, 31397, 155921, 360653, 370261, 1357201, 4652353, 2010733, ... but you can see that 4652353 > 2010733, so in any case it should be listed after, not before it. But above that, its merit is 10.03 < 10.20, the merit of 2010733, so it is not in a mistaken position: it shouldn't appear in the sequence. - Jose Brox, Dec 31 2005
The logarithmic (base 10) graph seems to be linearly asymptotic to n with slope ~ 1/log(10) which would imply that: log(prime with n-th record merit) ~ n as n goes to infinity. - N. J. A. Sloane, Aug 27 2010
The sequence b(n) = (prime(n+1)/prime(n))^n is increasing for terms prime(n) of this sequence. - Thomas Ordowski, May 04 2012
The smallest prime(n) such that (prime(n+1)/prime(n))^n is increasing: 2, 3, 7, 23, 113, 1129, 1327, ... (A205827). - Thomas Ordowski, May 04 2012
(prime(n+1)/prime(n))^n > 1 + merit(n) for n > 2, where merit(n) = (prime(n+1)-prime(n))/log(prime(n)). - Thomas Ordowski, May 14 2012
Merit(1) + merit(2) + ... + merit(n) =: S(n) ~ n, where merit(n) is as above. - Thomas Ordowski, Aug 03 2012
For the index of a(n), see the comment at A214935. - John W. Nicholson, Nov 21 2013

Examples

			The first few entries correspond to the following gaps. The table gives n, p, gap = q-p and the merit of the gap.
   1,       2,   1, 1.4427
   2,       3,   2, 1.82048
   3,       7,   4, 2.05559
   4,     113,  14, 2.96147
   5,    1129,  22, 3.12985
   6,    1327,  34, 4.72835
   7,   19609,  52, 5.26116
   8,   31397,  72, 6.95352
   9,  155921,  86, 7.19238
  10,  360653,  96, 7.50254
  11,  370261, 112, 8.73501
  12, 1357201, 132, 9.34782
		

References

  • Ed Pegg, Jr., Posting to Seq Fan mailing list, Nov 23 2005

Crossrefs

For the gaps, see A111871.

Programs

  • Mathematica
    With[{s = Map[(#2 - #1)/Log[#1] & @@ # &, Partition[Prime@ Range[10^6], 2, 1]]}, Map[Prime@ FirstPosition[s, #][[1]] &, Union@ FoldList[Max, s]]] (* Michael De Vlieger, Jul 19 2018 *)

Formula

a(n) = A277552(n) - A111871(n). - Bobby Jacobs, Nov 13 2016

Extensions

Corrected by Jose Brox, Dec 31 2005
Corrected and edited by Daniel Forgues, Oct 23 2009
Further edited by Daniel Forgues, Nov 01 2009, Nov 13 2009, Nov 24 2009

A205827 Primes prime(k) corresponding to the records in the sequence (prime(k+1)/prime(k))^k.

Original entry on oeis.org

2, 3, 7, 23, 113, 1129, 1327, 19609, 31397, 155921, 360653, 370261, 1357201, 2010733, 17051707, 20831323, 191912783, 436273009, 2300942549, 3842610773, 4302407359, 10726904659, 25056082087, 304599508537, 461690510011, 1346294310749, 1408695493609
Offset: 1

Views

Author

Thomas Ordowski, May 07 2012

Keywords

Comments

Probably A111870 is this sequence with the exception of the term a(4) = 23. - Farideh Firoozbakht, May 07 2012
For n from 5 to 28, a(n) = A111870(n-1). - Donovan Johnson, Oct 26 2012
The statement prime(k) > (prime(k+1)/prime(k))^k for k>=1 is a rewrite of the Firoozbakht conjecture (see link). - John W. Nicholson, Oct 27 2012
Values of k are in A214935.
The logarithmic (base 10) graph seems to be linearly asymptotic to n with slope ~ 1/log(10) which would imply that: log(prime(k)) ~ n as n goes to infinity. [Copy of comment by N. J. A. Sloane, Aug 27 2010 for A111870, copied and corrected for prime(k) by John W. Nicholson, Oct 29 2012]
(prime(k+1)/prime(k))^k ~ e^merit(k), where merit(k) = (prime(k+1)-prime(k))/log(prime(k)). - Thomas Ordowski, Mar 18 2013
Subset of A002386. - John W. Nicholson, Nov 19 2013
Copied comment from A111870 (modified variable to k): (prime(k+1)/prime(k))^k > 1 + merit(k) for k > 2, where merit(k) = (prime(k+1)-prime(k))/log(prime(k)). - Thomas Ordowski, May 14 2012 : Copied and modified by John W. Nicholson, Nov 20 2013

Examples

			The sequence (prime(k+1)/prime(k))^k for k=1,2,... starts with:
*1.500, *2.777, 2.744, *6.098, 2.305, 5.001, 2.178, 4.611, *8.054, 1.948, ...,
where records are marked with *. The corresponding primes are a(1)=prime(1)=2, a(2)=prime(2)=3, a(3)=prime(4)=7, a(4)=prime(9)=23, ...
		

Crossrefs

Programs

  • Mathematica
    t = {}; p = 2; best = 0; n = 0; While[n++; last = p; p = NextPrime[p]; p <= 100000, f = (p/last)^n; If[f > best, best = f; AppendTo[t, last]]]; t (* T. D. Noe, May 08 2012 *)
  • PARI
    record=0;for(n=1,75,current=(A000101[n]/A002386[n]*1.)^A005669[n];if(current>record,record=current;print1(A002386[n],", "))) \\ Each sequence is read in as a vector as to overcome PARI's primelimit. John W. Nicholson, Dec 01 2013

Formula

a(n) = A000040(A214935(n)).

Extensions

a(13)-a(25) from Donovan Johnson, May 08 2012
Definition corrected by Max Alekseyev, Oct 23 2012
Clarified definition with k as index of a(n)=prime(k) instead of index n, John W. Nicholson, Oct 24 2012
a(26)-a(28) from Donovan Johnson, Oct 26 2012
a(29)-a(38) from John W. Nicholson, Dec 01 2013

A134266 Primes associated with the prime gaps listed in A085237.

Original entry on oeis.org

2, 3, 5, 7, 13, 19, 23, 31, 47, 53, 61, 73, 83, 89, 113, 293, 317, 523, 887, 1129, 1327, 8467, 9551, 12853, 14107, 15683, 19609, 25471, 31397, 155921, 338033, 360653, 370261, 492113, 1349533, 1357201, 1561919, 2010733, 4652353, 11113933, 15203977, 17051707, 20831323, 47326693, 122164747, 189695659, 191912783
Offset: 1

Views

Author

David W. Wilson, Dec 31 2007

Keywords

Comments

The smallest prime p(n) such that p(n+1)-p(n) is nondecreasing. The smallest prime p(n) such that (p(n+1)/p(n))^p(n) is increasing. [Thomas Ordowski, May 26 2012]
a(n) is the last prime in the n-th sublist of prime numbers defined in A348178. - Ya-Ping Lu, Oct 19 2021

Crossrefs

See also A205827(n) = A000040(A214935(n)), A182514(n) = A000040(A241540(n)).
Cf. A348178.

Programs

  • Python
    from sympy import nextprime; p, r = 2, 0
    while p < 2*10**8:
        q = nextprime(p); g = q - p
        if g >= r: print(p, end = ', '); r = g
        p = q # Ya-Ping Lu, Jan 23 2024

Formula

a(n) = A000040(A085500(n)). - M. F. Hasler, Apr 26 2014

A241540 Indices of primes p in A182514, i.e., a(n) = primepi(p) = A000720(A182514(n)).

Original entry on oeis.org

1, 2, 4, 30, 217, 49749629143526
Offset: 1

Views

Author

M. F. Hasler, Apr 25 2014

Keywords

Comments

a(6) = A214935(33) = A000720(A205827(33)).

Crossrefs

A085500 Indices of primes where nondecreasing gaps occur.

Original entry on oeis.org

1, 2, 3, 4, 6, 8, 9, 11, 15, 16, 18, 21, 23, 24, 30, 62, 66, 99, 154, 189, 217, 1059, 1183, 1532, 1663, 1831, 2225, 2810, 3385, 14357, 29040, 30802, 31545, 40933, 103520, 104071, 118505, 149689, 325852, 733588, 983015, 1094421, 1319945, 2850174, 6957876, 10539432, 10655462
Offset: 1

Views

Author

Farideh Firoozbakht, Aug 15 2003

Keywords

Comments

A005669 is a subsequence of this sequence.

References

  • Richard K. Guy, Unsolved Problems in Number Theory, 3rd Edition, Springer, 2004, Section A8, pp. 31-39.

Crossrefs

Programs

  • Mathematica
    f[n_] := Prime[n+1]-Prime[n]; v1={}; v2={}; Do[If[f[n]>=If[n==1, 1, Last[v2]], v=n; v1=Append[v1, n]; v2=Append[v2, f[v]]; Print[v1]], {n, 105000000}]

Formula

a(n) = A000720(A134266(n)). - M. F. Hasler, Apr 26 2014

Extensions

a(45)-a(47) from Amiram Eldar, Sep 05 2024

A241623 Indices of the primes in A073861 (n-digit primes followed by a maximal gap).

Original entry on oeis.org

2, 24, 154, 1183, 3385, 40933, 325852, 2850174, 23163298, 203615628, 1820471368, 28106444830, 251265078335, 2921439731020, 6822667965940, 49749629143526, 2133658100875638, 20004097201301079
Offset: 1

Views

Author

M. F. Hasler, Apr 26 2014

Keywords

Comments

"Indices" here means the value of the function primepi = A000720, cf. formula.
Subsequence of A005669.

Crossrefs

See also A214935(n) = A000720(A205827(n)), A005669(n) = A000720(A002386(n)).

Formula

a(n) = A000720(A073861(n)).
a(n) = A005669(k) with k such that A002386(k) = A073861(n).
Showing 1-8 of 8 results.