mersenneforum.org Announcing a new Wagstaff PRP
 Register FAQ Search Today's Posts Mark Forums Read

 2021-06-28, 20:14 #1 ryanp     Jun 2012 Boulder, CO 5·61 Posts Announcing a new Wagstaff PRP (2^15135397+1)/3 is a Fermat Probable prime! (4556209 decimal digits) Also submitted to PRPTop. I am searching the range n=13M .. 17M currently, and nearly done. No other discoveries as of yet.
2021-06-28, 21:19   #2
paulunderwood

Sep 2002
Database er0rr

73418 Posts

Quote:
 Originally Posted by ryanp (2^15135397+1)/3 is a Fermat Probable prime! (4556209 decimal digits) Also submitted to PRPTop. I am searching the range n=13M .. 17M currently, and nearly done. No other discoveries as of yet.
Congrats

 2021-06-28, 23:38 #3 diep     Sep 2006 The Netherlands 32×83 Posts That's a very lucky find! Congrats on that one! Had you asked me i would've guessed next one might've lurked at 30M earliest and 70M latest.
2021-06-29, 11:29   #4
Dr Sardonicus

Feb 2017
Nowhere

22·3·401 Posts

Quote:
 Originally Posted by ryanp (2^15135397+1)/3 is a Fermat Probable prime! (4556209 decimal digits) Also submitted to PRPTop.

Congratulations! What base did you use for the PRP test? Three, perhaps?

2021-06-29, 13:24   #5
ryanp

Jun 2012
Boulder, CO

5·61 Posts

Quote:
 Originally Posted by Dr Sardonicus Congratulations! What base did you use for the PRP test? Three, perhaps?
I've run the Fermat PRP test with sllr64 using b=3, b=5, b=7 and b=11.

2021-06-29, 15:32   #6
Dr Sardonicus

Feb 2017
Nowhere

481210 Posts

Quote:
 Originally Posted by ryanp I've run the Fermat PRP test with sllr64 using b=3, b=5, b=7 and b=11.
Thanks!

Silly me, I failed to consider that you had tested multiple bases.

Of course, these numbers automatically "pass" the test to base 2. Paper and pencil suffices for this one.

If p > 3 is prime, N = (2^p + 1)/3, then (N-1)/2 = (2^(p-1) - 1)/3 is odd and divisible by p, so

N = (2^p + 1)/3 divides 2^p + 1, and 2^p + 1 divides 2^((N-1)/2) + 1, so N divides 2^((N-1)/2) + 1.

Now 2^((N-1)/2) + 1 divides 2^(N-1) - 1, so N divides 2^(N-1) -1, but does not divide 2^((N-1)/2) - 1.

2021-06-29, 16:59   #7
sweety439

"99(4^34019)99 palind"
Nov 2016
(P^81993)SZ base 36

52·112 Posts

Quote:
 Originally Posted by ryanp (2^15135397+1)/3 is a Fermat Probable prime! (4556209 decimal digits) Also submitted to PRPTop. I am searching the range n=13M .. 17M currently, and nearly done. No other discoveries as of yet.
Can you test all Wagstaff exponents below it? Currently only the Wagstaff exponents below 10 million are tested, see https://mersenneforum.org/showthread.php?t=24185

2021-06-29, 17:02   #8
sweety439

"99(4^34019)99 palind"
Nov 2016
(P^81993)SZ base 36

52·112 Posts

Quote:
 Originally Posted by ryanp I've run the Fermat PRP test with sllr64 using b=3, b=5, b=7 and b=11.
Only run Fermat test is dangerous, as there are many Carmichael numbers, see https://en.wikipedia.org/wiki/Bailli...Fermat%20tests, you should run either Miller-Rabin test or Baillie–PSW test.

Last fiddled with by sweety439 on 2021-06-29 at 17:02

2021-06-29, 20:23   #9
mathwiz

Mar 2019

C416 Posts

Quote:
 Originally Posted by sweety439 Only run Fermat test is dangerous, as there are many Carmichael numbers, see https://en.wikipedia.org/wiki/Bailli...Fermat%20tests, you should run either Miller-Rabin test or Baillie–PSW test.
As has been asked in numerous other threads: Is it that you lack the understanding, or the software/hardware, to do this yourself?

Maybe a mod can ban this guy until he stops flooding the forum with requests for other people to do things.

2021-06-29, 20:37   #10
paulunderwood

Sep 2002
Database er0rr

13·293 Posts

Quote:
 Originally Posted by sweety439 Only run Fermat test is dangerous, as there are many Carmichael numbers, see https://en.wikipedia.org/wiki/Bailli...Fermat%20tests, you should run either Miller-Rabin test or Baillie–PSW test.
The instructions on how to run sllr with switches to do a lucas test and a bpsw test are given in this post. If Ryan nor someone else does not step up to the plate in the meantime, I'll do it at the weekend.

Last fiddled with by paulunderwood on 2021-06-29 at 20:38

2021-06-29, 21:33   #11
VBCurtis

"Curtis"
Feb 2005
Riverside, CA

5·23·43 Posts

Quote:
 Originally Posted by mathwiz Maybe a mod can ban this guy until he stops flooding the forum with requests for other people to do things.
We gave him some time off to consider his behavior, and it hasn't changed much. I suppose your suggestion and this reply might be considered yet another warning to Sweety before the banhammer falls again.

 Similar Threads Thread Thread Starter Forum Replies Last Post GP2 Wagstaff PRP Search 414 2020-12-27 08:11 eight6225219 Miscellaneous Math 7 2019-11-21 14:59 chalsall Data 5 2012-12-31 15:57 chalsall Data 25 2011-01-25 07:00 zariss Linux 1 2004-08-23 21:42

All times are UTC. The time now is 19:44.

Sun Sep 19 19:44:03 UTC 2021 up 58 days, 14:13, 0 users, load averages: 1.15, 1.16, 1.16