mersenneforum.org  

Go Back   mersenneforum.org > Prime Search Projects > Wagstaff PRP Search

Reply
 
Thread Tools
Old 2010-02-21, 07:16   #12
10metreh
 
10metreh's Avatar
 
Nov 2008

2×33×43 Posts
Default

Congratulations!

Quote:
Originally Posted by CRGreathouse View Post
My apologies for being 'out of the loop'. But what is known about these so-called Vrba-Reix prps? All Wagstaff primes are VRprps, but is it known that composite VRprps are asymptotically rare?
I think we need to coin a new term: vrobable prime
10metreh is offline   Reply With Quote
Old 2010-02-21, 20:24   #13
T.Rex
 
T.Rex's Avatar
 
Feb 2004
France

38F16 Posts
Default

Quote:
Originally Posted by CRGreathouse View Post
Congrats, Tony!
Thanks !
Quote:
My apologies for being 'out of the loop'. But what is known about these so-called Vrba-Reix prps ?
Look at these explainations.
Quote:
All Wagstaff primes are VRprps, but is it known that composite VR prps are asymptotically rare?
We have no idea. However, based on my attempt to prove my conjecture for Mersennes, I think that some more conditions are required for a Wagstaff PRP to be a Prime: I think that the product of all S(n) (in: S(n)=S(n-1)^2-2) must be 1 (modulo Wq).
However, I guess that the probability for a Wagstaff composite to have the VR property is VERY VERY small, if non null. Since VR is based on a Cycle that is of length q (Wq=(2^q+1)/3), I think there is no way that a composite can have VR true. But, without a proof... it is only guesses. Any idea for building a proof of VR ?

Read the papers I've gathered on my Blog. It explains many things. If there is something you do not understand, just ask me.

Regards,

Tony
T.Rex is offline   Reply With Quote
Old 2010-02-21, 20:39   #14
ixfd64
Bemusing Prompter
 
ixfd64's Avatar
 
"Danny"
Dec 2002
California

22·3·191 Posts
Default

This is very impressive. I hope it gets proven soon!
ixfd64 is offline   Reply With Quote
Old 2010-02-23, 08:22   #15
ATH
Einyen
 
ATH's Avatar
 
Dec 2003
Denmark

1011010011002 Posts
Default

Congratulations on the new PRP!

Prime95 PRP-test concur :)

UID: athath/laptop, 2^4031399+1/3 is a probable prime! Wd4: 2044F973,00000000
ATH is offline   Reply With Quote
Old 2010-02-23, 08:58   #16
Cybertronic
 
Cybertronic's Avatar
 
Jan 2007
DEUTSCHLAND !

26710 Posts
Default

Congratulation for this nice number !!!!

10 years ago, I tried the Lucas Lehmer Test for this type of numbers.

