uva 11582 快速幂 Fibonacci循环节

题目：

紫书数论10.1 11582 P316

题解：

先找出Finonacci的循环节，这里注意a^b%n可以优化为(a%n)^b%n，
Unsigned long long 是比long long 在正数范围是前者大的，这里用unsigned long long 可以防止快速幂运算溢出
小于2^64的数要用unsigned long long

代码：

#include <iostream>#include <cstdio>using namespace std;const int maxn = 1000000+100;typedef long long LL;typedef unsigned long long ULL;LL f[maxn]={0,1};LL getT(LL n){    LL m = n*n+2;    int cnt=0;    for(LL i=2;i<=m;i++)    {        f[i]=(f[i-1]+f[i-2])%n;        cnt++;        if(f[i]==1&&f[i-1]==0)        {             return cnt;        }    }}int pow_mod(ULL a, ULL b, ULL m){    if(b==0) return 1;    int k=pow_mod(a, b/2, m);    k=k*k%m;    if(b%2) k=k*a%m;    return k;}int main(){    int T;    cin>>T;    ULL a,b,n;    while(T--)    {        scanf("%llu%llu%llu",&a,&b,&n);        LL round=getT(n);        cout<<round<<endl;        if(n!=1){        int tmp = pow_mod(a%round,b,round);        cout<<f[tmp]<<endl;        }        else            cout<<"0"<<endl;    }    return 0;}