OFFSET
1,3
COMMENTS
The old entry with this sequence number was a duplicate of A048298.
PROG
(PARI) { 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 }
/* Max Alekseyev, May 14 2008 */
CROSSREFS
KEYWORD
nonn
AUTHOR
Jack Brennen, May 13 2008
EXTENSIONS
a(6)-a(8) from Max Alekseyev, May 14 2008
Edited by Charles R Greathouse IV, Oct 28 2009
a(9) from Sean A. Irvine, Oct 29 2022
STATUS
approved