Luogu 3119 [USACO15JAN]草鉴定Grass Cownoisseur

思路很乱，写个博客理一理。

缩点 + dp。

首先发现把一个环上的边反向是意义不大的，这样子不但不好算，而且相当于浪费了一次反向的机会。反正一个强连通分量里的点绕一遍都可以走到，所以我们缩点之后把一个强连通分量放在一起处理。

设$st$表示缩点之后$1$所在的点，设$f_{x}$表示从$st$走到$x$的最长链，$g_{x}$表示从$x$走到$st$的最长链，因为把一个$DAG$上的边反向一下并不会走重复的点，那么我们最后枚举一下边$(x, y)$，把它反向，这样子$f_{x} + g_{y} - siz_{st}$就可以成为备选答案，更新$ans$即可。

注意到有可能整个图强连通，所以$ans$应初始化为$siz_{st}$。

考虑一下$f$和$g$怎么求，一种想法是$st$开始的最长路，我们可以在反图和正图上分别跑一遍$spfa$，这样子可以通过，但是我并不清楚在$DAG$上$spfa$的表现是不是稳定的，另一种想法就是$dp$，直接记搜搞一搞，但是直接记搜是错误的，因为有一些点是不可能走到的，所以在$dp$之前要先$dfs$一遍标记出所有的合法点，然后再进行记搜即可。

时间复杂度$O(n)$。

放上写得很丑很长还可能有锅的代码。

你谷的数据是真的水。

Code：

#include <cstdio>
#include <cstring>
#include <vector>
using namespace std;

const int N = 1e5 + ;

int n, m, tot = , head[N], dfsc = , dfn[N], low[N];
int scc = , f[N], g[N], inx[N], iny[N], top = , sta[N], bel[N], siz[N];
bool vis[N], ok[N];
vector <int> G1[N], G2[N];

struct Edge {
    int to, nxt;
} e[N];

inline void add(int from, int to) {
    e[++tot].to = to;
    e[tot].nxt = head[from];
    head[from] = tot;
}

inline void read(int &X) {
    X = ; char ch = ; int op = ;
    for(; ch > '' || ch < ''; ch = getchar())
        if(ch == '-') op = -;
    for(; ch >= '' && ch <= ''; ch = getchar())
        X = (X << ) + (X << ) + ch - ;
    X *= op;
}

inline int min(int x, int y) {
    return x > y ? y : x;
}

void tarjan(int x) {
    low[x] = dfn[x] = ++dfsc;
    vis[x] = , sta[++top] = x;
    for(int i = head[x]; i; i = e[i].nxt) {
        int y = e[i].to;
        if(!dfn[y]) {
            tarjan(y);
            low[x] = min(low[x], low[y]);
        } else if(vis[y])
            low[x] = min(low[x], dfn[y]);
    } 

    if(low[x] == dfn[x]) {
        ++scc;
        for(; sta[top + ] != x; --top) {
            vis[sta[top]] = ;
            siz[scc]++;
            bel[sta[top]] = scc;
        }
    }
}

inline void chkMax(int &x, int y) {
    if(y > x) x = y;
}

void dfs1(int x) {
    ok[x] = , vis[x] = ;
    for(unsigned int i = ; i < G1[x].size(); i++) {
        int y = G1[x][i];
        dfs1(y);
    }
}

int dp1(int x) {
    if(vis[x]) return f[x];
    vis[x] = ;
    int res = ;
    for(unsigned int i = ; i < G2[x].size(); i++) {
        int y = G2[x][i];
        if(ok[y]) chkMax(res, dp1(y));
    }
    f[x] = res + siz[x];
    return f[x];
}

void dfs2(int x) {
    ok[x] = , vis[x] = ;
    for(unsigned int i = ; i < G2[x].size(); i++) {
        int y = G2[x][i];
        dfs2(y);
    }
}

int dp2(int x) {
    if(vis[x]) return g[x];
    vis[x] = ;
    int res = ;
    for(unsigned int i = ; i < G1[x].size(); i++) {
        int y = G1[x][i];
        if(ok[y]) chkMax(res, dp2(y));
    }
    g[x] = res + siz[x];
    return g[x];
}

int main() {
//    freopen("testdata.in", "r", stdin);

    read(n), read(m);
    for(int i = ; i <= m; i++) {
        read(inx[i]), read(iny[i]);
        add(inx[i], iny[i]);
    } 

    for(int i = ; i <= n; i++)
        if(!dfn[i]) tarjan(i);

    for(int i = ; i <= m; i++) {
        if(bel[inx[i]] == bel[iny[i]]) continue;
        G1[bel[inx[i]]].push_back(bel[iny[i]]);
        G2[bel[iny[i]]].push_back(bel[inx[i]]);
    }

    memset(vis, , sizeof(vis));
    memset(ok, , sizeof(ok));
    dfs1(bel[]);

    memset(vis, , sizeof(vis));
    for(int i = ; i <= scc; i++) {
        if(ok[i]) dp1(i);
    }

    memset(vis, , sizeof(vis));
    memset(ok, , sizeof(ok));
    dfs2(bel[]);

    memset(vis, , sizeof(vis));
    for(int i = ; i <= scc; i++) {
        if(ok[i]) dp2(i);
    }

/*    for(int i = 1; i <= scc; i++)
        printf("%d ", g[i]);
    printf("\n");
    for(int i = 1; i <= scc; i++)
        printf("%d ", f[i]);
    printf("\n");   */

    int ans = siz[bel[]];
    for(int i = ; i <= m; i++) {
        int u = bel[iny[i]], v = bel[inx[i]];
        if(u == v) continue;
        if(f[u] && g[v]) chkMax(ans, f[u] + g[v] - siz[bel[]]);
    }

    printf("%d\n", ans);
    return ;
}