Atcoder Grand Contest 030 F - Permutation and Minimum(DP)
12 天以前做的题了,到现在才补/yun
做了一晚上+一早上终于 AC 了,写篇题解纪念一下
首先考虑如果全是 \(-1\) 怎么处理。由于我们不关心每个 pair
中的 max
是多少,并且显然每个 pair
的最小值都是两两不同的,因此我们可以考虑有多少个最小值组成的集合,然后答案乘上 \(n!\) 即可。而显然如果我们将“作为某个 pair
的最小值”的位置放上一个左括号,“不作为某个 pair
的最小值”的位置放上一个右括号,那么一个集合符合条件当且仅当这个集合对应的括号序列是一个合法括号序列,卡特兰数算算即可。
接下来思考存在某些数固定下来的情况怎么处理,显然我们可以将所有 pair
分为三类:\((-1,-1),(-1,x),(x,y)\),第三类的最小值显然是固定的,我们可以将这样的 \((x,y)\) 中的 \(x,y\) 设为访问过,然后只考虑那些没有访问过的元素。显然带上“某些值固定”这个限制条件后就无法直接套用组合数了,并且显然这个数据范围也不是让你纯组合数学就过的,因此考虑一个 DP。受到上面括号序列这个思想的启发,我们考虑这样一个过程:
- 我们考虑有一个长度为 \(m\) 的括号序列,你要在里面填上 \(m/2\) 个左括号和 \(m/2\) 个右括号,其中 \(m=2n-2\times\text{两个位置都已经确定了的 pair 的对数}\)
- 括号有方括号和圆括号之分,给定 \(k\) 个位置,填上的括号必须为方括号,另外 \(m-k\) 个位置填上的括号必须为圆括号。
- 填好左右括号之后,你要对于所有右方括号,在前面找一个左圆括号匹配,并且满足去掉这些已经匹配过的括号之后,剩下的括号在不分方括号和圆括号的情况下仍能组成一个合法括号序列
两个方案不同当且仅当括号序列不同或者存在某个右方括号,满足与其匹配的左圆括号位置不同。
够清晰的了吧
我们考虑 DP,事实证明正着 DP 是不好处理左圆括号与右方括号匹配这件事情的,因此考虑倒着 DP,\(dp_{i,j,k}\) 表示钦定了第 \(i\) 个位置到第 \(m\) 个位置上填括号的方案,还剩 \(j\) 个 ‘)’ 和 \(k\) 个 ‘]’ 没有匹配的方案数。那么考虑第 \(i-1\) 个位置:
- 如果它只能填方括号:
- 如果填左括号,那么它只能与圆括号匹配,\(dp_{i,j,k}\to dp_{i-1,j-1,k}\)
- 如果填右括号,那么未匹配的 ‘)’ 个数加一,\(dp_{i,j,k}\to dp_{i-1,j+1,k}\)
- 如果它只能填圆括号:
- 如果填左括号,那么有两种可能,与圆括号匹配,由于这种情况是不在乎它与哪个圆括号匹配的,因此贡献只有 \(1\),即 \(dp_{i,j,k}\to dp_{i-1,j-1,k}\),或者它也可以与方括号匹配,显然每个左圆括号与哪个右方括号匹配都是不同的,因此有 \(k\) 种可能,\(dp_{i,j,k}·k\to dp_{i-1,j,k-1}\)
- 如果填右括号,那么未匹配的 ‘)’ 个数加一,\(dp_{i,j,k}\to dp_{i-1,j,k+1}\)
随便推推即可,最终答案记得乘个 \(((n-m)/2)!\)。时空复杂度均 \(\mathcal O(n^3)\)
const int MAXN=600;
const int MOD=1e9+7;
int n,a[MAXN+5],vis[MAXN+5],id[MAXN+5],m=0,is[MAXN+5];
int dp[MAXN+5][MAXN+5][MAXN+5];
void add(int &x,int v){((x+=v)>=MOD)&&(x-=MOD);}
int main(){
scanf("%d",&n);
for(int i=1;i<=n<<1;i++) scanf("%d",&a[i]);
for(int i=1;i<=n<<1;i+=2) if(~a[i]&&~a[i+1]) vis[a[i]]=vis[a[i+1]]=1;
for(int i=1;i<=n<<1;i++) if(!vis[i]) id[i]=++m;
for(int i=1;i<=n<<1;i+=2){
if(~a[i]&&!~a[i+1]) is[id[a[i]]]=1;
if(!~a[i]&&~a[i+1]) is[id[a[i+1]]]=1;
} dp[m+1][0][0]=1;
// for(int i=1;i<=m;i++) printf("%d\n",is[i]);
for(int i=m;i;i--) for(int j=0;j<=m-i;j++) for(int k=0;k+j<=m-i;k++){
if(is[i]){
add(dp[i][j][k+1],dp[i+1][j][k]);
if(j) add(dp[i][j-1][k],dp[i+1][j][k]);
} else {
add(dp[i][j+1][k],dp[i+1][j][k]);
if(k) add(dp[i][j][k-1],1ll*dp[i+1][j][k]*k%MOD);
if(j) add(dp[i][j-1][k],dp[i+1][j][k]);
} //printf("%d %d %d %d\n",i,j,k,dp[i][j][k]);
} int res=dp[1][0][0];
int cc=0;for(int i=1;i<=(n<<1);i+=2) cc+=((!~a[i])&&(!~a[i+1]));
for(int i=1;i<=cc;i++) res=1ll*res*i%MOD;
printf("%d\n",res);
return 0;
}
/*
6
-1 -1 -1 -1 -1 -1 -1 -1 -1 -1 -1 -1
*/
Atcoder Grand Contest 030 F - Permutation and Minimum(DP)的更多相关文章
- Atcoder Grand Contest 034 F - RNG and XOR(FWT)
Atcoder 题面传送门 & 洛谷题面传送门 tsc 考试前 A 的题了,结果到现在才写这篇题解--为了 2mol 我已经一周没碰键盘了,现在 2mol 结束算是可以短暂的春天 短暂地卷一会 ...
- Atcoder Grand Contest 024 E - Sequence Growing Hard(dp+思维)
题目传送门 典型的 Atcoder 风格的计数 dp. 题目可以转化为每次在序列中插入一个 \([1,k]\) 的数,共操作 \(n\) 次,满足后一个序列的字典序严格大于前一个序列,问有多少种操作序 ...
- Atcoder Grand Contest 003 F - Fraction of Fractal(矩阵乘法)
Atcoder 题面传送门 & 洛谷题面传送门 Yet another AGC F,然鹅这次就没能自己想出来了-- 首先需注意到题目中有一个条件叫做"黑格子组成的连通块是四联通的&q ...
- JZOJ5405 & AtCoder Grand Contest 001 F. Permutation
题目大意 给出一个长度为\(n\)的排列\(P\)与一个正整数\(k\). 你需要进行如下操作任意次, 使得排列\(P\)的字典序尽量小. 对于两个满足\(|i-j|>=k\) 且\(|P_i- ...
- Atcoder Grand Contest 016 F - Games on DAG(状压 dp)
洛谷题面传送门 & Atcoder 题面传送门 如何看待 tzc 补他一个月前做的题目的题解 首先根据 SG 定理先手必输当且仅当 \(\text{SG}(1)=\text{SG}(2)\). ...
- Atcoder Grand Contest 010 C - Cleaning 树贪心(伪)
C - Cleaning 题目连接: http://agc010.contest.atcoder.jp/tasks/agc010_c Description There is a tree with ...
- Atcoder Grand Contest 001 D - Arrays and Palindrome(构造)
Atcoder 题面传送门 洛谷题面传送门 又是道思维题,又是道把我搞自闭的题. 首先考虑对于固定的 \(a_1,a_2,\dots,a_n;b_1,b_2,\dots,b_m\) 怎样判定是否合法, ...
- AtCoder Beginner Contest 122 D - We Like AGC (DP)
D - We Like AGC Time Limit: 2 sec / Memory Limit: 1024 MB Score : 400400 points Problem Statement Yo ...
- 【AGC030F】Permutation and Minimum(DP)
题目链接 题解 首先可以想到分组后,去掉两边都填了数的组. 然后就会剩下\((-1,-1)\)和\((-1,x)\)或\((x,-1)\)这两种情况 因为是最小值序列的情况数,我们可以考虑从大到小填数 ...
随机推荐
- better-scroll快速上手及封装(vue项目)
愿你有诗有梦,有坦荡荡的远方 本文声明:这是一篇学习coderwhy老师的vue2课程的一个笔记,所以本文章是在vue项目中实现,没学过vue的大佬们可以举一反三. 使用场景及介绍 BetterScr ...
- oo第四次博客-UML暨学期总结
一. 本单元两次作业架构设计 这两次作业实际上难度不大,不存在算法上的难题,大部分时间都是用在处理UML图中各个元素的关系上. 第一次UML主要处理UML类图.有UMLclass,UMLinterfa ...
- 2020BUAA软工个人项目作业
2020BUAA软工个人项目作业 17373010 杜博玮 项目 内容 这个作业属于哪个课程 2020春季计算机学院软件工程(罗杰 任健) 这个作业的要求在哪里 个人项目作业 我在这个课程的目标是 学 ...
- 2021.8.3考试总结[NOIP模拟29]
T1 最长不下降子序列 数据范围$1e18$很不妙,但模数$d$只有$150$,考虑从这里突破. 计算的式子是个二次函数,结果只与上一个值有关,而模$d$情况下值最多只有$150$个,就证明序列会出现 ...
- 攻防世界 杂项 1.base64÷4
666C61677B45333342374644384133423834314341393639394544444241323442363041417D 根据题目base64÷4得base16 在线工 ...
- 洛谷 P4587 [FJOI2016]神秘数
大鸽子 llmmkk 正在补8.3号咕掉的题 时隔两个月,再看到这道题,我又是一脸懵,这种思维的培养太重要了 链接: P4587 题意: 给出 \(n\) 个点的序列,\(m\) 次询问区间神秘数. ...
- 用C++实现的数独解题程序 SudokuSolver 2.7 及实例分析
引言:一个 bug 的发现 在 MobaXterm 上看到有内置的 Sudoku 游戏,于是拿 SudokuSolver 求解,随机出题,一上来是个 medium 级别的题: 073 000 060 ...
- sqlldr 导入有中文乱码问题
1.导入成功后,查看导入数据有乱码 2.查看字符集为uft8 select * from v$nls_parameters where PARAMETER like '%NLS_CHARACTERSE ...
- zabbix web管理页面 中文乱码问题
1.在自己电脑上找下图文件,C:\Windows\Fonts 2.上传到 /usr/share/zabbix/assets/fonts/ 目录下 可以看到 graphfont.ttf 是 /etc/a ...
- .Net Minimal Api 介绍
Minimal API是.Net 6中新增的模板,借助C# 10的一些特性以最少的代码运行一个Web服务.本文脱离VS通过VS Code,完成一个简单的Minimal Api项目的开发. 创建项目 随 ...