%I #23 Apr 24 2022 06:31:07
%S 1,1,13,55,1235,4615,55801,200343,8977475,36804235,367235363,
%T 1444888289,32062742231,120729974115,1205864254225,5201022002071,
%U 395884671433315,1603069490974835,15989295873680415,64312573140322525,1250332447587844829,5262481040435242585
%N Numerators in expansion of Euler transform of b(n) = 1/4.
%C Denominators of c(n) are 2^A004134(n).
%H Alois P. Heinz, <a href="/A061161/b061161.txt">Table of n, a(n) for n = 0..1000</a>
%H Geoffrey B. Campbell and A. Zujev, <a href="http://zujev.physics.ucdavis.edu/papers/Some%20almost%20partition%20theoretic%20identities.pdf">Some almost partition theoretic identities</a>, Preprint, 2016.
%H N. J. A. Sloane, <a href="/transforms.txt">Transforms</a>
%F Numerators of c(n), where c(n) = (1/(4*n))*Sum_{k=1..n} c(n-k)*sigma(k), n>0, c(0)=1.
%p b:= proc(n) option remember; `if`(n=0, 1, add(add(
%p d/4, d=numtheory[divisors](j))*b(n-j), j=1..n)/n)
%p end:
%p a:= n-> numer(b(n)):
%p seq(a(n), n=0..30); # _Alois P. Heinz_, Jul 28 2017
%t c[n_] := c[n] = If[n == 0, 1,
%t (1/(4n)) Sum[c[n-k] DivisorSigma[1, k], {k, 1, n}]];
%t a[n_] := Numerator[c[n]];
%t Table[a[n], {n, 0, 30}] (* _Jean-François Alcover_, Apr 24 2022 *)
%Y Cf. A000712, A061159, A061160.
%K easy,nonn,frac
%O 0,3
%A _Vladeta Jovovic_, Apr 17 2001
|