%I #10 Oct 19 2013 03:58:36
%S 2,2,9,40,144,1032,7824,48816,530712,5619456,50015232,707802624,
%T 9591229440,112586112000,1938416832000,32061464064000,469224244224000,
%U 9509523337728000,185772505605120000,3261008730405888000
%N Permanent of n X n matrix M with M(i,j) = (i+j) (mod 3).
%C Permanent of the n-th principal submatrix of A204263.
%t f[i_, j_] := Mod[i + j, 3];
%t m[n_] := Table[f[i, j], {i, 1, n}, {j, 1, n}]
%t TableForm[m[8]] (* 8x8 principal submatrix *)
%t Flatten[Table[f[i, n + 1 - i],
%t {n, 1, 14}, {i, 1, n}]] (* A204263 *)
%t Permanent[m_] :=
%t With[{a = Array[x, Length[m]]},
%t Coefficient[Times @@ (m.a), Times @@ a]];
%t Table[Permanent[m[n]], {n, 1, 22}] (* A204265 *)
%o (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
%Y Cf. A204263.
%K nonn
%O 1,1
%A _Clark Kimberling_, Jan 15 2012
|