数论二次剩余模板

最新推荐文章于 2019-09-02 15:31:01 发布

henucm

最新推荐文章于 2019-09-02 15:31:01 发布

阅读量441

点赞数

分类专栏：数论模板

本文链接：https://blog.csdn.net/henucm/article/details/100009964

版权

数论同时被 2 个专栏收录

51 篇文章 0 订阅

订阅专栏

模板

9 篇文章 0 订阅

订阅专栏

本文介绍了数论中的二次剩余概念，详细阐述了如何判断一个数在模P意义下是否为二次剩余，并附带相关证明过程。

摘要由CSDN通过智能技术生成

来自憨批的博客 ( 2k图片加载有点慢需要耐心等待 ) 附带证明

对于 P,n 如果存在x使得 x^2 ≡ n mod P 则称n是模P意义下的二次剩余

struct T{
    long long p, d;
};

long long Ksm(long long a, long long b, long long p) {
    long long res = 1;
    while(b) {
        if(b & 1) res = res * a % p;
        a = a * a % p;
        b >>= 1;
    }
    return res;
}

long long w;
//二次域乘法
T Mul_er(T a, T b, long long p) {
    T ans;
    ans.p = (a.p * b.p + a.d * b.d % p * w % p) % p;
    ans.d = (a.p * b.d % p + a.d * b.p % p) % p;
    return ans;
}
//二次域快速幂
T Ksm_er(T a, long long b, long long p) {
    T ans;
    ans.p = 1; ans.d = 0;
    while(b) {
        if(b & 1) ans = Mul_er(ans, a, p);
        a = Mul_er(a, a, p);
        b >>= 1;
    }
    return ans;
} 
//求勒让德符号
long long Legendre(long long a, long long p) {
    return Ksm(a, (p-1)>>1, p);
}

long long Recever(long long a, long long p) {
    a %= p;
    if(a < 0) a += p;
    return a;
}

long long solve(long long n, long long p) { //求解 x^2=n(mod p) 的x的值
    if(n % p == 0) return 0;
    if(p == 2) return 1;
    if(Legendre(n, p) + 1 == p) return -1;
    long long a = -1, t;
    while(1) {
        a = rand() % p;
        t = a * a - n;
        w = Recever(t, p);
        if(Legendre(w, p) + 1 == p) break;
    }
    T tmp;
    tmp.p = a; tmp.d = 1;
    T ans = Ksm_er(tmp, (p+1)>>1, p);
    return ans.p;
}
ll x = solve(n, p);
x==-1无解，否则
x和p-x都是合法解