UVA10870矩阵快速幂与递推公式的关系

/*UVA10870：
矩阵快速幂与递推公式的关系：关键就是这句话，以后做题能联想到即可
f(n)=a1f(n-1)+...adf(n-d)
*/
#include<iostream>
#include<stdio.h>
#include<string.h>
#include<algorithm>
#include<stdlib.h>
#include<math.h>
#include<queue>
#include<vector>
#include<map>
#define ULL unsigned long long
#define LL long long
#define Max_size 17

using namespace std;

int d,n,MOD;//
struct Mat
{
    LL mat[Max_size][Max_size];

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

void builtE()
{
    memset(E.mat,0,sizeof(E.mat));
    for (int i=0;i<Max_size;i++)
    {
        E.mat[i][i]=1;
    }
}
Mat operator^(Mat A,int x)
{
    Mat c;
    builtE();
    c=E;
    for ( ; x; x>>=1)
    {
        if ( x & 1)
            c = c * A;
        A = A * A;
    }
    return c;
}
Mat K,F;
int main()
{
    while(cin>>d>>n>>MOD)
    {
        memset(K.mat,0,sizeof(K.mat));
        memset(F.mat,0,sizeof(F.mat));
        for(int i=d-1;i>=0;i--) cin>>K.mat[i][d-1];
        for(int i=0;i<d;i++) K.mat[i+1][i]=1;
        K=K^(n-1);
        for(int i=0;i<d;i++) cin>>F.mat[0][i];
        F=F*K;
        cout<<F.mat[0][0]<<endl;
    }
    return 0;
}