JZOJ 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 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

• 对于一个点 $x$，区间加操作加的数为$$C_{x+k-l}^k=\frac{1}{k!}\ast (x-l+1)\ast (x-l+2)\ast \dots \ast (x-l+k)$$

• 我们发现可以暴力展开这个关于 $x$ 的多项式，并将其当做一个标记（长度为 $20$）。

• 之后在 $l$ 处加上这个标记，在 $r+1$ 处减去这个标记。

• 最后求答案时做前缀和，每次得到位置 $i$ 的一个多项式，将 $i$ 代入 $x$ 算即可。

• 时间复杂度 $O(n{k}^2)$ ，吸氧就能过啦。

• 如果用 $k$ 阶差分可以做到 $O(nk)$ 。

Code

#pragma GCC optimize(2)
#include<cstdio>
#include<cstring>
#include<cctype>
using namespace std;
typedef long long LL;
const int N=5e5+5,mo=1e9+7;
int f[21],g[21],h[N][21],d[21],c[21];
inline int read()
{
	int X=0,w=0; char ch=0;
	while(!isdigit(ch)) w|=ch=='-',ch=getchar();
	while(isdigit(ch)) X=(X<<3)+(X<<1)+(ch^48),ch=getchar();
	return w?-X:X;
}
void write(int x)
{
	if(x>9) write(x/10);
	putchar(x%10+'0');
}
inline int ksm(int x,int y)
{
	int s=1;
	while(y)
	{
		if(y&1) s=(LL)s*x%mo;
		x=(LL)x*x%mo;
		y>>=1;
	}
	return s;
}
int main()
{
	freopen("sequence.in","r",stdin);
	freopen("sequence.out","w",stdout);
	int n=read(),m=read();
	f[0]=g[0]=1;
	for(int i=1;i<=20;i++) f[i]=(LL)f[i-1]*i%mo;
	g[20]=ksm(f[20],mo-2);
	for(int i=19;i;i--) g[i]=(LL)g[i+1]*(i+1)%mo;
	while(m--)
	{
		int l=read(),r=read(),k=read();
		if(!k)
		{
			h[l][0]++;
			h[r+1][0]=(h[r+1][0]-1+mo)%mo;
			continue;
		}
		c[0]=(mo-l+1)%mo;
		c[1]=1;
		int len=1;
		for(int i=2;i<=k;i++)
		{
			c[++len]=0;
			for(int j=0;j<=len;j++) d[j]=c[j];
			for(int j=1;j<=len;j++) c[j]=d[j-1];
			for(int j=c[0]=0;j<len;j++) c[j]=(c[j]+(LL)d[j]*(i-l+mo))%mo;
		}
		for(int i=0;i<=len;i++)
		{
			h[l][i]=(h[l][i]+(LL)c[i]*g[k])%mo;
			h[r+1][i]=(h[r+1][i]-(LL)c[i]*g[k]%mo+mo)%mo;
		}
	}
	memset(c,0,sizeof(c));
	for(int i=1;i<=n;i++)
	{
		for(int j=0;j<=20;j++) c[j]=(c[j]+h[i][j])%mo;
		int ans=0,base=1;
		for(int j=0;j<=20;j++)
		{
			ans=(ans+(LL)c[j]*base)%mo;
			base=(LL)base*i%mo;
		}
		write(ans),putchar('\n');
	}
	return 0;
}