[BZOJ3622]已经没有什么好害怕的了 二项式反演

恰好$k$组不好求，先求至少$k$组，设恰好$k$组方案数为$G_k$，至少为$F_k$。
首先把$a$，$b$都排序，然后求出$t_i$表示$a_i$比$b$中多少个数大。设$f_{i,j}$为考虑$a_{1..i}$中选出$j$个比$b$大的组的方案数，有转移：

$$f_{i,j}=f_{i-1,j}+f_{i-1,j-1}*(t_i-j+1)$$
然后 $F_i$即为 $f_{n,i}*(n-i)!$，就是剩下的随意匹配。
那么有：
$$F_k=\sum_{i=k}^n {i \choose k}G_i$$
二项式反演就有：
$$G_k=\sum_{i=k}^n (-1)^{i-k}{i\choose k}F_i$$
$O(n^2)$求即可。
代码：

#include<iostream>
#include<cstdio>
#include<cstring>
#include<algorithm>
#define ll long long
using namespace std;
const int mod=1000000009;
int n,k,a[2010],b[2010],t[2010];
ll f[2010][2010],C[2010][2010],fac[2010];
int main()
{
    scanf("%d%d",&n,&k);
    k=n+k>>1;
    for(int i=1;i<=n;i++)
        scanf("%d",&a[i]);
    for(int i=1;i<=n;i++)
        scanf("%d",&b[i]);
    sort(a+1,a+n+1);
    sort(b+1,b+n+1);        
    for(int i=1;i<=n;i++)
        for(int j=(t[i]=t[i-1])+1;j<=n&&b[j]<a[i];j++,t[i]++);    
    for(int i=0;i<=n;C[i][0]=1,i++)
        for(int j=1;j<=i;j++)
            C[i][j]=(C[i-1][j-1]+C[i-1][j])%mod;
    fac[0]=1;
    for(int i=1;i<=n;i++)
        fac[i]=fac[i-1]*i%mod;  
    for(int i=0;i<=n;i++) f[i][0]=1; 
    for(int i=1;i<=n;i++)
        for(int j=1;j<=i;j++)        
            f[i][j]=(f[i-1][j-1]*max((t[i]-j+1),0)+f[i-1][j])%mod;
    ll ans=0;
    for(int i=k,sgn=1;i<=n;i++,sgn=-sgn)
        (ans+=(C[i][k]*f[n][i]%mod*fac[n-i]%mod*sgn+mod))%=mod;
    printf("%lld",ans);         
    return 0;
}