editing
approved
editing
approved
a(n) ~ c * d^n, where d = 1.3864622092472465020397266918102624708859968795203700659786636158522760956... and c = 0.15945087310540003725148530084775272562567007586487061850065597143186... - Vaclav Kotesovec, Aug 25 2017
approved
editing
reviewed
approved
proposed
reviewed
editing
proposed
R:= 1; :
C:= numtheory:-cfrac([0, [1, 1], seq([-x^ithprime(i), 1], i=1..15)]):
R:= 1;
for i from numtheory:-pi(50) to 1 by -1 do
S R:= series(C, 1-x^ithprime(i)/R, x, 51):;
od:
R:= series(1/R, x, 51):
seq(coeff(S, R, x, j), j=0..50); # Robert Israel, Apr 20 2017
Robert Israel, <a href="/A285407/b285407.txt">Table of n, a(n) for n = 0..5000</a>
C:= numtheory:-cfrac([0, [1, 1], seq([-x^ithprime(i), 1], i=1..15)]):
S:= series(C, x, 51):
seq(coeff(S, x, j), j=0..50); # Robert Israel, Apr 20 2017
approved
editing
proposed
approved