【AGC030F】Permutation and Minimum DP
题目大意
有一个长度为序列 \(a\),其中某些位置的值是 \(-1\)。
你要把 \(a\) 补成一个排列。
定义 \(b_i=\min(a_{2i-1},a_{2i})\),求有多少种可能的 \(b\)。
\(n\leq 300\)
题解
如果 \(a_{2i-1}\) 和 \(a_{2i}\) 都有值,就把这两个位置扔掉。
记 \(c_i\) 表示 \(i\) 这个值是否在初始的 \(a\) 中。
从后往前DP。记 \(f_{i,j,k}\) 表示已经处理完了 \(i\) 后面的数,有多少个 \(j>i,c_j=1\) 的数匹配的是 \(\leq i\) 的数,有多少个 \(j>i,c_j=0\) 的数匹配的是 \(\leq i\) 的数。
如果 \(c_i=1\) 且往后匹配的是 \(c_j=0\),那么方案数为 \(1\)。(因为 \(\min=i\))
如果 \(c_i=0\) 且往后匹配的是 \(c_j=0\),那么先暂定方案数为 \(1\)。(因为暂时不能确定 \(i\) 填在哪个位置。)记这种匹配对数为 \(cnt\)。
如果 \(c_i=0\) 且往后匹配的是 \(c_j=1\),那么方案数为 \(j\)。(因为可以确定 \(i\) 填在哪个位置。)
最后方案数要乘上 \(cnt!\),因为这些位置的 \(b\) 可以随便交换。
时间复杂度:\(O(n^3)\)
代码
#include<cstdio>
#include<cstring>
#include<algorithm>
#include<cstdlib>
#include<ctime>
#include<functional>
#include<cmath>
#include<vector>
#include<assert.h>
//using namespace std;
using std::min;
using std::max;
using std::swap;
using std::sort;
using std::reverse;
using std::random_shuffle;
using std::lower_bound;
using std::upper_bound;
using std::unique;
using std::vector;
typedef long long ll;
typedef unsigned long long ull;
typedef double db;
typedef std::pair<int,int> pii;
typedef std::pair<ll,ll> pll;
void open(const char *s){
#ifndef ONLINE_JUDGE
char str[100];sprintf(str,"%s.in",s);freopen(str,"r",stdin);sprintf(str,"%s.out",s);freopen(str,"w",stdout);
#endif
}
void open2(const char *s){
#ifdef DEBUG
char str[100];sprintf(str,"%s.in",s);freopen(str,"r",stdin);sprintf(str,"%s.out",s);freopen(str,"w",stdout);
#endif
}
int rd(){int s=0,c,b=0;while(((c=getchar())<'0'||c>'9')&&c!='-');if(c=='-'){c=getchar();b=1;}do{s=s*10+c-'0';}while((c=getchar())>='0'&&c<='9');return b?-s:s;}
void put(int x){if(!x){putchar('0');return;}static int c[20];int t=0;while(x){c[++t]=x%10;x/=10;}while(t)putchar(c[t--]+'0');}
int upmin(int &a,int b){if(b<a){a=b;return 1;}return 0;}
int upmax(int &a,int b){if(b>a){a=b;return 1;}return 0;}
const ll p=1000000007;
const int N=310;
void add(ll &a,ll b)
{
a=(a+b)%p;
}
int n;
int a[2*N];
ll f[2*N][N][N];
int b[2*N];
int c[2*N];
int t;
int main()
{
open2("f");
scanf("%d",&n);
for(int i=1;i<=2*n;i++)
scanf("%d",&a[i]);
int cnt=0;
for(int i=1;i<=2*n;i+=2)
if(a[i]!=-1&&a[i+1]!=-1)
b[a[i]]=b[a[i+1]]=2;
else if(a[i]!=-1)
b[a[i]]=1;
else if(a[i+1]!=-1)
b[a[i+1]]=1;
else
cnt++;
for(int i=1;i<=2*n;i++)
if(b[i]==1)
c[++t]=1;
else if(b[i]==0)
c[++t]=2;
f[t][0][0]=1;
for(int i=t;i>=1;i--)
for(int j=0;j<=t&&j<=n;j++)
for(int k=0;k<=t&&k<=n;k++)
if(c[i]==1)
{
add(f[i-1][j+1][k],f[i][j][k]);
if(k)
add(f[i-1][j][k-1],f[i][j][k]);
}
else
{
add(f[i-1][j][k+1],f[i][j][k]);
if(k)
add(f[i-1][j][k-1],f[i][j][k]);
if(j)
add(f[i-1][j-1][k],f[i][j][k]*j);
}
ll ans=f[0][0][0];
for(int i=1;i<=cnt;i++)
ans=ans*i%p;
printf("%lld\n",ans);
return 0;
}
【AGC030F】Permutation and Minimum DP的更多相关文章
- 【agc030f】Permutation and Minimum(动态规划)
[agc030f]Permutation and Minimum(动态规划) 题面 atcoder 给定一个长度为\(2n\)的残缺的排列\(A\),定义\(b_i=min\{A_{2i-1},A_{ ...
- 【AGC030F】Permutation and Minimum(DP)
题目链接 题解 首先可以想到分组后,去掉两边都填了数的组. 然后就会剩下\((-1,-1)\)和\((-1,x)\)或\((x,-1)\)这两种情况 因为是最小值序列的情况数,我们可以考虑从大到小填数 ...
- 【BZOJ4712】洪水(动态dp)
[BZOJ4712]洪水(动态dp) 题面 BZOJ 然而是权限题QwQ,所以粘过来算了. Description 小A走到一个山脚下,准备给自己造一个小屋.这时候,小A的朋友(op,又叫管理员)打开 ...
- 【题解】Jury Compromise(链表+DP)
[题解]Jury Compromise(链表+DP) 传送门 题目大意 给你\(n\le 200\)个元素,一个元素有两个特征值,\(c_i\)和\(d_i\),\(c,d \in [0,20]\), ...
- 【题解】Making The Grade(DP+结论)
[题解]Making The Grade(DP+结论) VJ:Making the Grade HNOI-D2-T3 原题,禁赛三年. 或许是我做过的最简单的DP题了吧(一遍过是什么东西) 之前做过关 ...
- 【题解】NOIP2017逛公园(DP)
[题解]NOIP2017逛公园(DP) 第一次交挂了27分...我是不是必将惨败了... 考虑这样一种做法,设\(d_i\)表示从该节点到n节点的最短路径,\(dp(i,k)\)表示从\(i\)节点 ...
- 【题解】284E. Coin Troubles(dp+图论建模)
[题解]284E. Coin Troubles(dp+图论建模) 题意就是要你跑一个完全背包,但是要求背包的方案中有个数相对大小的限制 考虑一个\(c_i<c_j\)的限制,就是一个\(c_i\ ...
- 【LeetCode】153. Find Minimum in Rotated Sorted Array 解题报告(Python)
[LeetCode]153. Find Minimum in Rotated Sorted Array 解题报告(Python) 标签: LeetCode 题目地址:https://leetcode. ...
- 【LeetCode】154. Find Minimum in Rotated Sorted Array II 解题报告(Python)
[LeetCode]154. Find Minimum in Rotated Sorted Array II 解题报告(Python) 标签: LeetCode 题目地址:https://leetco ...
随机推荐
- int16, int32, int64等类型说明
Int16 相当于 short 占2个字节 -32768 ~ 32767 Int32 相当于 int 占4个字节 -2147483648 ~ 2147483647 Int64 ...
- 在表单提交之前做校验-利用jQuery的submit方法
点击表单中的提交按钮,默认就会提交表单,如果要在表单提交之前做一些校验,那么就可以用jQuery的submit方法. 也就是jQuery的submit的方法执行顺序,在表单的提交之前.用法如下: $( ...
- Ext.isIterable
Ext.isIterable用于判断传入的参数是否为可迭代的 在这4种情况下,函数返回true 1:数组2:函数参数arguments3:HTML collections : NodeList4:HT ...
- h5与c3权威指南笔记--css3新属性选择器
[att*=val] 选择所有att属性值中包含val的.只要包含val值,不论val值在属性值的前面还是中间还是后面~ <style> div[class*=div]{ color: r ...
- 关于asp.net执行exe程序时权限不够的解决办法(2015.04.17更新)
一,本文背景 长话短说:asp.net项目中需要用到PDF转换成SWF文件,用户上传后自动调用pdf2swf.exe转换. 但有个问题,执行时权限不够,导致一直报错(滚动条一直在往下滚,刷屏中),见下 ...
- 前端入门25-福音 TypeScript
声明 本篇内容摘抄自以下来源: TypeScript 中文网 只梳理其中部分知识点,更多更详细内容参考官网. 正文-TypeScript 今天来讲讲有 Java 基础转 JavaScript 的福音: ...
- SAP HUM 拆包之后的HU号码依旧存在
比如HU 194811210666,已经被执行了Unpack操作. HUMO看这个HU号码, 执行, 可以发现,这个HU是空的,里面什么都没有. 双击HU号码,进入HU显示界面, 系统状态竟然是PHE ...
- 自定义修改Anaconda Jupyterlab Home目录
自定义修改Anaconda Jupyterlab Home目录 最近在使用Anaconda学习数据分析和机器学习,会使用到Jupyter,但是他默认目录是用户的目录,我并没有习惯将项目和资料放在C盘, ...
- Skyline Te Pro二次开发技能总结
前两天项目开发中,忽然一个Imagelabel的参数不会调了,但是前段时间可是很熟悉的.好吧,好记性不如烂笔头! 1. 模型弹出窗调试 这里的模型弹出框指涉及到模型操作的,比如监听模型选定事件.根据窗 ...
- python类的继承、封装和多态
摘自https://www.cnblogs.com/evablogs/p/6724477.html 继承 1 2 3 4 5 6 7 8 class Person(object): def _ ...