hdu 3491 Thieves

题意：
有一群盗贼要从s城出发去h城，现在告诉你一共有n个城，m条路，每一个城都要派一定的人去把守，现在问你一共要多少人才能堵住盗贼，使之不能从s城到h城，并前要求你不能在s城和h城拦截盗贼。
这是一个最小割的模板题，首先我们把每个城i拆成两个城i和i+n，权值为他需要把守的人数，s到s+n和h到h+n拆点后权值为inf，因为题目要求不能在s城和h城拦截盗贼，两个城之间的路连一条边，权值为正无穷，根据最小割=最大流，跑一遍最大流即可。


#include<stdio.h>
#include<string.h>
#include<iostream>
#define N 1000010   
#define inf 1<<30 
using namespace std;
struct Edge{  
int u,v,c,next;  
}edge[N];  
int head[N],pre[N],cur[N],dis[N],gap[N];  
int n,m,s,h,e,a,b,num,t,start,end;  

void add_adge(int u,int v,int c)  
{  
   edge[e].u=u; edge[e].v=v; edge[e].c=c;  edge[e].next=head[u]; head[u]=e++;  
   edge[e].u=v; edge[e].v=u; edge[e].c=0;  edge[e].next=head[v]; head[v]=e++;  
} 
 
int sap()
{
   int flow=0,aug=inf,u;
   bool flag;
   for(int i=0;i<=2*n+1;i++)
   {
      cur[i]=head[i];
      gap[i]=dis[i]=0;
   }
   gap[0]=2*n+2;
   u=pre[start]=start;
  while(dis[start]<2*n+2)
  {
      flag=0;
      for(int &j=cur[u];j!=-1;j=edge[j].next)
      {
          int v=edge[j].v;
          if(edge[j].c>0&&dis[u]==dis[v]+1)
          {
               flag=1;
               if(edge[j].c<aug) aug=edge[j].c;
               pre[v]=u;
               u=v;
               if(u==end)
               {
                  flow+=aug;
                  while(u!=start)
                  {
                     u=pre[u];
                     edge[cur[u]].c-=aug;
                     edge[cur[u]^1].c+=aug;
                  }
                  aug=inf;
               }
                break;
           }

     }
     if(flag)  continue;
     int mindis=2*n+2;
     for(int j=head[u];j!=-1;j=edge[j].next)
     {
          int v=edge[j].v;
          if(edge[j].c>0&&dis[v]<mindis)
          {
              mindis=dis[v];
              cur[u]=j;
          }
     }
    if((--gap[dis[u]])==0)  break;
    gap[dis[u]=mindis+1]++;
    u=pre[u];
   }
 return flow;
}

int main()
{
    scanf("%d",&t);
    while(t--)
    {
       scanf("%d%d%d%d",&n,&m,&s,&h);
       e=0; start=0; end=2*n+1;
       memset(head,-1,sizeof(head));
       for(int i=1;i<=n;i++)
       {
          scanf("%d",&num);
          if(i!=s&&i!=h) add_adge(i,i+n,num);
       }
       add_adge(0,s,inf);
       add_adge(s,s+n,inf);
       add_adge(h,h+n,inf);
       add_adge(h+n,2*n+1,inf);
       for(int i=1;i<=m;i++)
       {
          scanf("%d%d",&a,&b);
          add_adge(a+n,b,inf);
          add_adge(b+n,a,inf);
       }
       printf("%d\n",sap());
    }
    return 0;
}