BZOJ 3622 已经没有什么好害怕的了 DP 容斥原理

题目大意：给出长度为n的a,b两个序列，元素互不相同。求有多少种匹配满足$a_i>b_i$比$a_i< b_i$多m对。

将问题转化一下，即求有多少种匹配其中$a_i > b_i$恰好为$\frac{n+m}{2}$，其余均为$a_i < b_i$。
若设状态为$f_{i,j}$表示a序列中前$i$个元素匹配完恰好有j对为$a_i > b_i$，发现根本没法转移。

看到计数问题可以考虑容斥。

需要将状态设置得便于转移，将状态改为$f_{i,j}$表示a序列中前$i$个元素至少匹配了$j$对满足$a_i>b_i$，也就是说先找出$j$对，其他的暂时还不知道。
先将a,b序列排序，预处理出$maxx_i$表示$a_i < b_j$的最后一个位置在哪。转移方程为

$$f_{i,j}=f_{i-1,j}+f_{i-1,j-1}\times(maxx_i-j+1)$$

利用容斥，设$g_j$为a序列n个元素恰好匹配$j$组满足$a_i > b_i$的方案数，则

$$g_j=f_{n,j}\times fac_{n-j}-\sum_{k=j+1}^{k<=n}{g_k\times C_k^j}$$
至少匹配 $j$组$f_{n,j}$减去恰好匹配 $k$组，其中$i< k\le n$，即把大于 $j$的都减去。因为$k$组匹配要选出 j <script type="math/tex" id="MathJax-Element-243">j</script>个所以乘个组合数就好了。

#include <cstdio>
#include <algorithm>
#define N 2005
#define MOD 1000000009
using namespace std;
typedef long long LL;
int n,m,a[N],b[N],maxx[N];
LL f_pow(LL x,LL y) {
    LL tmp=1;
    while(y) {
        if(y&1) (tmp*=x)%=MOD;
        (x*=x)%=MOD;
        y>>=1;
    }
    return tmp;
}
LL fac[N],inv[N],f[N][N],g[N];
void init() {
    fac[1]=1;
    for(LL i=2;i<=n;i++) fac[i]=(fac[i-1]*i)%MOD;
    inv[n]=f_pow(fac[n],MOD-2);
    for(LL i=n-1;i;i--) inv[i]=inv[i+1]*(i+1)%MOD;
    return ;
}
LL C(int n,int m) { return fac[n]*inv[m]%MOD*inv[n-m]%MOD; }
int main() {
    scanf("%d%d",&n,&m);
    for(int i=1;i<=n;i++) scanf("%d",a+i);
    for(int i=1;i<=n;i++) scanf("%d",b+i);
    if(n+m&1) {
        printf("0\n");
        return 0;
    }
    (m+=n)/=2;
    sort(a+1,a+n+1), sort(b+1,b+n+1);
    for(int i=1;i<=n;i++) {
        maxx[i]=maxx[i-1];
        while(b[maxx[i]]<a[i] && maxx[i]<=n) maxx[i]++;
        maxx[i]--;
    }
    init();
    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]+f[i-1][j-1]*max(0,maxx[i]-j+1))%MOD;
    for(int i=n;i>=m;i--) g[i]=f[n][i]*fac[n-i]%MOD;
    for(int i=n;i>=m;i--)
        for(int j=i+1;j<=n;j++)
            (g[i]+=MOD-g[j]*C(j,i)%MOD)%=MOD;
    printf("%lld\n",g[m]);
    return 0;
}