login
Denominators of an Egyptian fraction for log(2), using only prime numbers and allowing repetitions.
2

%I #5 Jun 29 2015 02:08:51

%S 2,7,23,149,9967,27514363,217357018990141,

%T 1606507138380520696666019221,

%U 25278175600230973831063754126661264411663401021548343,2658247789993020763205588174060503375422404932907883834101966812060033051648561545028435381434723862291

%N Denominators of an Egyptian fraction for log(2), using only prime numbers and allowing repetitions.

%C Good approximation up to 800 decimal digits.

%p P:=proc(n) local a,i; a:=evalf(ln(2),100); for i from 1 by 1 to n do if 1/ithprime(i)<a then a:=a-1/ithprime(i); print(a); print(ithprime(i)); fi; od; end: P(100000);

%Y Cf. A139514-A139521, A139523.

%K easy,nonn,frac

%O 0,1

%A _Paolo P. Lava_ and _Giorgio Balzarotti_, Apr 28 2008