思路：

1. 首先 先普及一个性质： 循环矩阵*循环矩阵=循环矩阵

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

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

// by SiriusRen#include 
     
      #include 
      
       #define LL long longusing 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
       
        >=1;    }    for(int i=0;i
       
      
     

// by SiriusRen#include 
     
      #define LL long long#define f(Q,R) for(int Q=0;Q
      
       >=1;    }    f(i,n)printf("%lld ",f[i]);}
      
     

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