A136583 n such that floor(sqrt(10^(2*n-1))) is (probably) prime.
1, 2, 7, 18, 33, 206, 468, 1061, 6831, 40377
Offset: 1
Programs
-
Magma
for n in [1..10^6] do if IsPrime(Isqrt(10^(2*n-1))) then printf "%o, ", n; end if; end for; // Jason Kimberley, Sep 03 2011
-
Mathematica
rd = RealDigits[Sqrt[10], 10, 10^5][[1]]; Do[ If[ PrimeQ@ FromDigits@ Take[rd, n], Print@n], {n, 10^5}] (* Robert G. Wilson v, Jan 20 2008 *)
Extensions
a(6) - a(8) from Robert G. Wilson v, Jan 20 2008
Probable terms a(9) and a(10) from Jason Kimberley, Aug 19 and Sep 03 2011
Comments