原文链接https://www.cnblogs.com/zhouzhendong/p/9276479.html

题目传送门 - BZOJ3622

题意

  给定两个序列 $a,b$ ,各包含 $n$ 个数字。

  现在给 $a$ 中元素与 $b$ 中元素配对。问使得所有配对中 $a_?>b_?$ 的个数比 $a_?<b_?$ 的个数恰好多 $k$ 的方案总数。

  答案对 $10^9+9$ 取模,保证 $a$ 和 $b$ 中的所有数字互不相同。

  $n\leq 2000$

题解

  首先闭着眼睛排个序。

  然后,我们发现一个非常好的性质:如果要使得 $a_i$ 的配对 $b_?$ 满足 $a_i>b_?$ 那么,满足 $1\leq ?\leq c_i$ ,其中 $c_i$ 表示在 $1\leq j\leq n$ 中满足 $b_j<a_i$ 的最大 $j$ 值。

  我们先把问题转化一下,变成求恰好得到 $x$ 个 $a_?>b_?$ 的方案数。那么显然 $x=\cfrac{n+k}{2}$ 。如果 $x$ 不是整数,那么答案显然是 $0$ 。

  下面说的满足条件是指 $a_?>b_?$ 。

  于是我们考虑动态规划。

  记 $f_{i,j}$ 表示在 $a_1,a_2,...,a_i$ 中选择 $j$ 个,匹配比他小的 $b_i$ 的方案数。

  那么我们可以轻松确定转移方程:

$$f_{i,j}=f_{i-1,j}+f_{i-1,j-1}\times \max(c_i-j+1,0)$$

  我们假设 $f_i=f_{n,i}$ ,并设 $g_i$ 为配对所有的元素,恰好有 $i$ 个满足条件的方案数。

  于是下式成立:

$$f_i=\sum_{j=i}^n\binom{j}{i}g_j$$  

  回忆一下 $f$ 的定义,$f_i$ 表示在 $n$ 个里面选择 $i$ 个匹配 $i$ 个满足条件的方案数。

  考虑所有的最终情况。

  满足条件个数为 $j$ 的总共有 $g_j$ 种。每一个这样的方案中选择的 $j$ 个数中,任选 $i$ 个,所得到的方案,都会对 $f_i$ 有贡献。

  所以,$f_i=\sum_{j=0}^{n}\binom{j}{i}g_j=\sum_{j=i}^n\binom{j}{i}g_j$ 。

  于是我们来用 $f$ 表示 $g$ 。

$$\begin{eqnarray*}g_k&=&\sum_{i=k}^{n}g_i\binom{i}{k}\sum_{j=k}^{i}(-1)^{j-k}\binom{i-k}{j-k} \\&=&\sum_{i=k}^{n}g_i\sum_{j=k}^{i}(-1)^{j-k}\binom{i}{k}\binom{i-k}{j-k}\\&=&\sum_{i=k}^{n}g_i\sum_{j=k}^{i}(-1)^{j-k}\binom{i}{j}\binom{j}{k}\\&=&\sum_{i=k}^{n}(-1)^{i-k}\binom{i}{k}\sum_{j=i}^{n}g_j\binom{j}{i}\\&=&\sum_{i=k}^{n}f_i(-1)^{i-k}\binom{i}{k}\end{eqnarray*}$$

  然后按照这个公式算一算就可以了。

代码

#include <bits/stdc++.h>
using namespace std;
const int N=2005,mod=1e9+9;
int n,k,a[N],b[N],c[N];
int C[N][N],Fac[N];
int f[N][N];
int main(){
scanf("%d%d",&n,&k);
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);
for (int i=1;i<=n;i++)
for (c[i]=c[i-1];c[i]<n&&b[c[i]+1]<a[i];c[i]++);
for (int i=Fac[0]=1;i<=n;i++)
Fac[i]=1LL*Fac[i-1]*i%mod;
for (int i=0;i<=n;i++)
C[i][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])%mod;
for (int i=0;i<=n;i++)
f[i][0]=1;
for (int i=1;i<=n;i++)
for (int j=1;j<=n;j++)
f[i][j]=(1LL*f[i-1][j-1]*max(c[i]-j+1,0)+f[i-1][j])%mod;
int ans=0,x=(n+k)/2;
if (x*2!=n+k){
puts("0");
return 0;
}
for (int i=x;i<=n;i++)
ans=(1LL*f[n][i]*C[i][x]%mod*Fac[n-i]*((i-x+1)%2*2-1)+ans)%mod;
ans=(ans+mod)%mod;
printf("%d",ans);
return 0;
}

  

