P1155 双栈排序(二分图染色)

题目描述

Tom最近在研究一个有趣的排序问题。如图所示,通过2个栈S1和S2,Tom希望借助以下4种操作实现将输入序列升序排序。

操作a

如果输入序列不为空,将第一个元素压入栈S1

操作b

如果栈S1不为空,将S1栈顶元素弹出至输出序列

操作c

如果输入序列不为空,将第一个元素压入栈S2

操作d

如果栈S2不为空,将S2栈顶元素弹出至输出序列

如果一个1~n的排列P可以通过一系列操作使得输出序列为1,2,…,(n-1),n,Tom就称P是一个“可双栈排序排列”。例如(1,3,2,4)就是一个“可双栈排序序列”,而(2,3,4,1)不是。下图描述了一个将(1,3,2,4)排序的操作序列:<a,c,c,b,a,d,d,b>

当然,这样的操作序列有可能有几个,对于上例(1,3,2,4),<a,c,c,b,a,d,d,b>是另外一个可行的操作序列。Tom希望知道其中字典序最小的操作序列是什么。

输入输出格式

输入格式:

输入文件twostack.in的第一行是一个整数n。

第二行有n个用空格隔开的正整数,构成一个1~n的排列。

输出格式:

输出文件twostack.out共一行,如果输入的排列不是“可双栈排序排列”,输出数字0;否则输出字典序最小的操作序列,每两个操作之间用空格隔开,行尾没有空格。

输入输出样例

输入样例#1: 复制

4
1 3 2 4
输出样例#1: 复制

a b a a b b a b
输入样例#2: 复制

4
2 3 4 1
输出样例#2: 复制

0
输入样例#3: 复制

3
2 3 1
输出样例#3: 复制

a c a b b d

说明

30%的数据满足: n<=10

50%的数据满足: n<=50

100%的数据满足: n<=1000

题解一:

二分图染色+模拟 
1.首先考虑一个简单情况——单栈排序,显然有这样的一个事实:

a[i]和a[j] 不能压入同一个栈⇔存在一个k,使得i < j < k且a[k] < a[i] < a[j]

对应变量

I

J

K

 

下标

1

2

3

4

a[]

2

3

1

4

显然上面的i,j,k满足关系。

显然上面这个序列不能满足只用一个栈操作这个要求。
时间复杂度为O(n^3).对于n<=1000仍显吃力,对此可以用动态规划的思想,将上述复杂度降到O(n^2)。 
状态:f[i]=min(a[i],a[i+1], … ,a[n]) 
边界条件:f[n+1]=INF; 
状态转移方程:f[i]=min(f[i+1],a[i]); 
于是上述判断就转化为了f[j+1] < a[i] && a[i] < a[j] 
2.扩展到双栈排序: 
如果a[i]和a[j]不能在一个栈内,即连接一条i与j之间的无向边,接下来我们只需要判断这个图是否为二分图 
由于题目中说编号的字典序要尽可能的小,那么就把编号小的尽可能放到stack1 
判断二分图的方法可以采用黑白染色的方式,先从编号小的开始染,第一个顶点染成黑色,相邻的顶点染成不同的颜色,如果发现黑白冲突,那么说明这个图不是一个二分图,是不合法的,输出0. 
(DFS或BFS染色均可) 
3.染色后所有黑色的点进stack1,所有白色的点进stack2,最后模拟输出过程就可以了.

分析:首先,元素要么用一个栈排序,要么用两个栈排序,如果用一个栈排序,那么字典序可以保证最小,为什么要两个栈呢?因为会存在元素f(i),f(j)不能在一个栈里面排序.什么样的元素不能在同一个栈里面排序呢?当f(i) < f(j) f(i) > f(k),且i < j < k时不行,首先k必须要第一个弹出,因为j > i,在f(k)弹出之前f(i)和f(j)都在栈里面,而f(k)弹出之后f(j) > f(i),而f(j)在栈顶,所以不行.根据这个,我们可以把元素分到两个栈里去排序.可以把两个栈看作一个二分图,可以知道一个栈里面的点不能和栈里另一个点相连,如果满足二分图,那么就可以排序.怎么检测是不是二分图呢?把不能在一个栈里面排序的元素连边,给其中一个元素染色,另一个染不同的颜色,如果一个元素相连的颜色和自己相同则不是二分图.我们在染色的时候把最小的颜色染成最小的.这样在排序的时候就可以满足字典序.然后一个一个扫描,模拟排序即可.

 #include<cstdio>
