2019牛客暑期多校训练营(第九场)B:Quadratic equation(二次剩余)

本文探讨了在模p=1e9+7条件下,构造满足特定模运算的x和y值,使得(x+y)modp=b和(x*y)modp=c成立。通过分析二次剩余性质,给出了一种有效的求解策略,并提供了C++实现代码。

摘要生成于 C知道 ,由 DeepSeek-R1 满血版支持, 前往体验 >

【题解】

题意:给定p=1e9+7,构造x,y使其满足(x+y) mod p = b,(x*y) mod p = c .

思路:不考虑取模的情况下,(x-y)^2=(x+y)^2-4xy=b^2-4c 。在取模的意义下,(x-y)^2\equiv b^2-4c (mod p),因为a是模p的二次剩余的充分必要条件为 a^{(p-1)/2}\equiv 1(mod p) ,所以可以根据二次剩余求出x-y。

【代码】

#include<bits/stdc++.h>
using namespace std;
#define ll long long
const int mod=1e9+7,inv2=(mod+1)/2;
int Pow(int x,int n)
{
    int k=1;
    while(n){
        if(n&1) k=(ll)k*x%mod;
        x=(ll)x*x%mod;
        n>>=1;
    }
    return k;
}
int main()
{
    int T; scanf("%d",&T);
    while(T--){
        int a,b; scanf("%d%d",&a,&b);
        int k=(((ll)a*a-4ll*b)%mod+mod)%mod;
        int v=Pow(k,(mod+1)/4); //根据二次剩余求出x-y
        if((ll)v*v%mod!=k){
            puts("-1 -1");
            continue;
        }
        int x=(ll)(v+a)*inv2%mod;
        int y=(a-x+mod)%mod;
        if(x>y) printf("%d %d\n",y,x);
        else printf("%d %d\n",x,y);
    }
}

 

评论
添加红包

请填写红包祝福语或标题

红包个数最小为10个

红包金额最低5元

当前余额3.43前往充值 >
需支付:10.00
成就一亿技术人!
领取后你会自动成为博主和红包主的粉丝 规则
hope_wisdom
发出的红包
实付
使用余额支付
点击重新获取
扫码支付
钱包余额 0

抵扣说明:

1.余额是钱包充值的虚拟货币,按照1:1的比例进行支付金额的抵扣。
2.余额无法直接购买下载,可以购买VIP、付费专栏及课程。

余额充值