One-Way Conveyors 【Tarjan缩点+树链剖分】

One-Way Conveyors

 Aizu - 1408

---

  这道题上面推荐大家把maxN开的大一点，不然会RE，我倒是被这个坑了一下，还有一些细节，我将在后面徐徐道来。（目前最快还是比较的开心）

  先讲一下题意：有N个点，M条无向边，然后输入这M条无向边，再是一个K，表示我们有K个要求，每个u到v，我们现在要把原来的无向图变成有向图，也就是给每条边定向，还要满足所有K个u到v都要满足从u能抵达v。如果冲突就输出No，不然就是Yes加上每条边的方向，SPJ问题。

  接下去，讲一下我的想法，当然也是题解了，这个问题比较的繁琐，我们将它拆开来考虑。

环内

  环内的所有点都是满足两两之间可以任意的相互到达，所以环内的情况，环内的边是可以随便跑的。

环与环之间

  除了环，就是环与环之间的边了，如果说我们可以要从环A的u点到达环B的v点，那么说明环A到环B的边是由环A到环B的，我们可以确定任意两环的通道是唯一的，不然他们就在一个环内。所以，我们可以给环A到环B定向，因为他们的路径上有可能出现别的缩点的环，所以这些路径也会被定向。

如何解决环与环之间的割边的定向问题

  这里就可以用缩点之后的图来解决，不难发现，缩点之后的图是多棵有根树，我利用超级原点0来将他们都链接起来。所以，就变成了一棵树上的问题，现在我们要是的u到v定向，实际上就是u所属的环向v所属的环定向了，于是，我们可以利用树链剖分来解决这个问题，于是完成了这些割边的定向。

  如果两点本身不相互链接，那么还是要输出No，这里需要单独判断一下。

Yes时候输出的方法

  对于环内的，我们直接跑就可以了，对于环外的，也就是环与环链接的时候，我们可以利用树链剖分上维护的信息来确定，这里自己模拟一下很快就知道了。因为我们树链剖分操作的点也本身是环点，环与环链接的边又是可以唯一确定的，所以这个问题也就引刃而解了。

代码有点长，有些小细节上面要是没有提到的话，可以在评论区call me！

