A282497 - OEIS

%I #7 Feb 19 2017 14:04:56

%S 2,1,3,1,1,2,1,3,1,2,1,3,1,1,2,2,1,1,1,3,1,1,2,1,1,3,1,1,1,2,1,1,3,1,

%T 1,1,1,1,2,1,1,3,1,2,1,1,1,3,1,2,1,1,2,1,2,1,2,1,3,1,1,2,1,1,3,1,2,1,

%U 2,1,1,2,1,2,1,1,2,1,1,2,1,1,2,1,1,2,1,2,1,1,3,1,2,1,3,1,2,1,2,1,3,1,1,2,1,2,1,1,1,2,1,1,3,1,1,1,1,3,1,1,1,1,2,1,1,3,1,1,2,1

%N 'Somos expansion' of e: e=a(0)*sqrt(a(1)*sqrt(a(2)*sqrt(a(3)*sqrt(...)))). a(n)=floor(x(n)), x(n)=x(n-1)^2/a(n-1)^2, x(0)=e.

%C 1<=a(n)<=3 for all n. Reasoning: for x>1 it follows that 1<x/floor(x)<2.

%H Yuriy Sibirmovsky, <a href="/A282497/b282497.txt">Table of n, a(n) for n = 0..1999</a>

%F Product_{k>=0} a(k)^(1/2^k) = e.

%e Integer part of e is 2. Integer part of e^2/4 is 1.

%t $MaxExtraPrecision = 1000;

%t x0 = E;

%t Nm = 130;

%t j = 1;

%t Res = Table[1, {j, 1, Nm}];

%t While[j < Nm, Res[[j]] = Floor[x0]; x0 = N[(x0/ Res[[j]])^2, 20000];

%t   j++];

%Y Cf. A001113 (digits).

%K nonn

%O 0,1

%A _Yuriy Sibirmovsky_, Feb 16 2017