# Author: Manfred Scheucher # Date : 16.08.2015 sqrt2 = RealField(2^22)(sqrt(2)) A=dict() def a(n): if n not in A: A[n] = ceil(1/(sqrt2 - sum(1/a(j) for j in [0..n-1]))) return A[n] for n in [0..15]: print n,a(n)