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