OFFSET
1,1
COMMENTS
To shorthand the n-th - smallest n-digit prime it is convenient to subtract 10^(n-1) (n>1). Compare a(n) with A069100(n).
LINKS
David A. Corneth, Table of n, a(n) for n = 1..1000 (using b-file from A069100)
FORMULA
a(1)=2; at n>1 a(n)=prime(pi[10^(n-1)]+n)-10^(n-1)=A069100(n)-10^(n-1).
PROG
(PARI) a(n)=prime(primepi(10^(n-1))+n)-if(n==1, 0, 10^(n-1)) \\ Franklin T. Adams-Watters, Mar 07 2014
(PARI) a(n) = {if(n == 1, return(2)); my(t = 0); forprime(p = 10^(n-1), 10^n, t++; if(t==n, return(p - 10^(n-1))))} \\ David A. Corneth, Jun 16 2021
(Python)
from sympy import nextprime
def a(n): return nextprime(10**(n-1), ith=n) - 10**(n-1) * (n > 1)
print([a(n) for n in range(1, 51)]) # Michael S. Branicky, Jun 16 2021
CROSSREFS
KEYWORD
nonn,base
AUTHOR
Zak Seidov, May 12 2005
EXTENSIONS
a(3) corrected by Franklin T. Adams-Watters, Mar 07 2014
More terms from David A. Corneth, Jun 16 2021
STATUS
approved