#include<cstring>
#include<iostream>
#include<algorithm>
#define Ll long long
using namespace std;
int a[],v[],f[];
int q1[],top1,q2[],top2,now,l;
bool ok[][];
int n;
void gg(){printf("");exit();}
void dfs(int x,int y)
{
v[x]=y;
for(int i=;i<=n;i++)
if(ok[i][x])
if(v[i]==y)gg();else
if(v[i]==)dfs(i,y^);
}
int main()
{
scanf("%d",&n);
for(int i=;i<=n;i++)scanf("%d",&a[i]);
f[n+]=n+;
for(int i=n;i;i--)f[i]=min(f[i+],a[i]);
for(int i=;i<=n;i++)
for(int j=i+;j<=n;j++)
if(a[i]<a[j]&&a[i]>f[j])ok[i][j]=ok[j][i]=;
for(int i=;i<=n;i++)if(v[i]==)dfs(i,);
l=;now=;a[]=n+;
while(now<=n)
if(v[l]==&&a[q1[top1]]>a[l])printf("a "),q1[++top1]=l++;else
if(now==a[q1[top1]])printf("b "),top1--,now++;else
if(v[l]==&&a[q2[top2]]>a[l])printf("c "),q2[++top2]=l++;else
if(now==a[q2[top2]])printf("d "),top2--,now++;
}

可以用贪心或者二分图染色解决

下面讨论贪心的做法

首先我们来讨论单栈排序的性质.

很显然有如下三个性质:

  1. 栈顶元素一定要大于栈内元素

  2. 给定序列必须从左至右加入栈, 且只能加入一个栈(只有一个栈)

  3. 出栈时, 必须最小的先出栈, 然后大的才能出栈

假如我们不考虑当前已经出栈的元素, 根据3条件, 若当前最小值还未入栈, 那么这个栈一定不能弹, 只能一直加, 直到最小值入栈, 那么在这个过程中, 就有可能出现冲突(与1条件冲突). 用公式表达就是 a3 < a1 < a2a3<a1<a2, 如果出现这种情况就说明无解(单栈排序只有进栈和出栈两种选择, 所以没有多种方案).

那么再把单栈排序推广到双栈排序, 唯一的变数就是可以选择加入哪个栈了, 在单栈排序中会GG的 a3 < a1 < a2a3<a1<a2的情况也不会GG了, 考虑为什么不会GG了, 通过思考可以发现, 因为第二个栈可以容纳一个a2a2 所以可以多一条命, 那么如果出现a4 < a1 < a2 < a3a4<a1<a2<a3 的情况, 即使多一个栈也要GG.

不妨设s1s1为主栈, s2s2为辅栈(优先加入s1s1, 因为这样字典序小),

同过题目描述可知, 优先级关系如下:

  1. 加入s1s1

  2. 从s1s1弹出

  3. 加入s2s2

  4. 从s2s2弹出

我们按照这样的优先级进行每一个操作(即, 如果能加入s1s1就加入s1s1, 否则再考虑能否从s1s1弹出, 否则再考虑加入s2s2...), 弹栈操作无需过多判断, 能弹就弹, 这样一定是最优的, 入栈操作的判断稍微复杂一点, 我们需要判断是否可以把元素加入主栈, 假如当前需要考虑是否加入s1s1的元素为aa, 首先就需要判断它是否满足1条件, 如果满足1条件, 对应到上面的公式, 此时的aa应该是a2a2(为什么不是a1a1?,因为已经入栈的元素的序号一定比未入栈的元素的序号小, a2a2和a1a1刚好满足), 所以我们需要判断它后面是否存在一个a3a3比它大, 并且还有一个a4a4比它小, 如果存在, 那就说明这个元素应该进入辅栈(当然, 如果你发现它连辅栈都进入不了, 那就真的GG了).判断一个元素是否进入辅栈只需要判断这个元素是否满足条件1(比栈顶元素小), 不需要再像主栈那样判断后面是否有.., 为什么? 因为辅栈没有辅栈, 这不是三栈排序!(即使是也不能这么做), 你可能怀疑得这样不会错吗? 当然不会错, 因为我们有最强的1条件守护着辅栈, 在这时就像主栈那样判断.., 反而会错, 因为后面的元素不一定是加入辅栈的, 辅栈只能吃主栈剩下的东西, 假如去判断...,

