HDU 1757 A Simple Math Problem（矩阵快速幂模板）

题意：
题意很简单，不多说了。

思路：

|f(10) |       |a0 a1 a2 ...a8 a9|    |f(9)|
| f(9)  |       | 1   0   0 ... 0    0 |    |f(8)|
| .....  |   =  | ..    ...    ...   ...    |     | ..   |
| f(2) |        | 0   0   0 ... 0    0|     |f(1)|
| f(1) |　　 | 0   0   0 ... 1    0|     |f(0)|

为自己留个模板。

 #include<iostream>
 #include<algorithm>
 #include<cstring>
 #include<cstdio>
 #include<sstream>
 #include<vector>
 #include<stack>
 #include<queue>
 #include<cmath>
 #include<map>
 #include<set>
 using namespace std;
 typedef long long ll;
 typedef pair<int,int> pll;
 const int INF = 0x3f3f3f3f;
 const int maxn =  + ;

 int n,m;

 struct Matrix
 {
     int mat[][];
     Matrix operator *(Matrix a)
     {
         Matrix c;
         for(int i=;i<;i++)
         {
             for(int j=;j<;j++)
             {
                 c.mat[i][j]=;
                 for(int k=;k<;k++)
                 {
                     c.mat[i][j]+=(mat[i][k]*a.mat[k][j])%m;
                     c.mat[i][j]%=m;
                 }
             }
         }
         return c;
     }
 };

 Matrix base, ans;

 void init()
 {
     for(int i=;i<;i++)
     for(int j=;j<;j++)
     {
         if(i==j) ans.mat[i][j]=;
         else ans.mat[i][j]=;
         if(i-==j) base.mat[i][j]=;
         else base.mat[i][j]=;
     }
 }

 void pow(int n)
 {
     while(n)
     {
         if(n&)  ans=ans*base;
         base=base*base;
         n>>=;
     }
 }

 int main()
 {
     //freopen("in.txt","r",stdin);
     while(~scanf("%d%d",&n,&m))
     {
         init();
         for(int i=;i<;i++)  scanf("%d",&base.mat[][i]);
         if(n<)  {printf("%d\n",n%m);continue;}
         pow(n-);
         int sum=;
         for(int i=;i<;i++)  sum=(sum+ans.mat[][i]*(-i))%m;
         printf("%d\n",sum);
     }
     return ;
 }