Fast prime factorization in Python. Factors most 50-60 digit numbers within a minute or so (with PyPy).
The algorithm used depends on the size of the input
pollardPm1.pycontains an implementation of the large prime (two stage) variant of Pollard's p-1 algorithm.pollardRho.pycontains an implementation of Pollard's Rho algorithm with Brent's improvements.ecm.pycontains an implementation of Lenstra's elliptic curve factorization algorithm. It is inversionless (since it uses Montgomery coordinates), uses two stages, and uses Suyama's parametrization to generate random elliptic curves. It also contains an implementation of Montgomery's PRAC algorithm for scalar multiplication (thanks Paul Zimmerman!) but this turned out to be slower than the usual double-and-add algorithm weirdly.primeSieve.pycontains a bunch of prime sieves (Atkin, Eratosthenes, segmented Eratosthenes). Look at the file for specific benchmarks.
All you have to do is run the file factor.py, enter a number, and hit Enter. Here's an example in terminal:
python factor.py
Enter a number: 15
Factoring 15...
Number of digits: 2
Finding small prime factors...
Prime factors found: 3, 5
15 = 3^1 * 5^1
Time: 5.00679016113e-05 s
and another...
Enter number: 37897387397398739739826929827929827927927762729872987928
Factoring 37897387397398739739826929827929827927927762729872987928...
Number of digits: 56
Finding small prime factors...
Prime factors found: 2, 3
Factoring 1579057808224947489159455409497076163663656780411374497 with ECM...
Number of digits: 55
Bounds: 250000 128992510
Sieving primes...
Stage 2 found factor!
Found factor 67246307
Factoring 67246307...
Number of digits: 8
67246307 is prime!
Factoring 23481702991138940747474138758238071923617408171...
Number of digits: 47
Factoring 23481702991138940747474138758238071923617408171 with ECM...
Number of digits: 47
Bounds: 50000 12746592
Sieving primes...
Tried 40 random curves...
Tried 80 random curves...
Tried 120 random curves...
Tried 160 random curves...
Stage 2 found factor!
Found factor 4788272261623351
Factoring 4788272261623351...
Number of digits: 16
4788272261623351 is prime!
Factoring 4904003303934522319753958187821...
Number of digits: 31
4904003303934522319753958187821 is prime!
37897387397398739739826929827929827927927762729872987928 = 2^3 * 3^1 * 67246307^1 * 4788272261623351^1 * 4904003303934522319753958187821^1
Time: 24.7774269581 s
- A.O.L Atkin, D.J.Bernstein; Prime Sieves using Binary Quadratic Forms; Mathematics of Computation, 73-246: 1023-30
- Peter L Montgomery; Speeding the Pollard and Elliptical Methods of Factorization; Mathematics of Computation (Jan 1987), Issue 177: 243-264
- Montgomery, P.L.; Evaluating Recurrences of the form Xm+n = f(Xm, Xn, Xm-n) via Lucas Chains; Unpublished manuscript (Jan 1992)