|
%I #16 Nov 12 2023 13:00:36
%S 1,2,10,80,910,14180,294820,7898900,262166950,10401023300,
%T 480352145980,25403319156500,1524514142067820,103259722304891000,
%U 7855753426977689320,667222680003329539400,62774912079259300426030,6482982674137778455277720,728040412975153836997647580
%N Expansion of 1/(1-2x/(1-3x/(1-5x/(1-7x/(1-11x/...))))) where 2,3,5,7,11,... are the prime numbers.
%H Vaclav Kotesovec, <a href="/A269353/b269353.txt">Table of n, a(n) for n = 0..300</a>
%F G.f.: 1/(1-2x/(1-3x/(1-5x/(1-7x/(1-11x/(1-13x/(...))))))), by definition.
%p b:= proc(x, y) option remember; `if`(x=0 and y=0, 1,
%p `if`(x>0, b(x-1, y)*ithprime(y-x+1), 0)+`if`(y>x, b(x, y-1), 0))
%p end:
%p a:= n-> b(n$2):
%p seq(a(n), n=0..20); # _Alois P. Heinz_, Nov 12 2023
%t CoefficientList[Series[1/(1+ContinuedFractionK[-Prime[k]*x,1,{k,1,50}]),{x,0,50}],x]
%Y Cf. A000040, A285407.
%K nonn
%O 0,2
%A _Benedict W. J. Irwin_, Feb 25 2016
%E Offset corrected by _Vaclav Kotesovec_, Aug 26 2017
|
