|
MATHEMATICA
|
f[i_, j_] := Mod[2 i + 2 j, 3]; (* symmetric *)
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, 14}, {i, 1, n}]] (* A204429 *)
Permanent[m_] :=
With[{a = Array[x, Length[m]]},
Coefficient[Times @@ (m.a), Times @@ a]];
Table[Permanent[m[n]], {n, 1, 22}] (* A204430 *)
|