loj #6024. XLkxc

背景：

$\text{bzoj}$权限题？
没事$\text{loj}$有原题。

题目传送门：

https://loj.ac/problem/6024

题意：

求：${\sum }_{i=0}^n{\sum }_{j=1}^{a+i\cdot d}{\sum }_{l=1}^j{l}^k\mathrm{m}\mathrm{o}\mathrm{d}p$。

思路：

做法很妙，实在想不到。
一般人会去化简式子，可是根本化不动，我太菜了。
巨佬们给的做法是：
设$f_x={\sum }_{i=1}^x{i}^k$。
我们发现$f_x$是一个$k+1$次的多项式。
证明就是：
$$\begin{array}{rl}{f}_x& =\sum _{i=1}^x{i}^k\\ & =\frac{1}{k+1}\sum _{i=1}^{k+1}{C}_{k+1}^i\cdot B_{k+1-i}\cdot (x+1{)}^i\\ & =\frac{1}{k+1}\sum _{i=1}^{k+1}{C}_{k+1}^i\cdot B_{k+1-i}\cdot \sum _{j=0}^{i}C{(}_i^j)x^j{1}^{i-j}\\ & =\frac{1}{k+1}\sum _{i=1}^{k+1}{C}_{k+1}^i\cdot B_{k+1-i}\cdot \sum _{j=0}^{i}C{(}_i^j)x^j\end{array}$$
其中$B$是伯努利数。
$j$的上限是$i$，$i$的上限为$k+1$，故$x^j$的最大次数为$k+1$，因此得证。

设$g_x={\sum }_{i=1}^x{f}_i$。
我们发现$g$是一个$k+2$次的多项式。
证明就是：
$$\begin{array}{rl}{g}_x& =\sum _{i=1}^x\sum _{j=1}^j{j}^k\\ & =\sum _{i=1}^x\frac{1}{k+1}\sum _{j=1}^{k+1}{C}_{k+1}^j\cdot B_{k+1-j}\cdot (i+1{)}^j\\ & =\frac{1}{k+1}\sum _{i=1}^x\sum _{j=1}^{k+1}{C}_{k+1}^j\cdot B_{k+1-j}\cdot \sum _{l=0}^j{C}_j^l\cdot i^l\cdot 1^{j-l}\\ & =\frac{1}{k+1}\sum _{i=1}^x\sum _{j=1}^{k+1}{C}_{k+1}^j\cdot B_{k+1-j}\cdot \sum _{l=0}^j{C}_j^l\cdot i^l\\ & =\frac{1}{k+1}\sum _{j=1}^{k+1}{C}_{k+1}^j\cdot B_{k+1-j}\cdot \sum _{l=0}^j{C}_j^l\cdot \sum _{i=1}^x{i}^l\end{array}$$

上面我们已经证明了形如${\sum }_{i=1}^x{i}^l$是$l+1$次方的，$l$的上限是$j$，$j$的上限是$k+1$，当$l=j=k+1$时，为$k+2$次方，得证。

设$h_x={\sum }_{i=1}^x{g}_{a+i\cdot d}$。
我们发现$g$是一个$k+3$次的多项式。
证明就是：
在原来的基础上，再套一层即可。

剩下的按照次数拉格朗日插值即可。
时间复杂度好像有点大：$\Theta (n^3)$。
其实可以$\Theta (n^2)$（重心拉格朗日插值），当然我很懒，既然模板都没打，为什么要打呢。

代码：

#include<cstdio>
#include<cstring>
#include<algorithm>
#define LL long long
#define LD long double
#define mod 1234567891ll
using namespace std;
	int k;
	LL p,n,d;
	LL f[1010],g[1010],h[1010];
	struct node{LL x,y;} a[1010];
LL times(LL x,LL y)
{
	LD t=(LD)x*(LD)y-(LL)((LD)x*(LD)y/(LD)mod)*(LD)mod;
	LL tt=(LL)t;
	return tt>=0?tt%mod:(tt+mod)%mod;
}
LL dg(LL x,LL y)
{
	if(!y) return 1;
	LL op=dg(x,y>>1);
	if(y&1) return times(times(op,op),x); else return times(op,op);
}
LL inv(LL x)
{
	return dg(x,mod-2);
}
LL lagrange(LL x,int limit)
{
	LL sum=0;
	for(int i=0;i<=limit;i++)
	{
		LL op1=a[i].y,op2=1;
		for(int j=0;j<=limit;j++)
		{
			if(i==j) continue;
			op1=times(op1,(x-a[j].x+mod)%mod);
			op2=times(op2,(a[i].x-a[j].x+mod)%mod);
		}
		sum=(sum+times(op1,inv(op2)))%mod;
	}
	return sum;
}
int main()
{
	int T;
	scanf("%d",&T);
	while(T--)
	{
		scanf("%d %lld %lld %lld",&k,&p,&n,&d);
		for(int i=1;i<=k+3;i++)
			f[i]=(f[i-1]+dg(i,k))%mod;
		for(int i=1;i<=k+3;i++)
			g[i]=(g[i-1]+f[i])%mod;
		for(int i=0;i<=k+3;i++)
			a[i]=(node){(LL)i,g[i]};
		for(int i=0;i<=k+4;i++)
			h[i]=((!i?0:h[i-1])+lagrange((p+times(i,d))%mod,k+3))%mod;
		for(int i=0;i<=k+4;i++)
			a[i]=(node){(LL)i,h[i]};
		printf("%lld\n",lagrange(n,k+4));
	}
}