解题思路:
(原谅我盗图了,懒得弄这个)
矩阵快速幂直接干搞,最后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;
}