%I #9 Mar 08 2015 00:06:35
%S 2,3,13,179,8161,17319373,24288185721827,115965508643355137277722177,
%T 235520782452648150490328050427067627593114827790363,
%U 34927738441197833781363566983541975805101681093066970452740268721673910282017002081858804680098723909
%N Denominators of an Egyptian fraction for Catalan constant, using only prime numbers and allowing repetitions.
%C Good approximation up to 789 decimal digits.
%D Mohammad K. Azarian, Sylvester's Sequence and the Infinite Egyptian Fraction Decomposition of 1, Problem 958, College Mathematics Journal, Vol. 42, No. 4, September 2011, p. 330. Solution published in Vol. 43, No. 4, September 2012, pp. 340-342.
%p P:=proc(n) local a,i; a:=evalf(Catalan,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-A139519, A139521-A139523.
%K easy,nonn,frac
%O 0,1
%A _Paolo P. Lava_ and _Giorgio Balzarotti_, Apr 28 2008