在數論中,平方同餘是個經常被用於整數分解演算法的同餘關係。
給定一正整數 n,費馬因式分解法想找到兩數 x, y 滿足下列方程式:
![{\displaystyle x^{2}-y^{2}=n}](https://wikimedia.org/api/rest_v1/media/math/render/svg/6f0228f50232af30832df02e3ff40a35fd395198)
從上式我們便可以分解得 n = x2 - y2 = (x + y)(x - y)。 此演算法在實際用途上較慢因為我們需要試圖找出許多類似的數字,而只有一部分會滿足這個嚴格的式子。 然而, n 也可能可以被分解,如果我們能滿足以下較弱的平方同餘的情況:
![{\displaystyle x^{2}\equiv y^{2}{\pmod {n}}}](https://wikimedia.org/api/rest_v1/media/math/render/svg/bd665ed0e76232633eef21e39409409422637e4b)
![{\displaystyle x\not \equiv \pm y{\pmod {n}}}](https://wikimedia.org/api/rest_v1/media/math/render/svg/6e185c5b3c1d7c379b9a8a43e49b8b1735ca5c9d)
從上面二式我們可以輕易推得:
![{\displaystyle x^{2}-y^{2}\equiv 0{\pmod {n}}}](https://wikimedia.org/api/rest_v1/media/math/render/svg/0e861172a85faac7dfaa6c4a387e1a02ddfae7f4)
![{\displaystyle (x+y)(x-y)\equiv 0{\pmod {n}}}](https://wikimedia.org/api/rest_v1/media/math/render/svg/a656251c209ca1450dbec9abe9c76e40b6661cac)
這意味著 n 整除乘積 (x + y)(x - y)。 因此 (x + y) 以及 (x − y) 各自包含著 n 的因數, 但那些因數有可能是平凡的(等於 1 或是 n 本身)。 在這樣的情況下,我們需要找到其他的 x 和 y 。計算 (x + y, n) 和 (x - y, n) 的最大公因數將給我們這些因數;而這可以輾轉相除法快速得出。
平方同餘在整數分解演算法裡相當地有用、且廣為使用於,舉例來說,二次篩選法、 普通數域篩選法、連續分數因式分解法、 以及狄克森因式分解法。反過來說,因為要找到模一合數下的平方根概率上在多項式時間等價於分解該數, 任何整數分解演算法皆用於找到一個平方同餘數。
其可以使用因數基底去幫助更快找到平方同餘。 比起從頭開始找
,不如我們找到許多的
,其中 y 只有小的質因數, 並試著去將其中幾個 y 乘在一起去得到一個平方數(在等號的右側)。
我們取 n = 35 並得:
![{\displaystyle 6^{2}=36\equiv 1\equiv 1^{2}{\pmod {35}}}](https://wikimedia.org/api/rest_v1/media/math/render/svg/26e519ca3dcf62dffc95c523857c3047fcbcff1c)
因此分解成:
![{\displaystyle \gcd(6-1,35)\cdot \gcd(6+1,35)=5\cdot 7=35}](https://wikimedia.org/api/rest_v1/media/math/render/svg/8aa5d78871303a33b668520b880bc00e92bbea2e)
設 n = 1649, 並作為從一些非平方數的乘積建構出平方同餘的例子(參見 狄克森因式分解法), 首先我們得到幾個同餘式:
![{\displaystyle 41^{2}\equiv 32{\pmod {1649}}}](https://wikimedia.org/api/rest_v1/media/math/render/svg/6fb665037f4e6f670c357ba48eba31cf74b6137b)
![{\displaystyle 42^{2}\equiv 115{\pmod {1649}}}](https://wikimedia.org/api/rest_v1/media/math/render/svg/51b9ac8b604e4b6bbd22bfa5171434154745ff2f)
![{\displaystyle 43^{2}\equiv 200{\pmod {1649}}}](https://wikimedia.org/api/rest_v1/media/math/render/svg/64ec3ebf9bf91a7c47cb1aaea35e5990cea0445b)
在其之中,有兩個式子的數字只有小質數作為因數:
![{\displaystyle 32=2^{5}}](https://wikimedia.org/api/rest_v1/media/math/render/svg/8adac280c3894c00906bda191eaac08688cea4fc)
![{\displaystyle 200=2^{3}\cdot 5^{2}}](https://wikimedia.org/api/rest_v1/media/math/render/svg/15593c78d45aed8c9d949e9f92fe866f45eb895a)
且兩者的其中一個組合滿足質因數的次方皆為偶數,因此形成了一個平方數:
![{\displaystyle 32\cdot 200=2^{5+3}\cdot 5^{2}=2^{8}\cdot 5^{2}=(2^{4}\cdot 5)^{2}=80^{2}}](https://wikimedia.org/api/rest_v1/media/math/render/svg/b5b3aa5e847927fe629a17cd3bbc5a8a776e434e)
因此得出平方同餘的關係式:
![{\displaystyle 32\cdot 200=80^{2}\equiv 41^{2}\cdot 43^{2}\equiv 114^{2}{\pmod {1649}}}](https://wikimedia.org/api/rest_v1/media/math/render/svg/c40da692fcb13d3f7c4901bf04ad426158fa9443)
所以分別作為 x、y 值的 80 和 114 得出因數:
![{\displaystyle \gcd(114-80,1649)\cdot \gcd(114+80,1649)=17\cdot 97=1649.}](https://wikimedia.org/api/rest_v1/media/math/render/svg/788b84002a7e915588e0767878219d41c9073d5b)