Location via proxy:   [ UP ]  
[Report a bug]   [Manage cookies]                
A104851
Primes from merging of 10 successive digits in decimal expansion of e.
23
7427466391, 7413596629, 6059563073, 3490763233, 2988075319, 1573834187, 7021540891, 5408914993, 6480016847, 9920695517, 1838606261, 6062613313, 3845830007, 1692836819, 4425056953, 2505695369, 5490598793, 1782154249, 8215424999, 9229576351, 9519366803
OFFSET
1,1
COMMENTS
Scan decimal expansion of e from left to right, recording any 10-digit primes seen. - N. J. A. Sloane, Feb 05 2012
All the primes listed here must have 10 digits, i.e., "leading zeros are not allowed". Otherwise, one would also have some terms as 297606737 or 865746377 or 98793127 from A104850. - M. F. Hasler, Nov 01 2014
The original version read (1185790117, 1180978417, 1573834187, 1838606261, 1308008771, 1692836819, 1782154249, 1825288693, 1525971943, 1730123819, 1332069811, 1881593041, 1934580727, 1978623209, 1164218399, 1574862173, 1635834619, 1311914371, ...). These terms are obtained when using signed 32-bit integers, i.e., take the 10-digit numbers modulo 2^32, and select the primes between 10^9 and 2^31. - M. F. Hasler, Nov 01 2014
LINKS
MATHEMATICA
With[{de=FromDigits/@Partition[RealDigits[E, 10, 10000][[1]], 10, 1]}, Select[de, #>10^9&&PrimeQ[#]&]] (* Harvey P. Dale, Feb 05 2012 *)
PROG
(PARI) list_A104851(x=exp(1), m=10)=m=10^m; for(k=1, default(realprecision), isprime(p=x\.1^k%m)&&p*10>m&&print1(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 01 2014
CROSSREFS
KEYWORD
nonn,base
AUTHOR
Andrew G. West (WestA(AT)wlu.edu), Mar 27 2005
EXTENSIONS
Corrected by Harvey P. Dale, Feb 05 2012
Offset changed from 0 to 1 by Vincenzo Librandi, Apr 21 2013
STATUS
approved