OFFSET
0,3
LINKS
Alois P. Heinz, Table of n, a(n) for n = 0..750
FORMULA
a(n) ~ c * d^n / n^(3/2), where d = 18.843901825822305757579605844910623225182677164912157727486865738745..., c = 0.056432508184666421552433382560980970822794... . - Vaclav Kotesovec, Feb 24 2015
From Ilya Gutkovskiy, Apr 13 2019: (Start)
G.f. A(x) satisfies: A(x) = x*exp(7*Sum_{k>=1} (-1)^(k+1)*A(x^k)/k).
G.f.: A(x) = Sum_{n>=1} a(n)*x^n = x * Product_{n>=1} (1 + x^n)^(7*a(n)). (End)
MAPLE
with(numtheory):
a:= proc(n) option remember; `if`(n<2, n, -add(a(n-j)*add(
7*a(d)*d*(-1)^(j/d), d=divisors(j)), j=1..n-1)/(n-1))
end:
seq(a(n), n=0..30);
CROSSREFS
KEYWORD
nonn
AUTHOR
Alois P. Heinz, Feb 24 2015
STATUS
approved