mersenneforum.org  

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

Reply
 
Thread Tools
Old 2017-09-24, 20:43   #23
alpertron
 
alpertron's Avatar
 
Aug 2002
Buenos Aires, Argentina

2·683 Posts
Default

Due to Hasse theorem, the prime factor p can be found by ECM if the number p-1+q is smooth with respect to the selected bounds, where the absolute value of q is less than twice the square root of p. The value of q is unknown and depends on p and the elliptic curve being used.
alpertron is offline   Reply With Quote
Old 2017-09-25, 14:22   #24
LaurV
Romulan Interpreter
 
LaurV's Avatar
 
Jun 2011
Thailand

7·1,373 Posts
Default

Quote:
Originally Posted by alpertron View Post
Due to Hasse theorem, the prime factor p can be found by ECM if the number p-1+q is smooth with respect to the selected bounds, where the absolute value of q is less than twice the square root of p. The value of q is unknown and depends on p and the elliptic curve being used.
Correct.
LaurV is offline   Reply With Quote
Old 2017-09-25, 15:33   #25
alpertron
 
alpertron's Avatar
 
Aug 2002
Buenos Aires, Argentina

136610 Posts
Default

For example, in December 1995 Richard Brent found the prime factor p40 of the tenth Fermat number, completing its factorization.

p40 = 9409853205696664168149671432955079744397

The elliptic curve selected has sigma = 48998398, and the group order for the elliptic curve mod p40 can be computed as G = p40 - 1 + 149256697246264631024 that is smooth:

G = 22 * 3 * 5 * 17 * 312 * 53 * 67 * 233 * 739 * 5563 * 7901 * 20201 * 57163 * 309335137

Notice that p40 - 1 = 22 * 7715261 * 304910397901531269264567700073759 so this prime cannot be found using p-1 algorithm.
alpertron is offline   Reply With Quote
Old 2017-09-26, 05:50   #26
rcv
 
Dec 2011

14310 Posts
Default

Quote:
Originally Posted by LaurV View Post
Quote:
Originally Posted by RDS, few minutes ago

Enough!

Let n be composite. Assume that n is odd.
(if n is even then 2^n-1 splits as the difference of two squares: (2^(n/2) - 1)* (2^(n/2) + 1) )

There are three kinds of factors of 2^n-1.

