0,4

More precisely, consider 3 X 3 matrices with entries chosen from {0, 1, ..., n-1}, in which each row and column sums to n, where n >= 2. Then a(n) is the number of equivalence classes of such matrices under permutions of rows and columns and transpositions.

E. J. Morgan, Construction of Block Designs and Related Results. Ph.D. Dissertation, Univ. Queensland, 1978.

E. J. Morgan, On 3 X 3 matrices with constant row and column sum, Abstract 763-05-13, Notices Amer. Math. Soc., 26 (1979), page A-27.

M. F. Hasler, Table of n, a(n) for n=0,...,1000.

S. Plouffe, Approximations de S\'{e}ries G\'{e}n\'{e}ratrices et Quelques Conjectures}, Dissertation, Universit\'{e} du Qu\'{e}bec \`{a} Montr\'{e}al, 1992.

S. Plouffe, 1031 Generating Functions and Conjectures, Universit\'{e} du Qu\'{e}bec \`{a} Montr\'{e}al, 1992.

Let n = 3k, 3k-1 or 3k-2 according as n == 0, 2 or 1 mod 3, for n>=3.

Then a(n) = Sum_{ i=1,...,n-k } Sum_{ m=max(0,2i-n),...,floor(i/2) } Sum_{ r=0,...,floor(i/2)-m } c(i,m,r), where c(i,m,r) = n-2i+m+1 when m+r != i/2, or = floor((n-2i+m+2)/2) when m+r = i/2. [Typos corrected by Peter Pein, May 13 2008]

a(2) = 1:

110

101

011

a(3) = 3:

111 210 210

111 102 111

111 021 012

A005045:=-z**2*(-z**5+z**6-z**3+z+1)/((z**2+1)*(z**2+z+1)*(z+1)**2*(z-1)**5); [Conjectured by S. Plouffe in his 1992 dissertation. This matches at least 500 terms of the sequence, and is almost certainly correct, although no formal proof of its correctness has been given.]

Block[{k = Floor[(n + 2)/3]}, Sum[Sum[Sum[If[m + r == i/2, Floor[(n - 2*i + m + 2)/2], n - 2*i + m + 1], {r, 0, Floor[i/2 - m]}], {m, Max[2*i - n, 0], Floor[i/2]}], {i, 1, n - k}]]; Table[an, {n, 2, 100}] (from Peter Pein, May 13 2008)

(PARI from M. F. Hasler, Version 1, May 13 2008:)

A005045(n)={sum( i=1, n-(n+2)\3, sum( m=max(0, 2*i-n), i\2, sum( r=0, i\2-m, if( m+r!=i/2, n-2*i+m+1, (n-2*i+m+2)\2))))}

(PARI from M. F. Hasler, Version 2, much faster, May 13 2008:)

A005045(n)={sum( i=1, (2*n)\3, sum( m=max(0, 2*i-n), i\2, (n-2*i+m+1)*((i+1)\2-m)+(i%2==0)*(n-2*i+m+2)\2))}

Cf. A002817 for another version.

Sequence in context: A038505 A119971 A094272 this_sequence A069241 A092263 A076251

Adjacent sequences: A005042 A005043 A005044 this_sequence A005046 A005047 A005048

nonn,nice

njas

Edited by njas, May 12 2008, May 13 2008

More terms from Peter Pein, May 13 2008