%I #3 Aug 24 2012 10:50:01
%S 2,13,3433,770741,9157546933,68379695298643311959,
%T 177698848753108581519320740556520784033,
%U 191429752162010655176200572455622919752517437701651231352646653127630507601
%N Denominators of an Egyptian fraction for Euler-Mascheroni constant, using only prime numbers and allowing repetitions.
%C Good approximation up to 801 decimal digits
%p P:=proc(n) local a,i; a:=evalf(gamma,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-A139518, A139520-A139523.
%K easy,nonn,frac
%O 0,1
%A _Paolo P. Lava_ and _Giorgio Balzarotti_, Apr 28 2008