cp's OEIS Frontend

This is a front-end for the Online Encyclopedia of Integer Sequences, made by Christian Perfect. The idea is to provide OEIS entries in non-ancient HTML, and then to think about how they're presented visually. The source code is on GitHub.

A058389 Number of 3 X 3 matrices with nonnegative integer entries and all row sums equal to n, up to row and column permutation.

Original entry on oeis.org

1, 3, 14, 44, 129, 316, 714, 1452, 2775, 4963, 8478, 13838, 21827, 33306, 49504, 71754, 101871, 141807, 194128, 261570, 347633, 456026, 591384, 758596, 963657, 1212861, 1513806, 1874440, 2304225, 2813030, 3412466, 4114608, 4933519
Offset: 0

Views

Author

Vladeta Jovovic, Nov 24 2000

Keywords

Crossrefs

Programs

  • Mathematica
    a[n_] := (m = Mod[n, 6]; (n^3 + 9*n^2 + 39*n + 120)*n^3 + Which[m == 0, 12*(23*n^2 + 32*n + 24), m == 1 || m == 5, 249*n^2 + 303*n + 143, m == 2 || m == 4, 4*(69*n^2 + 96*n + 56), m == 3, 3*(83*n^2 + 101*n + 69)])/288; Table[a[n], {n, 0, 32}] (* Jean-François Alcover, Oct 12 2011, after Vladeta Jovovic *)
  • PARI
    \\ See A318951 for RowSumMats
    a(n)=RowSumMats(3, 3, n); \\ Andrew Howroyd, Sep 05 2018

Formula

a(n) = (1/6)*(C(C(n + 2, 2) + 2, 3) + 3/2*floor((n + 2)/2)*(C(n + 2, 2) - floor((n + 2)/2)) + 3*C(floor((n + 2)/2) + 2, 3) + 2*floor(C(n + 2, 2)/3) + 2*C(C(n + 2, 2) - 3*floor(C(n + 2, 2)/3) + 2, 3)).
Empirical G.f.: -(x^8 + 3*x^7 + 14*x^6 + 12*x^5 + 15*x^4 + 9*x^3 + 5*x^2 + 1) / ((x-1)^7*(x+1)^3*(x^2+x+1)). - Colin Barker, Dec 27 2012

Extensions

More terms from Marc LeBrun, Dec 11 2000