在数论中,平方同余是个经常被用于整数分解算法的同余关系。
给定一正整数 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 乘在一起去得到一个平方数(在等号的右侧)。
分解 35[编辑]
我们取 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)
分解 1649[编辑]
设 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)