HDU 4291 A Short problem

题意： 求 g(g(g(n))) mod 10⁹+ 7  ，其中g(n)=3*g(n-1)+g(n-2),  g(0)=0,g(1)=1;

分析： 当n 的3的时候内部就超精度了，需要本地暴力找出每层的循环节，最外面的对1000000007取模，找最外层循环节222222224， 在对222222224取模，

         找到次一层循环节183120，这样就可以通过三次快速幂求出结果了

         g(g(g(n)))%1000000007 = g(g(g(n)%183120)%222222224)%1000000007

#include<stdio.h>
#include<string.h>
void mul(__int64 a[][3],__int64 b[][3],__int64 mod)
{
    __int64 c[3][3];
    int i,j,k;
    for(i=1;i<=2;i++)
        for(j=1;j<=2;j++)
        {
            c[i][j]=0;
            for(k=1;k<=2;k++)
                c[i][j]+=((a[i][k]%mod)*(b[k][j]%mod))%mod;
        }
    for(i=1;i<=2;i++)
        for(j=1;j<=2;j++)
            b[i][j]=c[i][j]%mod;
}
__int64 cal(__int64 n,__int64 mod)
{
        if(n==0)
            return 0;
        __int64 a[3][3],b[3][3];
        a[1][1]=3; 
        a[1][2]=1;
        a[2][1]=1; 
        a[2][2]=0;
        b[1][1]=1;
        b[1][2]=b[2][1]=b[2][2]=0;
        n-=1;
        while(n)
        {
            if(n&1)
                mul(a,b,mod);
            mul(a,a,mod);
            n/=2;
        }
        return b[1][1];
}
int main()
{
    __int64 n;
    int i;
    while(scanf("%I64d",&n)!=EOF)
    {
        __int64 t1=cal(n,183120);
        __int64 t2=cal(t1,222222224);
        __int64 t3=cal(t2,1000000007);
        printf("%I64d\n",t3);
    }
    return 0;
}