【计数】cf223C. Partial Sums

考试时候遇到这种题只会找规律

You've got an array a, consisting of n integers. The array elements are indexed from 1 to n. Let's determine a two step operation like that:

1. First we build by the array a an array s of partial sums, consisting of n elements. Element number i (1 ≤ i ≤ n) of array s equals . The operation x mod y means that we take the remainder of the division of number x by number y.
2. Then we write the contents of the array s to the array a. Element number i (1 ≤ i ≤ n) of the array s becomes the i-th element of the array a (ai = si).

You task is to find array a after exactly k described operations are applied.

Input

The first line contains two space-separated integers n and k (1 ≤ n ≤ 2000, 0 ≤ k ≤ 109). The next line contains n space-separated integers a1, a2, ..., an — elements of the array a (0 ≤ ai ≤ 109).

Output

Print n integers  — elements of the array a after the operations are applied to it. Print the elements in the order of increasing of their indexes in the array a. Separate the printed numbers by spaces.

---

题目分析

可以从矩阵乘法开始想起，考虑转移矩阵，发现其主对角线下方全为0、元素按照次对角线对称。

或者就是找规律

3.16upd：

总觉得一个经典模型不应该用找规律这么假的方式随随便便搞掉吧……

网上的题解都是打表找规律 | 矩乘找规律 | dp找规律……

来自ZZK的新的理解方式：$a_i$的贡献也就是它走到$s^p_j$的方案数量。

 #include<bits/stdc++.h>
 #define MO 1000000007
 const int maxn = ;

 int n,k,a[maxn],b[maxn],inv[maxn],fac[maxn],pre[maxn];

 void init()
 {
     fac[] = fac[] = inv[] = inv[] = ;
     for (int i=; i<=n; i++)
         fac[i] = 1ll*fac[i-]*i%MO,
         inv[i] = MO-1ll*(MO/i)*inv[MO%i]%MO;
     pre[] = ;
     for (int i=; i<=n; i++)
         pre[i] = 1ll*pre[i-]*(k%MO+i-)%MO*inv[i]%MO;
 }
 int main()
 {
     scanf("%d%d",&n,&k);
     for (int i=; i<=n; i++) scanf("%d",&a[i]), b[i] = a[i];
     if (k){
         init();
         for (int i=; i<=n; i++)
         {
             b[i] = ;
             for (int j=; j<=i; j++)
                 b[i] = 1ll*(b[i]+1ll*pre[i-j]*a[j]%MO)%MO;
         }
     }
     for (int i=; i<=n; i++) printf("%d ",b[i]);
     return ;
 }

END