题意:有n个城市,每个城市有一定数量的警察,有一群小偷,从城市S,到T,问最少需要多少警察可以使小偷到不了T城市。
将每个城市的警察数量看成流量,那么问题就转化成求S - T的最小割。
将每个点拆成i , i + n ,流量是该点的值.然后跑一次最大流就可以了。
最大流最小割定理:任意一个流网络的最大流量等于该网络的最小的割的容量。
#include <iostream>
#include <cstdio>
#include <algorithm>
#include <string>
#include <cmath>
#include <cstring>
#include <queue>
#include <set>
#include <vector>
#include <stack>
#include <map>
#include <iomanip>
#define PI acos(-1.0)
#define Max 2005
#define inf 1<<28
#define LL(x) (x<<1)
#define RR(x) (x<<1|1)
#define REP(i,s,t) for(int i=(s);i<=(t);++i)
#define ll long long
#define mem(a,b) memset(a,b,sizeof(a))
#define mp(a,b) make_pair(a,b)
using namespace std;
inline void readint(int &ret)
{
char c;
do
{
c = getchar();
}
while(c < '0' || c > '9');
ret = c - '0';
while((c=getchar()) >= '0' && c <= '9')
ret = ret * 10 + ( c - '0' );
}
int city[111] ;
struct kdq
{
int s , e, l ,next ;
} ed[1111111] ;
int head[1111] ,num ;
void add(int s ,int e ,int l )
{
ed[num].s = s ;
ed[num].e = e ;
ed[num].l = l ;
ed[num].next = head[s] ;
head[s] = num ++ ;
ed[num].s = e ;
ed[num].e = s ;
ed[num].l = 0 ;
ed[num].next = head[e] ;
head[e] = num ++ ;
}
int deep[1111] ;
int qe[11111111] ;
void init()
{
mem(head,-1) ;
num = 0 ;
}
int S , T ;
int dinic_bfs()
{
mem(deep,-1) ;
deep[S] = 0 ;
int h = 0 , t = 0 ;
qe[h ++ ] = S ;
while( h > t )
{
int tt = qe[t ++ ] ;
for (int i = head[tt] ; ~i ; i = ed[i].next )
{
int e = ed[i].e ;
int l = ed[i].l ;
if(l > 0&& deep[e] == -1)
{
deep[e] = deep[tt] + 1 ;
qe[ h ++ ] = e ;
}
}
}
return deep[T] != -1 ;
}
int dinic_dfs(int now, int f)
{
if(now == T)return f ;
int flow = 0 ;
for (int i = head[now] ; ~i ; i = ed[i].next )
{
int e = ed[i].e ;
int l = ed[i].l ;
if(deep[e] == deep[now] + 1 && l > 0 && (f - flow ) > 0 )
{
int mm = min(l,f - flow) ;
int nn = dinic_dfs(e,mm) ;
flow += nn ;
ed[i].l -= nn ;
ed[i ^ 1].l += nn ;
}
}
if(!flow)deep[now] = -2 ;
return flow ;
}
int dinic()
{
int flow = 0 ;
while(dinic_bfs())
{
flow += dinic_dfs(S,inf) ;
}
return flow ;
}
int main()
{
int t ;
readint(t) ;
while( t -- )
{
int n , m , s, h ;
readint(n) ;
readint(m) ;
readint(s) ;
readint(h) ;
init() ;
S = 0 ,T = 2 * n + 1 ;
REP(i,1,n)
{
int a ;
readint(a) ;
if(i != s && i != h)
add(i , i + n , a) ;
else
add(i ,i + n ,inf) ;//起点出去的流量为inf,终点进来的流量也是inf
}
add(S,s,inf) ;
add(h + n ,T,inf) ;
REP(i,1,m)
{
int a , b ;
readint(a) ;
readint(b) ;
add(a + n ,b ,inf) ;//双向边
add(b + n , a ,inf) ;
}
printf("%d\n",dinic()) ;
}
return 0;
}