How to generate base10 representation of Mersenne-prime perfect numbers?

James Heinrich · 2012-03-03, 02:43

I found this site which has base10 representations of Mersenne-prime perfect numbers, but it's lacking the ones for M45-M47. How can I generate these?
Something that runs on Windows would be great, but if I do it with PARI I also have access to that on *nix.

Batalov · 2012-03-03, 02:56

Pari

Code:

allocatemem(800000000)
P=37156667;  
write("perf45",2^(P-1)*(2^P-1))
 
P=42643801;
write("perf46",2^(P-1)*(2^P-1))
 
P=43112609;
write("perf47",2^(P-1)*(2^P-1))
 
# in *nix
fold < perf47 > perf47.txt

kar_bon · 2012-03-03, 08:58

With pfgw:

Code:

pfgw -od -q"2^43112608*(2^43112609-1)" >M47.txt

James Heinrich · 2012-03-03, 13:17

Quote:

Originally Posted by kar_bon

With pfgw:

Code:

pfgw -od -q"2^43112608*(2^43112609-1)" >M47.txt

Perfect, thanks! Worked nicely.

I tried PARI, and failed. On *nix, all I could get was a Segmentation Fault. On Windows, it ran for 15 minutes and didn't seem to be doing anything so I aborted (by contrast, pfgw only took about 45 seconds).

science_man_88 · 2012-03-03, 13:28

Quote:

Originally Posted by James Heinrich

Perfect, thanks! Worked nicely.

I tried PARI, and failed. On *nix, all I could get was a Segmentation Fault. On Windows, it ran for 15 minutes and didn't seem to be doing anything so I aborted (by contrast, pfgw only took about 45 seconds).

I got a working script:

Code:

(09:26)>MeVec
%240 = [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, 1398269, 2976221, 3021377, 6972593, 13466917]
(09:26)>for(x=1,#MeVec, print((2^(MeVec[x]-1))*(2^MeVec[x]-1)))

just change MeVec to include the exponents you want.

Dubslow · 2012-03-03, 20:11

Maybe perhaps wait until we know whether or not we know M48?

Xitami · 2012-03-08, 00:20

http://bellard.org/mersenne.html

LaurV · 2012-03-08, 02:42

Quote:

Originally Posted by Xitami

http://bellard.org/mersenne.html

I love it! So, you are one of the IOCCC winners? Respect!

Xitami · 2012-03-08, 06:45

NO

LaurV · 2012-03-08, 07:05

Quote:

Originally Posted by Xitami

NO

Oh. Sorry. I don't know why I understood that is your web page the one you linked, maybe because of the first person speech there. Mea culpa. Anyhow, welcome to the forum.

Dubslow · 2012-03-08, 07:20

Well what makes you think it's his website? Also, he's been here longer than I have.

2012-03-03, 02:43	#1
James Heinrich "James Heinrich" May 2004 ex-Northern Ontario 2³·5·83 Posts	How to generate base10 representation of Mersenne-prime perfect numbers? I found this site which has base10 representations of Mersenne-prime perfect numbers, but it's lacking the ones for M45-M47. How can I generate these? Something that runs on Windows would be great, but if I do it with PARI I also have access to that on *nix.

2012-03-03, 02:56	#2
Batalov "Serge" Mar 2008 Phi(4,2^7658614+1)/2 10010010101100₂ Posts	Pari Code: allocatemem(800000000) P=37156667; write("perf45",2^(P-1)(2^P-1)) P=42643801; write("perf46",2^(P-1)(2^P-1)) P=43112609; write("perf47",2^(P-1)(2^P-1)) # in nix fold < perf47 > perf47.txt

2012-03-03, 08:58	#3
kar_bon Mar 2006 Germany 43·67 Posts	With pfgw: Code: pfgw -od -q"2^43112608*(2^43112609-1)" >M47.txt

Similar Threads
Thread	Thread Starter	Forum	Replies	Last Post
Right Perfect Prime Numbers	Housemouse	Math	34	2016-04-07 16:29
Have Found Principle to generate infinitive PRIME NUMBERS	Evgeny Dolgov	Miscellaneous Math	38	2010-09-05 17:45
Square numbers and binary representation	ET_	Miscellaneous Math	40	2010-06-06 12:55
Effective way to generate prime numbers (infinitive)	Evgeny Dolgov	Math	1	2003-12-08 09:25
Search for Mersenne primes by checking for perfect numbers	dsouza123	Miscellaneous Math	33	2003-09-02 16:18

2012-03-03, 20:11	#6
Dubslow Basketry That Evening! "Bunslow the Bold" Jun 2011 40<A<43 -89<O<-88 3×29×83 Posts	Maybe perhaps wait until we know whether or not we know M48?

2012-03-08, 00:20	#7
Xitami Apr 2010 2×7 Posts	http://bellard.org/mersenne.html

2012-03-08, 07:20	#11
Dubslow Basketry That Evening! "Bunslow the Bold" Jun 2011 40<A<43 -89<O<-88 3×29×83 Posts	Well what makes you think it's his website? Also, he's been here longer than I have.