hdu4289最小割

题意：一个小偷要从s市逃到d市，警察可以在任意一个城市逮捕他们，但是需要花费不同的钱，问在保证一定能够抓到小偷的情况下的最小花费是多少。

分析：每个城市都有一个花费，故进行拆点，在a,a+n之间连一套容量为cost[i]的线段。另外，两个相连通的城市，在(a+n,b),(b+n,a)之间连一条容量为无穷的边，最后求从s到d+n的最大流即可。

#include <iostream>
#include <cstdio>
using namespace std;
const int oo=1e9;
const int mm=111111;
const int mn=999;
int node,src,dest,edge;
int ver[mm],flow[mm],next[mm];
int head[mn],work[mn],dis[mn],q[mn];
void prepare(int _node,int _src,int _dest)//初始化
{
    node=_node,src=_src,dest=_dest;
    for(int i=0;i<=node;i++) head[i]=-1;
    edge=0;
}
void addedge(int u,int v,int c)//加边
{
    ver[edge]=v;flow[edge]=c;next[edge]=head[u];head[u]=edge++;
    ver[edge]=u;flow[edge]=0;next[edge]=head[v];head[v]=edge++;
}
bool Dinic_bfs()//将一个网络分层
{
    int i,u,v,l,r=0;
    for(i=0;i<=node;i++)
        dis[i]=-1;
    dis[q[r++]=src]=0;
    for(l=0;l<r;l++)
        for(i=head[u=q[l]];i>=0;i=next[i])
            if(flow[i]&&dis[v=ver[i]]<0)
            {
                dis[q[r++]=v]=dis[u]+1;
                if(v==dest) return 1;
            }
    return 0;
}
int Dinic_dfs(int u,int exp)//多路增广，从源点开始，用dfs从前一层向后一层反复寻找增广路
{
    if(u==dest) return exp;
    for(int &i=work[u],v,tmp;i>=0;i=next[i])
    {
        if(flow[i]&&dis[v=ver[i]]==dis[u]+1&&(tmp=Dinic_dfs(v,min(exp,flow[i])))>0)
        {
            flow[i]-=tmp;//如果该边不是最佳路径中的边，回溯时需要使用反向边
            flow[i^1]+=tmp;
            return tmp;
        }
    }
    return 0;
}
int Dinic_flow()//求最大流
{
    int i,ret=0,delta;
    while(Dinic_bfs())
    {
        for(i=0;i<=node;i++)
            work[i]=head[i];
        while(delta=Dinic_dfs(src,oo))  ret+=delta;
    }
    return ret;
}
int main()
{
    int n,m,src,dest,i,j,a,b;
    while(~scanf("%d%d",&n,&m))
    {
        scanf("%d%d",&src,&dest);
        prepare(2*n,src,dest+n);
        for(i=1;i<=n;i++)
        {
            scanf("%d",&a);
            addedge(i,i+n,a);
        }
        for(i=1;i<=m;i++)
        {
            scanf("%d%d",&a,&b);
            addedge(a+n,b,oo);
            addedge(b+n,a,oo);
        }
        printf("%d\n",Dinic_flow());
    }
    return 0;
}