ccf+tarjan+求有向图的强连通分量

点击打开链接

///tarjan的强连通分量
#include<iostream>
#include<cstdio>
#include<cstring>
#include<stack>
using namespace std;
const int maxn=1e4+5,maxm=1e5+5;
struct Side
{
    int to,next;
}side[maxm];
int node[maxn],low[maxn],dfn[maxn],belong[maxn];
int top,n,m,t,cur,ans=0,sum=0;
stack<int>s;
void add_side(int u,int v)
{
    side[top]=(Side){v,node[u]};
    node[u]=top++;
}
void dfs(int u)
{
    if(dfn[u]) return;
    dfn[u]=low[u]=t++;
    s.push(u);
    for(int i=node[u];i!=-1;i=side[i].next){
        int v=side[i].to;
        dfs(v);
        if(dfn[v]!=-1) low[u]=min(low[u],low[v]);
    }
    int temp;
    if(low[u]==dfn[u]){
    	sum=0;//初始化 
        do{
            temp=s.top();
            s.pop();
            belong[temp]=cur;
            dfn[temp]=-1;
            sum++;///一个强连通分量中点的个数 
        }while(temp!=u);
        cur++;///强连通分量的个数。 
        ans+=(sum*(sum-1))/2;
    }
}
int main()
{
    scanf("%d%d",&n,&m);
    top=0;
    memset(node,-1,sizeof(node));
    memset(dfn,0,sizeof(dfn));
    while(m--){
        int a,b;
        cin>>a>>b;
        add_side(a,b);
    }
    t=1,cur=1;///时间戳,强连通分量 
    for(int i=1;i<=n;i++){///从1开始 
    	if(!dfn[i]){
    		dfs(i);
    	}
    } 
    printf("%d\n",ans);
    return 0;
}