tarjan强连通分量

---

tarjan强连通分量

连通：无向图中，任意点i可以到任意点j

强连通：有向图中，任意点i可以到任意点j

vector<int’>c[ ]作为邻接表
dfn[ ]记录dfs顺序
low[ ]记录该点能去到的最早的dfs数

HDU1827

思路

tarjan找到闭环，缩点，只需联系闭环入度为0的

#include<iostream>
#include<cstring>
#include<stack>
#include<algorithm>
#include<vector>
#include<cstdio>
using namespace std;
int n,m,t,cnt,num,ans,vis[1005],dfn[1005],low[1005],cost[1005],belong[1005],lcost[1005],in[1005];
vector<int>c[1005];
stack<int>s;
void init(){
     for(int i=0;i<=1004;i++)c[i].clear();
        cnt=0;t=1;num=0;ans=0;
        memset(vis,0,sizeof(vis));
        memset(dfn,0,sizeof(dfn));
        memset(belong,0,sizeof(belong));
        memset(in,0,sizeof(in));
        for(int i=1;i<=n;i++)lcost[i]=1000000;
}
void tarjan(int x){
    s.push(x);
    vis[x]=1;
    dfn[x]=low[x]=t++;
    for(int i=0;i<c[x].size();i++){
        int temp=c[x][i];
        if(dfn[temp]==0){
            tarjan(temp);
            low[x]=min(low[x],low[temp]);
        }
        else if(vis[temp])low[x]=min(low[x],low[temp]);
    }
    if(low[x]==dfn[x]){
        int j;
        cnt++;
        while(1){j=s.top();
            s.pop();
            belong[j]=cnt;//标记，在同一个闭环的元素
            lcost[cnt]=min(lcost[cnt],cost[j]);//记录闭环中最小花费
            vis[j]=0;
            if(j==x)break;
        }
    }
}
void getdegree(){
    for(int i=1;i<=n;i++){//遍历有向的两点是不是在相同闭环里面，不在则说明闭环有入度
            for(int j=0;j<c[i].size();j++){
                int v=c[i][j];
                if(belong[v]!=belong[i])
                    in[belong[v]]++;//记录闭环入度数量
            }
        }
        
        for(int i=1;i<=cnt;i++){
            if(in[i]==0){
                num++;
                ans+=lcost[i];
            }
        }
}
int main(){
    while(scanf("%d%d",&n,&m)!=EOF){
        init();

        for(int i=1;i<=n;i++)cin>>cost[i];
        for(int i=1;i<=m;i++){
            int a,b;
            cin>>a>>b;
            c[a].push_back(b);
        }
        
        for(int i=1;i<=n;i++)
            if(dfn[i]==0)
                tarjan(i);

        getdegree();

        printf("%d %d\n",num,ans);//卡时间了，cout会TLE
    }
    return 0;
}

---

HDU3836

思路

求要加多少条边能全部强联通。
强连通说明每个缩点都有入度和出度
吧上题再记录一下出度即可；

#include<iostream>
#include<cstring>
#include<stack>
#include<algorithm>
#include<vector>
#include<cstdio>
using namespace std;
int n,m,t,cnt,num1,num2,vis[20005],dfn[20005],low[20005],cost[20005],belong[20005],lcost[20005],in[20005],out[20005];
vector<int>c[20005];
stack<int>s;
void init(){
     for(int i=0;i<=20004;i++)c[i].clear();
        cnt=0;t=1;num1=0;num2=0;
        memset(vis,0,sizeof(vis));
        memset(dfn,0,sizeof(dfn));
        memset(belong,0,sizeof(belong));
        memset(in,0,sizeof(in));
        memset(out,0,sizeof(out));
        for(int i=1;i<=n;i++)lcost[i]=1000000;
}
void tarjan(int x){
    s.push(x);
    vis[x]=1;
    dfn[x]=low[x]=t++;
    for(int i=0;i<c[x].size();i++){
        int temp=c[x][i];
        if(dfn[temp]==0){
            tarjan(temp);
            low[x]=min(low[x],low[temp]);
        }
        else if(vis[temp])low[x]=min(low[x],low[temp]);
    }
    if(low[x]==dfn[x]){
        int j;
        cnt++;
        while(1){j=s.top();
            s.pop();
            belong[j]=cnt;//标记，在同一个闭环的元素
            vis[j]=0;
            if(j==x)break;
        }
    }
}
void getdegree(){
    for(int i=1;i<=n;i++){//遍历有向的两点是不是在相同闭环里面，不在则说明闭环有入度
            for(int j=0;j<c[i].size();j++){
                int v=c[i][j];
                if(belong[v]!=belong[i]){
                    in[belong[v]]++;//记录闭环入度数量
                    out[belong[i]]++;
                }
            }
        }
        
        for(int i=1;i<=cnt;i++){
            if(in[i]==0)num1++;
            if(out[i]==0)num2++;
        }
}
int main(){
    while(scanf("%d%d",&n,&m)!=EOF){
        init();

        for(int i=1;i<=m;i++){
            int a,b;
            cin>>a>>b;
            c[a].push_back(b);
        }
        
        for(int i=1;i<=n;i++)
            if(dfn[i]==0)
                tarjan(i);

        getdegree();

        if(cnt==1)printf("0\n");
        else printf("%d\n",max(num2,num1));//卡时间了，cout会TLE
    }
    return 0;
}