HDU4612 Warm up 无向图的连通图缩点+树的直径

这道题乍一看一点思路都没有QAQ...我好菜啊55555...

直接搜了题解，发现并没有辣么难。。无向图的连通图缩点后，变成一棵树（就像有向图缩点后变成一个有向无环图），加一条边一定会构成一个环，为了使桥更少，那么这个环应该包含到更多的点或边，那就是树的直径了。这样的话，ans=缩点后的边数-直径边数。

一开始写了一发，RE了。。因为建树过程中v结点弄错了，应该是edge[i].to所属的联通子图。改过再一交，MLE。。发现我用了pair存储边的权值1，其实用不着，每次更新dis数组，直接+1就行，不用另外开空间（详见代码）。改过后再交发现WA了QAQ...怎么都不知道自己错哪儿了555...再一搜发现是有重边的。。。心情复杂.jpg

用了POJ3177的第二种方法处理重边，不过要注意去掉iscut这个变量，不然会MLE。。。

（详见我的另一篇博客：https://blog.csdn.net/Cc_Sonia/article/details/82694394）。

附上AC代码：

#include<cstdio>
#include<cstring>
#include<cstdlib>
#include<cmath>
#include<iostream>
#include<algorithm>
#include<vector>
#include<map>
using namespace std;
#define pb push_back
#define mkp make_pair
typedef pair<int,int>pp;

const int MAXN=200005;
const int MAXM=1000002*2;//因为是无向图,所以边数×2
struct Edge
{
    int to,next;
    //bool iscut;//是否为桥
}edge[MAXM];
int head[MAXN],tot;
int low[MAXN],dfn[MAXN],st[MAXN];
int be[MAXN];
int inde,top;
int block;//边双联通块数
bool inst[MAXN];
int bridge;

//vector<pp>tr[MAXN];
vector<int>tr[MAXN];
bool vis[MAXN];
int dis[MAXN];

void init()
{
    tot=0;
    memset(head,-1,sizeof(head));
}
void addedge(int u,int v)
{
    edge[tot].to=v;edge[tot].next=head[u];//edge[tot].iscut=false;
    head[u]=tot++;
}
void Tarjan(int u,int pre)
{
    int v;
    low[u]=dfn[u]=++inde;
    st[top++]=u;
    inst[u]=true;
    int sign=0;
    for(int i=head[u];i!=-1;i=edge[i].next)
    {
        v=edge[i].to;
        if(v==pre&&sign==0)//注意处理重边！
        {
            sign++;
            continue;
        }
        if(!dfn[v])
        {
            Tarjan(v,u);
            if(low[u]>low[v])
                low[u]=low[v];
            if(low[v]>dfn[u])
            {
                bridge++;
                //edge[i].iscut=true;
                //edge[i^1].iscut=true;
            }
        }
        else if(inst[v]&&low[u]>dfn[v])
             low[u]=dfn[v];
    }
    if(low[u]==dfn[u])
    {
        block++;
        do
        {
            v=st[--top];
            inst[v]=false;
            be[v]=block;
        }while(v!=u);
    }
}

void solve(int n)
{
    memset(dfn,0,sizeof(dfn));
    memset(inst,false,sizeof(inst));
    inde=top=0;
    bridge=0;block=0;
    Tarjan(1,0);
    //cout<<"bridge="<<bridge<<endl;

    for(int i=1;i<=n;i++)
        tr[i].clear();
    int u,v;
    for(int i=1;i<=n;i++)
        for(int j=head[i];j!=-1;j=edge[j].next)
        {
            //if(edge[j].iscut)//建树
            u=be[i];v=be[edge[j].to];//注意是edge[j].to!
            if(u!=v)
            {
                //tr[u].pb(mkp(v,1));tr[v].pb(mkp(u,1));
                tr[u].pb(v);tr[v].pb(u);
            }
        }
}

void dfs(int u)
{
    for(int i=0;i<tr[u].size();i++)
    {
        int v=tr[u][i];
        if(!vis[v])
        {
            vis[v]=true;
            //dis[v]=dis[u]+tr[u][i].second;
            dis[v]=dis[u]+1;
            dfs(v);
        }
    }
}

int main()
{
    int n,m;
    while(scanf("%d%d",&n,&m)==2)
    {
        if(n==0&&m==0)
            break;
        init();
        int u,v;
        for(int i=1;i<=m;i++)
        {
            scanf("%d%d",&u,&v);
            addedge(u,v);
            addedge(v,u);//无向图
        }
        solve(n);

        memset(vis,false,sizeof(vis));
        memset(dis,0,sizeof(dis));
        int st=1;
        vis[st]=true;
        dfs(st);
        int pos=-1,ma=-1;
        for(int i=1;i<=block;i++)//注意是block！
        {
            if(dis[i]>ma)
            {
                pos=i;
                ma=dis[i];
            }
        }

        memset(vis,false,sizeof(vis));
        memset(dis,0,sizeof(dis));
        st=pos;
        vis[st]=true;
        dfs(st);
        ma=-1;
        for(int i=1;i<=block;i++)
            ma=max(ma,dis[i]);
        //cout<<"ma="<<ma<<endl;

        printf("%d\n",bridge-ma);
    }
    return 0;
}