Is it a prime, so the result S=14 or S=194 ( S=2 is for Mp's ), but it is maybe not a proof.
Cybertronic is offline   Reply With Quote
Old 2010-02-23, 13:58   #17
T.Rex
 
T.Rex's Avatar
 
Feb 2004
France

911 Posts
Default

Quote:
Originally Posted by Cybertronic View Post
Congratulation for this nice number !!!!
Thanks !
Quote:
10 years ago, I tried the Lucas Lehmer Test for this type of numbers. Is it a prime, so the result S=14 or S=194 ( S=2 is for Mp's ), but it is maybe not a proof.
Hummm LLT for Mersenne makes use of the Tree of the DiGraph under x^2-2 . With Wagstaff numbers, there is no tree, only cycles. So LLT (meaning: ending with a 0, and then -2 and then 2) cannot work for Wagstaff.
Which seed did you use ? since you only talk of the last step.

Tony
T.Rex is offline   Reply With Quote
Old 2010-02-23, 13:59   #18
T.Rex
 
T.Rex's Avatar
 
Feb 2004
France

911 Posts
Default

Quote:
Originally Posted by ATH View Post
Congratulations on the new PRP!
Thanks !
Quote:
Prime95 PRP-test concur :)
UID: athath/laptop, 2^4031399+1/3 is a probable prime! Wd4: 2044F973,00000000
Thanks for double-checking !!!!

Tony
T.Rex is offline   Reply With Quote
Old 2010-02-23, 14:22   #19
Cybertronic
 
Cybertronic's Avatar
 
Jan 2007
DEUTSCHLAND !

3×89 Posts
Default

Quote:
Originally Posted by T.Rex View Post
Thanks !
Hummm LLT for Mersenne makes use of the Tree of the DiGraph under x^2-2 . With Wagstaff numbers, there is no tree, only cycles. So LLT (meaning: ending with a 0, and then -2 and then 2) cannot work for Wagstaff.
Which seed did you use ? since you only talk of the last step.

Tony
n=2^p+1 , not (2^p+1)/3
You start with S=14, after the LL procedure (S=S^2 mod n)-2 you get
S=14 or S=194 (14^2-2).

Here the UBASIC program:
10 P=3
20 N=2^P+1
30 T=1:S=14
40 S=(S^2 mod N)-2
50 inc T
60 if T=P then print P;S:P=nxtprm(P):goto 20
65 if P>100 then stop
70 goto 40


Here the plot:
3 5 -> 5=14-9
5 14
7 65 -> 65=194-129
11 14
13 194
17 14
19 194
23 14
29 76653758
31 194
37 42511304261
41 359734364441
43 194
47 17518790332040
53 6406117628425502
59 508334340273339569
61 194
67 21819313749230765027
71 1195598766420873147536
73 5630004848125169115977
79 194
83 5637954672539490707970155
89 20493566371156403602606220
97 68242359663609183220101332858
101 14
103 435037204291375683415946537180
107 32093918361200282030236449354443
109 414278313959563741717336012443446
113 5112917025335077160694569180088809
127 194
131 1973968232982423701206473030341408459198
137 47484206655905551903295319361957840172753
139 168667007687552910631152078108586066160726
149 18453965156596083542114781548941146664272227
151 262976984875073064821485878363989765425623023
157 49547112554957562135729742000244759963364760238
163 898773787877688180281701171238174776273408405979
167 14
173 6874700042006001666774484698531768324621714679913621
179 297595739783471837689716178237652964409196604508541930
181 1244797157598192014982163301734406121591070982512878662
191 14
193 246557495726160777169800308963987254019574345500178382517
197 154344194312622000819873725710721046832719797408900851567781
199 194

next p's are
313 194
347 14
701 14


The same like #7 List of Wagstaff Primes:
(2) 3 5 7 11 13 17 19 23 31 43 61 79 101 127 167 191 199 313 347 701

best

Last fiddled with by Cybertronic on 2010-02-23 at 14:31
Cybertronic is offline   Reply With Quote
Old 2010-02-23, 19:26   #20
Cybertronic
 
Cybertronic's Avatar
 
Jan 2007
DEUTSCHLAND !

3·89 Posts
Default

short addition :

I get with my program:

OK
list
10 P=3
20 N=2^P+1
30 T=1:S=14
40 S=(S^2 mod N)-2
50 inc T
60 if T=P then goto 100
70 goto 40
100 if or{S=(14 mod N),S=(194 mod N)} then print P;
120 P=nxtprm(P):goto 20
OK
run
3 5 7 11 13 17 19 23 31 43 61 79 101 127 167 191 199 313 347 701 1709 2617 3539 5807

exact all wagstaff prime exponents in logical order...

Last fiddled with by Cybertronic on 2010-02-23 at 19:30
Cybertronic is offline   Reply With Quote
Old 2010-02-24, 15:46   #21
Jeff Gilchrist
 
Jeff Gilchrist's Avatar
 
Jun 2003
Ottawa, Canada

7·167 Posts
Default

For fun I ran this exponent through my own GMP based Vrba-Reix PRP code on an Opteron and Core2 to verify and it comes out as prime:

W4031399 is prime! W4031399 = ...2388889642104198122875563 (1213572 digits)
Finished! Checked prime q values to: 4031399

Great work DUR!

Jeff.
Jeff Gilchrist is offline   Reply With Quote
Old 2010-02-24, 21:35   #22
T.Rex
 
T.Rex's Avatar
 
Feb 2004
France

911 Posts
Default

Quote:
Originally Posted by Jeff Gilchrist View Post
.... and it comes out as prime:
As a PRP. Sure it is a prime... but we lack a proof. So it is a PRP.
Thanks for the nice job, Jeff !!
Tony
T.Rex is offline   Reply With Quote
Reply

Thread Tools


Similar Threads
Thread Thread Starter Forum Replies Last Post
New Wagstaff PRP exponents ryanp Wagstaff PRP Search 26 2013-10-18 01:33
Hot tuna! -- a p75 and a p79 by Sam Wagstaff! Batalov GMP-ECM 9 2012-08-24 10:26
Wagstaff Conjecture davieddy Miscellaneous Math 209 2011-01-23 23:50
Best settings to factor Wagstaff p = (2^n +1) / 1 diep GMP-ECM 10 2010-07-26 21:33
30th Wagstaff prime T.Rex Math 0 2007-09-04 07:10

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

Thu Aug 13 14:59:13 UTC 2020 up 11:34, 1 user, load averages: 1.11, 1.22, 1.31

Powered by vBulletin® Version 3.8.11
Copyright ©2000 - 2020, 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.