#include <iostream>
#include <cstdio>
#include <cmath>
#include <string>
#include <cstring>
#include <algorithm>
#include <limits>
#include <vector>
#include <stack>
#include <queue>
#include <set>
#include <map>
#include <bitset>
//#include <unordered_map>
//#include <unordered_set>
#define lowbit(x) ( x&(-x) )
#define pi 3.141592653589793
#define e 2.718281828459045
#define INF 0x3f3f3f3f
#define HalF (l + r)>>1
#define lsn rt<<1
#define rsn rt<<1|1
#define Lson lsn, l, mid
#define Rson rsn, mid+1, r
#define QL Lson, ql, qr
#define QR Rson, ql, qr
#define myself rt, l, r
using namespace std;
typedef unsigned long long ull;
typedef unsigned int uit;
typedef long long ll;
const int maxN = 1e5 + 7, maxM = 1e6 + 7;
int N, M, Q, root[maxN];
inline int fid(int x) { return x == root[x] ? x : root[x] = fid(root[x]); }
struct Graph
{
    int head[maxN], cnt;
    struct Eddge
    {
        int nex, to;
        Eddge(int a=-1, int b=0):nex(a), to(b) {}
    } edge[maxM << 1];
    inline void addEddge(int u, int v)
    {
        edge[cnt] = Eddge(head[u], v);
        head[u] = cnt++;
    }
    inline void _add(int u, int v) { addEddge(u, v); addEddge(v, u); }
    inline void init()
    {
        cnt = 0;
        for(int i=0; i<=N; i++) head[i] = -1;
    }
} Old, Now;
bool vis[maxM << 1] = {false}, instack[maxN] = {false};
int dfn[maxN], tot, low[maxN], Stap[maxN], Stop, Belong[maxN], Bcnt;
void Tarjan(int u)
{
    dfn[u] = low[u] = ++tot;
    Stap[++Stop] = u;
    instack[u] = true;
    for(int i=Old.head[u], v; ~i; i=Old.edge[i].nex)
    {
        if(vis[i]) continue;
        vis[i] = vis[i ^ 1] = true;
        v = Old.edge[i].to;
        if(!dfn[v])
        {
            Tarjan(v);
            low[u] = min(low[u], low[v]);
        }
        else if(instack[v]) low[u] = min(low[u], dfn[v]);
    }
    if(low[u] == dfn[u])
    {
        int v; Bcnt++;
        do
        {
            v = Stap[Stop--];
            instack[v] =  false;
            Belong[v] = Bcnt;
        } while(v ^ u);
    }
}
bool read[maxM << 1] = {false};
inline void Build_New_Graph()
{
    for(int u=1; u<=N; u++)
    {
        for(int i=Old.head[u], v; ~i; i=Old.edge[i].nex)
        {
            if(read[i]) continue;
            read[i] = read[i ^ 1] = true;
            v = Old.edge[i].to;
            if(Belong[u] == Belong[v]) continue;
            Now._add(Belong[u], Belong[v]);
        }
    }
}
int siz[maxN], Wson[maxN], deep[maxN], fa[maxN];
void dfs_1(int u, int father)
{
    siz[u] = 1; deep[u] = deep[father] + 1; fa[u] = father;
    int maxx = 0;
    for(int i=Now.head[u], v; ~i; i=Now.edge[i].nex)
    {
        v = Now.edge[i].to;
        if(v == father) continue;
        dfs_1(v, u);
        siz[u] += siz[v];
        if(maxx < siz[v])
        {
            maxx = siz[v];
            Wson[u] = v;
        }
    }
}
int top[maxN], id[maxN], _Index;
void dfs_2(int u, int topy)
{
    top[u] = topy; id[u] = ++_Index;
    if(Wson[u]) dfs_2(Wson[u], topy);
    for(int i=Now.head[u], v; ~i; i=Now.edge[i].nex)
    {
        v = Now.edge[i].to;
        if(v == fa[u] || v == Wson[u]) continue;
        dfs_2(v, v);
    }
}
bool ok;
struct BIT_Tree
{
    int tree[maxN << 1] = {0}, lazy[maxN << 1] = {0};   //1 == down and 2 == up
    inline void pushdown(int rt)
    {
        if(lazy[rt])
        {
            tree[lsn] = lazy[rt]; tree[rsn] = lazy[rt];
            lazy[lsn] = lazy[rt]; lazy[rsn] = lazy[rt];
            lazy[rt]  = 0;
        }
    }
    inline void pushup(int rt) { tree[rt] = tree[lsn] | tree[rsn]; }
    void update(int rt, int l, int r, int ql, int qr, int state)    //state = 0(down) or 1(up)
    {
        if(ql <= l && qr >= r)
        {
            if(tree[rt] & (1 << (1 - state))) ok = false;
            tree[rt] = 1 << state; lazy[rt] = 1 << state;
            return;
        }
        int mid = HalF; pushdown(rt);
        if(qr <= mid) update(QL, state);
        else if(ql > mid) update(QR, state);
        else { update(QL, state); update(QR, state); }
        pushup(rt);
    }
    int query(int rt, int l, int r, int qx)
    {
        if(l == r) return tree[rt];
        int mid = HalF; pushdown(rt);
        if(qx <= mid) return query(Lson, qx);
        else return query(Rson, qx);
    }
    inline void Q_Range(int u, int v)   //from u to v
    {
        int op = 1;
        while(top[u] ^ top[v])
        {
            if(deep[top[u]] < deep[top[v]]) { op ^= 1; swap(u, v); }
            update(1, 1, _Index, id[top[u]], id[u], op);
            u = fa[top[u]];
        }
        if(deep[u] > deep[v])
        {
            update(1, 1, _Index, id[Wson[v]], id[u], op);
        }
        else if(deep[u] < deep[v])
        {
            update(1, 1, _Index, id[Wson[u]], id[v], op ^ 1);
        }
    }
} tree;
bool sured[maxM << 1] = {false}, seen[maxN] = {false}, Bas_point[maxN] = {false};
void OutPut_dfs(int u)
{
    for(int i=Old.head[u], v, bel_u, bel_v, op; ~i; i=Old.edge[i].nex)
    {
        if(sured[i]) continue;
        sured[i] = sured[i ^ 1] = true;
        v = Old.edge[i].to;
        if(Belong[u] ^ Belong[v])
        {
            bel_u = Belong[u]; bel_v = Belong[v];
            if(deep[bel_u] > deep[bel_v])
            {
                op = tree.query(1, 1, _Index, id[bel_u]);
                if(op == 1) printf("%d %d\n", v, u);
                else printf("%d %d\n", u, v);
            }
            else
            {
                op = tree.query(1, 1, _Index, id[bel_v]);
                if(op == 1) printf("%d %d\n", u, v);
                else printf("%d %d\n", v, u);
            }
        }
        else
        {
            printf("%d %d\n", u, v);
            if(Bas_point[v]) continue;
            Bas_point[v] = true;
            OutPut_dfs(v);
        }
    }
}
inline void init()
{
    tot = Stop = Bcnt = _Index = 0;
    Old.init(); Now.init(); ok = true;
    for(int i=1; i<=N; i++) root[i] = i;
}
int main()
{
    scanf("%d%d", &N, &M);
    init();
    for(int i=1, u, v, fu, fv; i<=M; i++)
    {
        scanf("%d%d", &u, &v);
        Old._add(u, v); fu = fid(u); fv = fid(v);
        if(fu ^ fv)
        {
            root[fu] = fv;
        }
    }
    for(int i=1; i<=N; i++) if(!dfn[i]) { Tarjan(i); Now._add(0, Belong[i]); }
    Build_New_Graph();
    dfs_1(0, 0);
    dfs_2(0, 0);
    scanf("%d", &Q);
    for(int i=1, x, y; i<=Q; i++)
    {
        scanf("%d%d", &x, &y);
        if(fid(x) ^ fid(y)) ok = false;
        if(!ok) continue;
        if(Belong[x] == Belong[y]) continue;
        tree.Q_Range(Belong[x], Belong[y]);
    }
    if(!ok) { printf("No\n"); return 0; }
    printf("Yes\n");
    for(int i=1; i<=N; i++)
    {
        if(seen[Belong[i]]) continue;
        seen[Belong[i]] = true;
        Bas_point[i] = true;
        OutPut_dfs(i);
    }
    return 0;
}