POJ1966 ZOJ2182<无向图点连通度 Isap版>

求顶点连通度的方法 ：
(1)把原图中每个结点拆成2个结点v1,v2,顶点v1到v2有一条弧<v1,v2>,容量为1
(2)对于原图中每条边e(u,v);在网络中有两条弧e1=<u2,v1> ,e2=<u1,v2>容量都为inf
(3)固定源点A2,枚举每个汇点(原图中不与A点相邻的),记录最小独立轨数目,该值就是顶点连通度

刚刚看到这位大神的博客，才知道固定一个源点的做法是错误的。。

地址：http://www.cnblogs.com/wuyiqi/archive/2011/10/09/2203516.html

固定一个源点的做法是错误的，因为如果刚好枚举的源点是割顶集里面的点就不行了
如如下数据5 6 (0,1) (0,2) (1,2) (0,3) (0,4) (3,4)
答案应该是1
但如果没有两两枚举源点汇点，答案就成了5了
所以，还是要两两枚举点的 

#include <iostream>
#include <cstdio>
#include <cstring>
#include <cstdlib>
#include <algorithm>
using namespace std;
#define inf 0x7fffffff
#define N 105
#define M 10005
int cur[N],gap[N],pre[N],dis[N];
int size,head[N];
bool s[N][N];
struct Edge
{
    int v,w,next;
    Edge(){}
    Edge(int V,int W,int NEXT):v(V),w(W),next(NEXT){}
}edge[M];

void InsertEdge(int u,int v,int w)
{
    edge[size] = Edge(v,w,head[u]);
    head[u] = size++;
    edge[size] = Edge(u,0,head[v]);
    head[v] = size++;
}
void Init(int n)
{
    size = 0;
    memset(head,-1,sizeof(head));
    for(int j=0; j<n; j++)
    InsertEdge(j,j+n,1);
    for(int j=0; j<n; j++)
    for(int k=j+1; k<n; k++)
    {
        if(s[j][k])
        {
            InsertEdge(j+n,k,inf);
            InsertEdge(k+n,j,inf);
        }
    }
}
int Isap(int st,int ed,int n) { //Isap模板 
    for(int i=0; i<=n; i++) {
        cur[i] = head[i];
        gap[i] = dis[i] = 0;
    }
    int u = pre[st] = st;
    int aug = inf ,maxflow = 0;
    while(dis[st] < n) {
loop:
        for(int &i=cur[u]; i!=-1; i=edge[i].next) {
            int v = edge[i].v;
            if(edge[i].w && dis[u] == dis[v] + 1) {
                aug = min(aug,edge[i].w);
                pre[v] = u;
                u = v;
                if(v==ed) {
                    maxflow += aug;
                    for(u=pre[u]; v!=st; v=u,u=pre[u]) {
                        edge[cur[u]].w -= aug;
                        edge[cur[ u]^1].w += aug;
                    }
                    aug = inf;
                }
                goto loop;
            }
        }
        int mindis = n;
        for(int i=head[u]; i!=-1; i=edge[i].next) {
            int v = edge[i].v;
            if(edge[i].w && dis[v]<mindis) {
                cur[u] = i;
                mindis = dis[v];
            }
        }
        if(--gap[dis[u]]==0) break;
        gap[dis[u] = mindis+1]++;
        u = pre[u];
    }
    return maxflow;
}

int main()
{
    int n,m,nv;
    while(scanf("%d%d",&n,&m)!=EOF)
    {
        nv = 2*n;
        for(int i=0; i<n; i++)
        for(int j=0; j<n; j++)
        s[i][j] = i==j;
        for(int i=0; i<m; i++)
        {
            int u,v;
            scanf(" (%d,%d)",&u,&v);
            s[u][v] = s[v][u] = 1;
        }
        int ans = inf;
        for(int i=0; i<n; i++) 
        {
            for(int j=0; j<n; j++)
            {
                if(s[i][j]) continue;
                Init(n);
                ans = min(ans,Isap(i+n,j,nv));
            }
        }
        if(ans == inf)ans = n;
        printf("%d
",ans);
    }
    return 0;
}