BZOJ3622 已经没有什么好害怕的了 动态规划 容斥原理 组合数学的更多相关文章

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

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

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

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

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

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

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

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

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

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

  6. BZOJ3622 已经没有什么好害怕的了(动态规划+容斥原理)

    显然可以转化为一个阶梯状01矩阵每行每列取一个使权值和为k的方案数.直接做不可做,考虑设f[i][j]为前i行权值和至少为j,即在其中固定了j行选1的方案数.设第i行从1~a[i]列都是1且a[i]+ ...

  7. BZOJ3622 已经没有什么好害怕的了

    Description Input Output Sample Input 4 2 5 35 15 45 40 20 10 30 Sample Output 4 HINT 输入的2*n个数字保证全不相 ...

  8. 【BZOJ3622】已经没什么好害怕的了 容斥原理+dp

    Description Input Output Sample Input 4 2 5 35 15 45 40 20 10 30 Sample Output 4 HINT 输入的2*n个数字保证全不相 ...

  9. 洛谷 P4859 && BZOJ3622: 已经没有什么好害怕的了

    题目描述 给出 \(n\) 个数 \(a_i\)​ ,以及 \(n\) 个数 \(b_i\)​ ,要求两两配对使得 \(a>b\) 的对数减去 \(a<b\) 的对数等于 \(k\) . ...

随机推荐

  1. [JavaScript]使用ArrayBuffer和Blob编辑二进制流

    Blob()构造方法返回一个新的Blob对象. 内容是包含参数array的二进制字节流. 语法 var aBlob = new Blob( array, options ); 参数 array is ...

  2. 程序包管理dpkg、apt-get、服务端openssh-server与客户端Xshell设置及lrzsz安装使用

    一.程序包管理器 dpkg.apt-get 1.dpkg 安装:sudo dpkg -i cmatrix_1.2a-5build3_amd64.deb 卸载:sudo dpkg -r cmatrix ...

  3. 彻底搞懂字符集编码:ASCII,Unicode 和 UTF-8

    一.ASCII 码 我们知道,计算机内部,所有信息最终都是一个二进制值.每一个二进制位(bit)有0和1两种状态,因此八个二进制位就可以组合出256种状态,这被称为一个字节(byte).也就是说,一个 ...

  4. centos6.5 有趣但是没有用的linux命令

    小火车 get http://dl.fedoraproject.org/pub/epel/6/x86_64/epel-release-6-8.noarch.rpm rpm -ivh epel-rele ...

  5. 弃 Java 而使用 Kotlin 的你后悔了吗?| kotlin将会是最好的开发语言

    自从 2011 年发布以来,Kotlin 凭借强大的功能在开发者中的欢迎程度与日俱增.且在一年前,Google 宣布 Kotlin 正式成为 Android 官方开发语言,由此引发了从 Java 迁移 ...

  6. 【MySql】delete用法

    delete 语句用于删除表中的数据, 基本用法为: delete from 表名称 where 删除条件; 以下是在表 students 中的实例: 删除 id 为 3 的行: delete fro ...

  7. Oracle imp exp 导入导出 执行脚本

    一:用命令 imp/exp 的方式进行数据的导入和导出 一:文件后缀名: 二:oracle  导出 exp 命令 echo 开始备份数据库 if not exist D:\oracle_bak\fil ...

  8. Confluence 6 发送 Confluence 通知到其他 Confluence 服务器

    你可以配置 Confluence 服务器向其他的 Confluence 服务器发送消息.在这种情况下,Confluence 服务器将不会显示 workbox. 希望发送消息到其他 Confluence ...

  9. Confluence 6 后台中的选择站点首页

    后台中的选择站点首页选择项. https://www.cwiki.us/display/CONFLUENCEWIKI/Configuring+the+Site+Home+Page

  10. Confluence 6 自定义主面板

    主面板(dashboard)是你 Confluence 站点的默认载入页面.这个页面能够给用户能够找到其他页面的所有必须的工具,重新进入未完成的工作或者快速导航到喜欢的空间和页面 站点的欢迎信息将会在 ...