OFFSET
0,3
LINKS
FORMULA
a(n) = 780*a(n-2) -194881*a(n-4) +22377420*a(n-6)
-1419219792*a(n-8) +55284715980*a(n-10)
-1410775106597*a(n-12) +24574215822780*a(n-14)
-300429297446885*a(n-16) +2629946465331120*a(n-18)
-16741727755133760*a(n-20) +78475174345180080*a(n-22)
-273689714665707178*a(n-24) +716370537293731320*a(n-26)
-1417056251105102122*a(n-28) +2129255507292156360*a(n-30)
-2437932520099475424*a(n-32) +2129255507292156360*a(n-34)
-1417056251105102122*a(n-36) +716370537293731320*a(n-38)
-273689714665707178*a(n-40) +78475174345180080*a(n-42)
-16741727755133760*a(n-44) +2629946465331120*a(n-46)
-300429297446885*a(n-48) +24574215822780*a(n-50)
-1410775106597*a(n-52) +55284715980*a(n-54)
-1419219792*a(n-56) +22377420*a(n-58)
-194881*a(n-60) +780*a(n-62) -a(n-64).
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[11, n];
CROSSREFS
KEYWORD
nonn,easy
AUTHOR
Alois P. Heinz, Mar 30 2012
STATUS
approved