[BZOJ1001][BeiJing2006]狼抓兔子
[BZOJ1001][BeiJing2006]狼抓兔子
试题描述
现在小朋友们最喜欢的"喜羊羊与灰太狼",话说灰太狼抓羊不到,但抓兔子还是比较在行的,
而且现在的兔子还比较笨,它们只有两个窝,现在你做为狼王,面对下面这样一个网格的地形:
左上角点为(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只狼,
才能完全封锁这条道路,你需要帮助狼王安排一个伏击方案,使得在将兔子一网打尽的前提下,参与的
狼的数量要最小。因为狼还要去找喜羊羊麻烦.
输入
第一行为N,M.表示网格的大小,N,M均小于等于1000.
接下来分三部分
第一部分共N行,每行M-1个数,表示横向道路的权值.
第二部分共N-1行,每行M个数,表示纵向道路的权值.
第三部分共N-1行,每行M-1个数,表示斜向道路的权值.
输入文件保证不超过10M
输出
输出一个整数,表示参与伏击的狼的最小数量.
输入示例
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
输出示例
14
数据规模及约定
见“输入”
题解
转化成对偶图跑最短路特别快。经典问题了。
这题强行出 n = 1 或 m = 1 的数据,真是 ****
#include <iostream> #include <cstdio> #include <cstdlib> #include <cstring> #include <cctype> #include <algorithm> #include <queue> using namespace std; int read() { int x = 0, f = 1; char c = getchar(); while(!isdigit(c)){ if(c == '-') f = -1; c = getchar(); } while(isdigit(c)){ x = x * 10 + c - '0'; c = getchar(); } return x * f; } #define maxn 1010 #define maxnode 2000010 #define maxm 6000010 #define oo 2147483647 int N, M, head[maxnode], nxt[maxm], to[maxm], dist[maxm]; void AddEdge(int a, int b, int c) { to[++M] = b; dist[M] = c; nxt[M] = head[a]; head[a] = M; swap(a, b); to[++M] = b; dist[M] = c; nxt[M] = head[a]; head[a] = M; return ; } int upt(int x, int y, int m) { return ((x - 1) * (m - 1) + y << 1) - 1; } int dnt(int x, int y, int m) { return upt(x, y, m) + 1; } struct Node { int u, d; Node() {} Node(int _, int __): u(_), d(__) {} bool operator < (const Node& t) const { return d > t.d; } }; priority_queue <Node> Q; int d[maxnode]; bool vis[maxnode]; void solve(int s, int t) { memset(vis, 0, sizeof(vis)); for(int i = 1; i <= N; i++) d[i] = oo; while(!Q.empty()) Q.pop(); d[s] = 0; Q.push(Node(s, 0)); while(!Q.empty()) { int u = Q.top().u; Q.pop(); if(vis[u]) continue; vis[u] = 1; if(u == t){ printf("%d\n", d[u]); break; } for(int e = head[u]; e; e = nxt[e]) if(d[to[e]] > d[u] + dist[e]) { d[to[e]] = d[u] + dist[e]; if(!vis[to[e]]) Q.push(Node(to[e], d[to[e]])); } } return ; } int main() { int n, m; scanf("%d%d", &n, &m); if(m == 1) { int mn = oo; for(int i = 1; i < n; i++) mn = min(mn, read()); return printf("%d\n", mn), 0; } if(n == 1) { int mn = oo; for(int i = 1; i < m; i++) mn = min(mn, read()); return printf("%d\n", mn), 0; } M = 0; memset(head, 0, sizeof(head)); N = ((n - 1) * (m - 1) << 1) + 2; int s = N - 1, t = N; for(int i = 1; i <= n; i++) for(int j = 1; j < m; j++) { int u, v; if(i == 1) { u = upt(i, j, m); AddEdge(s, u, read()); } if(1 < i && i < n) { u = upt(i, j, m); v = dnt(i - 1, j, m); AddEdge(u, v, read()); } if(i == n) { u = dnt(i - 1, j, m); AddEdge(u, t, read()); } } for(int i = 1; i < n; i++) for(int j = 1; j <= m; j++) { int u, v; if(j == 1) { u = dnt(i, j, m); AddEdge(u, t, read()); } if(1 < j && j < m) { u = upt(i, j - 1, m); v = dnt(i, j, m); AddEdge(u, v, read()); } if(j == m) { u = upt(i, j - 1, m); AddEdge(s, u, read()); } } for(int i = 1; i < n; i++) for(int j = 1; j < m; j++) { int u = upt(i, j, m), v = dnt(i, j, m); AddEdge(u, v, read()); } solve(s, t); return 0; }