【BZOJ3622】已经没有什么好害怕的了，两次DP

传送门
思路：
直接dp求答案很麻烦，因为限制条件较严格
不妨放宽限制条件
1.首先想到$(n-k) \mod 2=0$，不然就是0
2.升序排两遍，然后对于a[i]，我们记录b数组中所有比它小的元素的个数，记为pos[i]
3.f[i][j]指前i个元素中a数组有j个元素比b数组中元素大，其他i-j个元素位置不确定的方案数，显然我们有
$f[i][j]=f[i-1][j]+f[i-1][j-1]×(pos[i]-j+1)$
初始状态为$f[i][0]=1(i=1,2,...,n)$
4. $g[i]$表示到第n个元素，糖果>药片的组数为i，且剩下的为糖果<药片的方案数，显然$g[i]$是可以通过f[i][n]和之前算出的g减出来的
$\displaystyle g[i]=f[n][i]×(n-i)!-\sum^n_{j=i+1}g[j]×C_j^i$

#include<cstdio>
#include<algorithm>
#define N 2005
#define mo 1000000009
#define LL long long
using namespace std;
int n,k;
int fac[N],inv[N],a[N],b[N],pos[N],f[N][N],g[N];
int C(int m,int n){return (LL)fac[n]*inv[m]%mo*inv[n-m]%mo;}
main()
{
    scanf("%d%d",&n,&k);
    if (n-k&1) return puts("0"),0;
    inv[0]=fac[0]=1;
    for (int i=1;i<=n;++i) fac[i]=(LL)fac[i-1]*i%mo;
    inv[1]=1;
    for (int i=2;i<=n;++i) inv[i]=(LL)(mo-mo/i)*inv[mo%i]%mo;
    for (int i=2;i<=n;++i) inv[i]=(LL)inv[i-1]*inv[i]%mo;
    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);
    int now=1;
    for (int i=1;i<=n;++i)
    {
        while (now<=n&&b[now]<a[i]) ++now;
        --now;
        pos[i]=now;
    }
    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]+(LL)f[i-1][j-1]*(pos[i]-j+1)%mo)%mo;
    for (int i=n;i>=(n+k)/2;--i)
    {
        g[i]=(LL)f[n][i]*fac[n-i]%mo;
        for (int j=i+1;j<=n;++j)
            g[i]=(g[i]-(LL)g[j]*C(i,j))%mo;
    }
    printf("%d\n",(g[n+k>>1]+mo)%mo);
}