|
| |
|
|
A000348
|
|
Number of ways to pair up {1^2, 2^2, ..., (2n)^2 } so sum of each pair is prime.
|
|
1
| |
|
|
1, 1, 2, 4, 12, 9, 72, 160, 428, 2434, 3011, 10337, 126962, 264182, 783550, 5004266, 34340141, 176302123, 1188146567, 4457147441, 7845512385, 132253267889, 1004345333251, 3865703506342, 40719018858150, 213982561376958, 1266218151414286, 10976172953868304, 59767467676582641, 512279001476451101, 6189067229056357433
(list; graph; refs; listen; history; internal format)
|
|
|
|
OFFSET
| 1,3
|
|
|
LINKS
| L. E. Greenfield and S. J. Greenfield, Some Problems of Combinatorial Number Theory Related to Bertrand's Postulate, J. Integer Sequences, 1998, #98.1.2.
|
|
|
FORMULA
| a(n)=permanent(m), where the n-by-n matrix m is defined m(i,j) = 1 or 0, depending on whether (2i)^2+(2j-1)^2 is prime or composite, respectively. - T. D. Noe (noe(AT)sspectra.com), Feb 10 2007
|
|
|
PROG
| (PARI) permRWNb(a)=n=matsize(a)[1]; if(n==1, return(a[1, 1])); sg=1; nc=0; 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; nc+=z; x+=z*a[, j]; p+=prod(i=1, n, x[i], sg)); return(2*(2*(n%2)-1)*p) for(n=1, 24, a=matrix(n, n, i, j, isprime((2*i)^2+(2*j-1)^2)); print1(permRWNb(a)", ")) - Herman Jamke (hermanjamke(AT)fastmail.fm), May 13 2007
|
|
|
CROSSREFS
| Cf. A000341.
Sequence in context: A137369 A186118 A181817 * A168342 A141668 A087796
Adjacent sequences: A000345 A000346 A000347 * A000349 A000350 A000351
|
|
|
KEYWORD
| nonn,nice
|
|
|
AUTHOR
| greenfie(AT)math.rutgers.edu (S. J. Greenfield)
|
|
|
EXTENSIONS
| a(11)-a(16) from David W. Wilson (davidwwilson(AT)comcast.net).
a(17)-a(22) from T. D. Noe (noe(AT)sspectra.com), Feb 10 2007
a(23)-a(24) from Herman Jamke (hermanjamke(AT)fastmail.fm), May 13 2007
More terms from Sean A. Irvine (sairvin(AT)xtra.co.nz), Nov 14 2010
|
| |
|
|
