Magic squares of doubly even order: Difference between revisions

m
Line 329:
 
A magic one-liner:
<lang parigp>magicsquare(n)={matrix(n,n,i,j,k=i+j*n-n;if(bitand(38505,2^((j-1)%4*4+(i-1)%4)),k,n*n+1-k))}</lang>
 
Output:<pre>magicsquare(8)
Anonymous user