Ural 1277 cops ans thieves （最小割模型）

题目地址 ：http://acm.timus.ru/problem.aspx?space=1&num=1277

这里我们要拆点。把一个点拆成i,i‘ 。如何 i，j有边 ，在建边(i,j‘,inf),(j,i‘,inf)。 然后每个点点边(i‘,i,R[i])。这样建边以后，若要阻止 s到f的路径，那么必须破败一些边，那么我们为了是的边权最小，必须破坏边权小于inf的边，对应的就是图中拆点后的边(j’->j)  。实际上 这条边就代表了点j的点权。求最小割即是答案。 有一个需要注意的地方那个s==f的时候要特判。

VIEW CODE

#include<cstdio>
#include<algorithm>
#include<iostream>
#include<cmath>
#include<queue>
#include<stack>
#include<string>
#include<cstring>
#include<map>
#include<vector>
#include<set>
#include<ctime>
#include<stdlib.h>
using namespace std;
const int mmax= 210;
const int mod=1000000007;
const int inf=0x3fffffff;
using namespace std;

struct node
{
    int flow;
    int en;
    int next;
}E[40010+mmax];
int p[mmax];
int num;
void init()
{
    memset(p,-1,sizeof p);
    num=0;
}
void add(int st,int en,int flow)
{
    E[num].en=en;
    E[num].flow=flow;
    E[num].next=p[st];
    p[st]=num++;
    E[num].en=st;
    E[num].flow=0;
    E[num].next=p[en];
    p[en]=num++;
}

int d[mmax];
bool vis[mmax];
int qq[mmax];
int cur[mmax];
bool bfs(int st,int en)
{
    memset(vis,0,sizeof vis);
    int qcnt=0;
    qq[++qcnt]=st;
    d[st]=0;
    vis[st]=1;
    while(qcnt)
    {
        int x=qq[qcnt];
        qcnt--;
        for(int i=p[x]; i+1; i=E[i].next)
        {
            int v=E[i].en;
            if(!vis[v]&&E[i].flow)
            {
                vis[v]=1;
                qq[++qcnt]=v;
                d[v]=d[x]+1;
            }
        }
    }
    return vis[en];
}
int dfs(int st,int en,int  flow)
{
    if(st==en||flow==0)
        return flow;
    int f=0,dd;
    for(int &i=cur[st]; i+1;i=E[i].next)
    {
        int v=E[i].en;
        if(d[st]+1==d[v]&&(dd=dfs(v,en,min(flow,E[i].flow)))>0)
        {
            E[i].flow-=dd;
            E[i^1].flow+=dd;
            flow-=dd;
            f+=dd;
            if(flow==0)
                break;
        }
    }
    return f;
}
int dinic(int st,int en,int n)
{
    int flow=0;
    while(bfs(st,en))
    {
        for(int i=0;i<=n;i++)
            cur[i]=p[i];
        flow+=dfs(st,en,inf);
    }
    return flow;
}

int main()
{
    int k,n,m,s,f;
    while(cin>>k)
    {
        init();
        scanf("%d %d %d %d",&n,&m,&s,&f);
        for(int i=1;i<=n;i++)
        {
            int x;
            scanf("%d",&x);
            add(n+i,i,x);
        }
        for(int i=0;i<m;i++)
        {
            int u,v;
            scanf("%d %d",&u,&v);
            add(u,v+n,inf);
            add(v,u+n,inf);
        }
        if(s==f)
        {
            puts("NO");
            continue;
        }
        int ans=dinic(s,f+n,2*n);
        if(ans<=k)
            puts("YES");
        else
            puts("NO");
    }
    return 0;
}