bzoj3622

顯然蒟蒻是抄襲別人的。。。   
 #include <cstdio>  
    #include <cstring>  
    #include <iostream>  
    #include <algorithm>  
    #define M 2020  
    #define MOD 1000000009  
    using namespace std;  
    int n,k,s;  
    int a[M],b[M],next[M];  
    long long C[M][M],fac[M];  
    long long f[M][M],g[M];  
    //f[i][j]表示前i個糖果中已經確定有j組糖>藥的方案數  
    int main()  
    {  
        int i,j;  
        cin>>n>>k;  
        if(n+k&1)  
        {  
            cout<<0<<endl;  
            return 0;  
        }  
        s=n+k>>1;  
        for(i=1;i<=n;i++)  
            scanf("%d",&a[i]);  
        for(i=1;i<=n;i++)  
            scanf("%d",&b[i]);  
        sort(a+1,a+n+1);  
        sort(b+1,b+n+1);  
      
        for(j=1,i=1;i<=n;i++)  
        {  
            for(;j<=n&&b[j]<a[i];j++);  
            next[i]=j-1;  
        }  
      
        f[0][0]=1;  
        for(i=1;i<=n;i++)  
        {  
            f[i][0]=1;  
            for(j=1;j<=i;j++)  
                f[i][j]=(f[i-1][j]+f[i-1][j-1]*max(next[i]-(j-1),0))%MOD;  
        }  
    //  這個算出來顯然是用重複計數的，所以下面使用容斥 
        for(i=0;i<=n;i++)  
        {  
            C[i][0]=1;  
            for(j=1;j<=i;j++)  
                C[i][j]=(C[i-1][j]+C[i-1][j-1])%MOD;  
        }  
          
        for(fac[0]=1,i=1;i<=n;i++)  
            fac[i]=fac[i-1]*i%MOD;  
      
        for(i=n;i>=s;i--)  
        {  
            g[i]=f[n][i]*fac[n-i]%MOD;  //這裏顯然是有>=i種的方案數
            for(j=i+1;j<=n;j++)  
                (g[i]+=MOD-g[j]*C[j][i]%MOD)%=MOD;  //倒推重複的部分
//(n-i)! 是枚舉後面 n-i 可能的方案，個g[j]*C(j, i) 表示g[i] 中實際有 j 組藥片數 > 糖果數卻被計入 f[n][i] 的數量
 }  
      
        cout<<g[s]<<endl;  
        return 0;  
      
    }