kuangbin专题十九 UVA10870 构造矩阵

题意：
考虑递推关系式f(n)=a1*f(n-1)+a2*f(n-2)+….+ad*f(n-d)，计算f(n)%m
题解:
一般来说，看到一些公式，然后n比较大的时候，矩阵都会排上用场，我在这道题耗费了一下午还要看一下博客，之后发现，我的想法和思路是正确的，但是错就错在我构造错了，接下来让你们看看我具体错在哪里，希望跟我错一样的朋友以后不再犯这样的错误,ORZ。

数感是真的烂。。。那些数的对应关系弄错了,ORZ

#include<stdio.h>
#include<string.h>
#include<algorithm>
using namespace std;
#define LL long long int
const int MAXN=20;
struct node
{
    LL m[MAXN][MAXN];
    node()
    {
        memset(m,0,sizeof(m));
    }
};
int d,n,m,MOD;
LL f[20],e[MAXN];
node cla(node a,node b)
{
    node c;
    for(int i=0;i<d;i++)
        for(int j=0;j<d;j++)
            for(int k=0;k<d;k++)
            if(a.m[i][k]&&b.m[k][j])
            {
                c.m[i][j]+=a.m[i][k]*b.m[k][j];
                c.m[i][j]%=MOD; 
            } 
    return c;       
}
void POW(int k)
{
    node a,c;
    for(int i=0;i<d;i++) a.m[0][i]=e[i]; 
    for(int i=1;i<d;i++) a.m[i][i-1]=1;
    for(int i=0;i<d;i++) c.m[i][i]=1;
    while(k)
    {
        if(k&1) c=cla(c,a);
        a=cla(a,a);
        k>>=1;
    }
//  for(int i=0;i<d;i++)
//  {
//      for(int j=0;j<d;j++)
//      printf("%d ",c.m[i][j]);
//      printf("\n");
//  }

    LL sum=0;
    for(int i=0,j=d-1;i<d;i++,j--)
    {
        sum=(sum+c.m[0][i]*f[j])%MOD;
    }
    printf("%lld\n",sum);
}
int main()
{
    while(~scanf("%d%d%d",&d,&n,&m))
    {
        if(d==0&&n==0&&m==0)
        break;
        MOD=m;
        memset(f,0,sizeof(f));
        memset(e,0,sizeof(e));
        for(int i=0;i<d;i++)
        scanf("%lld",&e[i]),e[i]%=MOD;
        for(int i=0;i<d;i++)
        scanf("%lld",&f[i]),f[i]%=MOD;
        if(n>d)
        POW(n-d);
        else
        printf("%lld\n",f[n-1]);
    }   
}