Wolfram Function Repository
Instant-use add-on functions for the Wolfram Language
Function Repository Resource:
Factor an integer using Fermat's factorization algorithm
ResourceFunction["FactorIntegerFermat"][n] factors the integer n using Fermat's algorithm. |
A trivial factorization for Fermat's algorithm:
| In[1]:= |
| Out[1]= |
A slightly larger example:
| In[2]:= | ![]() |
| Out[2]= |
In the following example, we see that the factor returned by FactorIntegerFermat is close to
:
| In[3]:= | ![]() |
| Out[4]= |
The ratio of size difference between the square root and a factor to the size of the square root is around 1/2, which means the factor is removed from the square root of bign by approximately its fourth root.
| In[5]:= |
| Out[6]= |
| In[7]:= |
| Out[7]= |
| In[8]:= |
| Out[8]= |
FactorIntegerFermat scales well to larger semiprimes, providing the two primes are close to
:
| In[9]:= | ![]() |
| Out[12]= | ![]() |
Factoring the following semiprime requires 262143 iterations, which is more than the default number of iterations (216):
| In[13]:= | ![]() |
| Out[13]= |
Increasing MaxIterations obtains the factorization:
| In[14]:= | ![]() |
| Out[14]= | ![]() |
Lehman extended Fermat's algorithm to factor N=p q using a rational approximation to the ratio of the factors, p/q:
| In[15]:= |
| Out[15]= |
| In[16]:= | ![]() |
| Out[18]= |
This works as 10277/1199 is a good approximation to p/q:
| In[19]:= |
| Out[19]= |
Unlike FactorInteger, FactorIntegerFermat is not a general purpose factoring algorithm:
| In[20]:= |
| Out[20]= |
| In[21]:= |
| Out[21]= |
| In[22]:= |
| Out[22]= |
FactorInteger succeeds:
| In[23]:= |
| Out[23]= |
FactorIntegerFermat does not handle this case:
| In[24]:= |
| Out[24]= |
Wolfram Language 13.0 (December 2021) or above
This work is licensed under a Creative Commons Attribution 4.0 International License