OFFSET
1,1
COMMENTS
Permanent of the n-th principal submatrix of A204263.
MATHEMATICA
f[i_, j_] := Mod[i + j, 3];
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}]] (* A204263 *)
Permanent[m_] :=
With[{a = Array[x, Length[m]]},
Coefficient[Times @@ (m.a), Times @@ a]];
Table[Permanent[m[n]], {n, 1, 22}] (* A204265 *)
PROG
(PARI) matperm(M)=my(n=#M, t); sum(i=1, n!, t=numtoperm(n, i); prod(j=1, n, M[j, t[j]])); a(n)=matperm(matrix(n, n, i, j, (i+j)%3)) \\ Ralf Stephan, Oct 19 2013
CROSSREFS
KEYWORD
nonn
AUTHOR
Clark Kimberling, Jan 15 2012
STATUS
approved