bzoj 1066 : [SCOI2007]蜥蜴 网络流

题目链接

给一个n*m的图， 里面每一个点代表一个石柱， 石柱有一个高度。 初始时有些石柱上面有蜥蜴， 蜥蜴可以跳到距离他曼哈顿距离小于等于d的任意一个石柱上，跳完后， 他原来所在的石柱高度会减一， 如果高度变为0， 那么石柱消失， 无法在跳到这个位置上， 跳到的那个石柱高度不会发生改变， 同一时刻一个石柱无法站两个蜥蜴。问有多少蜥蜴无法跳出边界。

很裸的网络流， 如果一个石柱距离边界距离小于d， 那么向汇点连一条权值为inf的边， 如果一个石柱初始有蜥蜴， 那么源点向这个点连一条1的边， 每个点拆成两个点， u向u'连一条边， 权值为高度。 一个石柱向所有距离他曼哈顿距离小于d的石柱连边， 权值inf。 跑一遍网络流就可以...

#include<bits/stdc++.h>
using namespace std;
#define pb(x) push_back(x)
#define ll long long
#define mk(x, y) make_pair(x, y)
#define lson l, m, rt<<1
#define mem(a) memset(a, 0, sizeof(a))
#define rson m+1, r, rt<<1|1
#define mem1(a) memset(a, -1, sizeof(a))
#define mem2(a) memset(a, 0x3f, sizeof(a))
#define rep(i, a, n) for(int i = a; i<n; i++)
#define ull unsigned long long
typedef pair<int, int> pll;
const double PI = acos(-1.0);
const double eps = 1e-8;
const int mod = 1e9+7;
const int inf = 1061109567;
const int dir[][2] = { {-1, 0}, {1, 0}, {0, -1}, {0, 1} };
int a[25][25];
const int maxn = 2e5+5;
int q[maxn*2], head[maxn*2], dis[maxn/10], s, t, num;
struct node
{
    int to, nextt, c;
    node(){}
    node(int to, int nextt, int c):to(to), nextt(nextt), c(c){}
}e[maxn*2];
void init() {
    num = 0;
    mem1(head);
}
void add(int u, int v, int c) {
    e[num] = node(v, head[u], c); head[u] = num++;
    e[num] = node(u, head[v], 0); head[v] = num++;
}
int bfs() {
    mem(dis);
    dis[s] = 1;
    int st = 0, ed = 0;
    q[ed++] = s;
    while(st<ed) {
        int u = q[st++];
        for(int i = head[u]; ~i; i = e[i].nextt) {
            int v = e[i].to;
            if(!dis[v]&&e[i].c) {
                dis[v] = dis[u]+1;
                if(v == t)
                    return 1;
                q[ed++] = v;
            }
        }
    }
    return 0;
}
int dfs(int u, int limit) {
    if(u == t) {
        return limit;
    }
    int cost = 0;
    for(int i = head[u]; ~i; i = e[i].nextt) {
        int v = e[i].to;
        if(e[i].c&&dis[v] == dis[u]+1) {
            int tmp = dfs(v, min(limit-cost, e[i].c));
            if(tmp>0) {
                e[i].c -= tmp;
                e[i^1].c += tmp;
                cost += tmp;
                if(cost == limit)
                    break;
            } else {
                dis[v] = -1;
            }
        }
    }
    return cost;
}
int dinic() {
    int ans = 0;
    while(bfs()) {
        ans += dfs(s, inf);
    }
    return ans;
}
int main()
{
    int n, m, d;
    cin>>n>>m>>d;
    int mn = m*n, x;
    init();
    s = 2*mn, t = s+1;
    char c[25];
    for(int i = 0; i<n; i++) {
        scanf("%s", c);
        for(int j = 0; j<m; j++) {
            a[i][j] = c[j]-'0';
            if(!a[i][j])
                continue;
            if(i<d||j<d||i+d>=n||j+d>=m) {
                add(i*m+j+mn, t, inf);
            }
            add(i*m+j, i*m+j+mn, a[i][j]);
        }
    }
    int ans = 0;
    for(int i = 0; i<n; i++) {
        scanf("%s", c);
        for(int j = 0; j<m; j++) {
            if(c[j] == 'L') {
                add(s, i*m+j, 1);
                ans++;
            }
        }
    }
    for(int i = 0; i<n; i++) {
        for(int j = 0; j<m; j++) {
            if(!a[i][j])
                continue;
            for(int x = max(0, i-d); x<=min(n-1, i+d); x++) {
                for(int y = max(0, j-d); y<=min(m-1, j+d); y++) {
                    if(x == i&& y==j )
                        continue;
                    if(!a[x][y])
                        continue;
                    if(abs(x-i)+abs(y-j)>d)
                        continue;
                    add(i*m+j+mn, x*m+y, inf);
                }
            }
        }
    }
    ans -= dinic();
    cout<<ans<<endl;
    return 0;
}