【BZOJ3622】已经没有什么好害怕的了

Description

Input

Output

Sample Input

4 2
5 35 15 45
40 20 10 30

Sample Output

4

HINT

输入的2*n个数字保证全不相同。

还有输入应该是第二行是糖果,第三行是药片

题解:好吧这题不是神题,而是套路题,容斥+DP的套路在很多题中都用到过,不过我虽然知道套路,却被这题的第一步卡住了。

我们将两个序列从小到大排序。

好吧这步看起来可能很水,正常人看到无序的序列都会先想到排序,然而为什么要两个都从小到大排呢?排序的意义何在?一会讲。

题中给的是(糖果>药片)=(药片>糖果)+k,我们只需要用一点小学数学的知识就能算出(糖果>药片)的具体数量。下面是套路部分:我们设f[i][j]表示枚举到前i个,(糖果>药片)数至少为j的方案数。然后我们预处理出最大的k,满足药片k<糖果i,那么方程就是:

$f[i][j]=f[i-1][j]+f[i-1][j-1]*(k-j+1)$

现在知道排序的意义了吧?因为后面的i的决策区间一定包含前面的i的决策区间,所以只需要将其减去即可防止重复。

那么最后答案是什么呢?套路:ans=至少有k个的方案数-至少有k+1个的+至少有k+2个的方案数-。。。

所以答案=$\sum\limits_{j=k}^nf[n][j](-1)^{j-k}C_j^k$

P.S:这里解释一下容斥系数是$C_j^k$的原理吧,首先我们用f[i]代表上面的f[n][i],用g[i]代表恰好有i个符合条件的方案数。根据定义有:$f[i]=\sum\limits_{j=i}^ng[i]C_j^i$。所以我们上面得到的最终式子可以转化成:$\sum\limits_{j=k}^nf[j](-1)^{j-k}C_j^k=\sum\limits_{i=k}^ng[i]\sum\limits_{j=k}^i(-1)^{j-k}C_i^jC_j^k\\=\sum\limits_{i=k}^ng[i]\sum\limits_{j=k}^i(-1)^{j-k}C_i^kC_{i-k}^{j-k}\\=\sum\limits_{i=k}^ng[i]C_i^k\sum\limits_{j=k}^i(-1)^{j-k}C_{i-k}^{j-k}$

根据组合恒等式后面那个东西只在i=k的时候=1,其余时候=0,所以我们就得到了g[k]。

#include <cstdio>
#include <cstring>
#include <iostream>
#include <algorithm>
using namespace std;
typedef long long ll;
const ll P=1000000009;
int n,m;
ll ans;
ll f[2010][2010],c[2010][2010],jc[2010];
int a[2010],b[2010];
inline int rd()
{
int ret=0,f=1; char gc=getchar();
while(gc<'0'||gc>'9') {if(gc=='-') f=-f; gc=getchar();}
while(gc>='0'&&gc<='9') ret=ret*10+(gc^'0'),gc=getchar();
return ret*f;
}
int main()
{
n=rd(),m=rd();
if((n^m)&1)
{
printf("0");
return 0;
}
m=(n+m)>>1;
int i,j,k;
for(jc[0]=i=1;i<=n;i++) jc[i]=jc[i-1]*i%P;
for(i=1;i<=n;i++) a[i]=rd();
for(i=1;i<=n;i++) b[i]=rd();
for(i=0;i<=n;i++)
{
c[i][0]=1;
for(j=1;j<=i;j++) c[i][j]=(c[i-1][j-1]+c[i-1][j])%P;
}
sort(a+1,a+n+1),sort(b+1,b+n+1);
f[0][0]=1;
for(i=1;i<=n;i++)
{
for(k=1;k<=n&&b[k]<a[i];k++);
k--;
for(j=1;j<=i;j++) f[i][j]=(f[i-1][j]+f[i-1][j-1]*max(k-j+1,0))%P;
f[i][0]=f[i-1][0];
}
ll tmp=1;
for(i=m;i<=n;i++) f[n][i]=(f[n][i]*jc[n-i])%P,ans=(ans+tmp*f[n][i]*c[i][m]%P+P)%P,tmp=-tmp;
printf("%lld",(ans+P)%P);
return 0;
}

