Uva 11582(Fibonacci循环节)

点击打开链接

求fib[a^b]%n (a,b<2^64)

构造F[i]=fib[i]%n -> F[i]=(fib[i-2]+f[i-1])%n=(fib[i-2]%n+fib[i-1]%n)%n=(F[i-2]+F[i-1])%n

如果F[x-1]=F[i-1] &&F[x-2]=F[i-2] 则F[x]=F[i] F[x+1]=(F[x]+F[x-1)%n=F[i+1].....则找到循环节长度为loop=x-2-(i-1)+1

因为余数有n种 n*n个二元组中肯定有两个是相同的. 找到loop. fib[a^b]%n=fib[(a^b)%loop]%n

 

#include <iostream>
#include <cstring>
#include <algorithm>
using namespace std;
typedef unsigned long long ll;
const int N=1e6+100;
ll f[N];//fib[i]
ll F[N];// F[i]=(f[i])%n
ll pow_mod(ll a,ll b,ll mod)
{
	if(b==0) return 1;
	ll res=pow_mod(a,b/2,mod);
	res=res*res%mod;
	if(b%2)
	res=res*a%mod;
	
	return res;	
}
int main()
{
	ll a,b;
	int n;
	int t;
	cin>>t;
	while(t--)
	{
		cin>>a>>b>>n;//求f(a^b)%n
		//F[i]=f[i]%n
		// F[i]=(f[i-1]+f[i-2])%n=(f[i-1]%n+f[i-2]%n)%n=(F[i-1]+F[i-2])%n 
		// 当二元组(F[i-2],F[i-1])= (F[x-1],F[x-2])时 循环节产生: F[x]=F[i] F[x+1]=F[i+1].... 
		ll loop=1;
		F[0]=0;
		F[1]=1%n;
		for(int i=2;i<=n*n+100;i++)//余数有n种 则在n*n个二元组中肯的存在两个相等二元组->loop<=n*n 
		{
			F[i]=(F[i-1]+F[i-2])%n;	
			if(F[i]==F[1]&&F[i-1]==F[0])//找到循环周期为i-1 
			{
				loop=i-1;
				break;	
			}
		}
		//循环长度为loop  f[a^b]%n= f[(a^b)%loop] %n
		ll ans=pow_mod(a%loop,b,loop); 
		cout<<F[ans]<<endl;		
	}
	return 0;
}