OFFSET
0,2
LINKS
Alois P. Heinz, Table of n, a(n) for n = 0..10000
FORMULA
G.f.: 1/(1-x)^7 * 1/Product_{j>1} (1-x^j).
Euler transform of 7,1,1,1,... .
a(n) ~ 2 * 3^(5/2) * n^2 * exp(Pi*sqrt(2*n/3)) / Pi^6. - Vaclav Kotesovec, Oct 24 2018
MAPLE
a:= proc(n) option remember; `if`(n=0, 1, add(
(numtheory[sigma](j)+6)*a(n-j), j=1..n)/n)
end:
seq(a(n), n=0..40);
MATHEMATICA
nmax = 50; CoefficientList[Series[1/((1-x)^6 * Product[1-x^k, {k, 1, nmax}]), {x, 0, nmax}], x] (* Vaclav Kotesovec, Oct 24 2018 *)
PROG
(PARI) x='x+O('x^30); Vec(1/((1-x)^7*prod(j=2, 40, 1-x^j))) \\ G. C. Greubel, Oct 27 2018
(Magma) m:=30; R<x>:=PowerSeriesRing(Integers(), m); Coefficients(R!(1/((1-x)^7*(&*[1-x^j: j in [2..30]])))); // G. C. Greubel, Oct 27 2018
CROSSREFS
KEYWORD
nonn
AUTHOR
Alois P. Heinz, Oct 20 2018
STATUS
approved