OFFSET
0,2
COMMENTS
Euler transform of 3 * A046211(n).
LINKS
Seiichi Manyama, Table of n, a(n) for n = 0..1000
OEIS Wiki, Euler transform
FORMULA
G.f.: A(x) = 1 / ( Product_{k>=1} (1 - x^k)^A046211(k) )^3.
EXAMPLE
A(x)^3 = 1 + 9*x + 54*x^2 + 273*x^3 + 1242*x^4 + 5265*x^5 + 21231*x^6 + ... .
PROG
(PARI) b(n, k) = sumdiv(n, d, (gcd(d, k)==1)*(moebius(d)*k^(n/d)))/(k*n);
my(N=30, x='x+O('x^N)); Vec(1/prod(k=1, N, (1 - x^k)^b(k, 3))^3)
CROSSREFS
KEYWORD
nonn
AUTHOR
Seiichi Manyama, Jul 04 2024
STATUS
approved