hdu 3671 Boonie and Clyde

双连通分量

题意：给一个无向图，要求毁掉两个点，使图变得不连通，图一开始是连通的

因为要毁掉两个点，就不是简单的求割点，再看看数据范围，点数为1000，边数为10000，Tarjan的时间复杂度为O(E)，如果用枚举法，先枚举要毁掉的第一个点，再用Tarjan进行处理来找割点会不会超时呢？答案是不会，时间为O(v*E)，刚好是千万级别，不超

做法：先枚举要删除的第1个点，在原图中删除它，看看删除它后整个图的变化

        1.整个图变得不连通了（即这个点本身是割点），但是还没完要分类讨论一下

        （1）.整个图变为两部分，但是两部分刚好都是一个点，那么这两个点再毁掉哪个点都好，图的连通分支数都不会增加，这是一个特殊情况

　　　　　　例如，（1,2）（2,3）这种图，是无解的，任意毁掉两个点都无法增加图的连通分支，所以方案数为0

        （2）.整个图分为两部分，但是有一部分的点数为1，另一部分大于1，那么这时候只要在较大的那部分，任意毁掉一个点(无论是不是割点都行)，最后整个图都会至少被分为了两个部分

                （如果毁掉的是割点，将分成更多份），所以这样产生的方案数是V-2

        （3）.整个图分为了两个部分，两个部分的点数都大于1，那么任意在哪个部分毁掉那个点都可以（无论是不是割点都行），最后整个图都会至少分为两个部分，所以方案数为V-1

        （4）.整个图被分为了三个或更多的部分，那么也是在剩下的点中任意毁掉一个点都可以（无论那个点是不是割点），方案数为V-1

              （如果这个点刚好处于一个部分且这个部分只有它自己一个点，那么    毁掉后整个图的分支数减1；如果这个点在一个部分且这个部分不止它一个点且这个点不是割点，那么分支数  不会增加，如果是割点分支数为增加）

     2.删除第一个点后，整个图还是连通的（是连通，不是双连通）

        那么就在剩下的图中找割点，找到几个，方案数就是多少

 

 

最后注意一点，这样计算的结果，很容易想到是有重复的，但是不难想到，其实刚好重复了一次，因为对于一个图，方案是固定的，枚举了所有点，找出了所有已方案，相当于每个方案算了两次，最后答案除2即可

#include <iostream>

#include <cstdio>

#include <cstring>

using namespace std;

#define N 1010

#define M 10010



int n,tot,cancel,rts,num;

int e[N][M];

int dfn[N],low[N],cut[N],iscut[N],cnum,dcnt;



void dfs(int u ,int fa)

{

    dfn[u] = low[u] = ++dcnt;

    ++num;

    for(int i=1; i<=e[u][0]; i++)

    {

        int v = e[u][i];

        if(v == cancel || v == fa) continue;

        if(!dfn[v])

        {

            dfs(v,u);

            low[u] = min(low[u] , low[v]);

            if(fa == -1) {rts++; continue;}

            if(dfn[u] <= low[v])

            {

                if(!iscut[u])

                { iscut[u] = 1; cut[cnum++] = u;}

            }

        } 

        else low[u] = min(low[u] , dfn[v]);

    }

}



void solve()

{

    int res = 0,count,ok;

    for(cancel = 1; cancel<=n; cancel++) //枚举第1个删除的点

    {

        count = cnum = 0;

        memset(dfn,0,sizeof(dfn));

        memset(low,0,sizeof(low));

        memset(iscut,0,sizeof(iscut));

        ok = 0;

        for(int i=1; i<=n; i++)

            if(i!=cancel && !dfn[i])

            {

                count++;

                num = rts = 0;

                dfs(i , -1);

                if(num == 1) { ok++; continue; }

                if(rts >= 2) { iscut[i] = 1; cut[cnum++] = i; }

            }

        if(count >= 3) res += (n-1);

        else if(count == 2 && ok == 0)   res += (n-1);

        else if(count == 2 && ok == 1)   res += (n-2);

        else if(count == 1) res += cnum;

    }

    cout << res/2 << endl;

}



int main()

{

    int cas = 0;

    while(cin >> n >> tot)

    {

        if(!n && !tot) break;

        for(int i=1; i<=n; i++)

            e[i][0] = 0;

        while(tot--)

        {

            int u , v , c;

            cin >> u >> v;

            c = ++e[u][0];

            e[u][c] = v;

            c = ++e[v][0];

            e[v][c] = u;

        }

        printf("Case %d: ",++cas);

        solve();

    }

    return 0;

}