51nod 1161 Partial Sums

最新推荐文章于 2022-02-11 13:58:53 发布

算球？

最新推荐文章于 2022-02-11 13:58:53 发布

阅读量399

点赞数

分类专栏： ----组合数学----

本文链接：https://blog.csdn.net/gyhguoge01234/article/details/73194573

版权

----组合数学---- 专栏收录该内容

15 篇文章 0 订阅

订阅专栏

相关讨论里的没看懂，扯到矩阵就不会了，补线代。。。
参考：http://blog.csdn.net/qingshui23/article/details/52180926
http://blog.csdn.net/yitiaodacaidog/article/details/38071257
看着博客里的规律，自己画画，就出来规律了。可是真不知道怎么去找这规律。。。还要补数学。。
这个题还有V2，还要用到fft，想想就蛋疼。这里写图片描述

#include <iostream>
using namespace std;

typedef long long ll;
const int MAXN = 5555;
const ll mod = 1e9+7;
ll inv[MAXN],num[MAXN],res[MAXN];
ll N,K;

void init()
{
    inv[1] = 1;
    for(int i=2; i<= N; i++)
    {
        if(i >= mod) break;
        inv[i] = (mod - mod / i) * inv[mod % i] % mod;
    }

    res[0] = 1;
    for(int i = 1; i < N; ++i)
        res[i] = ((res[i-1]*(i+K-1)%mod)*inv[i])%mod;
}

int main()
{
    cin >> N >> K;
    init();
    for(int i = 1; i <= N; ++i)
        cin >> num[i];
    for(int i = 1; i <= N; ++i)
    {
        ll sum = 0;
        for(int j = 1; j <= i; ++j)
        {
            sum = (sum + res[i-j]*num[j])%mod;
        }
        cout << sum << endl;
    }
    return 0;
}