mersenneforum.org  

Go Back   mersenneforum.org > Great Internet Mersenne Prime Search > Math

Reply
 
Thread Tools
Old 2008-03-18, 00:55   #23
Mini-Geek
Account Deleted
 
Mini-Geek's Avatar
 
"Tim Sorbera"
Aug 2006
San Antonio, TX USA

17×251 Posts
Default

Quote:
Originally Posted by Anonymous View Post
Using the latest sr1sieve, I found that 2^3355584+1 is divisible by 5. I verified the factor with Pari/GP (which as far as I know uses plain old straight division, with no "shortcuts").
I did the same thing, for k=1-5 n=3355583 and checked all factors with Pari/GP.
Mini-Geek is offline   Reply With Quote
Old 2008-03-18, 00:58   #24
ColdFury
 
ColdFury's Avatar
 
Aug 2002

26×5 Posts
Default

Is this the same "friend" who had the radical compression technique?
ColdFury is offline   Reply With Quote
Old 2008-03-18, 01:01   #25
mdettweiler
A Sunny Moo
 
mdettweiler's Avatar
 
Aug 2007
USA (GMT-5)

186916 Posts
Default

Quote:
Originally Posted by Mini-Geek View Post
I did the same thing, for k=1-5 n=3355583 and checked all factors with Pari/GP.
I deleted my message, because I thought I put the output into Pari/GP wrong...but on second thought, though, it looks like I did do it correctly. :whoops: Can a mod please bring back my message?

Last fiddled with by xilman on 2008-03-18 at 07:50 Reason: Your wish is my command.
mdettweiler is offline   Reply With Quote
Old 2008-03-18, 01:07   #26
paulunderwood
 
paulunderwood's Avatar
 
Sep 2002
Database er0rr

22×941 Posts
Default

Code:
? n=3355584;forprime(p=3,1000,if(Mod(2,p)^n+1==0,print(p)))
769
The Windows version of NewPGen will not let me stop a test with divisor less that 1000 nor log divisors less than 2^32. I challenge jasong and his friends to reproduce their divisor of 2 by NewPGen

Last fiddled with by paulunderwood on 2008-03-18 at 01:24
paulunderwood is offline   Reply With Quote
Old 2008-03-18, 01:20   #27
R.D. Silverman
 
R.D. Silverman's Avatar
 
Nov 2003

22×5×373 Posts
Default

Quote:
Originally Posted by ewmayer View Post
Yes, Bob is dating himself, but Bob has an admirably reckless disregard for fashion and fads. The rest of us, slaves to hipster trends as we are, probably would prefer more up-to-date jargon such as "Dude, put down the crack pipe!" or even more trendily, the text-message-y "d00d, R U Hi?"
Crack gets one high. But quaaludes cause hallucinations similar to acid...
R.D. Silverman is offline   Reply With Quote
Old 2008-03-18, 01:49   #28
axn
 
axn's Avatar
 
Jun 2003

5,087 Posts
Default

Hmmm.... I might have run into an interesting bug in NewPGen. When I sieved k*b^n+1 form with (b=2,n=3355583,kmin=2,kmax=16), NewPGen declared all candidates as composite, instantly. However, if i give kmin=kmax=2, NewPGen just happily keeps going without finiding any factors!!

@Jasong, Try using the sieve option "b^n+k" with b=2, nmin=nmax=3355584 and k=1. NewPGen will instantly find the factor 769.
axn is offline   Reply With Quote
Old 2008-03-18, 02:29   #29
paulunderwood
 
paulunderwood's Avatar
 
Sep 2002
Database er0rr

22×941 Posts
Default

Quote:
Originally Posted by axn1 View Post
Hmmm.... I might have run into an interesting bug in NewPGen. When I sieved k*b^n+1 form with (b=2,n=3355583,kmin=2,kmax=16), NewPGen declared all candidates as composite, instantly. However, if i give kmin=kmax=2, NewPGen just happily keeps going without finiding any factors!!
It is a bug only if you not specify "include even k":

