|
|
PROG
| (PARI) NormaQ(a, n) = { local(aT) ; aT=mattranspose(a) ; if( a*aT == aT*a, 1, 0) ; } combMat(no, n) = { local(a, noshif) ; a = matrix(n, n) ; noshif=no ; for(co=1, n, for(ro=1, n, if( (noshif %2)== 1, a[ro, co] = 1, a[ro, co] = -1) ; noshif = floor(noshif/2) ; ) ) ; return(a) ; } { for (n = 1, 10, count = 0; a = matrix(n, n) ; for( no=0, 2^(n^2)-1, a = combMat(no, n) ; count += NormaQ(a, n) ; if(no%1000==0, print(n, " ", (no/2^(n^2)+0.), " ", count)) ; ) ; print(count) ; ) } - R. J. Mathar (mathar(AT)strw.leidenuniv.nl), Mar 14 2006
|
