JZOJ4876. 【NOIP2016提高A组集训第10场11.8】时空传送 拓扑序判断最长路是否合法

题目大意

给定一个$n$个点，$m$条边的有向图无环图，现在要求你删除一个点以及与它有关的边，使得图中的最长路最短。

$n\leq4*10^5$
$m\leq10^6$

解题思路

首先考虑当必须选某一条边时的最长路怎么算。假设这条边连接$u,v$。那么我们可以先预处理出$f_i$数组表示到$i$点的最长路是多长和$g_i$表示从$i$点往后走的最长路是多长，这两个都可以通过拓扑序$O(n)$得到。那么一条边的最长路就是$f_u+1+g_v$。但是现在要求删掉一个点。那么考虑这个点$p$会对包含哪些边的最长路没有影响。设一个点的拓扑序为$tp_i$，那么对于一条连接$u,v$的边，当$tp_u<p$并且$tp_v>p$时这个点对这条边肯定没有影响。但是有肯能不满足这个条件时。点$p$仍对这条边没有影响。但是仔细考虑一下就会发现，在这条最长路中肯定会有一条边符合要救，也就是说这条边的贡献肯定会被考虑上。另外为了防止类似一条链时断路的情况。比较最大值时也要把合法的$f_i$和$g_i$考虑上，判断是否合法的方法与上面类似。

程序

//YxuanwKeith
#include <cstring>
#include <cstdio>
#include <algorithm>
#include <set>

using namespace std;

const int MAXN = 75005, MAXM = 1e5 + 5;

multiset<int> d;
int n, m, deg[MAXN], ord[MAXN], f[MAXN], g[MAXN];
int tot, In[MAXN], Out[MAXN], Next[MAXM * 2], Go[MAXM * 2]; 

void link_in(int u, int v) {
    Next[++ tot] = In[u], In[u] = tot, Go[tot] = v;
}

void link_out(int u, int v) {
    Next[++ tot] = Out[u], Out[u] = tot, Go[tot] = v;
}

void add(int now) {
    for (int p = Out[now]; p; p = Next[p]) 
        d.insert(f[now] + g[Go[p]] + 1);
    d.insert(f[now]);
}

void del(int now) {
    for (int p = In[now]; p; p = Next[p]) 
        d.erase(d.find(g[now] + f[Go[p]] + 1));
    d.erase(d.find(g[now]));
}

void prepare() {
    static int d[MAXN];
    int l = 0, r = 0, cnt = 0;
    for (int i = 1; i <= n; i ++) 
        if (deg[i] == 0) d[++ r] = i, ord[++ cnt] = i;
    while (l <= r) {
        int now = d[++ l];
        for (int p = Out[now]; p; p = Next[p]) {
            int v = Go[p];
            deg[v] --;
            if (deg[v] == 0) {
                ord[++ cnt] = v;
                d[++ r] = v;
            }
        }
    }
    for (int i = 1; i <= n; i ++) {
        int now = d[i];
        for (int p = Out[now]; p; p = Next[p]) 
            f[Go[p]] = max(f[now] + 1, f[Go[p]]);
    }
    for (int i = n; i; i --) {
        int now = d[i];
        for (int p = In[now]; p; p = Next[p]) 
            g[Go[p]] = max(g[now] + 1, g[Go[p]]);
    }
}

int main() {
    scanf("%d%d", &n, &m);
    for (int i = 1; i <= m; i ++) {
        int u, v;
        scanf("%d%d", &u, &v);
        link_out(u, v);
        link_in(v, u);
        deg[v] ++;
    }
    prepare();
    for (int i = 1; i <= n; i ++) 
        d.insert(g[i]);
    int ans = n + 1, id;
    for (int i = 1; i <= n; i ++) {
        int now = ord[i];
        del(now);
        int s = *--d.end();
        if (s < ans || s == ans && now < id) ans = s, id = now;
        add(now);
    }
    printf("%d %d\n", id, ans);
}