矩阵快速幂uva10870

详见训练指南

#include<bits/stdc++.h>
using namespace std;
const int maxn=20;
int d,n,m;
int a[maxn],f[maxn];
struct Matrix
{
    long long mat[maxn][maxn];
    Matrix operator*(Matrix a)
    {
        Matrix res;
        for(int i=0;i<d;i++)
        {
            for(int j=0;j<d;j++)
            {
                res.mat[i][j]=0;
                for(int k=0;k<d;k++)
                {
                    res.mat[i][j]+=mat[i][k]*a.mat[k][j];
                    res.mat[i][j]%=m;
                }
            }
        }
        return res;
    }
};

Matrix matrix_pow(Matrix &A,int x)
{
    Matrix res;
    memset(res.mat,0,sizeof(res.mat));
    for(int i=0;i<d;i++)
        res.mat[i][i]=1;
    while(x)
    {
        if(x&1)res=res*A;
        A=A*A;
        x>>=1;
    }
    return res;
}
int main()
{
    while(scanf("%d%d%d",&d,&n,&m)!=EOF,d+n+m)
    {
        for(int i=d-1;i>=0;i--)scanf("%d",&a[i]);
        for(int i=0;i<d;i++)scanf("%d",&f[i]);
        Matrix A,B;
        memset(A.mat,0,sizeof(A.mat));
        for(int i=0;i<d;i++)A.mat[d-1][i]=a[i]%m;
        for(int i=0;i<d-1;i++)A.mat[i][i+1]=1;
        A=matrix_pow(A,n-d);
        for(int i=0;i<d;i++)B.mat[i][0]=f[i]%m;
        B=A*B;
        printf("%d\n",B.mat[d-1][0]%m);
    }
    return 0;
}