Telecowmunication_usaco5.4_最小点割集

Description

---

农夫约翰的奶牛们喜欢通过电邮保持联系，于是她们建立了一个奶牛电脑网络，以便互相交流。这些机器用如下的方式发送电邮：如果存在一个由c台电脑组成的序列a1,a2,…,a(c)，且a1与a2相连，a2与a3相连，等等，那么电脑a1和a(c)就可以互发电邮。
很不幸，有时候奶牛会不小心踩到电脑上，农夫约翰的车也可能碾过电脑，这台倒霉的电脑就会坏掉。这意味着这台电脑不能再发送电邮了，于是与这台电脑相关的连接也就不可用了。
有两头奶牛就想：如果我们两个不能互发电邮，至少需要坏掉多少台电脑呢？请编写一个程序为她们计算这个最小值和与之对应的坏掉的电脑集合。
以如下网络为例：

  1 
 / 
3-2 

这张图画的是有2条连接的3台电脑。我们想要在电脑1和2之间传送信息。电脑1与3、2与3直接连通。如果电脑3坏了，电脑1与2便不能互发信息了。

Input

---

第一行四个由空格分隔的整数：N,M,c1,c2.N是电脑总数(1<=N<=100)，电脑由1到N编号。M是电脑之间连接的总数(1<= M<=600)。最后的两个整数c1和c2是上述两头奶牛使用的电脑编号。连接没有重复且均为双向的(即如果c1与c2相连，那么c2与c1也相连)。两台电脑之间至多有一条连接。电脑c1和c2不会直接相连。 第2到M+1行 接下来的M行中，每行包含两台直接相连的电脑的编号。

Output

---

输出共有两行。第一行是使电脑c1和c2不能互相通信需要坏掉的电脑数目的最小值。第二行是排好序的坏掉的电脑的编号列表。注意c1和c2都不能坏掉。如果有多种可能情况，输出第一个数最小的一种，如果第一个数相同，则输出第二个数最小的一种，依此类推。

Analysis

---

想到拆点，但是具体仍然不知道怎么做，请教了一发幂萎
大概就是把每个点拆成i和i’，然后用i作为入点i’作为出点连权为INF的变，用权为1的边连接他们。这样就把点割变成了边割，分别删掉n条边跑n次最大流判断是否对最大流有影响就可以了

Code

---

/*
ID:wjp13241
PROG:telecow
LANG:C++
*/
#include <cstdio>
#include <cstring>
#include <cstdlib>
#include <algorithm>
#include <iostream>
#include <vector>
#include <cmath>
#include <stack>
#include <queue>
#define fo(i,a,b) for(int i=a;i<=b;i++)
#define dfo(i,a,b) for(int i=a;i>=b;i--)
#define fore(i,x,e) for(int i=ls[x];i;i=e[i].next)
#define fil(x,t) memset(x,t,sizeof(x))
#define STP system("pause")
#define min(x,y) x<y?x:y
#define max(x,y) x>y?x:y
#define PuB(v,x) v.push_back(x)
#define PoB(v) v.pop_back()
#define ld long double
#define ll long long
#define db double
#define INF 0x3f3f3f3f
#define LIM 100000000
#define EPS 1e-4
#define N 201
#define E N*20+1
using namespace std;
int c[N][N],bac[N][N],vis[N],p[E],cnt=0;
int add(int x,int y,int w){c[x][y]+=w;bac[x][y]+=w;}
int rec(int n)
{
    fo(i,1,n)
        fo(j,1,n)
            c[i][j]=bac[i][j];
}
int fnd(int x,int n,int ed,int mn)
{
    vis[x]=1;
    if (x==ed)
        return mn;
    fo(i,1,n)
        if (c[x][i]>0&&!vis[i])
        {
            int d=fnd(i,n,ed,min(mn,c[x][i]));
            if (d!=0)
            {
                c[x][i]-=d;
                c[i][x]+=d;
                return d;
            }
        }
    return 0;
}
int main()
{
//  freopen("telecow.in","r",stdin);
//  freopen("telecow.out","w",stdout);
    ios::sync_with_stdio(false);
    int n,m,st,ed;
    cin>>n>>m>>st>>ed;
    fil(c,0);
    fo(i,1,m)
    {
        int x,y;
        cin>>x>>y;
        add(x+n,y,INF);
        add(y+n,x,INF);
    }
    fo(i,1,n)
        add(i,i+n,1);
    rec(n*2);
    int mxFlow=0;
    fil(vis,0);
    while (int tmp=fnd(st+n,n*2,ed,INF*2))
    {
        fil(vis,0);
        mxFlow+=tmp;
    }
    cout<<mxFlow<<endl;
    bool first=true;
    fo(i,1,n)
        if (i^st&&i^ed)
        {
            rec(n*2);
            fo(j,1,n*2)
                c[i+n][j]=c[i][j]=c[j][i+n]=c[j][i]=0;
            int tot=0;
            fil(vis,0);
            while (int tmp=fnd(st+n,n*2,ed,INF*2))
            {
                fil(vis,0);
                tot+=tmp;
            }
            if (tot+1==mxFlow)
            {
                if (!first)
                    cout<<" ";
                else
                    first=false;
                cout<<i;
                p[i]=1;
                mxFlow--;
                fo(j,1,n*2)
                    bac[i+n][j]=bac[i][j]=bac[j][i+n]=bac[j][i]=0;
            }
        }
    cout<<endl;
    return 0;
}