OFFSET
0,3
COMMENTS
This sequence differs from A293110 first at n=8.
LINKS
Alois P. Heinz, Table of n, a(n) for n = 0..1000
FORMULA
G.f.: Product_{j>=1} 1/(1-x^j)^A007578(j).
a(n) ~ c * 7^n / n^(21/2), where c = 233774941.39802934196800791705821024006230754487492494942398064537776753785... - Vaclav Kotesovec, May 30 2019
MAPLE
g:= proc(n) option remember; `if`(n<4, [1, 1, 2, 4][n+1],
((4*n^3+78*n^2+424*n+495)*g(n-1) +(n-1)*(34*n^2+280*n+
305)*g(n-2) -2*(n-1)*(n-2)*(38*n+145)*g(n-3) -105*(n-1)
*(n-2)*(n-3)*g(n-4))/((n+6)*(n+10)*(n+12)))
end:
a:= proc(n) option remember; `if`(n=0, 1, add(add(g(d)
*d, d=numtheory[divisors](j))*a(n-j), j=1..n)/n)
end:
seq(a(n), n=0..35);
CROSSREFS
KEYWORD
nonn
AUTHOR
Alois P. Heinz, Oct 15 2017
STATUS
approved