这是一种方法：
使用连分数法生成 e的第 1 1000 位， @quantum 在Code to Generate e one Digit at a Time中回答，这是来自 @wnoise 在生成 2 的平方根的数字中的回答，这是“Haskell 代码的改编......一直在漂浮”：
def z(contfrac, a=1, b=0, c=0, d=1):
for x in contfrac:
while a > 0 and b > 0 and c > 0 and d > 0:
t = a // c
t2 = b // d
if not t == t2:
break
yield t
a = (10 * (a - c*t))
b = (10 * (b - d*t))
# continue with same fraction, don't pull new x
a, b = x*a+b, a
c, d = x*c+d, c
for digit in rdigits(a, c):
yield digit
def rdigits(p, q):
while p > 0:
if p > q:
d = p // q
p = p - q * d
else:
d = (10 * p) // q
p = 10 * p - q * d
yield d
def e_cf_expansion():
yield 1
k = 0
while True:
yield k
k += 2
yield 1
yield 1
def e_dec():
return z(e_cf_expansion())
gen = e_dec()
e = [str(gen.next()) for i in xrange(1000)]
e.insert(1, '.')
用于测试从Rosetta Code Primality_by_trial_division#Python中选择的整数的素数的函数：
def isprime(a):
if a < 2: return False
if a == 2 or a == 3: return True # manually test 2 and 3
if a % 2 == 0 or a % 3 == 0: return False # exclude multiples of 2 and 3
maxDivisor = a**0.5
d, i = 5, 2
while d <= maxDivisor:
if a % d == 0: return False
d += i
i = 6 - i # this modifies 2 into 4 and viceversa
return True
找出 e 中的前 10 位素数（我的贡献）：
for i in range(len(e[2:])-10):
x = int(reduce(operator.add,e[2:][i:i+10]))
if isprime(x):
print x
print i
break
这打印：
7427466391
98
这意味着 e 中的前 10 位素数出现在小数点后的第 98 位，这与“答案的位置”下的http://explorepdx.com/firsten.html一致。
生成 e 的数字的一种更简单的方法是使用欧拉级数展开，它可以使用从具有 100 位精度的欧拉数 (Python)改编的代码来完成，该代码使用 Python 的 Decimal 类以获得足够的精度：
import operator
import decimal as dc
def edigits(p):
dc.getcontext().prec = p
factorial = 1
euler = 2
for x in range(2, 150):
factorial *= x
euler += dc.Decimal(str(1.0))/dc.Decimal(str(factorial))
return euler
estring = edigits(150).to_eng_string()[2:]
for i in range(len(estring)-10):
x = int(reduce(operator.add,estring[i:i+10]))
if isprime(x):
print x
print i
break
这打印：
7427466391
98
正如@MarkDickinson 所指出的，一种更简单的方法是直接使用 decimal 模块以必要的精度生成 e 。例如：
import operator
import decimal
decimal.getcontext().prec = 150
e_from_decimal = decimal.Decimal(1).exp().to_eng_string()[2:]
for i in range(len(e_from_decimal)-10):
x = int(reduce(operator.add,e_from_decimal[i:i+10]))
if isprime(x):
print x
print i
break
这打印：
7427466391
98