Uva 11582 Colossal Fibonacci Numbers!
The i ‘th Fibonacci number f(i) is recursively defined in the following way:
• f(0) = 0 and f(1) = 1
• f(i + 2) = f(i + 1) + f(i) for every i ≥ 0
Your task is to compute some values of this sequence.
Input
Input begins with an integer t ≤ 10, 000, the number of test cases. Each test case consists of three integers a, b, n where 0 ≤ a, b < 2^64 (a and b will not both be zero) and 1 ≤ n ≤ 1000.
Output
For each test case, output a single line containing the remainder of f(a^b) upon division by n.
Sample Input
3
1 1 2
2 3 1000
18446744073709551615 18446744073709551615 1000
Sample Output
1
21
250
乍一看感觉此题没什么特别的难度,直到看到了本题样例中的
f(18446744073709551615^18446744073709551615),有种人生都悲剧了的感觉。
不过这样一看啊,出这么一个数据绝对是有什么鬼的,对于这种数列题目,我一般选择用Excel来试着计算一下。
A列是fibonacci数列,C列是对A取模的数,B项是取模的结果。
怎么样,是不是看出来了一点影子?B列的数是有循环的规律的!
可以试着修改一下C列的值(其实可以不需要C列,写mod的时候第二个参数直接填数即可,不过本人有意写成这样,便于观察和修改),发现不管是对几取模,结果总是以循环数的形式呈现的(当然,当C列数太大,就需要扩充一下A列中Fibonacci的项数咯)。
于是我们可以窥出为何题中n有1~1000的范围限制了,如果范围太大,我们在找循环节的时候(当然要找循环节),就有点尴尬了。
也就是说,这题可以分解为两个步骤:
1.找循环节
2.快速幂
通过寻找项数在循环节中的位置,便可解出答案^ _ ^
My Code
#include<iostream>
#include<cstdio>
using namespace std;
const int N=1001;
int f[N][N]; //f[i][j] 斐波那契数列第i向对j取模
int len[N],n,t; //k[i] n=i时循环节长度(从f[0]开始数)
unsigned long long a,b;
int ksm(unsigned long long a,unsigned long long b,int MODN) //快速幂
{
if (b==0) return 1;
int k=ksm(a,b/2,MODN);
return b%2==0?(k*k)%MODN:(k*k*a)%MODN;
}
int main()
{
f[0][1]=0,len[1]=1; //对n=1特殊处理
for (n=2;n<N;n++) //制表
{
f[0][n]=0,f[1][n]=1;
for (int i=2;i<N;i++)
{
f[i][n]=(f[i-1][n]+f[i-2][n])%n;
if (f[i][n]==0) {
len[n]=i;
break;
}
}
}
scanf("%d",&t);
while (t--)
{
scanf("%llu%llu%d",&a,&b,&n);
printf("%d\n",f[ksm(a%len[n],b,len[n])][n]);
}
return 0;
}
这里还有一个细节,a,b必须要用unsigned long long存起来。
运行的结果正确: