Usaco 5.4.5

这道题要求删去最少的点，使得网络中给定两点不连通。

这个最少的点在图论中有具体的定义，有名称为点连通度。

那么，如何求这个点连通度呢？

我们可以考虑从前面的知识进行转化，前面有一道题目是让我们求边连通度。

那么从边连通度转化为点连通度是否可行呢？答案是可行的。

可以把在网络中的每个点p(这里p代表其中一个点）,拆成两个点，p1,p2，在p1,p2之间连一条边，其容量为1.

那么对于原来就在网络中相连的两点（又因为该题为无向图），假设这两点为p,q，拆点后形成p1,p2,q1,q2,

我们可以在p2和q1之间连一条容量为INF（maxint),在q2和p1之间也连一条容量为INF的边。

做一遍最大流，就可以得出点连通度。

另外，题目还要求我们输出字典序最小的组成大小为点连通度的点割集。

这其实也可以借鉴前面那一道求边连通度的题目。

边连通度的题目，我们尝试按升序去删除每一条边，如果改变删除后，网络中的最大流减小，说名改边存在于边割集中。

同理，因为这道题目我们拆点了，并且拆点后的两个子点之间有边相连，所以我们可以尝试按升序去删除每一条拆点拆出来的边，

如果做最大流后流量减少了1（因为拆点出来的边容量为1），那么该点便属于点割集，一直做到最大流为0为止。

洋洋洒洒写了这么多，只是希望大家能懂。

/*
ID: volz.kz.g
PROB: telecow
LANG: C++
*/
#include <iostream>
#include <fstream>
#include <cstdio>
#include <cstdlib>
#include <cstring>
#include <cmath>
#define INF 0x7FFFFFFF
using namespace std;
ifstream fin("telecow.in");
ofstream fout("telecow.out");

int n,m,s,t,ans;
struct edge_node{
    int v,c,next;
}e[6000];
int edge_num;
int d[6000],c[6000],v[6000];
int g[501][501];
bool cut[501];

void init(){
    memset(e,0,sizeof(e));
    memset(v,0,sizeof(v));
    edge_num=0;
}
void insert(int a,int b,int c){
    //fout << a << " " << b << " " << c << endl;
    ++edge_num;
    e[edge_num].v=b;e[edge_num].c=c;e[edge_num].next=v[a];v[a]=edge_num;
    ++edge_num;
    e[edge_num].v=a;e[edge_num].c=0;e[edge_num].next=v[b];v[b]=edge_num;
}
void makegraph(){
    init();
    for (int i=1;i<=n;++i)
	for (int j=1;j<=n;++j)
	    if (g[i][j]){
		insert(i+n,j,INF);
		insert(j+n,i,INF);
	    }
    for (int i=1;i<=n;++i)
	if (cut[i])
	    insert(i,i+n,1);
	else 
	    insert(i,i+n,0);
}
int opp(int x){
    return (x%2==1)?x+1:x-1;
}
int maxflow(int x,int lim){
    int ret,min_num,max_flow,i,j;
    if (x==t) return lim;
    min_num=2*n;max_flow=0;
    j=v[x];
    while (j!=0){
	if (e[j].c>0){
	    if (d[e[j].v]+1==d[x]){
		ret=min(lim,e[j].c);
		ret=maxflow(e[j].v,ret);
		e[j].c-=ret;e[opp(j)].c+=ret;
		lim-=ret;max_flow+=ret;
		if (lim==0 || d[s]>=2*n) return max_flow;
	    }
	    min_num=min(min_num,d[e[j].v]+1);
	}
	j=e[j].next;
    }
    if (min_num==d[x]) return max_flow;
    --c[d[x]];
    if (c[d[x]]==0) d[s]=2*n;
    d[x]=min_num;c[d[x]]++;
    return max_flow;
}
int main(){
    memset(g,false,sizeof(g));
    memset(cut,true,sizeof(cut));
    fin >> n >> m >> s >>t;
    for (int i=1;i<=m;++i){
	int a,b;fin >> a >> b;
	g[a][b]=true;g[b][a]=true;
    }
    makegraph();
    c[0]=2*n;
    while (d[s]<2*n) ans+=maxflow(s+n,INF);
    fout << ans << endl;//输出点连通度
    int print[201],print_num=0;
    for (int i=1;i<=n;++i){
	if (i!=s && i!=t){
	    cut[i]=false;
	    makegraph();memset(c,0,sizeof(c));memset(d,0,sizeof(d));
	    int tmp=0;c[0]=2*n;
	    while (d[s]<2*n) tmp+=maxflow(s+n,INF);
	    if (tmp==ans-1){
		ans--;
		print[++print_num]=i;
	    }
	    else{
		cut[i]=true;
	    }
	    if (!ans) break;
	}
    }
    for (int i=1;i<print_num;++i) fout << print[i] << " ";fout << print[print_num] << endl;
    return 0;
}