|| 右下角 cnblogs 最小割 out cst con 分享 log

Description

現在小朋友們最喜歡的"喜羊羊與灰太狼",話說灰太狼抓羊不到，但抓兔子還是比較在行的， 而且現在的兔子還比較笨，它們只有兩個窩，現在你做為狼王，面對下面這樣一個網格的地形：

左上角點為(1,1),右下角點為(N,M)(上圖中N=4,M=5).有以下三種類型的道路 1:(x,y)<==>(x+1,y) 2:(x,y)<==>(x,y+1) 3:(x,y)<==>(x+1,y+1) 道路上的權值表示這條路上最多能夠通過的兔子數，道路是無向的. 左上角和右下角為兔子的兩個窩， 開始時所有的兔子都聚集在左上角(1,1)的窩裏，現在它們要跑到右下解(N,M)的窩中去，狼王開始伏擊 這些兔子.當然為了保險起見，如果一條道路上最多通過的兔子數為K，狼王需要安排同樣數量的K只狼， 才能完全封鎖這條道路，你需要幫助狼王安排一個伏擊方案，使得在將兔子一網打盡的前提下，參與的 狼的數量要最小。因為狼還要去找喜羊羊麻煩.

Input

第一行為N,M.表示網格的大小，N,M均小於等於1000. 接下來分三部分 第一部分共N行，每行M-1個數，表示橫向道路的權值. 第二部分共N-1行，每行M個數，表示縱向道路的權值. 第三部分共N-1行，每行M-1個數，表示斜向道路的權值. 輸入文件保證不超過10M

Output

輸出一個整數，表示參與伏擊的狼的最小數量.

Sample Input

3 4
5 6 4
4 3 1
7 5 3
5 6 7 8
8 7 6 5
5 5 5
6 6 6

Sample Output

14

HINT

2015.4.16新加數據一組，可能會卡掉從前可以過的程序。

平面圖求最小割，所以可以用spfa做。

 1 #include <iostream>
 2 #include <cstring>
 3 #include <cstdio>
 4 #include <queue>
 5 const int N = 1000 + 7, M = 2000000 + 7, inf = 1 << 30;
 6 using namespace std;
 7 int n,m,s,t,cnt,head[M],dis[M];
 8 bool used[M];
 9 struct edges
10 {
 
11     int nxt,to,wi;
12 } E[M * 3];
13  
14 int idid(int a,int b)
15 {
16     return (a - 1) * (m - 1) * 2 + b * 2;
17 }
18  
19 void add(int u,int v,int wi)
20 {
21     E[++cnt] = (edges){head[u],v,wi};head[u] = cnt;
22     E[++cnt] = (edges) {head[v],u,wi}; head[v] = cnt;
23 }
24  
25 void Init()
26 {
27     scanf("%d%d",&n,&m);
28     if(n == 1 || m == 1)
29     {
30         if(n == 1) swap(n,m); 
31         int ans =  inf, c;
32         for(int i = 1; i < n; ++i)
33         {
34             scanf("%d",&c);
35             ans = min(ans,c);
36         }
37         printf("%d\n",ans); return;
38     }
39     int c;s = 0, t = ((n-1) * (m-1) << 1) + 1;
40     for(int i = 1; i <= n; ++i)
41     for(int j = 1; j < m; ++j)
42     {
43         scanf("%d",&c); 
44         if(i == 1) add(s,j<<1,c);
45         else if(i == n) add(idid(i-1,j)-1,t,c);
46         else add(idid(i-1,j)-1,idid(i,j),c);
47     }
48  
49     for(int i = 1; i < n; ++i)
50     for(int j = 1; j <= m; ++j)
51     {
52         scanf("%d",&c);
53         if(j == 1) add(idid(i,j)-1,t,c);
54         else if(j == m) add(s,idid(i,j-1),c);
55         else add(idid(i,j-1),idid(i,j)-1,c);
56     }
57     for(int i = 1;i < n; ++i)
58     for(int j = 1 ; j < m; ++j)
59     {
60         scanf("%d",&c);
61         add(idid(i,j)-1,idid(i,j),c);
62     }
63 }
64  
65 void dij()
66 {
67     queue<int> Q;
68     memset(dis,127/2,sizeof(dis));
69     memset(used,false,sizeof(used)); used[s] = true;
70     Q.push(s); dis[s] = 0;
71     while(!Q.empty())
72     {
73         int u = Q.front(); Q.pop(); used[u] = false;
74         for(int i = head[u]; i; i = E[i].nxt)
75         {
76             int v = E[i].to;
77             if(dis[v] <= dis[u] + E[i].wi) continue;
78             dis[v] = dis[u] + E[i].wi;
79             if(used[v]) continue;
80             Q.push(v); used[v] = true;
81         }
82     }
83 }
84  
85 void Solve()
86 {
87     dij();
88     printf("%d\n",dis[t]);
89 }
90  
91 int main()
92 {
93 //  freopen("1.in","r",stdin);
94     Init();
95     if(n == 1 || m == 1) return 0;
96     Solve();
97     return 0;
98 }
99 

bzoj 1001 狼抓兔子