【提高组NOIP2008】双栈排序 (twostack.pas/c/cpp)

【题目描述】
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)排序的操作序列：

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

【输入】
输入文件twostack.in的第一行是一个整数n。
第二行有n个用空格隔开的正整数，构成一个1~n的排列。

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

样例输入	样例输出
4 1 3 2 4	a b a a b b a b
4 2 3 4 1	0
3 2 3 1	a c a b b d

【数据范围限制】

30%的数据满足： n<=10
50%的数据满足： n<=50
100%的数据满足： n<=1000

【题解】
这题有一点难度。

比赛时，我一看到这道题目，就想到了一种很WS的东西：暴力模拟！
然后就只拿了 27.3 分

其实这道题目要用到一种神奇的东西—— 二分图染色

大体的思路就是先处理好每一个数要进入的栈，再暴力模拟一遍就好了。

设 s[i] 表示输入的第i个数
设 f[i] 表示从第i个数起到第n个数的 最小值

我们会发现，总有一些数字是不能在同一个栈里的（由始至终都不可以哦~）。
如果 i < j < k，且 s[k] < s[i] < s[j]，那么s[i]和s[j]就一定不能出现在同一个栈里面。
论证方法就不写了。

但每次都枚举i,j,k，明显会时超。
所以我们就可以用到 f 数组了！
认真地再看一看这条语句：

如果 i < j < k，且 s[k] < s[i] < s[j]，那么s[i]和s[j]就一定不能出现在同一个栈里面。

也就是说，如果 $j+1到n之间的最小值<s[i]<s[j]$，i,j就不能同时出现在同一个栈中。

也就是$当 f[j+1] < s[i]<s[j]( i< j ) $时，i,j就不能同时出现在同一个栈中。

这样连边就能只用一个双重循环来做了。

再设 c[i] 表示第i个数要进去的栈（栈S1用1表示，栈S2用2表示，初始化为0）

我们就开始枚举1~N的所有数字i，只要c[i]为0，就把c[i]赋值（染色）为1（这样一定最优），然后开始枚举与 i 有连边的数 j，把c[j]赋值为和c[i]的值相反的数（即当c[i]=1,c[j]=2;c[i]=2,c[j]=1），并从j那里出发，继续染与j有连边的数……如果碰到c[j]已经被赋（染）过值了，并且c[i]=c[j]，那肯定是无法排序的了，输出0并结束程序就好了。

染色用DFS来做。

最后根据c[i]的值，暴力模拟一遍就可以了。

下面来说一下一些细节：

1. 连边时一定要连 双向边，譬如要在 x,y 之间连一条边，都要连两条单向的边 x,y 和 y,x，存边最好用 邻接表 ；
2. 连完边后，要把1到n的数都走一次 递归 ，并把与这个数相连的所有数染完色后，从这些数继续递归染色；
3. 每个数都只递归一次，因为重复地递归同一个数，是浪费时间而没有半点用处的；
4. 有一种特殊的情况，就是当栈S2的栈顶可以弹栈，而S1又可以从输入的栈加入一个数时，一定要先读入，后弹栈——因为如果先弹栈，后输入，那输出队列里存的就是 d a 了，这样明显不如 a d 那么优。

建议还不太明白的同学来模拟一些这个数据：

14
2 3 1 4 14 6 12 5 11 13 7 8 9 10

还有不懂的同学可以来找我。

下面附上标程，可以参考一下。

#include<cstdio>
#include<cstdlib>