# sum[n=1..inf]n^2/(n+1)! = e-1
# sum[n=1..inf]n^2/n! = 2e
import math
a = [1]
p = 1
for i in range(1, 11):
p = p * i
a.append(p)
print(a)
s = 0
t = 0
for i in range(1, 10):
s = s + (i*i) / a[i+1]
t = t + (i*i) / a[i]
print(s, t, math.e)
IyBzdW1bbj0xLi5pbmZdbl4yLyhuKzEpISA9IGUtMQojIHN1bVtuPTEuLmluZl1uXjIvbiEgPSAyZQoKaW1wb3J0IG1hdGgKCmEgPSBbMV0KcCA9IDEKZm9yIGkgaW4gcmFuZ2UoMSwgMTEpOgoJcCA9IHAgKiBpCglhLmFwcGVuZChwKQpwcmludChhKQoKcyA9IDAKdCA9IDAKZm9yIGkgaW4gcmFuZ2UoMSwgMTApOgoJcyA9IHMgKyAoaSppKSAvIGFbaSsxXQoJdCA9IHQgKyAoaSppKSAvIGFbaV0KcHJpbnQocywgdCwgbWF0aC5lKQo=
[1, 1, 2, 6, 24, 120, 720, 5040, 40320, 362880, 3628800]
1.718279045414462 5.436532738095238 2.718281828459045