OFFSET
0,3
LINKS
FORMULA
G.f.: -(x^30+x^29-154*x^28 +6777*x^26-1440*x^25-123961*x^24 +26752*x^23 +1132714*x^22-185889*x^21 -5684515*x^20+574750*x^19+16401668*x^18 -708928*x^17 -27757938*x^16+27757938*x^14+708928*x^13 -16401668*x^12 -574750*x^11+5684515*x^10 +185889*x^9-1132714*x^8-26752*x^7 +123961*x^6 +1440*x^5-6777*x^4+154*x^2-x-1) / (x^32-209*x^30+11936*x^28 -274208*x^26 +3112032*x^24-19456019*x^22 +70651107*x^20-152325888*x^18+196664896*x^16 -152325888*x^14+70651107*x^12 -19456019*x^10 +3112032*x^8-274208*x^6 +11936*x^4-209*x^2+1).
MATHEMATICA
A[1, 1] = 1;
A[m_, n_] := A[m, n] = Module[{i, j, s, t, M}, Which[m == 0 || n == 0, 1, m < n, A[n, m], True, s = Mod[n*m, 2]; M[i_, j_] /; j < i := -M[j, i]; M[_, _] = 0; For[i = 1, i <= n, i++, For[j = 1, j <= m, j++, t = (i - 1)*m + j - s; If[i > 1 || j > 1 || s == 0, If[j < m, M[t, t + 1] = 1]; If[i < n, M[t, t + m] = 1 - 2*Mod[j, 2]]]]]; Sqrt[Det[Array[M, {n*m - s, n*m - s}]] ]]];
a[n_] := A[9, n];
CROSSREFS
KEYWORD
nonn,easy
AUTHOR
Alois P. Heinz, Apr 15 2011
STATUS
approved