Ural 1277 - Cops and Thieves 无向图的最小点割

                 题意:

                          在一个无向图中小偷要偷东西..小偷从s点出发..要偷的东在点e...警察可以用一些警力封锁一些点让小偷无论如何都不能到达e....现在告诉能派出的最大警力以及封锁一个点所需要的警力..注意s,e不能封锁..问能否保证小偷无论如何都不能到达e....

                 题解:

                          很裸的最小点割了..无向图的..一样..做双边..拆点...拆点之间的容量为这个点所需的警察数量...其他边都是容量无穷大...跑最大流即为最小割...恶心的是数据里有起点=终点的..特判NO..

Program:

#include<iostream>
#include<algorithm>
#include<stdio.h>
#include<string.h>
#include<math.h>
#include<queue>
#define MAXN 505
#define MAXM 50005
#define oo 1000000007
#define ll long long
using namespace std; 
struct Dinic      
{      
       struct node    
       {     
             int c,u,v,next;    
       }edge[MAXM];    
       int ne,head[MAXN];    
       int cur[MAXN], ps[MAXN], dep[MAXN];    
       void initial()    
       {    
             ne=2;    
             memset(head,0,sizeof(head));     
       }    
       void addedge(int u, int v,int c)    
       {     
             edge[ne].u=u,edge[ne].v=v,edge[ne].c=c,edge[ne].next=head[u];    
             head[u]=ne++;    
             edge[ne].u=v,edge[ne].v=u,edge[ne].c=0,edge[ne].next=head[v];    
             head[v]=ne++;    
       }    
       int MaxFlow(int s,int t)    
       {                         
             int tr, res = 0;    
             int i,j,k,f,r,top;    
             while(1)    
             {    
                    memset(dep, -1, sizeof(dep));    
                    for(f=dep[ps[0]=s]=0,r=1;f!= r;)    
                       for(i=ps[f++],j=head[i];j;j=edge[j].next)    
                         if(edge[j].c&&dep[k=edge[j].v]==-1)    
                         {    
                               dep[k]=dep[i]+1;    
                               ps[r++]=k;    
                               if(k == t){  f=r; break;  }    
                         }    
                    if(dep[t]==-1) break;    
                    memcpy(cur,head,sizeof(cur));    
                    i=s,top=0;    
                    while(1)    
                    {    
                         if(i==t)    
                         {    
                               for(tr=oo,k=0;k<top;k++)    
                                  if(edge[ps[k]].c<tr)    
                                     tr=edge[ps[f=k]].c;    
                               for(k=0;k<top;k++)    
                               {    
                                     edge[ps[k]].c-=tr;    
                                     edge[ps[k]^1].c+=tr;    
                               }    
                               i=edge[ps[top=f]].u;    
                               res+= tr;    
                         }    
                         for(j=cur[i];cur[i];j=cur[i]=edge[cur[i]].next)     
                             if(edge[j].c && dep[i]+1==dep[edge[j].v]) break;     
                         if(cur[i])  ps[top++]=cur[i],i=edge[cur[i]].v;     
                         else    
                         {    
                                 if(!top) break;    
                                 dep[i]=-1;    
                                 i=edge[ps[--top]].u;    
                         }    
                   }    
             }    
             return res;    
      }    
}T;    
int main() 
{           
      int sum,n,m,s,e,i,x,y; 
      while (~scanf("%d",&sum))
      {
               scanf("%d%d%d%d",&n,&m,&s,&e);
               s=s<<1|1,e=e<<1,T.initial(); 
               for (i=1;i<=n;i++) scanf("%d",&x),
                                  T.addedge(i<<1,i<<1|1,x),T.addedge(i<<1|1,i<<1,x);
               while (m--)
               {
                        scanf("%d%d",&x,&y);
                        T.addedge(x<<1|1,y<<1,oo),T.addedge(y<<1|1,x<<1,oo);
               }              
               if (T.MaxFlow(s,e)<=sum && s!=e+1) printf("YES\n");
                                            else  printf("NO\n");
      }
      return 0;
}