bzoj4402: Claris的剑

第一眼——DP，然而本质不同。。。

若每种剑只取字典序最小的，可表示为：1 +  若干个（2 1）+ 2 + 若干个（3 2）+3 +……+k （+k-1）

然后，设剑最大宽度为k，必有序列1……2……3…………k…用组合数将剩下的“锯齿”放入“盒子”，加一个“不作为剑的部分盒子”。

然后线性求n!逆元，然后就没有然后了。

#include<iostream> 
#include<cstdio>
#define ll long long
#define N 2000005
#define p 1000000007
#define F(x,y) (ll)fac[(x)+(y)-1]*inv[(x)]%p*inv[(y)-1]%p //x个球，y个盒子 
using namespace std;
int inv[N],fac[N],n,m;ll Ans;
int main ()
{
	scanf("%d%d",&n,&m);
	fac[0]=1;inv[0]=1;
	fac[1]=1;inv[1]=1;
	for (int i=2;i<=n;i++)
	{
		fac[i]=(ll)fac[i-1]*i%p;
		inv[i]=p-(ll)inv[p%i]*(p/i)%p;
	}
	for (int i=2;i<=n;i++)
		inv[i]=(ll)inv[i]*inv[i-1]%p;
	Ans=1;
	for (int i=2;i<=min(m,n);i++)
	{
		Ans+=F((n-i)/2,i);
		if (i+1<=n) Ans+=F((n-i-1)/2,i);
		Ans%=p;
	}
	printf("%lld\n",Ans);
}