洛谷 1155 双栈排序 NOIP2008 图论 二分图染色

传送门

瞎扯(看题解的可以忽略这段)

• 估计不看一眼题解我是永远做不出来这道题目的，，
• 好像二分图染色的题目没有什么正常的，，
• 这道题藏的很深，，，
• 心机，，，=。=
• 看了眼题解之后思路还算清晰，
• 通过单栈拓展到双栈，进而归纳冲突发生条件，通过二分图避免冲突（其实这段还好理解）
• 然后遇到了一点点小麻烦，在染色时考虑对于多个联通块的染色方式， 在这里我是根据其在输入队列中的顺序来进行染色，即把出现早的点染成1，出现晚的点染成2，这样可以保证字典序
• 然后开始模拟，按顺序尝试b，a，d，c，保证字典序
• 最后按b，d的顺序弹出栈直到两栈全空

正经题解：

1. 题目中所给的双栈较为复杂，将其转化为单栈进行考虑

2. 对于一个栈中的i和j，当存在k使得$i < j < k，a[k] < a[i] < a[j]$时，i和j存在冲突，即不能用一个栈解决

   证明：
   充分性：即如果满足上述条件，那么q[i]和q[j]一定不能压入同一个栈。
   反证法:假设这两个数压入了同一个栈,那么压入q[k]，因为q[k]比q[i]和q[j]都小，所以很显然,当q[k]没有被弹出的时候，另两个数也都不能被弹出(否则输出序列的数字顺序就不是1,2,3,…,n了)。而之后,无论其它的数字在什么时候被弹出,q[j]总是会在q[i]之前弹出，而q[j]>q[i],这显然是不正确的.
   必要性：如果两个数不可以压入同一个栈，那么它们一定满足上述条件。
   证明逆否命题：也就是”如果不满足上述条件，那么这两个数一定可以压入同一个栈。”不满足上述条件有两种情况：
   情况1：对于任意$i<j<k且q[i]<q[j]，q[k]>q[i]$；（即对任意三个数，最小的总是在最前面）
   情况2：对于任意$i<j，q[i]>q[j]$。
   第一种情况：在q[k]被压入栈的时候，q[i]已经被弹出栈。那么,q[k]不会对q[j]产生任何影响(这里可能有点乱，因为看起来,$q[j]<q[k]$的时候是会有影响的，但实际上，这还需要另一个数r，满足$j<k<r且 q[r]<q[j]<q[k]$,也就是证明充分性的时候所说的情况。而事实上我们现在并不考虑这个r，所以说q[k]对q[j]没有影响)。
   第二种情况：可以发现这其实就是一个降序序列，所以所有数字都可以压入同一个栈。这样，原命题的逆否命题得证，所以原命题得证。

   以上证明来自sicily，传送门

3.因为以上判断方法的复杂度为$O(n^3)$，使用dp的思想，以$f[x]$表示从n到x的最小值，那么对于每一对(i,j)，可以通过检查$a[i]<a[j]，f[j+1] < a[i]$来判断i和j是否冲突，如果冲突即在i和j之间连边，表示i和j不能用同一个栈来解决
4. 按照输入的顺序，从每一个未被染色的点开始进行二分图染色，如果出现冲突直接输出0，exit
5. 按照b > a > d > c的顺序进行模拟，考虑到每次可能需要先弹栈再进栈，所以将弹出的优先级高于入栈，
6. 使用一个now变量，记录此时应该被弹出的数字，能简化思考及代码
7. 最后按顺序清空两个栈，即为结果

代码：

#include <cstdio>
#include <cstring>
#include <algorithm>

const int maxn = 1000 + 50;
const int maxm = 2000000 + 500;
const int inf = 0x7fffffff;
int n;
int a[maxn];
int f[maxn];
int last[maxn], pre[maxm], other[maxm];
int col[maxn];
int tot = 0;
bool flag = 0;
int sta1[maxn], sta2[maxn];
int statot2 = 0, statot1 = 0;
int anow = 0;

void add(int x, int y) {
    tot++;
    pre[tot] = last[x];
    last[x] = tot;
    other[tot] = y;
}

void pnt(int x) {
    int cur = 0;
    if (col[x] == 1) cur = 2;
    else cur = 1;
    if (flag) return;
    for (int p = last[x]; p; p = pre[p]) {
        int q = other[p];
        if (col[q] == 0) {
            col[q] = cur;
            pnt(q);
        } else if (col[q] != cur) {
            flag = 1;
            return;
        }
    }
}

int main () {
    scanf("%d", &n);
    for (int i = 1; i <= n; i++) scanf("%d", &a[i]);
    f[n+1] = inf;
    for (int i = n; i >= 1; i--) f[i] = std :: min(f[i+1], a[i]);
    for (int i = 1; i <= n; i++) 
        for (int j = i + 1; j <= n; j++) {
            if (a[i] < a[j] && f[j+1] < a[i]) {
                add(a[i], a[j]);
                add(a[j], a[i]);
            }
        }
    for (int i = 1; i <= n; i++) {
        if (flag) break;
        if (col[a[i]] == 0) {
            col[a[i]] = 1;
            pnt(a[i]);
        }
    }
    if (flag) {
        printf("0");
        exit(0);
    }
    anow = 1;
    for (int i = 1; i <= n; i++) {
        while (anow == sta1[statot1]) {
            anow++;
            statot1--;
            printf("b ");
        }
        if (col[a[i]] == 1) {
            printf("a ");
            statot1++;
            sta1[statot1] = a[i];
            continue;
        }
        while (anow == sta2[statot2]) {
            anow++;
            statot2--;
            printf("d ");
        }
        if (col[a[i]] == 2) {
            printf("c ");
            statot2++;
            sta2[statot2] = a[i];
        }
    }    
    while (anow <= n) {
        if (sta1[statot1] == anow) {
            printf("b ");
            anow++;
            statot1--;
        }
        if (sta2[statot2] == anow) {
            printf("d ");
            anow++;
            statot2--;
        }
    }
    return 0;
}