%I #12 Jan 12 2022 09:28:09
%S 1,1,1,2,2,2,2,2,2,2,2,2,2,2,2,2,2,2,2,2,2,3,3,3,3,3,3,3,3,3,3,3,3,3,
%T 3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,3,4,4,4,4,4,4,4,4,4,4,4,4,4,
%U 4,4,4,4,4,4,4,4,4,4,4,4,4,4,4,4,4,4,4,4,4,4,4,4,4,4,4,4,4,4,4
%N Symmetric matrix: f(i,j)=floor[(2i+2j+4)/4]-floor[(i+j+2)/4], by (constant) antidiagonals.
%C For n>=2, the number of occurrences of n is 16n-14. For a guide to related sequences and permanents, see A204551.
%e Northwest corner:
%e 1 1 2 2 2 2 3 3 3
%e 1 2 2 2 2 3 3 3 3
%e 2 2 2 2 3 3 3 3 4
%e 2 2 2 3 3 3 3 4 4
%e 2 3 3 3 3 4 4 4 4
%e 3 3 3 3 4 4 4 4 5
%e 3 3 3 4 4 4 4 5 5
%t f[i_, j_] := Floor[(2 i + 2 j + 4)/4] - Floor[(i + j + 2)/4];
%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}]] (* A204560 *)
%t Permanent[m_] :=
%t With[{a = Array[x, Length[m]]},
%t Coefficient[Times @@ (m.a), Times @@ a]];
%t Table[Permanent[m[n]], {n, 1, 15}] (* A204561 *)
%Y Cf. A204561, A204551.
%K nonn,tabl
%O 1,4
%A _Clark Kimberling_, Jan 16 2012
|