BZOJ 1066: [SCOI2007]蜥蜴( 最大流 )

结点容量..拆点然后随便写

---------------------------------------------------------------

#include<cstdio>#include<cstring>#include<algorithm> using namespace std; #define chk(x, y) (x >= 0 && x < R && y >= 0 && y < C) const int maxn = 10000;const int maxm = 2900; int h[maxn], cnt[maxn];int X[maxm][maxm], Y[maxm][maxm];int R, C, D, V, S, T, N;char s[maxm]; struct edge { int to, cap; edge *next, *rev;} E[1000000], *pt = E, *head[maxn], *p[maxn], *cur[maxn]; inline void Add(int u, int v, int w) { pt->to = v; pt->cap = w; pt->next = head[u]; head[u] = pt++;} inline void AddEdge(int u, int v, int w) { Add(u, v, w); Add(v, u, 0); head[u]->rev = head[v]; head[v]->rev = head[u];} void Init() { V = N = 0; scanf("%d%d%d", &R, &C, &D); memset(X, -1, sizeof X); for(int i = 0; i < R; i++) { scanf("%s", s); for(int j = 0; j < C; j++) if(s[j] > '0') { X[i][j] = V++; Y[i][j] = V++; AddEdge(X[i][j], Y[i][j], s[j] - '0'); } } S = V++; T = V++; for(int i = 0; i < R; i++) for(int j = 0; j < C; j++) if(~X[i][j]) for(int d = 1; d <= D; d++) for(int k = 0; k <= d; k++) { int x = i + k, y = j + d - k; if(chk(x, y) && ~X[x][y]) AddEdge(Y[i][j], X[x][y], maxn); x = i - k, y = j - d + k; if(chk(x, y) && ~X[x][y]) AddEdge(Y[i][j], X[x][y], maxn); x = i + k, y = j - d + k; if(chk(x, y) && ~X[x][y]) AddEdge(Y[i][j], X[x][y], maxn); x = i - k, y = j + d - k; if(chk(x, y) && ~X[x][y]) AddEdge(Y[i][j], X[x][y], maxn); } for(int i = 0; i < R; i++) { scanf("%s", s); for(int j = 0; j < C; j++) { if(~X[i][j] && (i < D || j < D || i + D >= R || j + D >= C)) AddEdge(Y[i][j], T, maxn); if(s[j] != '.') AddEdge(S, X[i][j], 1), N++; } } } void Solve() { memset(cnt, 0, sizeof cnt); memset(h, 0, sizeof h); cnt[0] = V; for(int i = 0; i < V; i++) cur[i] = head[i]; int Flow = 0; edge* e; for(int x = S, A = maxn; h[S] < V; ) { for(e = cur[x]; e; e = e->next) if(h[e->to] + 1 == h[x] && e->cap) break; if(e) { p[e->to] = cur[x] = e; A = min(A, e->cap); if((x = e->to) == T) { for(; x != S; x = p[x]->rev->to) { p[x]->cap -= A; p[x]->rev->cap += A; } Flow += A; A = maxn; } } else { if(!--cnt[h[x]]) break; h[x] = V; for(e = head[x]; e; e = e->next) if(h[e->to] + 1 < h[x] && e->cap) { h[x] = h[e->to] + 1; cur[x] = e; } cnt[h[x]]++; if(x != S) x = p[x]->rev->to; } } printf("%d\n", N - Flow); } int main() { Init(); Solve(); return 0;}---------------------------------------------------------------

1066: [SCOI2007]蜥蜴

Time Limit: 1 Sec Memory Limit: 162 MB
Submit: 2516 Solved: 1240
[Submit][Status][Discuss]

Description

在一个r行c列的网格地图中有一些高度不同的石柱，一些石柱上站着一些蜥蜴，你的任务是让尽量多的蜥蜴逃到边界外。 每行每列中相邻石柱的距离为1，蜥蜴的跳跃距离是d，即蜥蜴可以跳到平面距离不超过d的任何一个石柱上。石柱都不稳定，每次当蜥蜴跳跃时，所离开的石柱高度减1（如果仍然落在地图内部，则到达的石柱高度不变），如果该石柱原来高度为1，则蜥蜴离开后消失。以后其他蜥蜴不能落脚。任何时刻不能有两只蜥蜴在同一个石柱上。

Input

输入第一行为三个整数r，c，d，即地图的规模与最大跳跃距离。以下r行为石竹的初始状态，0表示没有石柱，1~3表示石柱的初始高度。以下r行为蜥蜴位置，“L”表示蜥蜴，“.”表示没有蜥蜴。

Output

输出仅一行，包含一个整数，即无法逃离的蜥蜴总数的最小值。

Sample Input

5 8 2
00000000
02000000
00321100
02000000
00000000
........
........
..LLLL..
........
........

Sample Output

1

HINT

100%的数据满足：1<=r, c<=20, 1<=d<=4

Source

Pku 2711 Leapin' Lizards