bzoj3622: 已经没有什么好害怕的了【动态规划+容斥原理】

解题思路：

当$n+k$是奇数时肯定无解。
然后$k=(n+k)/2$，直接算糖果大于药片的对数刚好有$k$个。

恰好为$k$ 的不好算, 那就来算 $≥ k$ 的.

首先对两个队的数组分别从小到大排序, 记 $cnt _i$ 表示 $a _i$ 比 B中多少元素大.

设$f(i,j)$ 表示考虑了 $a _1 ,...,a_i$ , 在其中选出$j$ 个, 且这$j$ 对都满足$a > b$，不考虑剩下的.

显然有
$f(i,j) = f(i − 1,j) + (cnt_ i− j + 1)f(i − 1,j − 1)$

可能有的同学觉得$f(n,i)=f(n,i)*(n-i)!$就是对数大于$i$的方案数，然后恰好$i$对就是$f(n,i+1)-f(n,i)$，其实并不是，因为$f(n,i)$中并不是每种满足的配对方式只算了一次，而是算了多次。

如：
糖果: $a_1，a_2，a_3$
药片: $b_1，b_2，b_3$
$a_1>b_1， a_2>b_2， a_3>b_3$

那么
$((a_1，b_1)，(a_2，b_2)，a_3不明) ((a_1，b_1)，(a_3，b_3)，a_2不明)$
就会被视为两种答案，其实是同一种。

考虑$g(j)$是恰好有$j$对的方案，那么它在$f(n,i)(i<j)$中被算了$C_j^i$次，所以正确答案为：
$g(i)=f(n,i)-\sum\limits_{j=i+1}^ng(j)$

#include<bits/stdc++.h>
#define ll long long
using namespace std;

int getint()
{
    int i=0,f=1;char c;
    for(c=getchar();(c!='-')&&(c<'0'||c>'9');c=getchar());
    if(c=='-')c=getchar(),f=-1;
    for(;c>='0'&&c<='9';c=getchar())i=(i<<3)+(i<<1)+c-'0';
    return i*f;
}

const int N=2005,p=1e9+9;
int n,k,a[N],b[N],cnt[N];
int f[N][N],c[N][N],fac[N];

int main()
{
    //freopen("lx.in","r",stdin);
    n=getint(),k=getint();
    if ((n+k)&1) {cout<<0;return 0;}
    k=(n+k)>>1;
    for(int i=0;i<=n;i++)c[i][0]=1;
    for(int i=1;i<=n;i++)
        for(int j=1;j<=i;j++)
            c[i][j]=(c[i-1][j]+c[i-1][j-1])%p;
    fac[0]=1;
    for(int i=1;i<=n;i++)fac[i]=(ll)fac[i-1]*i%p;
    for(int i=1;i<=n;i++)a[i]=getint();
    for(int i=1;i<=n;i++)b[i]=getint();
    sort(a+1,a+n+1);
    sort(b+1,b+n+1);
    for(int i=1;i<=n;i++)
        for(int j=1;j<=n;j++)
            if(a[i]>b[j])cnt[i]++;
            else break;
    f[0][0]=1;
    for(int i=1;i<=n;i++)
        for(int j=0;j<=i;j++)
        {
            f[i][j]=f[i-1][j];
            if(j&&cnt[i]-j+1>0)f[i][j]=(f[i][j]+(ll)(cnt[i]-j+1)*f[i-1][j-1]%p)%p;
        }
    for(int i=n;i>=k;i--)f[n][i]=(ll)f[n][i]*fac[n-i]%p;
    for(int i=n;i>=k;i--)
        for(int j=i+1;j<=n;j++)f[n][i]=(f[n][i]-(ll)c[j][i]*f[n][j]%p+p)%p;
    cout<<f[n][k];
    return 0;       
}