|
PROGRAM
|
(PARI) { a(n) = M=matrix(n, n, i, j, 2*(i==j)); r=0; b(1); r } { b(k) = local(z, t); if(k>n, z=t=0; for(i=1, n, for(j=1, i-1, if(M[ i, j ]==0, z++); if(abs(M[ i, j ])==1, t++); )); r+=3^z*2^t; return; ); forvec(x=vector(k-1, i, [ -1, 1 ]), for(i=1, k-1, M[ k, i ]=M[ i, k ]=x[ i ]); if( matdet(vecextract(M, 2^k-1, 2^k-1), 1)>0, b(k+1) ) ) } (Alekseyev)
|
|
EXTENSIONS
|
79505 from Wouter Meeussen (wouter.meeussen(AT)pandora.be), Sep 05 2003
Two more terms from Max Alekseyev (maxal(AT)cs.ucsd.edu), Dec 13 2005
15477341239385927 from Max Alekseyev (maxal(AT)cs.ucsd.edu), Nov 09 2006
|