OFFSET
0,3
LINKS
A. Knopfmacher and J. Knopfmacher, Counting irreducible factors of polynomials over a finite field, Discrete Math, 112 (1993), pp. 103-118.
FORMULA
G.f.: (1/(1-y*x)^2)*Product_{i>=2} 1/(1-x^i)^A001037(i).
Explicit formula given in Knopfmacher link.
EXAMPLE
1;
0, 2;
1, 0, 3;
2, 2, 0, 4;
4, 4, 3, 0, 5;
8, 8, 6, 4, 0, 6;
16, 16, 12, 8, 5, 0, 7;
32, 32, 24, 16, 10, 6, 0, 8;
MATHEMATICA
nn = 10; i[q_, r_] := 1/r Sum[MoebiusMu[r/d] q^d, {d, Divisors[r]}];
M[q_, n_, k_, r_] := If[0 <= k <= Floor[n/r - i[q, r]], Binomial[i[q, r] + k - 1, k]*q^(n - k*r) (1 - 1/q^r)^i[q, r], Binomial[i[q, r] + k - 1, k]* q^(n - k*r) Sum[(-1)^m Binomial[i[q, r], m] q^(-r m), {m, 0, Floor[n/r] - k}]]; Table[Table[M[2, n, k, 1], {k, 0, n}], {n, 0, nn}] // Grid
CROSSREFS
KEYWORD
nonn,tabl
AUTHOR
Geoffrey Critzer, Aug 13 2022
STATUS
approved