51nod 1161 Partial Sums(组合数学)

题目链接
把每次求前缀和操作在纸上画一下就会发现每个位置的值就等于这次的上一个位置的值加上上次的这个位置的值之和,就可以想到和机器人走方格有关,找规律之后我们发现求k次前缀和之后 a i a_i ai对位置 j ( j > = i ) j (j>=i) j(j>=i)的贡献的个数为下式,因为j-i最多n种,预处理逆元递推计算出这n个值就完事了。
C k − 1 + ( j − i ) j − i = ∏ p = k k − 1 + j − i p ( j − i ) ! C_{k - 1 + (j - i)}^{j - i}=\frac{\prod_{p=k}^{k-1+j-i}p}{(j-i)!} Ck1+(ji)ji=(ji)!p=kk1+jip

#include<bits/stdc++.h>
using namespace std;
typedef long long ll;
const int maxn = 5e3 + 5, mod = 1e9 + 7;
ll inv[maxn], f[maxn], a[maxn], ans[maxn];
int main()
{
    int n, k;
    cin >> n >> k;
    for (int i = 1; i <= n; ++i) cin >> a[i];

    inv[1] = 1;
    for (int i = 2; i <= n; ++i) inv[i] = (ll)(mod - mod / i) * inv[mod % i] % mod;

    f[0] = 1;
    for (int i = 1; i <= n; ++i) f[i] = f[i - 1] * (k + i - 1) % mod * inv[i] % mod;

    for (int i = 1; i <= n; ++i){
        for (int j = i; j <= n; ++j){
            ans[j] += f[j - i] * a[i] % mod;
            ans[j] %= mod;
        }
    }
    for (int i = 1; i <= n; ++i) cout << ans[i] << '\n';
    return 0;
}

  • 0
    点赞
  • 0
    收藏
    觉得还不错? 一键收藏
  • 0
    评论

“相关推荐”对你有帮助么?

  • 非常没帮助
  • 没帮助
  • 一般
  • 有帮助
  • 非常有帮助
提交
评论
添加红包

请填写红包祝福语或标题

红包个数最小为10个

红包金额最低5元

当前余额3.43前往充值 >
需支付:10.00
成就一亿技术人!
领取后你会自动成为博主和红包主的粉丝 规则
hope_wisdom
发出的红包
实付
使用余额支付
点击重新获取
扫码支付
钱包余额 0

抵扣说明:

1.余额是钱包充值的虚拟货币,按照1:1的比例进行支付金额的抵扣。
2.余额无法直接购买下载,可以购买VIP、付费专栏及课程。

余额充值