Description
求解 x2≡n(modp) 。 p 是一个奇质数。
Solution
由费马小定理
np−1≡1(modp)
所以
np−12≡±1(modp)
由欧拉准则
(np)≡np−12(modp)
其中
(np)
为勒让德符号。因为
x≡n12(modp)
,所以有
xp−1≡np−12≡1(modp)
。即
(np)=⎧⎩⎨⎪⎪1−10n在模p意义下是二次剩余n在模p意义下不是二次剩余n≡0(modp)
设
k=a2−n,ω=k√
则该方程的解为
x≡(a+ω)p+12(modp)
证明:若
k
是该模意义下的非二次剩余,则
(a+b)p≡ap+bp(modp)
则
x21≡≡≡≡≡≡≡≡(a+ω)p+1(a+ω)p(a+ω)(ap+ωp)(a+ω)(ap−1a+ωp−1ω)(a+ω)(a−ω)(a+ω)a2−ω2a2−(a2−n)n(modp)
另一解为
x2≡−x1(modp)
如何找到这个 k 呢。因为非二次剩余的个数大约有一半,随机几次即可。
后面去实现了一下,发现要重新定义一个复数域
struct Complex {
ll r, i;
Complex(ll _r = 0, ll _i = 0):r(_r), i(_i) {}
inline Complex operator +(ll x) {
return Complex((x + r) % P, i);
}
inline Complex operator *(Complex a) {
ll rr = (a.r * r % P + w * a.i % P * i % P + P * 3) % P,
ii = (a.i * r % P + a.r * i % P) % P;
return Complex(rr, ii);
}
};
inline ll Qr(ll x, ll p) {
if (Pow(x, (P - 1) / 2) == P - 1) return -1;
if (Pow(x, (P - 1) / 2) == 0) return 0;
ll a;
Complex k(0, 1);
while (true) {
a = rand();
if (Pow((a * a - x + P) % P, (P - 1) / 2) == P - 1) break;
}
w = (a * a - x + P) % P;
return Pow(k + a, (P + 1) / 2);
}