OFFSET
0,9
LINKS
Robert Israel, Table of n, a(n) for n = 0..10010 (rows 0 to 140, flattened)
EXAMPLE
The characteristic polynomials of M_0, M_1, ..., M_4 are 1, -x + 0, x^2 + 0x + 0, -x^3 + 0x^2 + 2x + 0 and x^4 +0 x^3 - 5x^2 - 2x + 0, so the first 15 terms are 1; -1,0; 1,0,0; -1,0,2,0; 1,0,-5,-2,0.
MAPLE
printf("[1]\n"):
for n from 1 to 10 do
M:= Matrix(n, n, (i, j) -> i mod j);
P:= (-1)^n*LinearAlgebra:-CharacteristicPolynomial(M, t);
C:= seq(coeff(P, t, i), i=n..0, -1);
printf("%a\n", [C]);
od: # Robert Israel, Dec 04 2024
MATHEMATICA
c[0]={1}; c[n_] := Reverse[CoefficientList[CharacteristicPolynomial[Table[Mod[i, j], {i, 1, n}, {j, 1, n}], x], x]]; Flatten[c/@Range[0, 10]]
CROSSREFS
KEYWORD
AUTHOR
Benoit Cloitre, Dec 15 2002
STATUS
approved