HDU5478 Can you find it【快速幂】

版权声明：转载什么的好说，附上友链就ojek了，同为咸鱼，一起学习。 https://blog.csdn.net/sodacoco/article/details/86503823

参看资料：

https://blog.csdn.net/lanshan1111/article/details/86499454

https://blog.csdn.net/queuelovestack/article/details/48754331

---

题目：

【害怕公式变形直接贴了图片】

题目大意：

       给定 k1,b1,k2,C，求满足公式的a,b值；

解题思路：

       来自大佬的讲解，

       得到最后的式子后，因为要将(a,b)按字典序输出，所以for循环将a从1开始循环至C，对于每一个a的值用快速幂求出

       a^(k1+b1)，然后由①式得到b，在用快速幂算出a^k1 和 b^k2以判断⑤式是否成立，成立输出即可，数据用long long类型。

实现代码：

#include<stdio.h>
#include<stdlib.h>
#include<math.h>
#include<cmath>
#include<algorithm>
#include<iostream>
#define exp 1e-10
using namespace std;
const int N = 200001;
const int inf = 1000000000;
const int mod = 2009;
__int64 Quick_Mod(int a, int b, int m){ //快速幂
    __int64 res = 1,term = a % m;
    while(b){
        if(b & 1) res = (res * term) % m;
        term = (term * term) % m;
        b >>= 1;
    }
    return res%m;
}
int main(){
    int C,k1,b1,k2,i,k=1;
    __int64 a,b,s;
    bool flag;
    while(~scanf("%d%d%d%d",&C,&k1,&b1,&k2)){
        flag=false;
        printf("Case #%d:\n",k++);
        for(i=1;i<C;i++){
            a=Quick_Mod(i,k1,C);
            b=C-Quick_Mod(i,k1+b1,C);
            s=Quick_Mod(b,k2,C);
            //printf("##%I64d %I64d\n",a,s);
            if(a==s)
                flag=true,printf("%d %I64d\n",i,b);
        }
        if(!flag)
            puts("-1");
    }
    return 0;
}