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

题意：

给两个长度为n，所有数互不相同的序列，然后两两配对

题解：

好题啊，几乎没什么思路。
首先假如n，k奇偶性不一样显然无解，要求有(n+k)/2对每组要求。
看到这种计数题考虑容斥，对两个序列排序后，$f[i][j]$表示前i个a列的配对，有至少j个符合要求。
预处理个$next[i]$，表示第i个数最多比b列第几个数大。
所以$f[i][j]=f[i-1][j]+f[i-1][j-1]*(next[i]-j+1)$
然后就要减去大于k对的情况，即求$f'[n][k]$
于是可以推出$f'[n][i]=f[n][i]*(n-i)!-\sum_{j=i+1}^n f'[n][j]*C_j^i$
这个挺好理解把，就是其他的随便匹配，然后减去>i的情况。
code：

#include<cstdio>
#include<cstdlib>
#include<cstring>
#include<iostream>
#include<algorithm>
#define LL long long
using namespace std;
const LL mod=1000000009;
LL n,k,a[2010],b[2010],next[2010];
LL g[2010],C[2010][2010];
LL f[2010][2010];
void pre()
{
    g[0]=1;
    for(LL i=1;i<=2000;i++) g[i]=g[i-1]*i%mod;
    for(LL i=0;i<=2000;i++)
    {
        C[i][0]=1;
        for(LL j=1;j<=i;j++) C[i][j]=(C[i-1][j]+C[i-1][j-1])%mod;
    }
}
int main()
{
    scanf("%lld %lld",&n,&k);
    if((n+k)%2!=0) {printf("0");return 0;}
    k=(n+k)>>1;
    for(LL i=1;i<=n;i++) scanf("%lld",&a[i]);
    for(LL i=1;i<=n;i++) scanf("%lld",&b[i]);
    sort(a+1,a+n+1);sort(b+1,b+n+1);
    for(LL i=1;i<=n;i++)
    {
        LL j=0;
        while(j<n&&b[j+1]<a[i]) j++;
        next[i]=j;
    }
    pre();
    memset(f,0,sizeof(f));
    f[0][0]=1;
    for(LL i=1;i<=n;i++)
        for(LL j=0;j<=i;j++)
        {
            f[i][j]=f[i-1][j];
            if(j>=1&&next[i]>=j) f[i][j]=(f[i][j]+f[i-1][j-1]*(next[i]-j+1)%mod)%mod;
        }
    for(LL i=n;i>=k;i--)
    {
        f[n][i]=f[n][i]*g[n-i]%mod;
        for(LL j=i+1;j<=n;j++) f[n][i]=(f[n][i]-f[n][j]*C[j][i]%mod)%mod;
        f[n][i]=(f[n][i]+mod)%mod;
    }
    printf("%lld",f[n][k]);
}