标签:缩点 init eof names code back for print main
#pragma comment(linker, "/STACK:1024000000,1024000000")
#include<cstdio>
#include<cstdlib>
#include<algorithm>
#include<iostream>
#include<cstring>
#include<vector>
using namespace std;
const int maxn=200010;
const int maxm=2000010;
int Laxt[maxn],Next[maxm],To[maxm],cnt,vis[maxn];
int dfn[maxn],low[maxn];
int times,ans,cut_cnt,n,m;
int scc[maxn],scc_cnt;
int dis[maxn],S;
int stk[maxn],top,Maxdis,Maxpos;
vector<int>G[maxn];
void _init()
{
memset(Laxt,0,sizeof(Laxt));
memset(dfn,0,sizeof(dfn));
memset(scc,0,sizeof(scc));
memset(vis,0,sizeof(vis));
ans=cut_cnt=top=scc_cnt=cnt=times=0;
}
void _add(int u,int v)
{
Next[++cnt]=Laxt[u];
Laxt[u]=cnt;
To[cnt]=v;
}
void _tarjan(int u,int v){
dfn[u]=low[u]=++times;
int num_v=0;
stk[top++]=u;
for(int i=Laxt[u];i;i=Next[i]){
if(To[i]==v){
num_v++;
if(num_v==1) continue;//保证重边
}
if(!dfn[To[i]]){
_tarjan(To[i],u);
if(low[u]>low[To[i]]) low[u]=low[To[i]];
if(dfn[u]<low[To[i]]) cut_cnt++;//割边,对应缩点后是边
}
else if(dfn[To[i]]<low[u]) low[u]=dfn[To[i]];
}
if(dfn[u]<=low[u]){//小于是个环,等于是个点,都要处理
G[++scc_cnt].clear();
for(;;){
int tmp=stk[--top];
scc[tmp]=scc_cnt;
if(tmp==u) break;
}
}
}
void _rebuild()//重建无环图
{
for(int i=1;i<=n;i++)
for(int j=Laxt[i];j;j=Next[j])
if(scc[i]!=scc[To[j]])
G[scc[i]].push_back(scc[To[j]]);
}
void _dfs(int u)
{
int i,L=G[u].size();
for(int i=0;i<L;i++)
if(!dis[G[u][i]]){
dis[G[u][i]]=dis[u]+1;
_dfs(G[u][i]);
}
}
void _findR()
{
memset(dis,0,sizeof(dis));
dis[1]=1;S=1;Maxdis=1;
_dfs(1);
for(int i=1;i<=scc_cnt;i++)
if(dis[i]>dis[S]) S=i;
memset(dis,0,sizeof(dis));
dis[S]=1;
_dfs(S);
for(int i=1;i<=scc_cnt;i++)
if(dis[i]>Maxdis) Maxdis=dis[i];
Maxdis--;
}
int main()
{
int i,j,k,u,v;
while(~scanf("%d%d",&n,&m)){
if(n==0&&m==0) return 0;
_init();
for(i=1;i<=m;i++){
scanf("%d%d",&u,&v);
_add(u,v);
_add(v,u);
}
for(i=1;i<=n;i++)
if(!dfn[i]) _tarjan(i,-1);
_rebuild();
_findR();
printf("%d\n",cut_cnt-Maxdis);
}
return 0;
}
标签:缩点 init eof names code back for print main
原文地址:http://www.cnblogs.com/hua-dong/p/7694069.html