[SDOI2008] 递归数列

题目描述：

雾。

题目分析：

矩阵快速幂。

搞出两个前缀和，相减即可

题目链接：

Luoogu 2461

Ac 代码：

// luogu-judger-enable-o2
#include <cstdio>
#include <iostream>
#include <algorithm>
#include <cstring>
#define int long long
int mod;
const int maxm=1010;
struct Matrix{
    int n,m;
    int a[21][21];
    Matrix()
    {
        memset(a,0,sizeof(a));
        n=m=0;
    }
};
Matrix f,sum1,sum2,t;
inline Matrix operator * (Matrix t1,Matrix t2)
{
    Matrix ans;
    ans.n=t1.n,ans.m=t2.m;
    for(int i=1;i<=t1.n;i++)
     for(int j=1;j<=t2.m;j++)
      for(int k=1;k<=t1.m;k++)
      {
        ans.a[i][j]=(ans.a[i][j]+(t1.a[i][k]*t2.a[k][j])%mod+mod)%mod;
        //printf("%d %d %d\n",t1.a[i][k],t2.a[k][j],t1.a[i][k]*t2.a[k][j]);
      }

    return ans;
}
int b[maxm],c[maxm],len;
inline Matrix operator ^ (Matrix x,int y)
{
    Matrix ans;
    ans.n=ans.m=len+1;
    for(int i=1;i<=len+1;i++) ans.a[i][i]=1;
    while(y)
    {
        if(y&1) ans=ans*x;
        x=x*x;
        y/=2; 
    }

    return ans;
}
main()
{
    scanf("%lld",&len);
    for(int i=1;i<=len;i++) scanf("%lld",&b[i]);
    for(int i=1;i<=len;i++) scanf("%lld",&c[i]);
    int n,m;
    scanf("%lld%lld%lld",&m,&n,&mod);
    for(int i=1;i<=len;i++) b[i]%=mod,c[i]%=mod;
    f.n=f.m=len+1;
    f.a[1][1]=1;
    for(int i=3;i<=len+1;i++) f.a[i][i-1]= 1;
    for(int i=2;i<=len+1;i++) f.a[1][i]=f.a[2][i]=c[i-1];
    t.n=len+1,t.m=1;
    for(int i=2;i<=len+1;i++)
    {
        t.a[i][1]=b[len-i+2];
        t.a[1][1]=(t.a[1][1]+b[i-1])%mod;
    }
    if(n-len>0) sum1=(f^(n-len))*t;
    else 
    for(int i=1;i<=n;i++)
      sum1.a[1][1]=(sum1.a[1][1]+b[i])%mod;
    if(m-len-1>0) sum2=(f^(m-len-1))*t;
    else
     for(int i=1;i<m;i++)
      sum2.a[1][1]=(sum2.a[1][1]+b[i])%mod;
    printf("%lld\n",((sum1.a[1][1]-sum2.a[1][1])%mod+mod)%mod);
    return 0;
}