简单的用矩阵快速幂求递推式
#include<iostream>
#include<cstdio>
#include<cstring>
#include<algorithm>
#include<cmath>
using namespace std;
typedef long long ll;
ll MOD;
int N;
struct node
{
ll a[20][20];
};
node shu,ans,mp;
//shu是输入的矩阵,ans是所求答案
node matrix(node x,node y)
{
for(int i=1;i<=N;i++)
for(int j=1;j<=N;j++){
mp.a[i][j]=0;
for(int p=1;p<=N;p++)
mp.a[i][j]=(mp.a[i][j]+x.a[i][p]*y.a[p][j]%MOD+MOD)%MOD;
//矩阵乘法
}
return mp;
}
void work(ll k)
{//矩阵快速幂
for(int i=1;i<=N;i++)
for(int j=1;j<=N;j++)
ans.a[i][j]=0;
for(int i=1;i<=N;i++) ans.a[i][i]=1;
node t=shu;
while(k){
if(k&1)
ans=matrix(ans,t);
k>>=1;
t=matrix(t,t);
}
}
ll a[20];
int main()
{
ll n,d,m;
while(~scanf("%lld%lld%lld",&d,&n,&m))
{
if(n==0&&m==0&&d==0) break;
N=d;MOD=m;
memset(shu.a,0,sizeof(shu.a));
for(int i=1;i<=d;i++)
scanf("%lld",&shu.a[1][i]);
for(int i=2;i<=d;i++)
shu.a[i][i-1]=1;
for(int i=1;i<=d;i++)
scanf("%lld",&a[i]);
if(n<=d)
{
printf("%lld\n",a[n]%MOD);
continue;
}
work(n-d);
ll ret=0;
for(int i=1;i<=d;i++)
ret=(ret+ans.a[1][i]*a[d-i+1])%MOD;
printf("%lld\n",ret);
}
return 0;
}