mersenneforum.org  

Go Back   mersenneforum.org > Extra Stuff > Blogorrhea > storflyt32

Reply
 
Thread Tools
Old 2014-11-23, 19:34   #67
Batalov
 
Batalov's Avatar
 
"Serge"
Mar 2008
Phi(4,2^7658614+1)/2

36·13 Posts
Default

Quote:
Originally Posted by danaj View Post
What am I missing that makes that important enough to make a post about? It takes only a few milliseconds to determine that result.
You are not missing anything. The poster is an attention-seeker or a lame bot.
Batalov is offline   Reply With Quote
Old 2014-11-23, 19:51   #68
TheMawn
 
TheMawn's Avatar
 
May 2013
East. Always East.

11·157 Posts
Default

Tell us more about Patterns.
TheMawn is offline   Reply With Quote
Old 2014-11-23, 21:07   #69
danaj
 
"Dana Jacobsen"
Feb 2011
Bangkok, TH

22·227 Posts
Default

Quote:
Originally Posted by CRGreathouse View Post
Is there a special form for the divisors of a number of the form a*b^c + d, in cases like this where c is relatively large?
Short answer: PFGW claims it can optimize this with d "up to 42-bits" but I know nothing at all about the details. My pretests are just using basic trial division (gcd for tiny primes, then simple treesieve for larger).

In the case of the composite mentioned, there is a small factor so easily caught -- it takes under 0.04s for my is_prime to find it just using simple trial division tests based on the input size. OpenFPGW and Pari/GP don't do much trial division, but still don't take very long for the first PRP test to indicate composite (~30s and 4.5 minutes respectively). If I skip the pretests and force just BPSW or base-2 SPRP it's 2.5min for ntheory (as expected, PFGW is much faster at this size).

I don't think LLR would cover this, especially not with the /3. Ignoring the /3, PFGW-3.7.7 indicates it can optimize k * b^n +/- c for c to 42-bits. I haven't looked into it at all however -- the LLR code I've added to my software is braindead simple, just doing b=2, c = -1, k = not-div-by-3. c=+1 is easily done via BLS75-T5. Riesel discusses the k-div-3 case for k*2^n-1 Bosma (1993) has some info, and there's a good post from a few years ago on this forum that has a lot of details, but I haven't gone into it (and it doesn't answer your question anyway).

Berrizbeitia et al. have a generalization to A*m^s+w where w has restrictions.
danaj is offline   Reply With Quote
Old 2014-11-23, 21:09   #70
science_man_88
 
science_man_88's Avatar
 
"Forget I exist"
Jul 2009
Dumbassville

838410 Posts
Default

Quote:
Originally Posted by CRGreathouse View Post
Is there a special form for the divisors of a number of the form a*b^c + d, in cases like this where c is relatively large?
I hope that's rhetorical if(gcd(a,d)!=1 or gcd(b,d)!=1, then gcd(a,d) or gcd(b,d) is a divsor) or a typo.

Last fiddled with by science_man_88 on 2014-11-23 at 21:10
science_man_88 is offline   Reply With Quote
Old 2014-11-23, 21:23   #71
CRGreathouse
 
CRGreathouse's Avatar
 
Aug 2006

3·1,993 Posts
Default

Quote:
Originally Posted by science_man_88 View Post
I hope that's rhetorical if(gcd(a,d)!=1 or gcd(b,d)!=1, then gcd(a,d) or gcd(b,d) is a divsor) or a typo.
In this case (a = 23, b = 10, c = 46554, d = -11) both GCDs are 1. But I'm not talking about finding a single factor but a class of factors, just like prime factors of Mersenne numbers have a special form. In this case I checked for divisibility by all small primes directly (modular exponentiation) but I wondered if there was a faster way. (I was pretty sure that there was one and that people here would know more about this than I do.)
CRGreathouse is offline   Reply With Quote
Old 2014-11-23, 21:38   #72
CRGreathouse
 
CRGreathouse's Avatar
 
Aug 2006

3×1,993 Posts
Default

Quote:
Originally Posted by danaj View Post
In the case of the composite mentioned, there is a small factor so easily caught -- it takes under 0.04s for my is_prime to find it just using simple trial division tests based on the input size.
My experience was similar -- 57ms with modular trial division in PARI/GP*, or 92ms using a big GCD.

* Of course if I did it in PARI directly it would be a good bit faster but I just did it in gp.

Quote:
Originally Posted by danaj View Post
OpenFPGW and Pari/GP don't do much trial division, but still don't take very long for the first PRP test to indicate composite (~30s and 4.5 minutes respectively). If I skip the pretests and force just BPSW or base-2 SPRP it's 2.5min for ntheory (as expected, PFGW is much faster at this size).

