OFFSET
3,1
LINKS
Simone Severini, www-users.york.ac.uk/~ss54
PROG
(PARI) permRWNb(a)=n=matsize(a)[1]; if(n==1, return(a[1, 1])); sg=1; in=vectorv(n); x=in; x=a[, n]-sum(j=1, n, a[, j])/2; p=prod(i=1, n, x[i]); for(k=1, 2^(n-1)-1, sg=-sg; j=valuation(k, 2)+1; z=1-2*in[j]; in[j]+=z; x+=z*a[, j]; p+=prod(i=1, n, x[i], sg)); return(2*(2*(n%2)-1)*p) for(n=1, 23, a=matrix(n, n, i, j, if(i==1||j==1||i==n||j==n, 1, -1+2*isprime((i-1)*n+j))); print1(permRWNb(a)", ")) - Herman Jamke (hermanjamke(AT)fastmail.fm), May 15 2007
CROSSREFS
KEYWORD
sign
AUTHOR
Simone Severini, Feb 15 2006
EXTENSIONS
More terms (and corrected definition) from Herman Jamke (hermanjamke(AT)fastmail.fm), May 15 2007
STATUS
approved