JZOJ-senior-5922. 【NOIP2018模拟10.23】sequence

Time Limits: 1000 ms Memory Limits: 524288 KB Detailed Limits

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 2

见选手目录下的sequence/sequence2.in与sequence/sequence2.ans。
该组样例的数据范围同第 1 个测试点。

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。除此之外，对于每个数据点，还满足以下限制。

Solution

我们发现k很小，考虑一下 $C(n,k)$ 是什么呢？
首先我们知道 $C$ 这个组合数是跟杨辉三角形有关的，于是我们就在往这方面去想
对于每一次增加的便利值，相当于将一段数的每一个位置加上杨辉三角形中对应的某一列的某一段
$C(n,k)$ 是一个第 $0$ 列全为 $1$ 的杨辉三角形，每一列就是前一列的前缀和，而这里的 $k$ 就相当于做 $k$ 次前缀和，于是我们就用差分来做这题
设 $F[k][i]$ 表示当前差分数组还要进行 $K$ 次累加才可以得到正确答案
对于每次优化，将 $F[k][l]$++， $F[k][r+1]$ --，打好标记
最后从下往上进行求前缀操作，一切看起来都是那么的顺利

但是，我们发现如果一个差分数组长这样
(line k) 1 0 0 0 -1 0 0
它在进行第一次操作时
(line k-1) 1 1 1 1 0 0 0
十分的正常
但是它在第二次操作时
(line k-2) 1 2 3 4 4 4 4
它就不是我们想要的模样
我们希望能在第五个位置打上-4标记

于是，对于每次优化，我们都需要打上修正标记
至于标记是什么，随手推推就好了，详见代码

Code

#include<algorithm>
#include<cstdio>

#define fo(i,a,b) for(int i=a;i<=b;++i)
#define fd(i,a,b) for(int i=a;i>=b;--i)
#define ll long long

using namespace std;

const int N=5e5+5,K=21,P=1e9+7;
int n,m,l,r,k;
ll jc[N],ny[N],f[K][N];

void dec(ll &x,ll y)
{
	x-=y; if(x<0) x=((x%P)+P)%P;
}

void inc(ll &x,ll y)
{
	x+=y; if(x>=P) x%=P;
}

ll ksm(ll a,ll b)
{
	ll s=1;
	for(;b;a=a*a%P,b>>=1)
		if(b&1) s=s*a%P;
	return s;
}

ll C(int x,int y)
{
	return jc[x]*ny[y]%P*ny[x-y]%P;
}

int main()
{
	freopen("sequence.in","r",stdin);
	freopen("sequence.out","w",stdout);
	scanf("%d%d",&n,&m);
	jc[0]=1;
	fo(i,1,N-5) jc[i]=jc[i-1]*i%P;
	ny[N-5]=ksm(jc[N-5],P-2);
	fd(i,N-6,0) ny[i]=ny[i+1]*(i+1)%P;
	fo(i,1,m)
	{
		scanf("%d%d%d",&l,&r,&k);
		++f[k][l],--f[k][r+1]; int p=1;
		fd(j,k-1,0) dec(f[j][r+1],C(r-l+p,p)),++p;//差分修正标记
	}
	fd(i,20,1)
		fo(j,1,n)
		{
			inc(f[i][j],f[i][j-1]);
			inc(f[i-1][j],f[i][j]);
		}
	fo(i,1,n) inc(f[0][i],f[0][i-1]);
	fo(i,1,n) printf("%lld\n",f[0][i]);
}