题意:输入两个非负整数a,b(0 <= a,b < 2^64)和 n(1<= n < 1000),输出f(a^b) % n,f为斐波那契数列。
思路:所有计算都是对n取模的,不妨设F(i) = f(i) % n,当f(i - 1),f(i)出现重复时,整个序列就开始重复。余数最多有n种,所以最多n*n项就会出现重复,所以先计算出周期,再用快速幂求出a^b%M即可。
#include <stdio.h>
#include <bits/stdc++.h>
using namespace std;
typedef unsigned long long ull;
int f[1007 * 1007];
int q_pow(ull a,ull b,int Mod)
{
ull res,base;
res = 1;
base = a;
while(b) {
if(b & 1) res = (res * base) % Mod;
base = (base * base) % Mod;
b >>= 1;
}
return res % Mod;
}
int main(void)
{
int T,n,M,index,i;
ull a,b;
scanf("%d",&T);
while(T--) {
scanf("%llu %llu %d",&a,&b,&n);
if(a == 0 || n == 1) printf("%d\n",0);
else {
f[0] = 0,f[1] = 1;
for(i = 2; i <= n * n + 10; i++) {
f[i] = (f[i - 2] + f[i - 1]) % n;
if(f[i - 1] == f[0] && f[i] == f[1]) {
M = i - 1;
break;
}
}
index = q_pow(a % M,b,M);
printf("%d\n",f[index]);
}
}
return 0;
}