mersenneforum.org

mersenneforum.org (https://www.mersenneforum.org/index.php)
-   Information & Answers (https://www.mersenneforum.org/forumdisplay.php?f=38)
-   -   Very (large) PRPs? (https://www.mersenneforum.org/showthread.php?t=21294)

PawnProver44 2016-05-11 06:37

Very (large) PRPs?
 
I was reading this [URL="http://www.mersenne.ca/prp.php?show=1&min_exponent=1009&max_exponent=6972593"]list[/URL], when I found out it wasn't too hard to find all those PRPs. I looked at the end of the list, and found that "random" PRPs are discovered with 50,000 digits. What should I best use to find "random" PRPs say between 50k and 200k digits? Not that I would want to ever submit them on the [URL="http://www.primenumbers.net/prptop/prptop.php"]PRPtop[/URL] because the site simply does not accept "random" PRPs and there are no such PRPs there either. I would expect this to be done in about a week too, (200k digits almost), but correct me if that is not possible, please.

LaurV 2016-05-11 06:53

The "list" [U]will[/U] accept "random" PRPs. Find one with more than 46k digits and send it there. Just find one.

axn 2016-05-11 08:58

The [URL="http://www.primenumbers.net/prptop/submit.php"]submission page[/URL] lists 20,000 digit as the minimum. But to actually enter Top 10000, it should be 42k+.

The submission page also says that
[QUOTE]please type one number by line using the following format :
number,its number of digits in base 10
and don't type their decimal expansion[/QUOTE]

So it is possible that if your number doesn't have a compact form, it might not be accepted.

PawnProver44 2016-05-11 16:52

Can you find any examples of "random" PRPs on the site? I was looking for a specific program that I could input the large number in a .txt file (about 200,000 digits) and call nextprime(x) or previousprime(x), Pfgw is too slow is my guess?

Batalov 2016-05-11 18:45

[QUOTE=PawnProver44;433640]...and call nextprime(x) or previousprime(x), Pfgw is too slow is my guess?[/QUOTE]
PRPtop is full of them as it is:
[URL="http://www.primenumbers.net/prptop/searchform.php?form=10^x%2Bn&action=Search"]nextprime(10^x)[/URL] (and next and next and next)
[URL="http://www.primenumbers.net/prptop/searchform.php?form=12%5Ex%2Bn&action=Search"]nextprime(12^x)[/URL]
[URL="http://www.primenumbers.net/prptop/searchform.php?form=%282^67588%2Ba%29%3F&action=Search"]some other random example[/URL]
Finding one is very easy, if that's your goal or if you find it fun thing to do. Most contributors there don't stop until they submit tens of thousands.

PawnProver44 2016-05-11 23:17

First, with what text editor can carry 200k+digits (need at least 400k digits) so I can input into pfgw (or some other program which can give me a direct result regardless the size for nextprime(x) = ?):

input.txt:
------------------------------------------
nextprime(597612129456797651389670218204770656.........(199,928 more digits).........107345799455301271768055018207593331)-597612129456797651389670218204770656.........(199,928 more digits).........107345799455301271768055018207593331
------------------------------------------

And get a result like this

-f -i input.txt

nextprime(x)-x is a 3-PRP!

-t -q input.txt

nextprime(597612129456797651389670218204770656.........(199,928 more digits).........107345799455301271768055018207593331)

is a Fermat and Lucus PRP!

I don't want to risk notepad, but someone who has tried that please let me know.

Pfgw output I do nextprime(x)-x, but when i input nextprime(x) it gives me:

nextprime(x) is a 3-PRP!

So I do not know what nextprime(x) is, but I can ignore that using nextprime(x)-x.

Batalov 2016-05-11 23:26

[QUOTE=PawnProver44;433691]So I do not know what nextprime(x) is, but I can ignore that.[/QUOTE]
Dude, I can tell you that nextprime(10^1000000000) is a billion-digit prime (even without calculating it; I can ignore that, too).

Where can I get my EFF money!?? :max:

alpertron 2016-05-12 11:45

[QUOTE=Batalov;433694]Dude, I can tell you that nextprime(10^1000000000) is a billion-digit prime (even without calculating it; I can ignore that, too).

Where can I get my EFF money!?? :max:[/QUOTE]

prevprime(10^1000000000) is a billion-digit prime but nextprime(10^1000000000) is not.

Batalov 2016-05-12 16:53

Good point! I lost a "+" in the message: it is a billion[COLOR=DarkRed]+[/COLOR] digit number and qualifies for the EFF prize [SPOILER]not[/SPOILER]!

PawnProver44 2016-05-12 20:53

Here are my running times with Wolfram | Alpha? (random 10, 20, 40, 80 digit numbers) Can anyone compare to Pfgw to see if it is faster?

nextprime(1615940397) = 1615940399 (1.2 s)

nextprime(14068924638939642180) = 14068924638939642181 (1.3 s)

nextprime(3718587170008038384873784221963098409550) = 3718587170008038384873784221963098409563 (1.3 s)

nextprime(97783010505674166402735696905325823805408510851148150244960782744188145048695585) = 97783010505674166402735696905325823805408510851148150244960782744188145048695591 (1.3 s)

I was expecting the time to be at lest 16 times as much for doubling the number of digits? Is this true for Pfgw? Thanks to anyone who can give approximate running times for Pfgw for (1,000; 2,000; 4,000; 8,000; digits PRPs please) as a larger version of my test.

Edit: I forgot to call nextprime(x) for random 80 digit number.

alpertron 2016-05-12 21:05

[QUOTE=PawnProver44;433770]Here are my running times with Wolfram | Alpha? (random 10, 20, 40 digit numbers) Can anyone compare to Pfgw to see if it is faster?

nextprime(1615940397) = 1615940399 (1.2 s)

nextprime(14068924638939642180) = 14068924638939642181 (1.3 s)

nextprime(3718587170008038384873784221963098409550) = 3718587170008038384873784221963098409563 (1.3 s)

nextprime(97783010505674166402735696905325823805408510851148150244960782744188145048695585) = 97783010505674166402735696905325823805408510851148150244960782744188145048695591 (1.3 s)

[/QUOTE]

From those timings, it is clear that the real time for nextprime is 0.0 s for all instances, and 1.2 s - 1.3 s is lost in the Web service. It can even include an artificial delay where no computation is done, so you cannot send too many requests to the server.


All times are UTC. The time now is 14:39.

Powered by vBulletin® Version 3.8.11
Copyright ©2000 - 2021, Jelsoft Enterprises Ltd.