A139053 - OEIS

%I #18 May 07 2021 08:45:17

%S 2,3,5,11,13,17,101,103,107,1009,1013,1019,10007,10009,10037,100003,

%T 100019,100043,1000003,1000033,1000037,10000019,10000079,10000103,

%U 100000007,100000037,100000039,1000000007,1000000009,1000000021,10000000019,10000000033,10000000061

%N Array read by rows: row n lists the first 3 primes with n digits.

%H Michael S. Branicky, <a href="/A139053/b139053.txt">Table of n, a(n) for n = 1..3000</a>

%t np[n_]:=Module[{a=NextPrime[10^n]},{a,NextPrime[a], NextPrime[NextPrime[ a]]}]; Flatten[Array[np,12,0]] (* _Harvey P. Dale_, Dec 14 2011 *)

%t Flatten@Array[NextPrime[10^#,{1,2,3}]&,12,0] (* _Giorgos Kalogeropoulos_, May 06 2019 *)

%o (Python)

%o from sympy import nextprime

%o def auptodigs(maxdigits):

%o   alst = []

%o   for n in range(1, maxdigits+1):

%o     p1 = nextprime(10**(n-1))

%o     p2 = nextprime(p1)

%o     p3 = nextprime(p2)

%o     alst.extend([p1, p2, p3])

%o   return alst

%o print(auptodigs(11)) # _Michael S. Branicky_, May 07 2021

%Y Cf. A000040, A003617, A073914, A139052, A139054.

%K base,nonn,tabf

%O 1,1

%A _Omar E. Pol_, Apr 08 2008

%E More terms from _Max Alekseyev_, Dec 12 2011

%E a(32) and beyond from _Michael S. Branicky_, May 07 2021