login
A248804
Prime sieve of e.
6
18, 8284, 90452, 60287, 526624, 24, 36999, 669, 6, 24076630, 5945, 382, 85, 6, 427, 9, 3, 0, 2, 7, 66, 0, 57, 429526059, 0, 813, 2862, 4, 763, 5, 9525, 90, 573, 4, 21540, 934884, 46066808, 6847, 85, 4234544, 53, 77, 9, 551702, 8386062, 8, 7520, 38265, 760, 7
OFFSET
1,1
LINKS
Manfred Scheucher, Sage Script (Note: should also run in pure python with a few modifications).
PROG
(Python)
from decimal import *
getcontext().prec = 300
def primes(n):
""" Returns a list of primes < n """
sieve = [True] * n
for i in range(3, int(n**0.5)+1, 2):
if sieve[i]:
sieve[i*i::2*i]=[False]*((n-i*i-1)//(2*i)+1)
return [2] + [i for i in range(3, n, 2) if sieve[i]]
b = primes(429526060)
x = (Decimal(1).exp())
y = str(x).replace(".", "")
for x in b:
y = y.replace(str(x), " ", 1) #replace first occurrence only
f = [int(x) for x in y.split()[:-1]]
print(f)
# David Consiglio, Jr., Jan 03 2015
CROSSREFS
Cf. A001113, A245770 (prime sieve of Pi), A248831 (prime sieve of sqrt(2)).
Sequence in context: A159405 A265450 A213402 * A060617 A222202 A201986
KEYWORD
nonn,base
AUTHOR
Jared Kish, Oct 14 2014
EXTENSIONS
More terms from David Consiglio, Jr., Jan 03 2015
More terms from Manfred Scheucher, May 25 2015
STATUS
approved