mersenneforum.org

mersenneforum.org (https://www.mersenneforum.org/index.php)
-   Data (https://www.mersenneforum.org/forumdisplay.php?f=21)
-   -   Mersenne number factored (disbelievers are biting elbows) (https://www.mersenneforum.org/showthread.php?t=19407)

LaurV 2014-08-01 06:48

[QUOTE=Batalov;379476]GIMPS actually has all LL residues. Iirc from EWMayer, LL residues can also be used, but the interested party will have to dig up all details.[/QUOTE]
I don't believe that. Who reported full residues to GIMPS? What gimps have is the last hex digits of each residue, which can't be of any use here.

Batalov 2014-08-01 07:39

[QUOTE=ewmayer;342684]Testing the cofactor for rigorous primality is expensive, but PRP-testing (starting with the full LL residue) can be done similarly cheaply as cofactor-PRP testing is done for Fermat numbers (starting with a Pépin-test residue).

In the late 1990s I adapted the procedure used for Fermats to PRP-test a bunch of M(p) cofactors (and discovered numerous PRP cofactors up to my search limit of a then-lofty 20k digits or so). My version of the test needed a Pépin-style test residue to be generated for the target M(p), but Peter Montgomery [URL="http://ndatech.com/mersenne/archives/digest/v01_0368.txt"]showed me how to use an LL-test residue[/URL], thus eliminating the redundant LL-test/Pépin-style-mod-powering effort.[/QUOTE]
Scroll the linked discussion until "Saving storage when testing cofactors of Mp". I remembered correctly: the storage space [I]can [/I]be saved, but it is not the lowest 64bit that should have been saved but the few hundred [I]highest [/I]bits. So, as GIMPS "free-of-charge" residues go now, they probably won't help. Because LL test is a test in a weird base [TEX]\omega,\ \bar{\omega}[/TEX] (in diguise), the result is more twisted than in simple base 3. (Peter Montgomery showed how to use upper bits, ... I wonder if lower bits could be an equivalent fingerprint. Tonight I won't think about it.)

Anyway, Dario most likely has all Pépin-test residue RES64s now in his logs (up to a certain limit), so he has everything he needs for very fast new incoming factor tests within his range.

axn 2014-08-01 08:19

[QUOTE=Batalov;379476]All one needs to do is [U]keep all residues[/U]: 3^Mp (mod Mp) or even just their last 64 bits. [/QUOTE]

I don't see how the last 64 bits help.

Lets say F divides Mp. By Fermat, 3^(Mp/F) = 3 (mod (Mp/F)). Raising both sides to F, we get 3^Mp = 3^F (mod Mp/F). Now if we have the full residue 3^Mp (mod Mp), we can easily calculate 3^Mp (mod Mp/F). But last 64-bits? Are you saying 3^Mp = 3^F (mod Mp/F) ==> 3^Mp = 3^F (mod Mp). I just tried it with M2311 (77567729423209*4514379640917651135021865565129IPRP-652). LHS holds, but RHS doesn't.

alpertron 2014-08-01 11:28

[QUOTE=axn;379483]I don't see how the last 64 bits help.

Lets say F divides Mp. By Fermat, 3^(Mp/F) = 3 (mod (Mp/F)). Raising both sides to F, we get 3^Mp = 3^F (mod Mp/F). Now if we have the full residue 3^Mp (mod Mp), we can easily calculate 3^Mp (mod Mp/F). But last 64-bits? Are you saying 3^Mp = 3^F (mod Mp/F) ==> 3^Mp = 3^F (mod Mp). I just tried it with M2311 (77567729423209*4514379640917651135021865565129IPRP-652). LHS holds, but RHS doesn't.[/QUOTE]

Anyway, for numbers in this range, it is easier to perform a PRP test than to find a new prime factor of the Mersenne number. There is no double check of PRP, so there is still the doubt of having incorrect residues. This is why I would prefer to perform a PRP test on the new cofactor.

flagrantflowers 2014-08-04 19:48

Did some rudimentary math on the PRP data @mersenne.ca for the odds of finding a positive PRP result given the number of factors found (three or less versus four or more).

Below M(2M) there are 3384 exponents with 4 or more factors, of these 78 have a positive PRP result.
~2.3% chance of a positive prp result given four or more factors.

For three factors or less: total of ~110801 exponents factored, subtract exponents with more than four factors leaves us with ~107417.

There are 203 positive PRP results for 3 factors or less, gives:
~0.19% chance of a positive PRP result given three or less factors.

[B]So the odds of finding a positive PRP with more than three factors is about ten times higher than for exponents with three factors or less.
[/B]
This analysis does not differentiate between one, two and three factors. The data is not readily available but I get the feeling that this would further lower the odds of a positive PRP for a single factor.

Either way it seems that running PRPs on exponents with numerous factors is more fruitful than exponents with a single, or couple, of factors.

EDIT: there seem to only be ~75 positive PRP tests for a single factor. This is only ~25% of the positive PRP tests found to date. Yet the large majority of exponents only have one factor.

R.D. Silverman 2014-08-04 22:24

[QUOTE=flagrantflowers;379686]Did some rudimentary math on the PRP data @mersenne.ca for the odds of finding a positive PRP result given the number of factors found (three or less versus four or more).
[/QUOTE]

No, you did NOT do any math. What you did was mindless numerology
on a small sample set and then jump to some clearly erroneous
conclusions. Your conclusions, if correct would violate Erdos-Kac.

Stop babbling. Stop the pseudo-math. Stop pretending that what you
are doing is math. Stop equating mindless numerology with "performing
rudimentary math".

If, in fact, you think that you did some math, please show it to us.
What you have presented is not math and does not constitute meaningful
analysis of the data.

alpertron 2014-08-05 00:36

[QUOTE=flagrantflowers;379686]Did some rudimentary math on the PRP data @mersenne.ca for the odds of finding a positive PRP result given the number of factors found (three or less versus four or more).

Below M(2M) there are 3384 exponents with 4 or more factors, of these 78 have a positive PRP result.
~2.3% chance of a positive prp result given four or more factors.

For three factors or less: total of ~110801 exponents factored, subtract exponents with more than four factors leaves us with ~107417.

There are 203 positive PRP results for 3 factors or less, gives:
~0.19% chance of a positive PRP result given three or less factors.

[B]So the odds of finding a positive PRP with more than three factors is about ten times higher than for exponents with three factors or less.
[/B]
This analysis does not differentiate between one, two and three factors. The data is not readily available but I get the feeling that this would further lower the odds of a positive PRP for a single factor.

Either way it seems that running PRPs on exponents with numerous factors is more fruitful than exponents with a single, or couple, of factors.

EDIT: there seem to only be ~75 positive PRP tests for a single factor. This is only ~25% of the positive PRP tests found to date. Yet the large majority of exponents only have one factor.[/QUOTE]

Notice that in the range you are considering, the lowest numbers have much more factoring effort done, so I think the statistics are not correct. It is like comparing apples and oranges. For example, if you consider the range 500K to 2M you will find a completely different picture.

R.D. Silverman 2014-08-05 00:53

[QUOTE=alpertron;379712]Notice that in the range you are considering, the lowest numbers have much more factoring effort done, so I think the statistics are not correct. It is like comparing apples and oranges. For example, if you consider the range 500K to 2M you will find a completely different picture.[/QUOTE]

It isn't even that! He babbles about numbers with 3 factors, vs. numbers
with 4 factors, when all he is looking at is numbers with 3 SMALL factors;
i.e. almost all of the numbers with three factors are not even included in
his data because most of these have factors whose size lie beyond his
search range. What he is doing is MEANINGLESS.

alpertron 2014-08-05 12:03

There is a formula on section 4.5.4 (Factoring into Primes) of Knuth's TAOCP that computes the probability of having the second largest prime factor less than some bound (expressed as log B / log N, B is the bound, N is the composite number).

That can be used here, but it appears that the formula is not easy to compute. For example, we could fix the bound to 10[SUP]30[/SUP] and then for all prime exponents between 1M and 2M, add all probabilities. That sum could be (approximately) the expected number of PRP factors in the mentioned range that we could found in a reasonable amount of time.

flagrantflowers 2014-08-05 16:21

[QUOTE=alpertron;379712]Notice that in the range you are considering, the lowest numbers have much more factoring effort done, so I think the statistics are not correct. It is like comparing apples and oranges. For example, if you consider the range 500K to 2M you will find a completely different picture.[/QUOTE]


I considered ignoring <M(10K) as you're correct there are few exponents with more than one factor and the picture does drastically change. But if you ignore that range the odds of positive PRP results as a function of factors is more disparate between 1-3 and more than four.

I find it interesting that there are obviously much fewer (as a percent) of exponents with less than three factors versus more than four.

You're correct I didn't do any "math". I calculated some uncorrected percentages. The order of magnitude difference, despite my "pseudomath", seems to be too large to simply be dumb "numerology".

flagrantflowers 2014-08-05 16:28

[QUOTE=R.D. Silverman;379714]It isn't even that! He babbles about numbers…[/QUOTE]


Unfortunately we only have PRP tests completed up to M(2M) so going beyond that range doesn't make sense as I would be including mostly candidates without a PRP result (positive or negative).

All I'm saying is the percent of exponents with a positive PRP result and having more than three factors is a very large portion of positive tests. Despite the fact that there are much more numbers with one factor.

I'm not trying to say that this means anything. I'm trying to say that this is functionally what the odds are on the tests done to date. This is obviously from this sample set and doesn't, ideally, represent the population of all <M(2M).

I don't pretend to understand the underlying factor theorems.


All times are UTC. The time now is 06:47.

Powered by vBulletin® Version 3.8.11
Copyright ©2000 - 2021, Jelsoft Enterprises Ltd.