题目大意
给出T组数据,输入两个非负整数a, b 和正整数n(0 <= a, b < 2^64,1 <= n <= 1000) , 求f(a ^ b) % n, 其中钦定f(0) = 0, f(1) = 1, f(i + 2) = f(i) + f(i + 1) (i 为非负整数), %为mod。
解析
一道数论题:
会发现整个数列取模之后有循环节
于是我们通过找到那个循环节来简化运算
然后熟悉一下“幂取模”的模板还有模运算的一些基本性质
规律最大为n*n,实际没那么大,时间大约为100ms
#include<queue>
#include<cstdio>
#include<cstring>
#include<iostream>
#include<algorithm>
typedef unsigned long long ll;
using namespace std;
const int maxn=1000+10;
ll a,b;
int f[maxn*maxn],n,m;
int power(ll a,ll p)
{
if(p==1)
return a%m;
if(p==0)
return 1;
if(p==2){
return a*a%m;
}
return p%2==1?a*power(a*a%m,p/2)%m:power(a*a%m,p/2)%m;
}
inline void solve()
{
cin>>a>>b>>n;
if(n==1||!a){printf("0\n");return ;}
f[1]=1,f[2]=1;
for(int i=3;i<=n*n+10;i++)
{
f[i]=f[i-1]+f[i-2];
f[i]%=n;
if(f[i]==f[2]&&f[i-1]==f[1])
{
m=i-2;
break;
}
}
int k=power(a%m,b);
printf("%d\n",f[k]);
}
int main()
{
int n;
cin>>n;
while(n--)
solve();
return 0;
}