mersenneforum.org  

Go Back   mersenneforum.org > Extra Stuff > Miscellaneous Math

Reply
 
Thread Tools
Old 2010-09-08, 11:15   #430
science_man_88
 
science_man_88's Avatar
 
"Forget I exist"
Jul 2009
Dumbassville

26×131 Posts
Default

Quote:
Originally Posted by CRGreathouse View Post
I'm not sure what you're saying. Do you have an example of a prime larger than 10^1999 which does not qualify for #19?*

* Or should I say, #19 as of 07 Sep 10 08:41 PM, since these definitions are fairly malleable.
see (2.5*10^1999)*2^2+c = 10^2000 + c so just find all c that give primes and then check for the special forms he doesn't allow.


forstep(n=1,100,2,if(isprime(10^2000+n),print(n))) is what I've been using up to you what n you search.
science_man_88 is offline   Reply With Quote
Old 2010-09-08, 12:45   #431
CRGreathouse
 
CRGreathouse's Avatar
 
Aug 2006

3×1,993 Posts
Default

OK, so #19 as stated definitely covers precisely those primes greater than 10^1999. It suffices to use c\in\{101,103\} and b = n = 2.
CRGreathouse is offline   Reply With Quote
Old 2010-09-08, 13:04   #432
science_man_88
 
science_man_88's Avatar
 
"Forget I exist"
Jul 2009
Dumbassville

26·131 Posts
Default

Quote:
Originally Posted by CRGreathouse View Post
OK, so #19 as stated definitely covers precisely those primes greater than 10^1999. It suffices to use c\in\{101,103\} and b = n = 2.
I never knew that it had a covering set and my code doesn't give 101,103 because i start at 10^2000 not 10^1999 but pretty much.
science_man_88 is offline   Reply With Quote
Old 2010-09-08, 13:38   #433
CRGreathouse
 
CRGreathouse's Avatar
 
Aug 2006

135338 Posts
Default

Quote:
Originally Posted by science_man_88 View Post
I never knew that it had a covering set and my code doesn't give 101,103 because i start at 10^2000 not 10^1999 but pretty much.
4N + 1 and 4N + 3 cover all primes > 2, right? So do that with 4N + 101 and 4N + 103.

10^1999 is the first 2000-digit number, just like 10^9 is the first 10-digit number.

Last fiddled with by CRGreathouse on 2010-09-08 at 13:39
CRGreathouse is offline   Reply With Quote
Old 2010-09-08, 19:39   #434
CRGreathouse
 
CRGreathouse's Avatar
 
Aug 2006

3·1,993 Posts
Default

I have verified the primality of 7984559573504259856359124657, a p28, through trial division.
CRGreathouse is offline   Reply With Quote
Old 2010-09-09, 02:19   #435
3.14159
 
3.14159's Avatar
 
May 2010
Prime hunting commission.

69016 Posts
Default

Quote:
Originally Posted by Charles
I have verified the primality of 7984559573504259856359124657, a p28, through trial division.
I'm inclined to call bullshit, but, whatever, I won't bother. I'm busy looking for a 166350-digit prime anyway.

Odds: 1 in 7294.

So, 1 n 7294 candidates should be prime.

Each test should take about 10 minutes.

10 * 7294 = 72940 minutes ≈ 52 days.

Actually, it is 8.5 minutes. 8.5 * 7294 ≈ 43 days.

Last fiddled with by 3.14159 on 2010-09-09 at 02:30
3.14159 is offline   Reply With Quote
Old 2010-09-09, 04:10   #436
CRGreathouse
 
CRGreathouse's Avatar
 
Aug 2006

3·1,993 Posts
Default

Quote:
Originally Posted by 3.14159 View Post
I'm inclined to call bullshit, but, whatever, I won't bother.
That's your prerogative. I spent many hours trial-dividing (though, happily, my prediction on its timing was accurate -- I guessed the completion time within 5 minutes!).

Of course, your skepticism bolsters one of my points on this category: it's essentially unverifiable. By contrast, if we were using modern algorithms, we could provide certificates:
Code:
[2 5 1]

[3 3 1]

[11 3 1]

[103 3 1]

[149 2 1]

[9376643 3 1]
Someone proposed that you could ask for residues mod all the primes less than the square root of the number, but asymptotically that takes something like 2\sqrt n/\log 2 bits, or 32 TB in this case... using a mixed-radix method allows for smaller 'certificates', but I doubt the improvement is great.

Last fiddled with by CRGreathouse on 2010-09-09 at 04:12
CRGreathouse is offline   Reply With Quote
Old 2010-09-09, 04:14   #437
CRGreathouse
 
CRGreathouse's Avatar
 
Aug 2006

3×1,993 Posts
Cool

Quote:
Originally Posted by 3.14159 View Post
I'm busy looking for a 166350-digit prime anyway.

Odds: 1 in 7294.
If your base has no small prime factors, you must have sieved up to about 6 trillion to get odds that good.
CRGreathouse is offline   Reply With Quote
Old 2010-09-09, 04:26   #438
3.14159
 
3.14159's Avatar
 
May 2010
Prime hunting commission.

24×3×5×7 Posts
Default

Quote:
Originally Posted by Charles
If your base has no small prime factors, you must have sieved up to about 6 trillion to get odds that good.
I used b = 2, and indeed sieved to about 6.2 trillion.
3.14159 is offline   Reply With Quote
Old 2010-09-09, 04:28   #439
3.14159
 
3.14159's Avatar
 
May 2010
Prime hunting commission.

32208 Posts
Default

Code:
[2 5 1]

[3 3 1]

[11 3 1]

[103 3 1]

[149 2 1]

[9376643 3 1]
Quit using PARI, guy. It doesn't work for anything larger than 2000 digits.

Last fiddled with by 3.14159 on 2010-09-09 at 04:28
3.14159 is offline   Reply With Quote
Old 2010-09-09, 04:36   #440
CRGreathouse
 
CRGreathouse's Avatar
 
Aug 2006

3×1,993 Posts
Default

Quote:
Originally Posted by 3.14159 View Post
Quit using PARI, guy. It doesn't work for anything larger than 2000 digits.
Considering that the number has 28 digits, and 28 <= 2000, I think I was using the appropriate tool.

Basically, I'm interested in things that will expand my mind rather than simply test the speed of my CPU. When there's a good reason I'll use powerful tools outside of Pari, but for the most part if I need other tools I'm in the range of brute CPU calculations that I don't care about.
CRGreathouse is offline   Reply With Quote
Reply



Similar Threads
Thread Thread Starter Forum Replies Last Post
Prime posting thread, part 2. (With a catch.) 3.14159 Miscellaneous Math 55 2010-11-19 23:55
Tiny range request .... 555.1M petrw1 LMH > 100M 1 2010-07-13 15:35
Other primes thread nuggetprime No Prime Left Behind 32 2009-10-21 21:48
Error: tiny factoring failed 10metreh Msieve 26 2009-03-08 23:28
Tiny error on nfsnet pages. antiroach NFSNET Discussion 1 2003-07-08 00:27

All times are UTC. The time now is 22:51.


Fri Aug 6 22:51:09 UTC 2021 up 14 days, 17:20, 1 user, load averages: 3.71, 4.06, 3.90

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