A061644 "Right perfect numbers": primes of the form 1 + a perfect number.
7, 29, 33550337, 137438691329
Offset: 1
Links
- Mersenne Forum, Thread 10336
- C. Rivera, Puzzle 203
Programs
-
Mathematica
pn={6, 28, 496, 8128, 33550336, 8589869056, 137438691328, 2305843008139952128, 2658455991569831744654692615953842176, 191561942608236107294793378084303638130997321548169216};lst={};Do[p=pn[[n]]+1;If[PrimeQ[p],AppendTo[lst,p]],{n,Length[pn]}];lst... and/or...PerfectNum[n_]:=Plus@@Divisors[n]/2;lst={};Do[p=PerfectNum[n];If[p==n&&PrimeQ[p+1],AppendTo[lst,p+1]],{n,10!}];lst (* Vladimir Joseph Stephan Orlovsky, Jan 27 2009 *)
-
PARI
forprime(p=2,100,P=2^p-1;Q=P*(P+1)/2+1;if(isprime(P)&&isprime(Q),print1(Q,","))) \\ Edward Jiang, Sep 10 2014
Formula
P(p)*[P(p)+1]/2 + 1 is prime, where P(p) is a Mersenne prime. (corrected by Lekraj Beedassy, May 01 2009)
Comments