洛谷P1345 奶牛的电信Telecowmunication

题目描述

农夫约翰的奶牛们喜欢通过电邮保持联系，于是她们建立了一个奶牛电脑网络，以便互相交流。这些机器用如下的方式发送电邮：如果存在一个由 cc 台电脑组成的序列a_1,a_2,\cdots ,a_ca1​,a2​,⋯,ac​，且 a_1a1​ 与 a_2a2​ 相连，a_2a2​ 与 a_3a3​ 相连，等等。那么电脑 a_1a1​ 和 a_cac​ 就可以互发电邮。

很不幸，有时候奶牛会不小心踩到电脑上，农夫约翰的车也可能碾过电脑，这台倒霉的电脑就会坏掉。这意味着这台电脑不能再发送电邮了，于是与这台电脑相关的连接也就不可用了。

有两头奶牛就想：如果我们两个不能互发电邮，至少需要坏掉多少台电脑呢？请编写一个程序为她们计算这个最小值。

以如下网络为例：

   1*
  /
 3 - 2*

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

输入格式

第一行：四个由空格分隔的整数：N,M,c_1,c_2N,M,c1​,c2​。NN 是电脑总数，电脑由 11 到 NN 编号。MM 是电脑之间连接的总数。后面的两个整数 c_1c1​ 和 c_2c2​ 是上述两头奶牛使用的电脑编号。连接没有重复且均为双向的（即如果 c_1c1​ 与 c_2c2​ 相连，那么 c_2c2​ 与 c_1c1​ 也相连）。两台电脑之间至多有一条连接。电脑 c_1c1​ 和 c_2c2​ 不会直接相连。

第 22 到 M+1M+1 行：接下来的 MM 行中，每行包含两台直接相连的电脑的编号。

输出格式

一行，一个整数，表示使电脑 c_1c1​ 和 c_2c2​ 不能互相通信需要坏掉的电脑数目的最小值。

输入输出样例

输入 #1复制

3 2 1 2
1 3
2 3

输出 #1复制

1

说明/提示

对于 100\%100% 的数据：1\le N \le 1001≤N≤100，1\le M \le 6001≤M≤600。

上代码：

#include <bits/stdc++.h>
    #define vc vector
    #define INF ((int)(1e9))
    #define LINF ((ll)(1e18))
    #define pb push_back
    #define mp make_pair
    #define ll long long
    #define _tp template
    #define _tyn typename
    #define sint short int
    #define ull unsigned ll
    #define pii pair<int,int>
    #define uint unsigned int
    #define ms(_data) memset(_data,0,sizeof(_data))
    #define fin(_filename) freopen(_filename,"r",stdin)
    #define fout(_filename) freopen(_filename,"w",stdout)
    #define msn(_data,_num) memset(_data,_num,sizeof(_data))
    using namespace std;
    _tp<_tyn T>void mymax( T &_a , T _b ){ if( _a < _b ) _a = _b; }
    _tp<_tyn T>void mymin( T &_a , T _b ){ if( _a > _b ) _a = _b; }
    void print(int _x){printf("%d\n",_x);}
    void print(ll _x){printf("%I64d ",_x);}
    #define il inline
    il int in(){
        char c = getchar();
        int ret = 0;
        while( c < '0'  ||  c > '9' ) c = getchar();
        while( c >= '0'  &&  c <= '9' ){
            ret *= 10;
            ret += c-'0';
            c = getchar();
        }
        return ret;
    }
    il void read( int &x ){
        x = in();
    }
    il void read( int &x, int &y ){
        x = in(); y = in();
    }
    il void read( int &x1 , int &x2 , int &x3 ){
        x1 = in(); x2 = in(); x3 = in();
    }
    il void read( int &x1 , int &x2 , int &x3 , int &x4 ){
        x1 = in(); x2 = in(); x3 = in(); x4 = in();
    }
    /
    /
    #define MAXN 110
    struct Edge{
        int to,cap,rev;
        Edge(){}
        Edge( int tt , int cc , int rr ){
            to = tt;
            cap = cc;
            rev = rr;
        }
    };
    int n,m,s,t;
    vc<Edge> e[MAXN<<1];
    vc<Edge> inp[MAXN<<1];
    il void addedge( int f , int t , int c ){
        inp[f].pb( Edge(t,c,inp[t].size()) );
        inp[t].pb( Edge(f,0,inp[f].size()-1) );
        e[f].pb(Edge());
        e[t].pb(Edge());
    }
    il void f5(){
        for( int i = 1 ; i <= n+n+1 ; i++ )
            for( uint j = 0 ; j < inp[i].size() ; j++ )
                e[i][j] = inp[i][j];
    }
    int lev[MAXN<<1];    //Level
    il void bfs( int st ){
        ms(lev);
        queue<int> q;
        q.push(st);
        lev[st] = 1;
        while( !q.empty() ){
            int now = q.front();
            q.pop();
            for( uint i = 0 ; i < e[now].size() ; i++ ){
                if( e[now][i].cap <= 0 ) continue;
                int v = e[now][i].to;
                if( !lev[v] ){
                    lev[v] = lev[now]+1;
                    q.push(v);
                }
            }
        }
    }
    bool vis[MAXN<<1];
    il int dfs( int pos , int flow ){
        if( pos == t ) return flow;
        for( uint i = 0 ; i < e[pos].size() ; i++ ){
            Edge x = e[pos][i];
            int u = x.to;
            if( vis[u] ) continue;
            if( x.cap <= 0 ) continue;
            if( lev[u] != lev[pos]+1 ) continue;
            vis[u] = 1;
            int tans = dfs( u , min(flow,x.cap) );
            vis[u] = 0;
            if( tans > 0 ){
                e[pos][i].cap -= tans;
                e[u][x.rev].cap += tans;
                return tans;
            }
        }
        return 0;
    }
    int ans = 0;
    int main(){
        read(n,m,s,t);
        for( int i = 1 ; i <= n ; i++ ){
            addedge(i,n+i,1);
        }
        for( int i = 1 ; i <= m ; i++ ){
            int a,b;
            read(a,b);
            addedge(a+n,b,INF);
            addedge(b+n,a,INF);
        }
        f5();
        while(1){
            bfs(s+n);
            if( !lev[t] ) break;
            ms(vis);
            int tans;
            while( ( tans = dfs(s+n,INF) ) > 0 ){
                ans += tans;
                ms(vis);
            }
            //cout << "[Tans] " << tans << endl;
        }
        printf("%d\n",ans);
        return 0;
    }