BZOJ 2655 calc(拉格朗日插值)

题目
推荐一篇博客

这种因为某一个东西范围很大的，多半用插值多项式来解决。

#include<bits/stdc++.h>
#define maxn 1505
using namespace std;

int A,n,mod,f[maxn][maxn],a[maxn],sums[maxn],invf[maxn]={1,1},inv[maxn]={1,1};

int main(){
	scanf("%d%d%d",&A,&n,&mod);
	for(int i=0;i<=3*n;i++) f[0][i] = 1;
	for(int i=1;i<=n;i++)
		for(int j=1;j<=3*n;j++)
			f[i][j] = (1ll * f[i-1][j-1] * i * j + f[i][j-1]) % mod;
	for(int i=0;i<=2*n;i++) a[i] = f[n][i+n];
	A -= n;
	n = 2 * n;
	sums[n+1] = 1;
	for(int i=n;i>=0;i--) sums[i] = 1ll * sums[i+1] * (A-i) % mod;
	for(int i=2;i<=n;i++) inv[i] = 1ll * (mod - mod / i) * inv[mod % i] % mod , 
		invf[i] = 1ll * invf[i-1] * inv[i] % mod;
	int ans = 0;
	for(int i=0,ps=1;i<=n;i++){
		ans = (ans + ((n-i&1)?-1ll:1ll) * a[i] * invf[n-i] % mod * invf[i] % mod * ps % mod * sums[i+1]) % mod;
		ps = 1ll * ps * (A-i) % mod;
	}
	printf("%d\n",(ans+mod)%mod);
}

但是和这道题的容斥解法一比，插值好像没有任何技术含量。。。。。。
dalao博客

$f[i]$表示长度为$i$的没有重复元素的序列和。
$g[i]$表示${\sum }_{j=1}^A{j}^i$
然后我们来尝试直接
$f[i]=f[i-1]g[1]$
发现$f[i-1]$中（最多只）可能有一个和$g[1]$中造成贡献（被加入到最尾端）的那个数相同。
那么
$f[i]=f[i-1]g[1]-\left(\genfrac{}{}{0ex}{}{i-1}{1}\right)f[i-2]g[2]$
那个组合数代表我们强制让一个数等于最尾端的数并且插入了前面$i-2$个数中。
但是发现$f[i-2]$中(最多只)可能有一个和$g[2]$中造成贡献的那个数相同,减多了。
那么
$f[i]=f[i-1]g[1]-\left(\genfrac{}{}{0ex}{}{i-1}{1}\right)f[i-2]g[2]+\left(\genfrac{}{}{0ex}{}{i-1}{2}\right)2!f[i-3]g[3]$
钻出来的$2!$是因为如果有3个相同的数字,$\left(\genfrac{}{}{0ex}{}{i-1}{1}\right)f[i-2]g[2]$是固定了最后一个，另外一个在剩余$2$个相同的中选一个来贡献，所以一个方案会被它统计两次。
那么去重就需要去两倍。
那么不难发现实际上.
$f[i]={\sum }_{j=1}^i(-1{)}^{i-j-1}(\genfrac{}{}{0ex}{}{i-1}{i-j-1})(i-j-1)!f[i-j]g[j]$
$g$用伯努利数计算即可。
这种前一个只会算多后一个的容斥模型比较奇怪。。。。。。