HDU 6470 【矩阵快速幂】

题目链接：http://acm.hdu.edu.cn/showproblem.php?pid=6470

写这道题是为了让自己不要忘记矩阵快速幂如何推出矩阵式子的。

注意 代码是TLE的！！

 #include<stdio.h>
 #include<string.h>
 #define mem(a, b) memset(a, b, sizeof(a))
 typedef long long ll;
 const int mod = ;

 ll n;

 struct Matrix
 {
     ll a[][];
 }A, res, temp;

 Matrix Multiply(Matrix a, Matrix b)
 {
     Matrix ans;
     mem(ans.a, );
     for(int i = ; i <= ; i ++)
         for(int j = ; j <= ; j ++)
             for(int k = ; k <= ; k ++)
             {
                 ans.a[i][j] += a.a[i][k] * b.a[k][j] % mod;
                 ans.a[i][j] %= mod;
             }
     return ans;
 }

 int main()
 {
     mem(A.a, );  //构造矩阵
     A.a[][] = A.a[][] = A.a[][] = A.a[][] = A.a[][] = A.a[][] = A.a[][] = A.a[][] = A.a[][] = A.a[][] = ;
     A.a[][] = A.a[][] = ;
     A.a[][] = A.a[][] = ;
     int T;
     scanf("%d", &T);
     while(T --)
     {
         temp = A;  //初始化temp矩阵
         mem(res.a, ); //初始化为单位矩阵
         for(int i = ; i <= ; i ++)
             res.a[i][i] = ;
         scanf("%lld", &n);
         n -= ;
         while(n)
         {
             if(n % )
                 res = Multiply(res, temp);
             temp = Multiply(temp, temp);
             n /= ;
         }
         printf("%lld\n", ( * res.a[][] % mod +  * res.a[][] % mod +  * res.a[][] % mod +  * res.a[][] % mod +  * res.a[][] % mod + res.a[][] % mod) % mod);
     }
     return ;
 }