(1) Algebraic. These are factors of 2^m-1 for m | n.
(2) Intrinsic. If n = p*m and 2^m-1 is divisible by p, then
2^pm -1 will be divisible by p^2 (i.e. you pick up another power of p via Hensel's Lemma)
(3) Primitive. All primitive prime factors of 2^n-1 are
primes that do not divide 2^j - 1 for j < n. All primitive prime factors are 1 mod 2n.
Thank you, RDS, for the marvelously simple condition for when the intrinsic factor occurs!

Also thank you for the information which leads to a clarification to my post #8 of this thread.

I said formula (1) was the "primitive". To clarify, I believe formula (1) is often called "phi". Where "phi" is usually the full "primitive", but "phi" may be the product of a small prime "intrinsic" times the "primitive". So, correcting formula (1) from my previous post:

The "phi" of that form is \frac{(x^{p\cdot q}-1) (x^1-1)}{(x^p-1)(x^q-1)}=\text{<primitive> or <instrinsic>}\times\text{<primitive>}\qquad\qquad(1)

My personal database (of slightly extended Cunningham forms) does have separate phi and primitive columns, so I have no excuse for being so careless with my explanation. An SQL query of my database shows 33 official Cunningham numbers with base 2 (exponent limits of 1300, 1300, and 2600) that possess an intrinsic factor in my database. That may or may not match the official count, but it should give the OP a rough (order of magnitude) estimate of how often intrinsic factors must be considered. (With the RDS-provided condition telling the OP exactly when intrinsic factors must be considered!)
rcv is offline   Reply With Quote
Old 2017-09-26, 05:54   #27
LaurV
Romulan Interpreter
 
LaurV's Avatar
 
Jun 2011
Thailand

7×1,373 Posts
Default

Quote:
Originally Posted by alpertron View Post
G = p40 - 1 + 149256697246264631024
Ok, and now you tell, for normal people like us who have no idea, where the 149256697246264631024 comes from?
Hehe...
LaurV is offline   Reply With Quote
Old 2017-09-26, 12:17   #28
alpertron
 
alpertron's Avatar
 
Aug 2002
Buenos Aires, Argentina

101010101102 Posts
Default

For prime numbers of that size, you can compute the group order by using Schoof's algorithm. There are faster algorithms for prime numbers greater than 100 digits, for example Schoof–Elkies–Atkin algorithm.

Last fiddled with by alpertron on 2017-09-26 at 12:17
alpertron is offline   Reply With Quote
Old 2017-09-26, 14:33   #29
Dr Sardonicus
 
Dr Sardonicus's Avatar
 
Feb 2017
Nowhere

10010001000112 Posts
Default

Quote:
Originally Posted by alpertron View Post
For example, in December 1995 Richard Brent found the prime factor p40 of the tenth Fermat number, completing its factorization.

p40 = 9409853205696664168149671432955079744397

The elliptic curve selected has sigma = 48998398, and the group order for the elliptic curve mod p40 can be computed as G = p40 - 1 + 149256697246264631024 that is smooth:

G = 22 * 3 * 5 * 17 * 312 * 53 * 67 * 233 * 739 * 5563 * 7901 * 20201 * 57163 * 309335137

Notice that p40 - 1 = 22 * 7715261 * 304910397901531269264567700073759 so this prime cannot be found using p-1 algorithm.
I don't know much at all about ECM. Hmm, I can probably look up the factorization of F10 (google google) Ahh, there it is! (read read)

Two things: In

FACTORIZATION OF THE TENTH FERMAT NUMBER

they refer to the expression p + 1 - q, not p - 1 + q.
Quote:
Let g = |G| be the order of G. We have g = p + 1 − t, where t satisfies Hasse’s inequality t^2 < 4p. The number of curves with given t can be expressed in terms of the Kronecker class number of t^2 − 4p (see [46]).
Second, the values you cite are for a different number:
Quote:
In December 1995, using program C with B1 = 370,000, D = 255, e = 6, we found the 40-digit factor

p'40 = 9409853205696664168149671432955079744397

of

p252 − 1,

where p252 is the largest prime factor of F10 . See §9 for the application to proving primality of p252 . The curve is defined as in §2.3 with

σ = 48998398,

and the group order is

g = 2^2*3*5*17*31^2*53*67*233*739*5563*7901*20201*57163*309335137.
You can check that for this p and g, p + 1 - g = 149256697246264631022.

As to the elliptic curve that found the p40 prime factor of F10,
Quote:
The group order is

g = p40 + 1 − 3674872259129499038 = 2^2*3^2*5*149*163*197*7187*18311*123677*226133*314263*4677853
Dr Sardonicus is offline   Reply With Quote
Old 2018-04-11, 20:59   #30
bhelmes
 
bhelmes's Avatar
 
Mar 2016

1010110012 Posts
Default

Dear Charles,

the mersenne numbers are embedded in the quadratic function f(n)=2n²-1
where for all n equal powers of two the f(n) is a mersenne number.

All primes p sieves double periodically the function values of f(n)
p|f(n) -> p|f(n+kp) and p|f(n-kp), k element of N

The mersenne numbers build a folding of this quadratic function.

I think you can use the prime factors of the mersenne numbers in order to make a sieve if you consider the differences between the mersenne numbers based on the regarding n of the quadratic function.

Is it possible to calculate all factors of mersenne numbers
if p of Mp is not a prime ?

It might not be the fastest way.

The result of the factorisation of mersenne numbers
are part of the sieving construction of that quadratic function f(n)=2n^2-1

By the way the function f(n)=2n^2-1 is embedded in the
function f(u,v)=u^2 - v^2 +2uv which seems to be a metric ?

Was it proven that there are infinitive primes p=2n^2-1 ?

Greetings from the primes
Bernhard

Last fiddled with by bhelmes on 2018-04-11 at 21:40
bhelmes is offline   Reply With Quote
Old 2018-04-11, 22:46   #31
JeppeSN
 
JeppeSN's Avatar
 
"Jeppe"
Jan 2016
Denmark

23·3·7 Posts
Question

Quote:
Originally Posted by LaurV View Post
This reminds me that long ago I requested eliminating the "exponent must be prime" criteria from mfaktX, or giving the possibility to the user to have a command line switch to disable the primality test - there was the time when I was trying to factor 2^18xx-1 (I forgot the xx, sorry, but it was a composite exponent in 1800's, and right now the FDB looks down too, so I can not look for it) which had (and maybe still has) a composite cofactor of about 300 digits left.
(A long time after.)

I am thinking 2^1891 - 1 because 1891 = 31*61, so that famous primes M31 and M61 are the algebraic factors?

/JeppeSN
JeppeSN is offline   Reply With Quote
Old 2018-04-12, 04:27   #32
LaurV
Romulan Interpreter
 
LaurV's Avatar
 
Jun 2011
Thailand

7·1,373 Posts
Default

That has a C5xx left, it was not so high, I may be wrong about the "18" part. Grr... I am getting older... Anyhow, the discussion is here on the forum somewhere, but it does not matter anymore.
LaurV is offline   Reply With Quote
Old 2018-04-12, 14:25   #33
CRGreathouse
 
CRGreathouse's Avatar
 
Aug 2006

3×1,993 Posts
Default

Quote:
Originally Posted by bhelmes View Post
Is it possible to calculate all factors of mersenne numbers
if p of Mp is not a prime ?
In general this is very hard.

Quote:
Originally Posted by bhelmes View Post
Was it proven that there are infinitive primes p=2n^2-1 ?
No, this is beyond current capabilities.
CRGreathouse is offline   Reply With Quote
Reply



Similar Threads
Thread Thread Starter Forum Replies Last Post
Composite P-1 factors not showing up under recently cleared exponents? ixfd64 PrimeNet 2 2018-02-28 07:54
PrimeNet has composite factors recorded James Heinrich PrimeNet 4 2011-09-16 14:40
is M21934219 composite ? S485122 Data 50 2011-01-10 10:28
Missing factors at the 'Known Factors' page MatWur-S530113 PrimeNet 11 2009-01-21 19:08
F10,21=10^(2^21)+1 is composite Shaopu Lin Factoring 2 2004-10-31 13:48

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


Fri Jul 16 18:47:45 UTC 2021 up 49 days, 16:35, 1 user, load averages: 2.78, 4.33, 4.50

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.