mersenneforum.org  

Go Back   mersenneforum.org > Prime Search Projects > Five or Bust - The Dual Sierpinski Problem

Reply
 
Thread Tools
Old 2009-09-04, 22:59   #23
Zuzu
 
Sep 2009

11 Posts
Default

Quote:
Originally Posted by philmoore View Post
Well, Zuzu, I hope that you're right and I'm wrong, as your figures certainly look more optimistic! I will take another look at it this weekend, when I'm not doing school work. That said, I admit that statistics is not my strong suit, although I felt like I had a good "seat of the pants" feel for it when I was doing graduate work in physics. I estimated the 1.32153 expected number as follows: each candidate Nhas probability of being prime equal to {e^\gamma\ln F}\over{\ln N}, where F is the factoring depth. There are approximately 328,300 candidates left in the sieve for all exponents n between 0 and 50 million, so I assumed that the candidates are equally distributed and replaced the sum by an integral for n from 4.65 million to 50 million. Next, I assumed that since the probabilities for each candidate are independent, that if p is the probability of at least one candidate being prime, then p^2 is the probability of at least two, p^3 is the probability of at least three, and so forth. Summing these probabilities gives the expected number of primes as p\over{1-p}, so if we set this equal to the expected number 1.32, we get p as about 1.32/2.32, or about 56.9%. Note that this gives a 32.3% chance of at least two primes if we search all the candidates, but note that there is only about a 50% chance that the two primes will then belong to different k values. I also translated this into the language of Poisson processes, where e^{-w\log_2(n_f/n_i)} was the probability of finding no primes for all exponents between n_i and n_f. My weights for w were .17700 for 41693 and .17566 for 40291. (Divide these weights by 2 and you will see that they are a little smaller than those on Payam Samidoost's webpage, but those weights were computed by Brennen's Proth weight applet with a fairly low cutoff.)

I will look up my notes from the last discovery and see if I can resolve this. Any comments are welcome, as I may be making some wrong assumptions here!
Phil,

Thank you for having provided the relevant information: the number of surviving candidates for a given factoring depth (sieving depth ?). Then the probability of an individual candidate being prime is approximately independent of the Proth weight. I'll make the calculations and confirm or infirm my previous considerations...

Last fiddled with by Zuzu on 2009-09-04 at 23:01 Reason: misleading formulation, needing precision
Zuzu is offline   Reply With Quote
Old 2009-09-05, 18:13   #24
Jeff Gilchrist
 
Jeff Gilchrist's Avatar
 
Jun 2003
Ottawa, Canada

7×167 Posts
Default

Quote:
Originally Posted by paleseptember View Post
As suggested, I went and bought a lottery ticket. Strangely enough, no joy.
Did you make sure all the numbers you picked on the ticket were PRP as well? If not...

Jeff.
Jeff Gilchrist is offline   Reply With Quote
Old 2009-09-06, 00:44   #25
philmoore
 
philmoore's Avatar
 
"Phil"
Sep 2002
Tracktown, U.S.A.

21358 Posts
Default

We are now listed on Henri and Renaud Lifchitz' page PRP Records with first, second, and fourth places! Our "production score" climbed from 45.70 to 47.76, and the record probable prime I found a year ago is now in seventh place. Go team!
philmoore is offline   Reply With Quote
Old 2009-09-09, 06:48   #26
geoff
 
geoff's Avatar
 
Mar 2003
New Zealand

22058 Posts
Default

So that is why they call Australia the lucky country. Not just gold mines and sunshine, they have all the big PRPs too :-)

Well done again Ben!
geoff is offline   Reply With Quote
Reply

Thread Tools


Similar Threads
Thread Thread Starter Forum Replies Last Post
probable largest prime. sudaprime Miscellaneous Math 11 2018-02-05 08:10
Gigantic Probable Prime Triplet found Cybertronic Twin Prime Search 18 2011-08-20 13:36
Fourth probable prime found, one to go! philmoore Five or Bust - The Dual Sierpinski Problem 22 2010-01-01 00:23
Another record probable prime found! philmoore Five or Bust - The Dual Sierpinski Problem 15 2009-02-08 19:43
Record probable prime found! philmoore Five or Bust - The Dual Sierpinski Problem 18 2009-01-28 19:47

All times are UTC. The time now is 07:24.

Tue Sep 29 07:24:24 UTC 2020 up 19 days, 4:35, 0 users, load averages: 1.31, 1.42, 1.54

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.