【tarjan缩点+分层图】草鉴定Grass Cownoisseur

P3119 [USACO15JAN] 草鉴定Grass Cownoisseur

tarjan缩点+建分层图

缩点：单向边+可重复走一个牧场

缩完点后， 对于不在原图中的一条边(u, v),  if(sd[u] != sd[v]) 连边， 建立新图编号为1~tot//sd[] 该点属于哪个强连通分量， tot是强连通的数量

分层图：可逆向走一次，分层图顶点编号为 tot+1 ~ 2 * n；对于缩点后新图中的一条边(u, v), 将该边终点与分层图的起点连边， 达到逆向走一次的目的add(v, u+tot)

分层图自身也要连边， 对于缩点后的新图(u, v) ， 分层图连相应的边add(u+tot, v+tot)

记得分层图上的点权值与缩点后新图权值相同！！

因为最后会走回1号点，会加上1号点所在强连通的权值， 所以开始初始dis只可以为0

自信代码姹紫嫣红竟因忘删注释？！数组开的正好re两个点， 开大十倍就好了QAQ

#include<queue>
#include<cstdio>
#include<iostream>
using namespace std;
const int sz = 1000010;
int n, m, tim = 0, top = 0, num = 0, num2 = 0, tot = 0;
int h[sz], dfn[sz], low[sz], sta[sz], x[sz], y[sz], sd[sz];
int h2[sz], dis[sz ], node[sz];
bool vis[sz], book[sz ];
struct Egde {
    int nxt, to, dis;
}e[sz], e2[sz << 1];
int read() {
    char ch = getchar(); int x = 0, f = 1;
    while(ch < '0' || ch > '9') {if(ch == '-') f = -1; ch = getchar();}
    while(ch >= '0' && ch <= '9') {x = x * 10 + ch - '0'; ch = getchar();}
    return x * f;
}
void add(int from, int to) {
    e[++num].nxt = h[from];
    e[num].to = to;
    h[from] = num;
}
void add2(int from, int to) {
    e2[++num2].nxt = h2[from];
    e2[num2].to = to;
     h2[from] = num2;
}
void tarjan(int u) {
    low[u] = dfn[u] = ++tim;
    sta[++top] = u;
    vis[u] = 1;
    for(int i = h[u]; i; i = e[i].nxt) {
        int v = e[i].to;
        if(!dfn[v]) {
            tarjan(v);
            low[u] = min(low[v], low[u]);
        }
        else if(vis[v]) low[u] = min(low[u], dfn[v]);
    }
    if(dfn[u] == low[u]) {
            tot++;
        while(u != sta[top+1]) {
            sd[sta[top]] = tot;
            node[tot]++;
            vis[sta[top]] = 0;
            top--;
        }
    }
}
void spfa(int s) {
    queue<int>q;
    for(int i = 1; i <= (n<<1); i++) {
        dis[i] = 0;
        book[i] = 0;
    }
    book[s] = 1;
    q.push(s);
    while(!q.empty()) {
        int u = q.front();
        q.pop();
        book[u] = 0;
        for(int i = h2[u]; i; i = e2[i].nxt) {
            int v = e2[i].to;
            if(dis[v] < dis[u] + node[v]) {
                dis[v] = dis[u] + node[v];
                if(!book[v]) {
                    book[v] = 1;
                    q.push(v);
                }
            }
        }
    }
}
int main() {
    scanf("%d%d", &n, &m);
    for(int i = 1; i <= m; i++) {
        x[i] = read(), y[i] = read();
        add(x[i], y[i]);
    }    
    for(int i = 1; i <= n; i++)
        if(!dfn[i]) tarjan(i);
    for(int i = 1; i <= m; i++) {
        int u = sd[x[i]], v = sd[y[i]];
        if(u != v) {
            add2(u, v);
            add2(v, u+tot);
            add2(u+tot, v+tot);
        }    
    }
    for(int i=1;i<=tot;i++) node[i+tot]=node[i];
    spfa(sd[1]);
    printf("%d", dis[sd[1]+tot]);
    return 0;
}

和aq重学了个缩点QAQ