1. 程式人生 > >1001: [BeiJing2006]狼抓兔子

1001: [BeiJing2006]狼抓兔子

online emp etc front algo bsp 兩個 pty ons

1001: [BeiJing2006]狼抓兔子

Time Limit: 15 Sec Memory Limit: 162 MB
Submit: 23595 Solved: 5940

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新加數據一組,可能會卡掉從前可以過的程序。

code

 1 #include<cstdio>
 2 #include<queue>
 3 #include<algorithm>
 4 #include<cstring>
 5 
 6 using namespace std;
 7 const int MAXN = 2000100;
 8 struct Edge{
 9     int to,w,nxt;
10     Edge(){}
11     Edge(int a,int b,int c){to = a,w = b,nxt = c;}
12 }e[MAXN<<2
]; 13 int head[MAXN],dis[MAXN]; 14 bool vis[MAXN]; 15 queue<int>q; 16 int n,m,tot; 17 18 int read() 19 { 20 int x = 0,f = 1;char ch = getchar(); 21 while (ch<0||ch>9) {if (ch==-)f=-1; ch = getchar();} 22 while (ch>=0&&ch<=9) x = x*10+ch-0, ch = getchar(); 23 return x; 24 } 25 void add_edge(int u,int v,int w) 26 { 27 e[++tot] = Edge(v,w,head[u]); 28 head[u] = tot; 29 e[++tot] = Edge(u,w,head[v]); 30 head[v] = tot; 31 } 32 void spfa() 33 { 34 memset(dis,0x3f,sizeof(dis)); 35 q.push(0); 36 vis[0] = true; 37 dis[0] = 0; 38 while (!q.empty()) 39 { 40 int u = q.front(); 41 q.pop(); 42 for (int i=head[u]; i; i=e[i].nxt) 43 { 44 int v = e[i].to,w = e[i].w; 45 if (dis[v]>dis[u]+w) 46 { 47 dis[v] = dis[u]+w; 48 if (!vis[v]) 49 { 50 q.push(v); 51 vis[v] = true; 52 } 53 } 54 } 55 vis[u] = false; 56 } 57 } 58 int main() 59 { 60 n = read(), m = read(); 61 int x,y,z,tmp = (2*(n-1)*(m-1))+1; 62 for (int i=1; i<=n; ++i) 63 for (int j=1; j<m; ++j) 64 { 65 z = read(); 66 x = i==1? 0:(2*(i-1)-1)*(m-1)+j; 67 y = i==n? tmp:(2*(i-1))*(m-1)+j; 68 add_edge(x,y,z); 69 } 70 for (int i=1; i<n; ++i) 71 for (int j=1; j<=m; ++j) 72 { 73 z = read(); 74 x = j==1? tmp:(2*(i-1))*(m-1)+j-1; 75 y = j==m? 0:(2*(i-1))*(m-1)+j-1+m; 76 add_edge(x,y,z); 77 } 78 for (int i=1; i<n; ++i) 79 for (int j=1; j<m; ++j) 80 { 81 z = read(); 82 x = (2*(i-1))*(m-1)+j; 83 y = (2*(i-1)+1)*(m-1)+j; 84 add_edge(x,y,z); 85 } 86 spfa(); 87 printf("%d",dis[tmp]); 88 return 0; 89 }

1001: [BeiJing2006]狼抓兔子