hdoj 5478 Can you find it

最新推荐文章于 2021-06-08 10:58:09 发布

squee_spoon

最新推荐文章于 2021-06-08 10:58:09 发布

阅读量427

点赞数

CC 4.0 BY-SA版权

分类专栏：数学相关文章标签： 2015ICPC

本文链接：https://blog.youkuaiyun.com/squee_spoon/article/details/48754495

数学相关专栏收录该内容

69 篇文章

订阅专栏

本文介绍了一种针对特定模方程求解的算法实现，该算法通过枚举参数a并利用快速幂运算来寻找满足特定模方程的整数解对(a, b)，适用于给定的大素数C及三个整数k1, b1, k2的情况下。

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

Problem Description

Given a prime number

C(1≤C≤2×105), and three integers k1, b1, k2

(1≤k1,k2,b1≤109). Please find all pairs (a, b) which satisfied the equation

ak1⋅n+b1+

bk2⋅n−k2+1 = 0 (mod C)(n = 1, 2, 3, ...).

枚举a，把n=1代入，算出b。根据a，b，代入n=2，如果等式成立，就输出。

#include <iostream>     
using namespace std;     

#define ll long long

const int maxn = 100010;

ll c,k1,B1,k2;

ll p;

ll Quick_Pow(ll a,ll n) {
    ll ret=1;
    ll temp=a%p;
    while (n){
        if (n&1) ret=ret*temp%p;
        temp=temp*temp%p;
        n>>=1;
    }
    return ret;
}

int main(){
    int cas = 0;
    while(cin>>c>>k1>>B1>>k2){
        cas++;
        printf("Case #%d:\n",cas);
        bool hassol=0;
        for(int a=1;a<c;a++){
            //a^(k1*n+b1)
            //b^(k2*n-k2+1)
            p=c;
            ll a1 = Quick_Pow(a,k1+B1);
            ll b = c - a1;
            
            bool ok=1;
            ll an=a1;
            ll bn=b;
            ll atmp = Quick_Pow(a,k1);
            ll btmp = Quick_Pow(b,k2);
            an*=atmp;	bn*=btmp;
            an%=c;		bn%=c;
            if( (an+bn)%c==0 ){
                hassol=1;
                printf("%I64d %I64d\n",a,b);
            }
        }
        if(!hassol){
            printf("-1\n");
        }
    }
    return 0;
}