Code:
Base = 2
Include even k's ? (y) : n
n = 3355583
kmin = 2
kmax = 2
Type : k.b^n+1
The output file newout already exists - continue? (y) : y
Are you ready to start sieving ? (y) : y
CPU capabilities: CMOV: Supported. SSE2: Supported.
Using bitmap : allocating 0.0Kb of RAM...
...succeeded
0 candidates
Sieving numbers with 1010132 decimal digits
Initialising sieve...
Done.
02:26:06 Starting sieve (reduced screen output while p < 1 billion)
Stopping...o 33554467 (33.5 million)
02:26:10 Saving results to newout...
02:26:10 ...done.
Output file 'newout' generated, contains 0 k's
Paul ought to iron this out

I have emailed him

Last fiddled with by paulunderwood on 2008-03-18 at 02:39
paulunderwood is offline   Reply With Quote
Old 2008-03-18, 02:35   #30
Jay
 
Jay's Avatar
 
Dec 2007

3410 Posts
Default

Quote:
Originally Posted by jasong View Post
If you tell NewPGen to sieve k=1 to 5 for the equation k*2^3355583+1, it will say everything has been sieved and all have been found composite, and it will do that in less than a second.
Quote:
Originally Posted by jasong View Post
I'll repeat my challenge. Find a factor for 2^3355584+1 and test it without shortcuts. If my friends(or tolerators, possibly) at Riesel Sieve tell me that they've duplicated the result, then I'll shut up. Otherwise, you should seriously consider showing so-called cranks more respect.
Not certain which you want factors for, but I'll assume you mean the second.

The following factors were found via ECM using Prime95.

2^3355584+1: 769*274177*961934081

Using GMP I then obtained 2^3355584+1 Mod(f) for each factor obtained via GMP. Each resulted in a mod result of 0.

My understanding is that Prime95 uses an FFT computation in the interests of speed, while GMP uses integer math. However I could very easily be incorrect regarding these statements. However both packages are open source and nobody has ever found any errors in the computation methods of either.

Two different methods used to obtain the factors and to test the validity of factors. Both methods open source and reviewed by a large number of people.

I doubt that any conspiracy exists with the goal to ignore potential primes. It appears that 2^3355584+1 not only tests as being a composite, but also has (some) easily identifiable factors.

Was there anything else your friend would like to state, for the record?
Jay is offline   Reply With Quote
Old 2008-03-18, 03:19   #31
axn
 
axn's Avatar
 
Jun 2003

5,087 Posts
Default

Quote:
Originally Posted by paulunderwood View Post
It is a bug only if you not specify "include even k":
D'oh! Of course. Still a bug, just not the one I imagined.
axn is offline   Reply With Quote
Old 2008-03-18, 03:47   #32
jasong
 
jasong's Avatar
 
"Jason Goatcher"
Mar 2005

3·7·167 Posts
Default

comment deleted.

Last fiddled with by jasong on 2008-03-18 at 04:25
jasong is offline   Reply With Quote
Old 2008-03-18, 03:51   #33
Jay
 
Jay's Avatar
 
Dec 2007

2×17 Posts
Default

Quote:
Originally Posted by jasong View Post
bc is running the equation ((2^3355584)+1)/789 right now.

Oops, the file gave out a number that was cut off at the top after it was printed. I want to make sure I'm doing this right, so after I solve one, I can turn around and plug in the answer in the place 789 is above. If 789 is a factor, than switching the numbers should get the result 789.
If you're attempting to triple check the factors I found, you may want to use the correct factors. I said

2^3355584+1: 769*274177*961934081
Jay is offline   Reply With Quote
Reply

Thread Tools


Similar Threads
Thread Thread Starter Forum Replies Last Post
Proof of Legendre's conjecture, that there is always a prime between n^2 and (n+1)^2 MarcinLesniak Miscellaneous Math 41 2018-03-29 16:30
Legendre's prime counting function pbewig Information & Answers 0 2011-07-14 00:47
What is Legendre Symbol? slowing down sr2sieve? cipher Software 3 2009-05-20 13:35
Computing n-th power residue symbols geoff Sierpinski/Riesel Base 5 2 2006-10-24 00:09
defective memory chip? ixfd64 Hardware 2 2004-11-28 05:45

All times are UTC. The time now is 23:20.


Fri Aug 6 23:20:58 UTC 2021 up 14 days, 17:49, 1 user, load averages: 4.08, 4.11, 4.06

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.