jzoj5922. 【NOIP2018模拟10.23】sequence（组合数）

5922. 【NOIP2018模拟10.23】sequence

Description
小 F 是一位 Hack 国的居民，他生活在一条长度为 n 的街道上，这个街道上总共有 n 个商店。每个商店里售卖着不同的 Hack 技能包，每个商店本身也会有个便利值。初始时，每个商店的便利值均为 0。每一天，街道上都会有一些商店优化改造。
具体来说，对于每一天，优化改造的商店都是一个连续的区间 l ∼ r，每次优化改造也会有一个优化参数 k。对于所有 l ≤ i ≤ r ，第 i 个商店的便利值会增加 。
小 F 想知道，m 天之后，每个商店的便利值分别是多少。由于小 F 并不喜欢高精度，因此你只需要输出便利值对 10^9 + 7 取模的结果。

Input
从文件sequence.in中读入数据。
第 1 行，两个整数 n, m 表示街道的长度与天数。
接下来的 m 行，每行三个整数 l, r, k，表示第 i 天优化改造的商店区间和优化参数。

Output
输出到文件sequence.out中，共 n 行。
每行 1 个整数，表示第 i 个商店的便利值对 109 + 7 取模的结果。

Sample Input
5 3
1 4 3
2 5 0
3 4 2

Sample Output
1
5
12
24
1
第 1 次操作之后，每个商店的便利值分别为 1, 4, 10, 20, 0。
第 2 次操作之后，每个商店的便利值分别为 1, 5, 11, 21, 1。
第 3 次操作之后，每个商店的便利值分别为 1, 5, 12, 24, 1。

Data Constraint
对于 100% 的数据，满足 1 ≤ n, m ≤ 5 × 10^5, 0 ≤ k ≤ 20。除此之外，对于每个数据点，还满足以下限制。

分析：可以发现组合数有一个简单的性质，即C(n,m)=C(n-1,m)+C(n-1,m-1)，我们可以从这个式子中获得启发。考虑一个下标从 0 开始的数列，这个数列的每个数均为 1。
我们对这个数列做 k 阶前缀和。通过这个简单的式子，容易发现这个数列的第 i项即为C(k+i,k)。于是我们可以得出一个这样的做法。
对于这个数列，维护这个数列的 k 阶差分。对于一次修改操作，我们只需要在各阶差分数组上修改，最后做一遍 k 阶的前缀和即可。注意差分数组上修改时，区间边界要相应地减掉一些值。

代码

#include <cstdio>
#include <algorithm>
#define N 600000
#define ll long long
#define mo 1000000007
using namespace std;

ll c[N][30],f[N][30];
int n,m;

int main()
{
	freopen("sequence.in", "r", stdin);
	freopen("sequence.out","w", stdout);
	scanf("%d%d", &n, &m);
	c[0][0] = 1;
	for (int i = 1; i <= n + 30; i++)
	{
		c[i][0] = 1;
		for (int j = 1; j <= 20; j++) 
			c[i][j] = (c[i - 1][j - 1] + c[i - 1][j]) % mo;
	}
	while (m--)
	{
		int l, r, k;
		scanf("%d%d%d", &l, &r, &k);
		for (int i = 0; i <= k; i++)
		{
			f[l][i] = (f[l][i] + c[k][k - i]) % mo;
			f[r + 1][i] = (f[r + 1][i] - c[r + 1 - l + k][k - i] + mo) % mo;
		}
	}
	for (int i = 1; i <= n; i++)
		for (int j = 0; j <= 20; j++)
			f[i + 1][j] = ((f[i + 1][j] + f[i][j]) % mo + f[i][j + 1]) % mo;
	for (int i = 1; i <= n; i++)
		printf("%lld\n", f[i][0]);
}