mersenneforum.org Announcing a new Wagstaff PRP
 User Name Remember Me? Password
 Register FAQ Search Today's Posts Mark Forums Read

 2021-06-28, 20:14 #1 ryanp     Jun 2012 Boulder, CO 17516 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

4,129 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 2·17·23 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

53×109 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

373 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

53·109 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

2×1,693 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

338610 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

2·131 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

4,129 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

122378 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 14:38.

Sun May 22 14:38:00 UTC 2022 up 38 days, 12:39, 0 users, load averages: 1.57, 1.04, 1.08

Copyright ©2000 - 2022, Jelsoft Enterprises Ltd.

This forum has received and complied with 0 (zero) government requests for information.

Permission is granted to copy, distribute and/or modify this document under the terms of the GNU Free Documentation License, Version 1.2 or any later version published by the Free Software Foundation.
A copy of the license is included in the FAQ.

≠ ± ∓ ÷ × · − √ ‰ ⊗ ⊕ ⊖ ⊘ ⊙ ≤ ≥ ≦ ≧ ≨ ≩ ≺ ≻ ≼ ≽ ⊏ ⊐ ⊑ ⊒ ² ³ °
∠ ∟ ° ≅ ~ ‖ ⟂ ⫛
≡ ≜ ≈ ∝ ∞ ≪ ≫ ⌊⌋ ⌈⌉ ∘ ∏ ∐ ∑ ∧ ∨ ∩ ∪ ⨀ ⊕ ⊗ 𝖕 𝖖 𝖗 ⊲ ⊳
∅ ∖ ∁ ↦ ↣ ∩ ∪ ⊆ ⊂ ⊄ ⊊ ⊇ ⊃ ⊅ ⊋ ⊖ ∈ ∉ ∋ ∌ ℕ ℤ ℚ ℝ ℂ ℵ ℶ ℷ ℸ 𝓟
¬ ∨ ∧ ⊕ → ← ⇒ ⇐ ⇔ ∀ ∃ ∄ ∴ ∵ ⊤ ⊥ ⊢ ⊨ ⫤ ⊣ … ⋯ ⋮ ⋰ ⋱
∫ ∬ ∭ ∮ ∯ ∰ ∇ ∆ δ ∂ ℱ ℒ ℓ
𝛢𝛼 𝛣𝛽 𝛤𝛾 𝛥𝛿 𝛦𝜀𝜖 𝛧𝜁 𝛨𝜂 𝛩𝜃𝜗 𝛪𝜄 𝛫𝜅 𝛬𝜆 𝛭𝜇 𝛮𝜈 𝛯𝜉 𝛰𝜊 𝛱𝜋 𝛲𝜌 𝛴𝜎𝜍 𝛵𝜏 𝛶𝜐 𝛷𝜙𝜑 𝛸𝜒 𝛹𝜓 𝛺𝜔