【BZOJ3622】已经没有什么好害怕的了 容斥+DP的更多相关文章

  1. [BZOJ3622]已经没有什么好害怕的了(容斥DP)

    给定两个数组a[n]与b[n](数全不相等),两两配对,求“a比b大”的数对比“b比a大”的数对个数多k的配对方案数. 据说做了这题就没什么题好害怕的了,但感觉实际上这是一个套路题,只是很难想到. 首 ...

  2. BZOJ3622 已经没有什么好害怕的了 二项式反演+DP

    题目传送门 https://lydsy.com/JudgeOnline/problem.php?id=3622 题解 首先显然如果 \(n - k\) 为奇数那么就是无解.否则的话,"糖果& ...

  3. BZOJ3622(容斥+dp)

    思路 "恰k个"考虑求至少k.k+1--个容斥 题面说所有数字都不同,可以将所求转化为糖比药多的组数恰为\((n+k)/2\)的方案数 \(f[i][j]\)数组我觉得更好的理解方 ...

  4. P4859-已经没有什么好害怕的了【容斥,dp】

    正题 题目链接:https://www.luogu.com.cn/problem/P4859 题目大意 两个长度为\(n\)的序列\(a,b\)两两匹配,求\(a_i>b_i\)的组数比\(a_ ...

  5. 【BZOJ3622】已经没有什么好害怕的了(动态规划,容斥)

    [BZOJ3622]已经没有什么好害怕的了(动态规划,容斥) 题面 BZOJ 题解 很明显的,这类问题是要从至少变成恰好的过程,直接容斥即可. 首先我们要求的是(糖果>药片)=(药片>糖果 ...

  6. BZOJ3622 已经没有什么好害怕的了 【dp + 二项式反演】

    题目链接 BZOJ3622 题解 既已开题 那就已经没有什么好害怕的了 由题目中奇怪的条件我们可以特判掉\(n - k\)为奇数时答案为\(0\) 否则我们要求的就是糖果大于药片恰好有\(\frac{ ...

  7. [bzoj3622]已经没有什么好害怕的了_动态规划_容斥原理

    bzoj-3622 已经没有什么好害怕的了 题目大意: 数据范围:$1\le n \le 2000$ , $0\le k\le n$. 想法: 首先,不难求出药片比糖果小的组数. 紧接着,我开始的想法 ...

  8. bzoj3622已经没有什么好害怕的了

    bzoj3622已经没有什么好害怕的了 题意: 给n个数Ai,n个数Bi,将Ai中的数与Bi中的数配对,求配对Ai比Bi大的比Bi比Ai大的恰好有k组的方案数.n,k≤2000 题解: 蒟蒻太弱了只能 ...

  9. bzoj3622已经没有什么好害怕的了 dp+组合+容斥(?)

    3622: 已经没有什么好害怕的了 Time Limit: 10 Sec  Memory Limit: 256 MBSubmit: 1033  Solved: 480[Submit][Status][ ...

随机推荐

  1. 【服务器防护】centos iptables 防火墙设置 mac过滤

    1.阻止MAC地址为XX:XX:XX:XX:XX:XX主机的所有通信: iptables -A INPUT -s 192.168.1.21 -m mac --mac-source XX:XX:XX:X ...

  2. centos7下mysqldump+crontab自动备份数据库

    1.创建文件夹(存放备份数据) mkdir /bak mkdir /bak/mysqldata 2.编写脚本 vi /usr/sbin/bakmysql.sh 脚本内容如下 DATE=`date +% ...

  3. 转 redis使用场景 简介

    Redis实战(五) 聊聊Redis使用场景 发表于 2016-11-21 | 数据存储 | Redis 文章目录 1. 使用场景说明 1.1. 计数器 1.2. 排行榜 1.3. 用于存储时间戳 1 ...

  4. unity, Root Motion

    (引自:http://tieba.baidu.com/p/4323644080) 然后详细看了下这个文档:http://docs.unity3d.com/Manual/RootMotion.html ...

  5. Memcached真的过时了吗?

    Memcached真的过时了吗? 这两年Redis火得可以,Redis也常常被当作Memcached的挑战者被提到桌面上来.关于Redis与Memcached的比较更是比比皆是.然而,Redis真的在 ...

  6. Linux 内存中的Cache,真的能被回收么?

    您真的了解Linux的free命令么? 在Linux系统中,我们经常用free命令来查看系统内存的使用状态.在一个RHEL6的系统上,free命令的显示内容大概是这样一个状态: 这里的默认显示单位是k ...

  7. keepalived virtual_router_id 44

    在同一局或网内如果有多个keepalived 的话 virtuall_router_id 44 (不能相同,但同一对,是一定相同)

  8. struts2设置加载非默认路径的struts.xml文件解决方案

    方案一:   首先我们要明白struts2加载配置文件都是从它自己的jar包和\WEB-INF\classes两个默认的位置加载的,如果你想改变strusts2中的文件的默认加载路径,可以在web项目 ...

  9. 从零开始,跟我一起做jblog项目(二)Maven

    从零开始,跟我一起做jblog项目(一)引言 从零开始,跟我一起做jblog项目(二)Maven maven是一个项目管理工具,尤其适用于JAVA世界 在jblog的开发前期,还没有系统使用过mave ...

  10. windows7环境下使用pip安装MySQLdb

    1.首先,需要确定你已经安装了pip.在Python2.7的安装包中,easy_install.py和pip都是默认安装的.可以在Python的安装目录先确认,如果\Python27\Scripts里 ...