UVA 12167 Proving Equivalences(强连通)

注意最后需要枚举所有不是强连通内部边的边。

#include<iostream>
#include<cstdio>
#include<cstring>
#include<algorithm>
#include<vector>
#include<stack>
using namespace std;
const int maxn=20000+5;
vector<int> G[maxn];
int pre[maxn],lowlink[maxn],sccno[maxn],dfs_clock,scc_cnt;
stack<int> s;
int n,m;
struct edge{int u,v;}e[50000+5];
void dfs(int u)
{
    pre[u]=lowlink[u]=++dfs_clock;
    s.push(u);
    for(int i=0;i<G[u].size();i++){
        int v=G[u][i];
        if(!pre[v]){
            dfs(v);
            lowlink[u]=min(lowlink[u],lowlink[v]);
        }
        else if(!sccno[v]){
            lowlink[u]=min(lowlink[u],pre[v]);
        }
    }
    if(lowlink[u]==pre[u]){
        scc_cnt++;
        for(;;){
            int x=s.top();s.pop();
            sccno[x]=scc_cnt;
            if(x==u) break;
        }
    }
}
void find_scc(){
    dfs_clock=0;scc_cnt=0;
    memset(pre,0,sizeof(pre));
    memset(sccno,0,sizeof(sccno));
    for(int i=0;i<n;i++)
        if(!pre[i]) dfs(i);
}
int in[maxn],out[maxn];
int main()
{
    int t;
    scanf("%d",&t);
    while(t--)
    {
        scanf("%d%d",&n,&m);
        for(int i=0;i<n;i++) G[i].clear();
        for(int i=0;i<m;i++)
        {
            scanf("%d%d",&e[i].u,&e[i].v);
            e[i].u--;e[i].v--;
            G[e[i].u].push_back(e[i].v);
        }
        find_scc();
        memset(in,0,sizeof(in));
        memset(out,0,sizeof(out));
        for(int i=0;i<m;i++){
            int pu=sccno[e[i].u];
            int pv=sccno[e[i].v];
            if(pu==pv) continue;
            out[pu]=1;
            in[pv]=1;
        }
        int a=0,b=0;
        for(int i=1;i<=scc_cnt;i++){
            if(in[i]==0) a++;
            if(out[i]==0) b++;
        }
        int ans=max(a,b);
        if(scc_cnt==1) ans=0;
        printf("%d\n",ans);
    }
    return 0;
}