I don't think LLR would cover this, especially not with the /3. Ignoring the /3, PFGW-3.7.7 indicates it can optimize k * b^n +/- c for c to 42-bits. I haven't looked into it at all however -- the LLR code I've added to my software is braindead simple, just doing b=2, c = -1, k = not-div-by-3. c=+1 is easily done via BLS75-T5. Riesel discusses the k-div-3 case for k*2^n-1 Bosma (1993) has some info, and there's a good post from a few years ago on this forum that has a lot of details, but I haven't gone into it (and it doesn't answer your question anyway).

Berrizbeitia et al. have a generalization to A*m^s+w where w has restrictions.
That's interesting, but you're talking about primality/compositeness testing, while I'm talking about (partial) factorization. That is: I want the equivalent of the 2kp + 1 for Mersenne numbers, not the equivalent of LLR/BLS.

Last fiddled with by CRGreathouse on 2014-11-23 at 21:38
CRGreathouse is offline   Reply With Quote
Old 2014-11-23, 22:33   #73
danaj
 
"Dana Jacobsen"
Feb 2011
Bangkok, TH

38C16 Posts
Default

Quote:
Originally Posted by CRGreathouse View Post
That's interesting, but you're talking about primality/compositeness testing, while I'm talking about (partial) factorization. That is: I want the equivalent of the 2kp + 1 for Mersenne numbers, not the equivalent of LLR/BLS.
I realized that partway through writing my reply, but kept on with it anyway . I claim ignorance of special factor forms, but would enjoy seeing any positive responses.
danaj is offline   Reply With Quote
Old 2014-11-24, 00:56   #74
R.D. Silverman
 
R.D. Silverman's Avatar
 
Nov 2003

22×5×373 Posts
Default

Quote:
Originally Posted by danaj View Post
What am I missing that makes that important enough to make a post about? It takes only a few milliseconds to determine that result.
The last time that I dared to question the purpose of posting a result, I was (to summarize)
accused of being an asshole by Gordon in a fairly lengthy diatribe.

Allow me to point out that there has been no similar reply here.......
R.D. Silverman is offline   Reply With Quote
Old 2014-11-24, 02:47   #75
LaurV
Romulan Interpreter
 
LaurV's Avatar
 
Jun 2011
Thailand

7·1,373 Posts
Default

To be back on track: just pointing out that, according with factorDB, this number is known to be composite since (about) March 2011. As Batalov said, the guy is either innocent (in the bad sense) or a bot.

[edit: I am also interested in the discussion started by CRG's question, for some "general" case]

Last fiddled with by LaurV on 2014-11-24 at 02:48
LaurV is offline   Reply With Quote
Old 2014-11-25, 06:05   #76
storflyt32
 
Feb 2013

2·229 Posts
Default

Quoting here.

Batalov said:

I am moving this discussion to "Homework help".
It has nothing to do with new Fermat Factors.

If you take a hammer and a screw, hammer the screw into the wall and declare: "the hammer is not working and the screw is all bent", -- this whole thing only demonstrates that you use a hammer where you need a screwdriver.

-

Oh, definitely we missed some factors between the P252 and the P564 (Fermat).

Or, what the heck else.

Getting back to it, of course.

Last fiddled with by storflyt32 on 2014-11-25 at 06:06
storflyt32 is offline   Reply With Quote
Old 2014-11-25, 22:53   #77
TheMawn
 
TheMawn's Avatar
 
May 2013
East. Always East.

11·157 Posts
Default

Quote:
Originally Posted by storflyt32 View Post
Quoting here.

Batalov said:

I am moving this discussion to "Homework help".
It has nothing to do with new Fermat Factors.

If you take a hammer and a screw, hammer the screw into the wall and declare: "the hammer is not working and the screw is all bent", -- this whole thing only demonstrates that you use a hammer where you need a screwdriver.

-

Oh, definitely we missed some factors between the P252 and the P564 (Fermat).

Or, what the heck else.

Getting back to it, of course.

And that pretty much seals it. Bot.

TheMawn is offline   Reply With Quote
Reply



Similar Threads
Thread Thread Starter Forum Replies Last Post
Composite being Prime kruoli FactorDB 5 2018-02-16 16:54
How can I prove this PRP prime? siegert81 Math 2 2014-11-19 10:24
How do I prove a 4000 digit number is prime?? VJS Lounge 4 2005-05-09 20:56
How do you prove a number is prime? Alien Math 12 2004-01-07 11:36
why not stop when composite is prove? mdjvz Software 4 2003-09-28 17:13

All times are UTC. The time now is 19:56.


Fri Jul 16 19:56:57 UTC 2021 up 49 days, 17:44, 1 user, load averages: 1.81, 2.02, 2.30

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.