OFFSET
0,3
LINKS
Seiichi Manyama, Table of n, a(n) for n = 0..426
FORMULA
a(0) = 1 and a(n) = (-1) * (n-1)! * Sum_{k=1..n} k*A001157(k)*a(n-k)/(n-k)! for n > 0.
E.g.f.: Product_{k>=1} (1 - x^k)^f(k), where f(k) = (1/k) * Sum_{j=1..k} gcd(k,j)^3. - Ilya Gutkovskiy, Aug 17 2021
PROG
(PARI) N=66; x='x+O('x^N); Vec(serlaplace(exp(-sum(k=1, N, sigma(k, 2)*x^k))))
CROSSREFS
KEYWORD
sign
AUTHOR
Seiichi Manyama, Oct 30 2017
STATUS
approved