UVA - 11582 Colossal Fibonacci Numbers!（取模）

题意：输入两个非负整数a，b（0 <= a,b < 2^64）和 n（1<= n < 1000)，输出f（a^b) % n，f为斐波那契数列。

思路：所有计算都是对n取模的，不妨设F(i) = f(i) % n，当f(i - 1)，f(i)出现重复时，整个序列就开始重复。余数最多有n种，所以最多n*n项就会出现重复，所以先计算出周期，再用快速幂求出a^b%M即可。

#include <stdio.h>
#include <bits/stdc++.h>
using namespace std;
typedef unsigned long long ull;
int f[1007 * 1007];
int q_pow(ull a,ull b,int Mod)
{
    ull res,base;
    res = 1;
    base = a;
    while(b) {
        if(b & 1) res = (res * base) % Mod;
        base = (base * base) % Mod;
        b >>= 1;
    }
    return res % Mod;
}
int main(void)
{
    int T,n,M,index,i;
    ull a,b;
    scanf("%d",&T);
    while(T--) {
        scanf("%llu %llu %d",&a,&b,&n);
        if(a == 0 || n == 1) printf("%d\n",0);
        else {
            f[0] = 0,f[1] = 1;
            for(i = 2; i <= n * n + 10; i++) {
                f[i] = (f[i - 2] + f[i - 1]) % n;
                if(f[i - 1] == f[0] && f[i] == f[1]) {
                    M = i - 1;
                    break;
                }
            }
            index = q_pow(a % M,b,M);
            printf("%d\n",f[index]);
        }
    }
    return 0;
}