OFFSET
1,1
COMMENTS
Next term has n > 113821. - Eric W. Weisstein, Nov 04 2015
Also: concatenation of A007376(1 .. a(n)) is prime. - M. F. Hasler, Oct 23 2019
LINKS
Eric Weisstein's World of Mathematics, Champernowne Constant Digits
Eric Weisstein's World of Mathematics, Consecutive Number Sequences
Eric Weisstein's World of Mathematics, Constant Primes
Eric Weisstein's World of Mathematics, Integer Sequence Primes
Eric Weisstein's World of Mathematics, Smarandache Prime
MATHEMATICA
f[0] = 0; f[n_Integer] := 10^(Floor[Log[10, n]] + 1)*f[n - 1] + n; Do[If[PrimeQ[FromDigits[Take[IntegerDigits[f[n]], n]]], Print[n]], {n, 1, 3000}]
Cases[FromDigits /@ Rest[FoldList[Append, {}, RealDigits[N[ChampernowneNumber[], 1000]][[1]]]], p_?PrimeQ :> IntegerLength[p]] (* Eric W. Weisstein, Nov 04 2015 *)
PROG
(Python)
from itertools import count, islice
from sympy import isprime
def A071620_gen(): # generator of terms
c, l = 0, 0
for n in count(1):
for d in str(n):
c = 10*c+int(d)
l += 1
if isprime(c):
yield l
CROSSREFS
KEYWORD
nonn,base,hard,more
AUTHOR
Robert G. Wilson v, Jun 21 2002
EXTENSIONS
Edited by Charles R Greathouse IV, Apr 28 2010
a(6) = 4347 from Eric W. Weisstein, Jul 14 2013
a(7) = 37735 from Eric W. Weisstein, Jul 15 2013
STATUS
approved