那么这就是在强制认为这后面的东西是属于辅栈的, 那肯定是不行的!

 #include <cstdio>
#include <cstring>
#include <stack>
using namespace std;
const int N = ;
const char s[] = {'r', 'a', 'b', 'c', 'd'};
int data[N], n;
stack<int> s1;
stack<int> s2;
int ans[N * ], top = ;
int Max[N]; bool check(int pos) {
if (s2.empty()) return true;
int i;
for (i = pos + ; i <= n && (data[i] < data[pos] || data[i] < s2.top()); i++);
for (int j = i + ; j <= n; j++) if (data[j] < data[pos]) return ;
return ;
} int main() {
// freopen("testdata (7).in", "r", stdin);
scanf("%d", &n);
for (int i = ; i <= n; i++)
scanf("%d", &data[i]); for (int i = n; i >= ; i--) {
Max[i] = max(Max[i + ], data[i]);
} bool flag = true;
int pos = , now = ;
while(pos <= n) {
switch() {
case :
if ( (s1.empty() || s1.top() > data[pos]) && check(pos)) {
s1.push(data[pos]);
ans[++top] = ;
pos++;
break;
}
case :
if (s1.empty() == && s1.top() == now) {
ans[++top] = ;
now++;
s1.pop();
break;
}
case :
if ((s2.empty() || s2.top() > data[pos])) {
s2.push(data[pos]);
ans[++top] = ;
pos++;
break;
}
case :
if (s2.empty() == && s2.top() == now) {
ans[++top] = ;
now++;
s2.pop();
break;
}
flag = false;
}
if (flag == false) break;
}
while ((s1.empty() == || s2.empty() == ) && flag) {
if (s1.empty() == && s1.top() == now) {
ans[++top] = ;
now++;
s1.pop();
}
if (s2.empty() == && s2.top() == now) {
ans[++top] = ;
now++;
s2.pop();
}
} if (flag == ) printf("0\n");
else {
for (int i = ; i <= top; i++) {
putchar(s[ans[i]]);
putchar(' ');
}
}
return ;
}

