HDU3836 Equivalent Sets :Tarjan缩点
Description
一个N个点与M条边的有向图,求出需要加多少使其变为一个强联通分量.
Input
若干组测试数据.
每组数据第一行为N,M
后面M行为两个正整数A,B,代表A到B有一条有向边.
Output
对应每组数据输出一个整数,及要加的最小边数.
Sample Input
4 0
3 2
1 2
1 3
Sample Output
4
2
HINT
N <= 20000,M <= 50000
题解
可以先进行缩点.然后就有一个无环有向图。
我们可以尝试发现只有然所有连通块的入度与出度不为0就可以变成一个连通块。
那么就让加的每一条边从一个出度为0的连通块连接一个入度为0的连通块即可。
所以答案就算max(出度为0的连通块的个数,入度为0的连通块的个数)
#include <cstdio>
#include <iostream>
#include <cmath>
#include <stack>
#include <algorithm>
#include <cstring>
#include <climits>
#define MAXN 40000+10
#define MAXM 70000+10
using namespace std;
int head[MAXN],num,n,m;
int dfn[MAXN],low[MAXN],vis[MAXN],dfnum,col[MAXN],co;
int ru[MAXN],rans,cans,cu[MAXN],t;
stack<int> st;
struct Edge{
int from,to,next;
}edge[MAXN<<1];
void add(int from,int to)
{
edge[++num].next=head[from];
edge[num].from=from;
edge[num].to=to;
head[from]=num;
}
void tarjian(int x)
{
dfn[x]=low[x]=++dfnum;
vis[x]=1;st.push(x);
for(int i=head[x];i;i=edge[i].next)
{
if(!dfn[edge[i].to])
{
tarjian(edge[i].to);
low[x]=min(low[x],low[edge[i].to]);
}else if(vis[edge[i].to]) low[x]=min(low[x],dfn[edge[i].to]);
}
if(dfn[x]==low[x])
{
col[x]=++co;
vis[x]=0;
while(st.top()!=x){col[st.top()]=co;vis[st.top()]=0;st.pop();}
st.pop();
}
}
int main()
{
while(scanf("%d%d",&n,&m)!=EOF)
{
num=0;rans=cans=co=dfnum=0;
memset(ru,0,sizeof(ru));
memset(cu,0,sizeof(cu));
memset(edge,0,sizeof(edge));
memset(head,0,sizeof(head));
memset(dfn,0,sizeof(dfn));
memset(low,0,sizeof(low));
memset(col,0,sizeof(col));
memset(vis,0,sizeof(vis));
while( !st.empty() )
st.pop();
for(int i=1;i<=m;i++)
{
int x,y;
scanf("%d%d",&x,&y);
add(x,y);
}
for(int i=1;i<=n;i++)
if(!col[i])
tarjian(i);
for(int i=1;i<=num;i++)
{
if(col[edge[i].from]!=col[edge[i].to])
ru[col[edge[i].to]]++,cu[col[edge[i].from]]++;
}
for(int i=1;i<=co;i++)
{
if(ru[i]==0) rans++;
if(cu[i]==0) cans++;
}
if(co!=1)
printf("%d\n",max(rans,cans));
else printf("0\n");
}
return 0;
}