|
|
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
(list;
graph;
refs;
listen;
history;
text;
internal format)
|
|
|
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)
|
|
CROSSREFS
|
|
|
KEYWORD
|
nonn,base
|
|
AUTHOR
|
|
|
EXTENSIONS
|
|
|
STATUS
|
approved
|
|
|
|