A204157 Array: row n shows the coefficients of the characteristic polynomial of the n-th principal submatrix of max(3i-j, 3j-i), as in A204156.

1, -1, -13, -4, 1, 88, 78, 9, -1, -496, -704, -260, -16, 1, 2560, 4960, 3080, 650, 25, -1, -12544, -30720, -26784, -9856, -1365, -36, 1, 59392, 175616, 197120, 104160, 25872, 2548, 49, -1, -274432, -950272, -1304576, -901120, -327360, -59136, -4368, -64, 1, 1245184, 4939776, 8017920, 6849024, 3294720

Offset: 1

Author

Clark Kimberling, Jan 12 2012

Keywords

• tabl
• sign

Comments

Let p(n)=p(n,x) be the characteristic polynomial of the n-th principal submatrix. The zeros of p(n) are real, and they interlace the zeros of p(n+1). See A202605 and A204016 for guides to related sequences.

Examples

			Top of the array:
 1....-1
-13...-4.....1
 88....78....9.....-1
-496..-704..-260...-16...1
		

References

• (For references regarding interlacing roots, see A202605.)

Crossrefs

Cf. A204156, A202605, A204016.

Programs

Mathematica

f[i_, j_] := -1 + Max[3 i - j, 3 j - i];
m[n_] := Table[f[i, j], {i, 1, n}, {j, 1, n}]
TableForm[m[8]] (* 8x8 principal submatrix *)
Flatten[Table[f[i, n + 1 - i],
  {n, 1, 15}, {i, 1, n}]]    (* A204156 *)
p[n_] := CharacteristicPolynomial[m[n], x];
c[n_] := CoefficientList[p[n], x]
TableForm[Flatten[Table[p[n], {n, 1, 10}]]]
Table[c[n], {n, 1, 12}]
Flatten[%]                   (* A204157 *)
TableForm[Table[c[n], {n, 1, 10}]]