|
|
PROG
| (PARI) /* from Max Alekseyev */
{ npairs(n) = loca(r, q, z); r=0;
forvec(p=vector(n, i, [1, 2*n]),
q = eval( setminus( Set(vector(2*n, i, i)), Set(p) ) );
for(j=1, n!,
z=numtoperm(n, j);
if(type( sum(j=1, #p, p[j]/q[z[j]]) )=="t_INT", r++); ); , 2); r }
|