| DATA
|
13, 29, 37, 53, 61, 73, 97, 101, 109, 127, 137, 149, 157, 173, 181, 193, 197, 229, 257, 269, 277, 293, 317, 337, 349, 373, 389, 397, 421, 457, 461, 487, 509, 541, 557, 569, 577, 613, 617, 641, 653, 661, 673, 677, 701, 709, 733, 757, 773, 797, 809, 811, 821
|
| PROG
|
(PARI) mx2pmp(n, m) = { v=vector(n+1); c=0; for(a=1, n, for(x=1, n, y = a*x^2+1; if(isprime(y), c++; v[c]=y; break); ); ); v1=vecsort(v); for(j=1, n, if(v1[j]<>0 && v1[j]==v1[j+1], print1(v1[j]", "))) }
|