[bzoj3622]已经没有什么好害怕的了

最新推荐文章于 2018-10-07 08:37:00 发布

FZHvampire

最新推荐文章于 2018-10-07 08:37:00 发布

阅读量947

点赞数 1

分类专栏：容斥动态规划

本文链接：https://blog.csdn.net/FZHvampire/article/details/50771291

版权

动态规划同时被 2 个专栏收录

35 篇文章 0 订阅

订阅专栏

容斥

4 篇文章 0 订阅

订阅专栏

3622: 已经没有什么好害怕的了

Time Limit: 10 Sec Memory Limit: 256 MB
Submit: 197 Solved: 113
[Submit][Status][Discuss]
Description

这里写图片描述

Input

这里写图片描述

Output

这里写图片描述

Sample Input

4 2

5 35 15 45

40 20 10 30

Sample Output

4
HINT

这里写图片描述

输入的2*n个数字保证全不相同。
还有输入应该是第二行是糖果，第三行是药片

这里写图片描述

g[i]在求的时候不能直接把后面的g[j]全加起来的原因是在计算f[][]*()!的这个式子的时候，有些数列会被算很多遍，所以要容斥一下。

#include<iostream>
#include<cstdio>
#include<cstring>
#include<algorithm>
using namespace std;
#define LL long long
#define D 1000000009
const int N=2010;
int n,m,a[N],b[N];
LL p[N],f[N][N],g[N],jc[N],c[N][N];
int main(){
    int i,j;
    scanf("%d%d",&n,&m);
    if((m+n)%2){printf("0\n");return 0;}
    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=0,i=1;i<=n;++i){
        while(j+1<=n&&b[j+1]<a[i]) j+=1;
        p[i]=j;
    }
    f[0][0]=1;
    for(i=1;i<=n;++i)
      for(f[i][0]=f[i-1][0],j=1;j<=i;++j) f[i][j]=(f[i-1][j]+(f[i-1][j-1]*max(0LL,(p[i]-(LL)j+1)))%D)%D;
    for(jc[1]=1,i=2;i<=n;++i) jc[i]=jc[i-1]*(LL)i%D;
    for(i=0;i<=n;++i){
        c[i][0]=c[i][i]=1;
        for(j=1;j<i;++j) c[i][j]=(c[i-1][j-1]+c[i-1][j])%D;
    }
    for(i=n;i;--i){
        g[i]=f[n][i]*jc[n-i]%D;
        for(j=i+1;j<=n;++j){
            g[i]-=g[j]*c[j][i];
            g[i]=(g[i]%D+D)%D;
        }
    }
    printf("%lld\n",g[(n+m)/2]);
}