参见我上一篇博客。
没有缩点的tarjin
传送门!!
下面是有缩点的
#include<cstdio>
int vet[2000],head[2000],next[2000];
int edgenum,index,top,cnt; //不做说明,这4个数初值都是0;
int dfn[2000],low[2000],instack[2000],stack[2000],belong[2000]; //dfn low belong 记得每次清零。(如果要做多次的话)
int n,m;
inline int add(int x,int y) {
vet[++edgenum]=y;
next[edgenum]=head[x];
head[x]=edgenum;
}
int tarjin(int u) {
low[u]=dfn[u]=++index;
instack[u]=true;
stack[++top]=u;
for(int e=head[u]; e; e=next[e]) {
int v=vet[e];
if(!dfn[v]) {
tarjin(v);
if(low[u]>low[v])
low[u]=low[v];
} else if(instack[v]&& low[u]>dfn[v])
low[u]=dfn[v];
}
if(dfn[u]==low[u]) {
cnt++;
int j=0;
while(j!=u) {
j=stack[top--];
instack[j]=false;
belong[j]=cnt;
}
}
}
void work() {
for (int i=1; i<=n; i++)
if (!dfn[i])
tarjin(i);
}
void zip(){
for(int i=1;i<=n;i++)
for(int e=head[i];e;e=next[e])
if(belong[i]!=belong[vet[e]])
add(belong[i]+n,belong[vet[e]]+n);
}
int main() {
scanf("%d%d",&n,&m);
int x,y;
for(int i=1; i<=m; i++) {
scanf("%d%d",&x,&y);
add(x,y);
}
work();
zip();
}