A198161 Primes from merging of 10 successive digits in decimal expansion of sqrt(2).
4142135623, 8872420969, 9698078569, 7537694807, 7973799073, 7846210703, 2644121497, 9935831413, 6592750559, 7010955997, 1472851741, 5251407989, 2533965463, 5339654633, 6152583523, 1525835239, 3950547457, 5750287759, 5996172983, 4084988471, 6668713013
Offset: 1
Links
- Vincenzo Librandi, Table of n, a(n) for n = 1..1000
Crossrefs
Programs
-
Mathematica
With[{len=10},Select[FromDigits/@Partition[RealDigits[Sqrt[2],10,1000][[1]],len,1],IntegerLength[#]==len&&PrimeQ[#]&]]
-
PARI
A198161(n, x=sqrt(2), m=10, silent=0)={m=10^m; for(k=1, default(realprecision), (isprime(p=x\.1^k%m)&&p*10>m)||next; silent||print1(p", "); n--||return(p))} \\ The optional arguments can be used to produce other sequences of this series (cf. Crossrefs). Use e.g. \p999 to set precision to 999 digits. - M. F. Hasler, Nov 02 2014
Comments