|
PROG
|
(PARI)
{e=[2, 3, 5, 7, 13, 17, 19, 31, 61, 89, 107, 127, 521, 607, 1279, 2203, 2281, 3217, 4253, 4423, 9689, 9941, 11213, 19937, 21701, 23209, 44497, 86243, 110503, 132049, 216091, 756839, 859433, 1257787]; } /* exponents of Mersenne primes */
for(n=1, #e, p=(2^e[n]-1)*(2^(e[n]-1)); if(ispseudoprime(p+1), print1(p, ", ")));
|