%I #28 Dec 12 2023 07:38:46
%S 0,0,0,8,16,8,24,24,24,32,40,32,48,48,48,56,64,56,72,72,72,80,88,80,
%T 96,96,96,104,112,104,120,120,120,128,136,128,144,144,144,152,160,152,
%U 168,168,168,176,184,176,192,192,192,200,208,200,216,216,216,224,232,224
%N Number of 3 X 3 reduced magic squares with distinct values and maximum value 2n; also, with magic sum 3n.
%C In a reduced magic square the row, column, and two diagonal sums must all be equal (the "magic sum") and the minimum entry is 0. The maximum entry is necessarily even and = (2/3)*(magic sum).
%C a(n) is a quasipolynomial with period 6.
%C The second differences of A108576 are a(n/2) for even n and 0 for odd n. The first differences of A108578 are a(n).
%H Thomas Zaslavsky, <a href="/A174256/b174256.txt">Table of n, a(n) for n = 1..10000</a>.
%H Matthias Beck and Thomas Zaslavsky, <a href="https://cs.uwaterloo.ca/journals/JIS/VOL13/Zaslavsky/sls.html">Six Little Squares and How Their Numbers Grow </a>, J. Int. Seq. 13 (2010), 10.6.2.
%H Matthias Beck and Thomas Zaslavsky, <a href="https://people.math.binghamton.edu/zaslav/Tpapers/SLSfiles/">"Six Little Squares and How their Numbers Grow" Web Site</a>: Maple worksheets and supporting documentation.
%H <a href="/index/Rec#order_05">Index entries for linear recurrences with constant coefficients</a>, signature (0,1,1,0,-1).
%F a(n) = 8*A174257(n).
%F G.f.: 8*x^4 * (2*x+1) / ((x^2-1) * (x^3-1)). [amended by _Georg Fischer_, Apr 17 2020]
%F a(n) = 2*(6*n - 13 - 8*cos(2*n*Pi/3) - 3*cos(n*Pi))/3. - _Wesley Ivan Hurt_, Oct 04 2018
%t Take[CoefficientList[Series[(8x^8 (2x^2+1))/((x^4-1)(x^6-1)),{x,0,120}],x],{1,-1,2}] (* _Harvey P. Dale_, Aug 07 2017 *)
%Y Cf. A108576, A108577, A174257.
%K nonn,easy
%O 1,4
%A _Thomas Zaslavsky_, Mar 14 2010