BZOJ1163: [Baltic2008]Mafia

最新推荐文章于 2019-05-22 22:29:49 发布

原创最新推荐文章于 2019-05-22 22:29:49 发布 · 412 阅读

0 ·

本内容遵循CC 4.0 BY-SA版权协议

标签

#最大流最小割 #BZOJ

蒟蒻zxl的Blog专栏同时被 3 个专栏收录

460 篇文章

订阅专栏

BZOJ

281 篇文章

订阅专栏

图论---网络流最小割

6 篇文章

订阅专栏

本文介绍了一种解决最大流最小割问题的经典算法实现。通过拆点技巧优化网络流问题求解过程，并提供了完整的C++代码示例。该算法适用于寻找网络中最大流及其对应的最小割情况。

最大流最小割

题目传送门

仍然是套路拆点，就不细讲了。

注意要从s的入点跑到t的出点。

代码：

#include<cstdio>
#include<cstring>
#include<algorithm>
#define MAXN 400
#define MAXM 20000
#define inf 0x7fffffff
using namespace std;
struct edge{
    int next,to,v,flow;
}ed[MAXM*6+5];
int n,m,k,s,t;
int h[MAXN+5],cop[MAXN+5],dis[MAXN+5],que[MAXN+5];
bool f[MAXN+5];
inline char readc(){
    static char buf[100000],*l=buf,*r=buf;
    if (l==r) r=(l=buf)+fread(buf,1,100000,stdin);
    if (l==r) return EOF; return *l++;
}
inline int _read(){
    int num=0; char ch=readc();
    while (ch<'0'||ch>'9') ch=readc();
    while (ch>='0'&&ch<='9') { num=num*10+ch-48; ch=readc(); }
    return num;
}
void addedge(int x,int y,int z){
    ed[k].next=h[x]; ed[k].to=y; ed[k].v=z; h[x]=k++;
    ed[k].next=h[y]; ed[k].to=x; ed[k].v=0; h[y]=k++;
}
bool bfs(){
    memset(f,false,sizeof(f));
    int r=0,w=1; que[1]=s; f[s]=true; dis[s]=0;
    while (r<w){
        int x=que[++r];
        for (int i=h[x];~i;i=ed[i].next)
            if (!f[ed[i].to]&&ed[i].v>ed[i].flow){
                int v=ed[i].to;
                que[++w]=v; f[v]=true; dis[v]=dis[x]+1;
            }
    }
    return f[t];
}
int dfs(int x,int rem){
    if (x==t||rem==0) return rem; int sum=0,p;
    for (int &i=cop[x];~i;i=ed[i].next)
        if (dis[ed[i].to]==dis[x]+1&&(p=dfs(ed[i].to,min(rem,ed[i].v-ed[i].flow)))){
            sum+=p,ed[i].flow+=p,ed[i^1].flow-=p,rem-=p;
            if (!rem) break;
        }
    return sum;
}
int maxflow(){
    int ans=0;
    while (bfs()){
        memcpy(cop,h,sizeof(cop));
        ans+=dfs(s,inf);
    }
    return ans;
}
int main(){
    memset(h,-1,sizeof(h));
    n=_read(); m=_read(); s=_read(); t=_read();
    for (int i=1;i<=n;i++){
        int d=_read(); addedge(i,i+n,d);
    }
    for (int i=1;i<=m;i++){
        int u=_read(),v=_read();
        addedge(u+n,v,inf); addedge(v+n,u,inf);
    }
    t+=n;
    return printf("%d\n",maxflow()),0;
}