hdu1757 A Simple Math Problem（矩阵快速幂）

题目点我点我点我

解题思路：

（原谅我盗图了，懒得弄这个）

矩阵快速幂直接干搞，最后Fn就是第一个矩阵的第一行乘以第二个矩阵。

最后提醒自己注意矩阵的初始化。

/* ***********************************************
┆  ┏┓　　　┏┓ ┆
┆┏┛┻━━━┛┻┓ ┆
┆┃　　　　　　　┃ ┆
┆┃　　　━　　　┃ ┆
┆┃　┳┛　┗┳　┃ ┆
┆┃　　　　　　　┃ ┆
┆┃　　　┻　　　┃ ┆
┆┗━┓　马　┏━┛ ┆
┆　　┃　勒　┃　　┆　　　　　　
┆　　┃　戈　┗━━━┓ ┆
┆　　┃　壁　　　　　┣┓┆
┆　　┃　的草泥马　　┏┛┆
┆　　┗┓┓┏━┳┓┏┛ ┆
┆　　　┃┫┫　┃┫┫ ┆
┆　　　┗┻┛　┗┻┛ ┆
************************************************ */

#include <stdio.h>
#include <string.h>
#include <iostream>
#include <algorithm>
#include <vector>
#include <queue>
#include <stack>
#include <set>
#include <map>
#include <string>
#include <math.h>
#include <stdlib.h>
using namespace std;

#define rep(i,a,b) for (int i=(a),_ed=(b);i<=_ed;i++)
#define per(i,a,b) for (int i=(b),_ed=(a);i>=_ed;i--)
const int inf_int = 2e9;
const long long inf_ll = 2e18;
#define ll long long
#define ull unsigned long long

#define maxn 15
int n = 10;
int mod;

struct Mat
{
    int mat[maxn][maxn];
};

Mat operator * (Mat a,Mat b)
{
    Mat c;
    memset(c.mat,0,sizeof c.mat);
    for(int k=0;k<n;k++)
    {
        for(int i=0;i<n;i++)
        {
            for(int j=0;j<n;j++)
            {
                c.mat[i][j] += ((a.mat[i][k] * b.mat[k][j]) % mod);
                c.mat[i][j] %= mod;
            }
        }
    }
    return c;
}

Mat operator ^ (Mat a,int k)
{
    Mat c;
    memset(c.mat,0,sizeof c.mat);
    for(int i=0;i<n;i++)c.mat[i][i] = 1;
    while(k)
    {
        if(k&1)c = c * a;
        a = a * a;
        k >>= 1;
    }
    return c;
}


int main()
{
    //freopen("in.txt","r",stdin);
    //freopen("out.txt","w",stdout);
	int k;
	while(~scanf("%d%d",&k,&mod))
    {
        Mat a;
        memset(a.mat,0,sizeof a.mat);
        for(int i=0;i<n;i++)
            scanf("%d",&a.mat[0][i]);
        if(k<10)
        {
            cout<<k%mod<<endl;
            continue;
        }
        for(int i=1;i<n;i++)
        {
            for(int j=0;j<n;j++)
            {
                if(i==(j+1))
                {
                    a.mat[i][j] = 1;
                }
            }
        }
        Mat res;
        res = a^(k-9);
        int ans=0;
        for(int i=0;i<n;i++)
        {
            ans += ((res.mat[0][i]*(n-i-1))%mod);
            ans %= mod;
        }
        cout<<ans<<endl;

    }
    return 0;
}