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 = 1000 + 5;

int n,m;

struct Matrix
{
    int mat[12][12];
    Matrix operator *(Matrix a)
    {
        Matrix c;
        for(int i=0;i<10;i++)
        {
            for(int j=0;j<10;j++)
            {
                c.mat[i][j]=0;
                for(int k=0;k<10;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=0;i<10;i++)
    for(int j=0;j<10;j++)
    {
        if(i==j) ans.mat[i][j]=1;
        else ans.mat[i][j]=0;
        if(i-1==j) base.mat[i][j]=1;
        else base.mat[i][j]=0;
    }
}

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

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