login
The OEIS is supported by the many generous donors to the OEIS Foundation.

 

Logo
Hints
(Greetings from The On-Line Encyclopedia of Integer Sequences!)
A065582 Smallest prime ending in exactly n 9's. 11
19, 199, 1999, 49999, 199999, 2999999, 19999999, 799999999, 10999999999, 59999999999, 1099999999999, 34999999999999, 59999999999999, 499999999999999, 14999999999999999, 139999999999999999 (list; graph; refs; listen; history; text; internal format)
OFFSET
1,1
COMMENTS
From Jeppe Stig Nielsen, Jul 30 2022: (Start)
Can decrease, for example a(25) < a(24). So not the same as Smallest prime ending in n or more 9s.
a(n) can contain other 9s as well, for example a(46), a(118), a(156). (End)
LINKS
Hugo Pfoertner, Table of n, a(n) for n = 1..500 (first 100 terms from Harry J. Smith).
MATHEMATICA
Do[a = Table[9, {n} ]; k = 0; While[ b = FromDigits[ Join[ IntegerDigits[k], a]]; Mod[k, 10] == 9 || !PrimeQ[b], k++ ]; Print[b], {n, 1, 17} ]
pe9[n_]:=Module[{k=1, rh=FromDigits[PadRight[{}, n, 9]]}, While[!PrimeQ[ k 10^n+rh], k++]; k 10^n+rh]; Array[pe9, 20] (* Harvey P. Dale, Mar 26 2012 *)
PROG
(PARI) { for (n=1, 100, t=10^n; b=t - 1; d=0; while(!isprime(b + t*d), d++; if (d%10==9, d++)); write("b065582.txt", n, " ", b + t*d) ) } \\ Harry J. Smith, Oct 23 2009
(Python)
from sympy import isprime
def a(n):
pow, end, i = 10**n, 10**n-1, 1
while i%10 == 9 or not isprime(i*pow + end): i += 1
return i*pow + end
print([a(n) for n in range(1, 20)]) # Michael S. Branicky, Jul 30 2022
CROSSREFS
Sequence in context: A135162 A185687 A067272 * A241021 A055558 A160452
KEYWORD
nonn,base
AUTHOR
Robert G. Wilson v, Nov 28 2001
STATUS
approved

Lookup | Welcome | Wiki | Register | Music | Plot 2 | Demos | Index | Browse | More | WebCam
Contribute new seq. or comment | Format | Style Sheet | Transforms | Superseeker | Recents
The OEIS Community | Maintained by The OEIS Foundation Inc.

License Agreements, Terms of Use, Privacy Policy. .

Last modified April 25 10:22 EDT 2024. Contains 371967 sequences. (Running on oeis4.)