BZOJ-1339/1163:[Baltic2008]Mafia

[文章目录]

Description

给你一张无向图,n个点,m条边。给出每个点的代价,询问最少花多少的代价控制点能够使得a-b不通过被控制的点不能够互相到达。n<=200,m<=20000

最大流最小割。
将每个点拆点,流量为代价。那么原图的最小割即为使得a-b不能连通的最小代价。
一开始脑残点数组开了210.殊不知拆点需要开两倍。

#include <queue>
#include <cstdio>
#include <cstring>
#include <algorithm>
using namespace std;
int n,m,s,t,ans;
const int inf=2000000000;
int head[410],to[101000],nxt[101000],f[101000],cnt=1;
inline void add(int x,int y,int z)
{
    to[++cnt]=y;
    nxt[cnt]=head[x];
    head[x]=cnt;
    f[cnt]=z;
    to[++cnt]=x;
    nxt[cnt]=head[y];
    head[y]=cnt;
    f[cnt]=0;
}
int dis[410];
queue<int>q;
bool bfs()
{
    memset(dis,-1,sizeof dis);
    while(!q.empty()) q.pop();
    dis[s]=0; q.push(s);int x;
    while(!q.empty())
    {
        x=q.front(); q.pop();
        for(int i=head[x];i;i=nxt[i])
            if(f[i]>0&&dis[to[i]]<0)
            {
                dis[to[i]]=dis[x]+1;
                if(to[i]==t) return true;
                q.push(to[i]);
            }
    }
    return false;
}
int dinic(int x,int flow)
{
    if(x==t) return flow;
    int xx,tmp=flow;
    for(int i=head[x];i;i=nxt[i])
        if(f[i]>0&&dis[to[i]]==dis[x]+1)
        {
            xx=dinic(to[i],min(f[i],tmp));
            if(!xx) dis[to[i]]=-1;
            f[i]-=xx; f[i^1]+=xx; tmp-=xx;
            if(!tmp) break;
        }
    return flow-tmp;
}
int main()
{
    scanf("%d%d%d%d",&n,&m,&s,&t); t=t+n;
    int x,y;
    for(int i=1;i<=n;i++)
    {
        scanf("%d",&x);
        add(i,i+n,x);
    }
    while(m--)
    {
        scanf("%d%d",&x,&y);
        add(x+n,y,inf); add(y+n,x,inf);
    }
    while(bfs()) ans+=dinic(s,inf);
    printf("%d",ans);
    return 0;
}

发表评论

邮箱地址不会被公开。 必填项已用*标注