1339 / 1163: [Baltic2008]Mafia

1163: [Baltic2008]Mafia

Time Limit: 10 Sec  Memory Limit: 162 MB
Submit: 96  Solved: 60
[Submit][Status][Discuss]

Description

匪徒准备从一个车站转移毒品到另一个车站,警方准备进行布控. 对于每个车站进行布控都需要一定的代价,现在警方希望使用最小的代价控制一些车站,使得去掉这些车站后，匪徒无法从原定的初始点到达目标点

Input

第一行输入N,M代表车站的总个数,及有多少条双向边连接它们. 2<=n<=200 , 1 <=m<=20000. 第二行给出两个数a,b,代表匪徒的出发点及目标点.1<=a,b<=N,a<>b. 再下来有N行,给出对第i个车站进行布控所需要的Money,其不超过10 000 000 再下来M行,用于描述图的结构.

Output

最少需要多少Money

Sample Input

5 6
5 3
2
4
8
3
10
1 5
1 2
2 4
4 5
2 3
3 4

Sample Output

5

HINT

 

Source

题解：一个呵呵呵的最小割，将每个题目中的车站拆分成两个点，然后连接上一条边权为该站点费用的边，然后别的题目中所述的边直接连双向边，注意别弄乱了就是了

 /**************************************************************
     Problem:
     User: HansBug
     Language: Pascal
     Result: Accepted
     Time: ms
     Memory: kb
 ****************************************************************/

 type
     point=^node;
     node=record
                g,w:longint;
                next,anti:point;
     end;
 var
    i,j,k,l,m,n,s,t,ans:longint;
    a:array[..] of point;
    d,dv:array[..] of longint;
 function min(x,y:longint):longint;
          begin
               if x<y then min:=x else min:=y;
          end;
 procedure add(x,y,z:longint);
           var p:point;
           begin
                new(p);p^.g:=y;p^.w:=z;p^.next:=a[x];a[x]:=p;
                new(p);p^.g:=x;p^.w:=;p^.next:=a[y];a[y]:=p;
                a[x]^.anti:=a[y];a[y]^.anti:=a[x];
           end;
 function dfs(x,flow:longint):longint;
          var p:point;k:longint;
          begin
               if x=t then exit(flow);
               dfs:=;p:=a[x];
               while p<>nil do
                     begin
                          if (p^.w<>) and (d[x]=(d[p^.g]+)) then
                             begin
                                  k:=dfs(p^.g,min(flow-dfs,p^.w));
                                  if p^.w<>maxlongint then dec(p^.w,k);
                                  if p^.anti^.w<>maxlongint then inc(p^.anti^.w,k);
                                  inc(dfs,k);if dfs=flow then exit;
                             end;
                          p:=p^.next;
                     end;
               if d[s]=n then exit;
               dec(dv[d[x]]);
               if dv[d[x]]= then d[s]:=n;
               inc(d[x]);inc(dv[d[x]]);
          end;
 begin
      readln(n,m);
      for i:= to n* do a[i]:=nil;
      readln(s,t);s:=s*-;t:=t*;
      for i:= to n do
          begin
               readln(j);
               add(i*-,i*,j);
          end;
      for i:= to m do
          begin
               readln(j,k);
               add(j*,k*-,maxlongint);
               add(k*,j*-,maxlongint);
          end;
      fillchar(dv,sizeof(dv),);
      fillchar(d,sizeof(d),);
      ans:=;n:=n*;
      while d[s]<n do inc(ans,dfs(s,maxlongint));
      writeln(ans);
      readln;
 end.