思路：

1. 首先 先普及一个性质： 循环矩阵*循环矩阵=循环矩阵
   
   由于此题是距离小于d的都加上一个数。
   
   那么 构造矩阵的时候 我们发现 诶呦 这是个循环矩阵
   
   看看数据范围 n^2log(k)可以过。
   
   那就把这个矩阵改一改。
   
   因为这是个循环矩阵， 所以呢 只用保存一行就可以了。
   
   每回做乘法的时候只做第一行的乘法。
   
   for(i) for(j) temp[i]+=a[j]*b[(i+j)%n];
   
   就这么着 搞搞就能过了。

2. （好像可以用FFT？ 表示并不会）

// by SiriusRen

#include <cstdio>

#include <cstring>

#define LL long long

using namespace std;

long long f[666],co[666],c[666];

int n,m,d,k;

void mul(LL *a,LL *b){

    memset(c,0,sizeof(c));

    for(int i=0;i<n;i++)

        for(int j=0;j<n;j++)

            c[i]+=a[j]*b[(i+j)%n];

    for(int i=0;i<n;i++)b[i]=c[i]%m;

}

int main(){

    scanf("%d%d%d%d",&n,&m,&d,&k);

    for(int i=0;i<n;i++)scanf("%lld",&f[i]);

    for(int i=0;i<=d;i++)co[i]=co[n-i]=1;

    while(k){

        if(k&1)mul(co,f);

        mul(co,co),k>>=1;

    }

    for(int i=0;i<n;i++)printf("%lld ",f[i]);

}

// by SiriusRen

#include <cstdio>

#define LL long long

#define f(Q,R) for(int Q=0;Q<R;Q++)

using namespace std;

LL f[666],co[666],c[666];

int n,m,d,k;

void mul(LL *a,LL *b){

    f(i,n)c[i]=0;

    f(i,n)f(j,n)c[(i+j)%n]+=a[i]*b[j];

    f(i,n)b[i]=c[i]%m;

}

int main(){

    scanf("%d%d%d%d",&n,&m,&d,&k);

    f(i,n)scanf("%lld",&f[i]);

    f(i,d+1)co[i]=co[n-i]=1;

    while(k){

        if(k&1)mul(co,f);

        mul(co,co),k>>=1;

    }

    f(i,n)printf("%lld ",f[i]);

}

Code length能进前三的存在哈哈哈