A338438 a(n) = n! * [x^n] Product_{k=1..n, gcd(n,k) = 1} (1 + x^k/k!).

1, 1, 0, 3, 4, 15, 6, 168, 64, 171, 130, 44418, 804, 802750, 2380, 683595, 5782144, 840363295, 40410, 24358246734, 221953840, 3114081516, 5372831860, 62269802986834, 25703183400, 548547778815000, 3028000483716, 9604167089258628, 30673543523224, 13242158988496348746

Offset: 0

Author

Ilya Gutkovskiy, Oct 27 2020

Keywords

• nonn

Crossrefs

Cf. A007837, A036998, A338436.

Programs

Mathematica

Table[n! SeriesCoefficient[Product[(1 + Boole[GCD[n, k] == 1] x^k/k!), {k, 1, n}], {x, 0, n}], {n, 0, 29}]