点击打开链接
///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;
}