%I #25 Jul 05 2023 17:02:02
%S 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,
%T 126,735,120,126,105,1,763,616,2730,336,5320,960,1260,1008,840,336,1,
%U 2619,2884,15498,756,41580,4320,11340,6720,6804,7560,4320,3024,2268
%N Triangle T(n,k) giving number of labeled cyclic subgroups of order k in symmetric group S_n, n >= 1, 1 <= k <= g(n), where g(n) = A000793(n) is Landau's function.
%C A057731 contains zeros. This sequence contains only positive values of A057731(n,k)/A000010(k). - _Alois P. Heinz_, Feb 16 2013
%H Alois P. Heinz, <a href="/A074881/b074881.txt">Rows n = 1..42, flattened</a>
%F T(n,k) = A057731(n,k)/A000010(k).
%e Triangle begins:
%e 1;
%e 1, 1;
%e 1, 3, 1;
%e 1, 9, 4, 3;
%e 1, 25, 10, 15, 6, 10;
%e 1, 75, 40, 90, 36, 120;
%e 1, 231, 175, 420, 126, 735, 120, 126, 105;
%e ...
%t nmax = 10;
%t 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];
%t 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_ *)
%o (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
%Y Row sums give A051625.
%Y Cf. A000010, A181949.
%K nonn,tabf
%O 1,5
%A _Vladeta Jovovic_, Sep 30 2002