BZOJ3622: 已经没有什么好害怕的了

题意：给定两个长度为N的互不相同的序列，将其两两配对，求有多少种配对使得a[i]>b[i]的对数减去其他对数恰好为k。
一道经典的计数题，但思路理清后发现并不算难。
首先设x为a[i]小于b[i]的对数。
x+x+k=n
2*x=n-k
所以当n-k为奇数，答案为0。
先令ne[i]表示小于a[i]的最后一个b[j]的位置。首先考虑设f[i][j]表示到前i个糖果，有j个位置满足配对,且不考虑后续情况。
转移显然为
$$f[i][j] = f[i-1][j]+f[i-1][j-1]*(ne[i]-j+1)$$
如果我们要考虑后续呢？
我们可以设g[i]表示只有i组糖果大于药片
此时方程为
$$g[i] = f[n][i]*(n-i)!-\sum_{j=i+1}^n g[j]*\dbinom{j}{i}$$
 前者枚举对数形态而后者表示实际有j组数对但被统计到答案中去。

 倒推即可。复杂度o(n^2)。

实际上笔者在写代码时g,f调换了顺序。

#include<bits/stdc++.h>
using namespace std;
typedef long long ll;
int read(){
	int x=0,f=1; char ch=getchar();
	while (ch<'0'||ch>'9'){if (ch=='-') f=-1; ch=getchar();}
	while (ch>='0'&&ch<='9'){x=x*10+ch-'0'; ch=getchar();}
	return x*f;
}
const int mod=1e9+9;
int a[2005],b[2005],ne[2005];
ll f[2005],g[2005][2005],c[2005][2005],fac[2005];
int main(){
	int n=read(),k=read();
	if ((n-k)%2==1){printf("0\n"); return 0;}
	for (int i=1;i<=n;i++)a[i]=read();
	for (int i=1;i<=n;i++)b[i]=read();
	sort(a+1,a+1+n); sort(b+1,b+1+n);
	ne[0]=1;
	k=(k+n)>>1;
	for (int i=1;i<=n;i++){
		for (int j=ne[i-1];j<=n;j++){
			if (b[j]>a[i]) break; ne[i]=j;
		}
	}
	for (int i=0;i<=n;i++){
		c[i][0]=1;
		for (int j=1;j<=i;j++)
		c[i][j]=(c[i-1][j]+c[i-1][j-1])%mod;
	}
	g[0][0]=1;
	for (int i=1;i<=n;i++)
	for (int j=0;j<=i;j++){
	g[i][j]=g[i-1][j];
	if (j&&ne[i]-j+1>=1) 
	(g[i][j]+=g[i-1][j-1]*(ne[i]-j+1))%=mod;}
	fac[0]=1;
	for (int i=1;i<=n;i++)
	fac[i]=(fac[i-1]*i)%mod;
	for (int i=n;i>=k;i--){
		(f[i]=g[n][i]*fac[n-i])%=mod;
		for (int j=i+1;j<=n;j++)
		f[i]=(f[i]-(f[j]*c[j][i])+mod)%mod;
	}
	printf("%lld\n",f[k]);
	return 0;
}