OFFSET
1,5
COMMENTS
A057731 contains zeros. This sequence contains only positive values of A057731(n,k)/A000010(k). - Alois P. Heinz, Feb 16 2013
LINKS
Alois P. Heinz, Rows n = 1..42, flattened
EXAMPLE
Triangle begins:
1;
1, 1;
1, 3, 1;
1, 9, 4, 3;
1, 25, 10, 15, 6, 10;
1, 75, 40, 90, 36, 120;
1, 231, 175, 420, 126, 735, 120, 126, 105;
...
MATHEMATICA
nmax = 10;
T[n_, k_] := n! SeriesCoefficient[O[x]^(n+1) + Sum[MoebiusMu[k/i]*Exp[ Sum[x^j/j, {j, Divisors[i]}]], {i, Divisors[k]}], {x, 0, n}]/ EulerPhi[k];
Table[DeleteCases[Table[T[n, k], {k, 1, 2 nmax}], 0], {n, 1, nmax}] // Flatten (* Jean-François Alcover, Sep 16 2019, after Andrew Howroyd *)
PROG
(PARI) T(n, k)={n!*polcoeff(sumdiv(k, i, moebius(k/i)*exp(sumdiv(i, j, x^j/j) + O(x*x^n))), n)/eulerphi(k)} \\ Andrew Howroyd, Jul 02 2018
CROSSREFS
KEYWORD
nonn,tabf
AUTHOR
Vladeta Jovovic, Sep 30 2002
STATUS
approved