UVA-11582 巨大的斐波那契数！(快速幂取模+打表)

最新推荐文章于 2021-10-20 01:34:15 发布

原创最新推荐文章于 2021-10-20 01:34:15 发布 · 307 阅读

0 ·

本内容遵循CC 4.0 BY-SA版权协议

标签

#数论

数论专栏收录该内容

9 篇文章

订阅专栏

博客围绕计算f(i)除以n的余数展开。先给出输入要求，即两个非负整数a、b和正整数n。解题思路是对n取模，利用二元组重复确定周期M。同时提醒注意使用unsigned long long类型、打表和快速幂计算、考虑n=1或a=0的特殊情况，最后给出代码。

题意：

输入两个非负整数a、b和正整数n(0<=a,b< $2^{64}$ ,1<=n<=1000),你的任务是计算f( $a^{b}$ )除以n的余数。其中f(0)=f(1),且对于所有非负整数i，f(i+2)=f(i+1)+f(i)。

思路：

所有计算都是对n取模的，不妨设F(i)=f(i) mod n。不难发现，当二元组(F(i),F(i+1))出现重复时，整个序列就开始重复。例如，n=3，序列F(i)的前10项为1，1，2，0，2，2，1，0，1，1，第9，10项和前两项完全一样。

那么，多久会出现重复呢？因为余数最多n种，所以最多 $n^{2}$ 项就会出现重复。设周期为M，则只需计算F(0)~F( $n^{2}$ )直到找到周期M。

另外，有几点要注意！！！

1、题目要求a,b在0~ $2^{64}$ -1，而long long的范围是- $2^{63}$ ~ $2^{63}$ -1，unsigned long long的范围是0~ $2^{64}$ -1，所以此题要用unsigned long long！！！(对应的printf要用%llu)

2、时间复杂度的要求：此题要打表，先计算出n从1到1000时的F(i)和对应的周期M;计算 $a^{b}$ 时用快速幂。

3、考虑特殊情况，即n=1或a=0时直接输出答案0。

代码：

#include<bits/stdc++.h>
using namespace std;
typedef unsigned long long ll;
int t;
ll a,b,m[1010],n,f[1010][50010];
ll qsm(ll m1,ll m2,ll mod){
	ll res=1;
	while(m2){
		if(m2&1)
			res=(res*m1)%mod;
		m2>>=1;
		m1=(m1*m1)%mod;
	}
	return res%mod;
}
int main(){
	for(int k=2;k<=1000;k++){ 
		f[k][0]=0;
		f[k][1]=1%k;
		m[k]=0;
		for(int i=2;i<=k*k;i++){
			f[k][i]=(f[k][i-1]+f[k][i-2])%k;
			if(f[k][i-1]==f[k][0]&&f[k][i]==f[k][1]){
				m[k]=i-1;
				break;
			}
		}
	} 
	scanf("%d",&t);
	while(t--){
		scanf("%llu%llu%llu",&a,&b,&n);
		if(n==1||a==0)
			printf("0\n");
		else{
			ll c=qsm(a%m[n],b,m[n]);
			printf("%llu\n",f[n][c]);
		}
	}
	return 0;
}