poj1815(最小割点集)

【题目大意】
现代社会人们都靠电话通信。A 与 B 能通信当且仅当 A 知道 B 的电话号或者 A
知道 C 的电话号且 C 与 B 能通信。若 A 知道 B 的电话号，那么 B 也知道 A 的电
话号。然而不好的事情总是会发生在某些人身上，比如他的电话本丢了，同时他
又换了电话号，导致他跟所有人失去了联系。现在给定 N 个人之间的通信关系
以及特定的两个人 S 和 T，问最少几个人发生不好的事情可以导致 S 与 T 无法通
信并输出这些人。如果存在多组解，输出字典序最小的一组。（2 <= N <= 200）
【建模方法】
此题是求一个最小点割集，但是要求输出方案且要字典序最小。首先仍是拆点构
图，求一次最小割记为 ans。之后我们采取贪心的策略，从 1 至 N 枚举删点，如
果最小割不变，说明该点不可能在最小割中，我们再把该点加入到网络中；否则
最小割一定变小了，用这个较小值更新 ans，记录该点是一个解并不再将其放回。

重复这个过程，把所有点都扫描一遍后结果就

#include<cstdio>
#include<cstring>
#include<vector>
#include<queue>
using namespace std;
int bmap[205][205];
typedef int cap_type;
#define MAX_V 800 + 30 + 16
struct edge
{
	int to, rev;
	cap_type cap;

	edge(int to, cap_type cap, int rev) : to(to), cap(cap), rev(rev)
	{}
};

vector <edge> G[MAX_V];
int level[MAX_V];
int iter[MAX_V];
void add_edge(int from, int to, int cap)
{
	G[from].push_back(edge(to, cap, G[to].size()));
	G[to].push_back(edge(from, 0, G[from].size() - 1));
}
void bfs(int s)
{
	memset(level, -1, sizeof(level));
	queue<int> que;
	level[s] = 0;
	que.push(s);
	while (!que.empty())
	{
		int v = que.front();
		que.pop();
		for (int i = 0; i < G[v].size(); ++i)
		{
			edge &e = G[v][i];
			if (e.cap > 0 && level[e.to] < 0)
			{
				level[e.to] = level[v] + 1;
				que.push(e.to);
			}
		}
	}
}
cap_type dfs(int v, int t, cap_type f)
{
	if (v == t)
	{
		return f;
	}
	for (int &i = iter[v]; i < G[v].size(); ++i)
	{
		edge &e = G[v][i];
		if (e.cap > 0 && level[v] < level[e.to])
		{
			cap_type d = dfs(e.to, t, min(f, e.cap));
			if (d > 0)
			{
				e.cap -= d;
				G[e.to][e.rev].cap += d;
				return d;
			}
		}
	}

	return 0;
}
cap_type max_flow(int s, int t)
{
	cap_type flow = 0;
	for (;;)
	{
		bfs(s);
		if (level[t] < 0)
		{
			return flow;
		}
		memset(iter, 0, sizeof(iter));
		cap_type f;
		while ((f = dfs(s, t, 0x3f3f3f3f)) > 0)
		{
			flow += f;
		}
	}
}
int ans[205];
int vis[205];
int tot = 0;
int main()
{
    int n,s,t;
    while(scanf("%d%d%d",&n,&s,&t)!=EOF){
    for(int j=0;j<2*n+10;j++) G[j].clear(); tot = 0;
    for(int i=1;i<=n;i++) for(int j=1;j<=n;j++) scanf("%d",&bmap[i][j]);
    if(bmap[s][t]==1) {printf("NO ANSWER!\n"); return 0; }
    for(int i=1;i<=n;i++)
    {
        if(i==s||i==t) add_edge(i,i+n,0x3f3f3f3f);
        else add_edge(i,i+n,1);
    }
    for(int i=1;i<=n;i++)
        for(int j=i+1;j<=n;j++)
        {
            if(bmap[i][j]==0) continue;
            add_edge(i+n,j,0x3f3f3f3f);
            add_edge(j+n,i,0x3f3f3f3f);
        }
    int res = max_flow(s,t+n);
    memset(vis,0,sizeof(vis));
    printf("%d\n",res);
    for(int i=1;i<=n;i++)
    {
        if(i==s||i==t) continue;
        vis[i] = 1;
        for(int j=0;j<2*n+10;j++) G[j].clear();
        for(int j=1;j<=n;j++)
        {
            if(vis[j]) continue;
            if(j==s||j==t) add_edge(j,j+n,0x3f3f3f3f);
            else add_edge(j,j+n,1);
        }
        for(int j=1;j<=n;j++)
        {
            if(vis[j]) continue;
            for(int k=1;k<=n;k++)
            {
                if(vis[k]||!bmap[j][k]||j==k) continue;
                add_edge(j+n,k,0x3f3f3f3f);
            }
        }
        int tmp = max_flow(s,t+n);
        //printf("%d %d %d\n",i,tmp,res);
        if(tmp<res) ans[tot++] = i,res = tmp;
        else vis[i] = 0;
    }
    if(tot>0){
    for(int i=0;i<tot;i++)
    {
        if(i!=0) printf(" ");
        printf("%d",ans[i]);
    }
    printf("\n");}
    }
    return 0;
}