Recessional velocity

Shanks' square forms factorization is a method for integer factorization devised by Daniel Shanks as an improvement on Fermat's factorization method.

The success of Fermat's method depends on finding integers $x$ and $y$ such that $x^{2}-y^{2}=N$, where $N$ is the integer to be factored. An improvement (noticed by Kraitchik) is to look for integers $x$ and $y$ such that $x^{2}\equiv y^{2}{\pmod {N}}$. Finding a suitable pair $(x,y)$ does not guarantee a factorization of $N$, but it implies that $N$ is a factor of $x^{2}-y^{2}=(x-y)(x+y)$, and there is a good chance that the prime divisors of $N$ are distributed between these two factors, so that calculation of the greatest common divisor of $N$ and $x-y$ will give a non-trivial factor of $N$.

A practical algorithm for finding pairs $(x,y)$ which satisfy $x^{2}\equiv y^{2}{\pmod {N}}$ was developed by Shanks, who named it Square Forms Factorization or SQUFOF. The algorithm can be expressed in terms of continued fractions, or in terms of quadratic forms. Although there are now much more efficient factorization methods available, SQUFOF has the advantage that it is small enough to be implemented on a programmable calculator.

Algorithm

Input: $N$, the integer to be factored, which must be neither a prime number nor a perfect square, and a small multiplier $k$.

Output: a non-trivial factor of $N$.

The algorithm:

Initialize $P_{0}=\lfloor {\sqrt {kN}}\rfloor ,Q_{0}=1,Q_{1}=kN-P_{0}^{2}.$

Repeat

$b_{i}=\left\lfloor {\frac {\lfloor {\sqrt {kN}}\rfloor +P_{i-1}}{Q_{i}}}\right\rfloor ,P_{i}=b_{i}Q_{i}-P_{i-1},Q_{i+1}=Q_{i-1}+b_{i}(P_{i-1}-P_{i})$

until $Q_{i}$ is a perfect square.

Initialize $b_{0}=\left\lfloor {\frac {\lfloor {\sqrt {kN}}\rfloor -P_{i-1}}{\sqrt {Q_{i}}}}\right\rfloor ,P_{0}=b_{0}{\sqrt {Q_{i}}}+P_{i-1},Q_{0}={\sqrt {Q_{i}}},Q_{1}={\frac {kN-P_{0}^{2}}{Q_{0}}}$

Repeat

$b_{i}=\left\lfloor {\frac {\lfloor {\sqrt {kN}}\rfloor +P_{i-1}}{Q_{i}}}\right\rfloor ,P_{i}=b_{i}Q_{i}-P_{i-1},Q_{i+1}=Q_{i-1}+b_{i}(P_{i-1}-P_{i})$

until $P_{i+1}=P_{i}.$

Then if $f=\gcd(N,P_{i})$ is not equal to $1$ and not equal to $N$, then $f$ is a non-trivial factor of $N$. Otherwise try another value of $k$.

Shanks' method has time complexity $O({\sqrt[{4}]{N}})$.

Stephen S. McMasters (see link in External Link section) wrote a more detailed discussion of the mathematics of Shanks' method, together with a proof of its correctness.

Example

N = 11111, k = 1

P₀ = 105 Q₀ = 1 Q₁ = 86

P₁ = 67 Q₁ = 86 Q₂ = 77

P₂ = 87 Q₂ = 77 Q₃ = 46

P₃ = 97 Q₃ = 46 Q₄ = 37

P₄ = 88 Q₄ = 37 Q₅ = 91

P₅ = 94 Q₅ = 91 Q₆ = 25

Here Q₆ is a perfect square

P₀ = 104 Q₀ = 5 Q₁ = 59

P₁ = 73 Q₁ = 59 Q₂ = 98

P₂ = 25 Q₂ = 98 Q₃ = 107

P₃ = 82 Q₃ = 107 Q₄ = 41

P₄ = 82

Here P₃ = P₄

gcd(11111, 82) = 41, which is a factor of 11111.

References

• 20 year-old Real Estate Agent Rusty from Saint-Paul, has hobbies and interests which includes monopoly, property developers in singapore and poker. Will soon undertake a contiki trip that may include going to the Lower Valley of the Omo.
  
  My blog: http://www.primaboinca.com/view_profile.php?userid=5889534
• 20 year-old Real Estate Agent Rusty from Saint-Paul, has hobbies and interests which includes monopoly, property developers in singapore and poker. Will soon undertake a contiki trip that may include going to the Lower Valley of the Omo.
  
  My blog: http://www.primaboinca.com/view_profile.php?userid=5889534
• 20 year-old Real Estate Agent Rusty from Saint-Paul, has hobbies and interests which includes monopoly, property developers in singapore and poker. Will soon undertake a contiki trip that may include going to the Lower Valley of the Omo.
  
  My blog: http://www.primaboinca.com/view_profile.php?userid=5889534

External links

• 2005, Stephen S. McMathTemplate:Dead link
• Analysis and Improvement of the Continued Fraction Method of Factorization, Daniel Shanks (transcribed by S.McMath 2004)
• SQUFOF Notes, Daniel Shanks (transcribed by S. McMath 2004)
• Square Form Factorisation, Jason Gower and Samuel Wagstaff