login
The OEIS is supported by the many generous donors to the OEIS Foundation.

 

Logo
Hints
(Greetings from The On-Line Encyclopedia of Integer Sequences!)
A172675 Number of 3*n X 3 0..2 arrays with row sums 2 and column sums 2*n. 1
21, 2385, 352128, 57775905, 10060071021, 1820016119376, 338183208699840, 64089909936535329, 12331175198408791725, 2401214665364782652385, 472159936393091112404160, 93594776429965445731933200 (list; graph; refs; listen; history; text; internal format)
OFFSET
1,1
COMMENTS
Apparently the sequence gives the even-order terms in the diagonal of the rational function R(x,y,z) = 1/(1 - (x^2 + y^2 + z^2 - x*y - y*z + x*z)), i.e., a(n) = [(x*y*z)^(2*n)] R(x,y,z), n >= 1. - Gheorghe Coserea, Aug 09 2018
This is because those even-order terms would be the same if y were replaced by -y, giving S(x,y,z) = 1/(1 - (x^2 + y^2 + z^2 + x*y + y*z + x*z)) = Sum_{i>=0} (x^2 + y^2 + z^2 + x*y + y*z + x*z)^i, and the terms for (x*y*z)^(2*k) come from i=3*k in this sum. - Robert Israel, Jan 15 2023
LINKS
MAPLE
f:= proc(n) coeftayl((x^2+y^2+z^2+x*y+y*z+x*z)^(3*n), [x, y, z]=[0, 0, 0], [2*n, 2*n, 2*n]) end proc:
map(f, [$1..30]); # Robert Israel, Jan 15 2023
MATHEMATICA
a[n_] := SeriesCoefficient[(x^2 + y^2 + z^2 + x*y + y*z + x*z)^(3n), {x, 0, 2n}, {y, 0, 2n}, {z, 0, 2n}];
Table[a[n], {n, 1, 12}] (* Jean-François Alcover, May 17 2023, after Robert Israel *)
CROSSREFS
Sequence in context: A033510 A123844 A172622 * A068254 A144853 A131314
KEYWORD
nonn
AUTHOR
R. H. Hardin, Feb 06 2010
STATUS
approved

Lookup | Welcome | Wiki | Register | Music | Plot 2 | Demos | Index | Browse | More | WebCam
Contribute new seq. or comment | Format | Style Sheet | Transforms | Superseeker | Recents
The OEIS Community | Maintained by The OEIS Foundation Inc.

License Agreements, Terms of Use, Privacy Policy. .

Last modified July 15 16:36 EDT 2024. Contains 374333 sequences. (Running on oeis4.)