二次剩余（简单理解）

二次剩余

定义

$p$是一个奇素数。
$$x^2\equiv n(modp)$$
判断n是否有解勒让德符号 $\equiv$ 欧拉判别准则。

勒让德符号
$$\left(\frac{n}{p}\right)=\{\begin{array}{ll}1& p\nmid n\wedge n是p的二次剩余\\ -1& p\nmid n\wedge n不是p的二次剩余\\ 0& p\mid n\end{array}$$
在区间$[1,p-1]$中存在解的数称为$p$的二次剩余，有$\frac{p-1}{2}$个。

欧拉判别式：
$$\left(\frac{n}{p}\right)\equiv n^{\frac{(p-1)}{2}}(modp)$$
简单理解：
$$n^{\frac{(p-1)}{2}}\equiv (n^{(p-1)}{)}^{\frac{1}{2}}(modp)\equiv \sqrt{1}\equiv \pm 1(modp)$$
等于-1时无解，所以此时n不是p的二次剩余

Cipolla 算法(求解算法)

理解类似复数域的东西。

$i=\sqrt{-1}$，而$n$可能在模$p$的情况下得不到整数解，所以我们也弄一个复数域。

就像是当$n$不是$p$的二次剩余的时候，我们设$w^2\equiv n$。

当n是p的二次剩余时，我们找到一个数$a$使$(a^2-n)$不是p的二次剩余，$w^2\equiv (a^2-n)$。则$x^2\equiv n(modp)$的解是$(a-w{)}^{\frac{(p+1)}{2}}$

证明：

先证明一些小定理：

定理1：$(a+b{)}^p\equiv a^p+b^p$
$$\begin{array}{rl}(a+b{)}^p& \equiv \sum _{i=0}^p{C}_p^i{a}^{p-i}{b}^i\\ & \equiv \sum _{i=0}^p\frac{p!}{i!(p-i!)}{a}^{p-i}{b}^i\\ & \equiv a^p+b^p(modp)\end{array}$$
​ 定理2：$w^p\equiv -w(modp)$
$$\begin{array}{rl}{w}^p& \equiv w^{p-1}w\\ & \equiv (w^2{)}^{\frac{p-1}{2}}w\\ & \equiv (n{)}^{\frac{p-1}{2}}w\\ & \equiv -w\end{array}$$
证明：$x^2\equiv n(modp)$的解是$(a-w{)}^{\frac{(p+1)}{2}}$，设$w^2=(a^2-n)$
$$\begin{array}{rl}x& \equiv (a-w{)}^{\frac{p+1}{2}}\\ & \equiv ((a-w{)}^{p+1}{)}^{\frac{1}{2}}\\ & \equiv ((a-w{)}^p(a-w){)}^{\frac{1}{2}}\\ & \equiv ((a^p-w^p)(a-w){)}^{\frac{1}{2}}\\ & \equiv ((a+w)(a-w){)}^{\frac{1}{2}}\\ & \equiv (a^2-w^2{)}^{\frac{1}{2}}\\ & \equiv (a^2-(a^2-n){)}^{\frac{1}{2}}\\ & \equiv n^{\frac{1}{2}}(modp)\end{array}$$
$\therefore x\equiv (a+w{)}^{\frac{(p+1)}{2}}\equiv (modp)$

模板题

【模板】二次剩余

#include<bits/stdc++.h>

using namespace std;
mt19937_64 gen(time(0));
typedef long long ll;

ll n, mod, w;

ll qpow(ll x, ll y) {
    ll ans = 1;
    while(y) {
        if(y & 1) ans = ans * x % mod;
        x = x * x % mod;
        y >>= 1;
    }
    return ans;
}

struct Y {
    ll x, y;
    Y operator * (const Y & a) {
        Y tmp;
        tmp.x = (x * a.x % mod + y * a.y % mod * w % mod) % mod;
        tmp.y = (x * a.y % mod + y * a.x % mod) % mod;
        return tmp; 
    }
    Y operator + (const Y & a) {
        Y tmp;
        tmp.x = (x + a.x) % mod;
        tmp.y = (y + a.y) % mod;
        return tmp;
    }
};

Y ksm(Y a, ll x) {
    Y ans = {1, 0};
    while(x) {
        if(x & 1) ans = ans * a;
        a = a * a;
        x >>= 1;
    }
    return ans;
}

int main() {
#ifndef ONLINE_JUDGE
    freopen("in.txt", "r", stdin);
    freopen("out.txt", "w", stdout);
#endif
    int t;
    scanf("%d", &t);
    while(t--) {
        scanf("%lld%lld", &n, &mod);
        n %= mod;
        if(n == 0) printf("0\n");
        else if(qpow(n, (mod-1)/2) == mod-1) puts("Hola!");
        else {
            ll tmp;
            while(1) {
                tmp = gen() % mod;
                w = (tmp * tmp % mod - n + mod) % mod;
                if(qpow(w, (mod-1)/2) == mod-1) break;
            }
            Y a = {tmp, 1};
            a = ksm(a, (mod+1)/2);
            tmp = mod - a.x;
            if(tmp == a.x) cout << a.x << endl;
            else cout << min(tmp, a.x) << ' ' << max(tmp, a.x) << endl;
        }
    }
    return 0;
}