P1155 双栈排序(二分图染色)的更多相关文章

  1. LOJ P1155 双栈排序 二分图染色 图论

    https://www.luogu.org/problem/show?pid=P1155 题解: https://www.byvoid.com/zhs/blog/noip2008-twostack 开 ...

  2. NOIP2008双栈排序[二分图染色|栈|DP]

    题目描述 Tom最近在研究一个有趣的排序问题.如图所示,通过2个栈S1和S2,Tom希望借助以下4种操作实现将输入序列升序排序. 操作a 如果输入序列不为空,将第一个元素压入栈S1 操作b 如果栈S1 ...

  3. [luogu1155 NOIP2008] 双栈排序 (二分图染色)

    传送门 Description Input 第一行是一个整数 n . 第二行有 n 个用空格隔开的正整数,构成一个 1−n 的排列. Output 共一行,如果输入的排列不是"可双栈排序排列 ...

  4. [NOIp2008] 双栈排序 (二分图染色 + 贪心)

    题意 给你一个长为 \(n\) 的序列 \(p\) ,问是否能够通过对于两个栈进行 push, pop(print) 操作使得最后输出序列单调递增(即为 \(1 \cdots n\) ),如果无解输出 ...

  5. 洛谷P1155 双栈排序题解(图论模型转换+二分图染色+栈)

    洛谷P1155 双栈排序题解(图论模型转换+二分图染色+栈) 标签:题解 阅读体验:https://zybuluo.com/Junlier/note/1311990 原题地址:洛谷P1155 双栈排序 ...

  6. [NOIP2008] 提高组 洛谷P1155 双栈排序

    题目描述 Tom最近在研究一个有趣的排序问题.如图所示,通过2个栈S1和S2,Tom希望借助以下4种操作实现将输入序列升序排序. 操作a 如果输入序列不为空,将第一个元素压入栈S1 操作b 如果栈S1 ...

  7. P1155 双栈排序

    题目描述 Tom最近在研究一个有趣的排序问题.如图所示,通过2个栈S1和S2,Tom希望借助以下4种操作实现将输入序列升序排序. 操作aaa 如果输入序列不为空,将第一个元素压入栈S1​ 操作b 如果 ...

  8. 洛谷——P1155 双栈排序

    题目描述 Tom最近在研究一个有趣的排序问题.如图所示,通过2个栈S1和S2,Tom希望借助以下4种操作实现将输入序列升序排序. 操作a 如果输入序列不为空,将第一个元素压入栈S1 操作b 如果栈S1 ...

  9. 洛谷$P1155$ 双栈排序 贪心+二分图匹配

    正解:贪心+二分图匹配 解题报告: 传送门$QwQ$ 跪了,,,我本来以为我$NOIp$做得差不多了,,,然后康了一眼发现没做多少啊其实$QAQ$ 然后来康题趴$QwQ$ 首先考虑如果只有一个栈的情况 ...

随机推荐

  1. Leetcode0006--ZigZag Conversion

    [转载请注明]https://www.cnblogs.com/igoslly/p/9017638.html 来看一下题目: The string "PAYPALISHIRING" ...

  2. [Android]有关外部链接唤醒App需要注意的坑

    移动互联网发展到今天,一个移动app需要和各种各样的外部链接关联,它不再仅仅从手机的桌面启动,更多的将会从其他的应用.浏览器链接.短信.二维码或者微信分享等渠道启动,这里涉及到的是各种各样的营销渠道和 ...

  3. 服务器端 CentOS 下配置 JDK 和 Tonmcat 踩坑合集

    一.配置 JDK 时,在 /etc/profile 文件下配置环境变量,添加   #java environment export JAVA_HOME=/usr/java/jdk- export CL ...

  4. 白盒-CNN纹理深度可视化: 使用MIT Place 场景预训练模型

    MIT发文:深度视觉的量化表示................ Places2 是一个场景图像数据集,包含 1千万张 图片,400多个不同类型的场景环境,可用于以场景和环境为应用内容的视觉认知任务. ...

  5. linux 汇编 - 函数调用

    Linux 汇编-函数调用 *:first-child { margin-top: 0 !important; } body>*:last-child { margin-bottom: 0 !i ...

  6. JAVA趣味逻辑算法

    /**已知4位同学中的一位数学考了100分,当小李询问这4位是谁考了100分时,4个人的回答如下: A说:不是我. B说:是C C说:是D. D说:他胡说. 已知三个人说的是真话,一个人说的是假话.现 ...

  7. js的基础运用

    总结: 1.定义:分为隐式定义和显式定义可以先定义后赋值. 2.+:当两边都是数值则运行加法运算,若一遍是字符型则进行拼接. 3.数值变字符:数值变量.toString()方法. 字符变数值:通过加一 ...

  8. zabbix实现163邮件报警

    Zabbix 邮件报警 电脑登录网易邮箱配置,把自己的授权码看一下,并写入配置文件 server端安装配置邮件服务器 [root@server ~]# yum -y install mailx dos ...

  9. models中,对于(Small)IntegerField类型字段的choices参数在前端的展示

    # models.py class UserInfo(models.Model): gender_choices = ( (1, "男"), (2, "女"), ...

  10. Linux之加密(基于key认证、建立私有云CA)

    对称加密: 一般的加密是用一个密码加密文件,解密用同样的密码,加密解密用一把密钥 非对称加密: 一个密码加密文件,解密却用另外一组密码,意思就是加密解密的密码不一样,其结果就是用这一组密钥中的一个来加 ...