OFFSET
0,3
LINKS
Vaclav Kotesovec, Table of n, a(n) for n = 0..440
FORMULA
a(n)/n! ~ 3^(1/4) * exp(2*sqrt(6*n)/Pi) / (Pi * 2^(3/4) * n^(3/4)).
E.g.f.: Product_{k>=1} 1 / (1 - x^k)^f(k), where f(k) = (1/k) * Sum_{j=1..k} mu(gcd(k,j)). - Ilya Gutkovskiy, Aug 17 2021
a(0) = 1; a(n) = (n-1)! * Sum_{k=1..n} phi(k) * a(n-k)/(n-k)!. - Seiichi Manyama, Apr 29 2022
E.g.f.: exp( Sum_{n>=1} (mu(n)/n) * x^n/(1 - x^n) ), where mu(n) = A008683(n). - Paul D. Hanna, Jun 24 2023
MATHEMATICA
nmax = 20; CoefficientList[Series[Exp[Sum[EulerPhi[k]*x^k/k, {k, 1, nmax}]], {x, 0, nmax}], x] * Range[0, nmax]!
a[n_] := a[n] = If[n == 0, 1, Sum[EulerPhi[k]* a[n-k], {k, 1, n}]/n]; Table[n! a[n], {n, 0, 20}]
PROG
(PARI) a(n) = if(n==0, 1, (n-1)!*sum(k=1, n, eulerphi(k)*a(n-k)/(n-k)!)); \\ Seiichi Manyama, Apr 29 2022
CROSSREFS
KEYWORD
nonn
AUTHOR
Vaclav Kotesovec, Sep 05 2018
STATUS
approved