Uva 11582 Colossal Fibonacci Numbers!

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存起来。
运行的结果正确：