蜥蜴_bzoj1066_最大流

31 篇文章 0 订阅

Description


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

Input


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

Output


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

Analysis


  1. 网络流
  2. 每根柱子有限制所以拆点,连边是它的高度
  3. 每根柱子只有一只蜥蜴所以源点和蜥蜴连边是1
  4. 柱子随便走所以能相互到达的柱子连边INF
  5. 所有能一步跳出的柱子连汇点
  6. 事后写出了dinic+边集(尴尬脸

Code


#include <cstdio>
#include <cstring>
#include <queue>
#define rep(i, a, b) for (int i = a; i <= b; i++)
#define fill(x, t) memset(x, t, sizeof(x))
#define N 1001
#define E N * N * 8 + 1
#define INF 0x3f3f3f3f
using namespace std;
struct edge{int y, w, next;}e[E];
int map[N][N], cur[N], dis[N], ls[N], edgecnt;
int dist(int x1, int y1, int x2, int y2){
    return (x1 - x2) * (x1 - x2) + (y1 - y2) * (y1 - y2);
}
int min(int x, int y){
    return x<y?x:y;
}
int getPoint(int m, int x, int y){
    return (x - 1) * m + y;
}
int add(int x, int y, int w){
    e[++edgecnt] = (edge){y, w, ls[x]}; ls[x] = edgecnt;
    e[++edgecnt] = (edge){x, 0, ls[y]}; ls[y] = edgecnt;
    //printf("x %d->y %d, v:%d\n", x, y, w);
    return 0;
}
int bfs(int st, int ed){
    queue<int>q;
    fill(dis, 0);
    q.push(st);
    dis[st] = 1;
    while (!q.empty()){
        int now = q.front(); q.pop();
        int tmp = (now + 1) / 2;
        // printf("%d\n",tmp);
        for (int i = ls[now]; i; i = e[i].next){
            if (e[i].w > 0 && !dis[e[i].y]){
                dis[e[i].y] = dis[now] + 1;
                q.push(e[i].y);
                if (e[i].y == ed){
                    return 1;
                }
            }
        }
    }
    return 0;
}
int find(int now, int ed, int mn){
    if (now == ed || !mn){
        return mn;
    }
    int ret = 0;
    for (int &i = cur[now]; i; i = e[i].next){
        if (e[i].w > 0 && dis[now] + 1 == dis[e[i].y]){
            int d = find(e[i].y, ed, min(mn - ret, e[i].w));
            e[i].w -= d;
            e[i ^ 1].w += d;
            ret += d;
            if (ret == mn){
                break;
            }
        }
    }
    return ret;
}
int dinic(int n, int m, int st, int ed){
    int mxFlow = 0;
    while (bfs(st, ed)){
        rep(i, 0, ed){
            cur[i] = ls[i];
        }
        mxFlow += find(st, ed, INF);
    }
    return mxFlow;
}
int main(void){
    int n, m, d;
    scanf("%d%d%d",&n,&m,&d);
    edgecnt = 1;
    rep(i, 1, n){
        char s[31];
        scanf("%s", s + 1);
        rep(j, 1, m){
            map[i][j] = s[j] - '0';
            if (map[i][j]){
                add(getPoint(m, i, j) * 2 - 1, getPoint(m, i, j) * 2, map[i][j]);
            }
        }
    }
    int st = getPoint(m, n, m) * 2 + 1, ed = getPoint(m, n, m) * 2 + 2;
    int tot = 0;
    rep(i, 1, n){
        char s[31];
        scanf("%s", s + 1);
        rep(j, 1, m){
            if (s[j] == 'L'){
                add(st, getPoint(m, i, j) * 2 - 1, 1);
                tot ++;
            }
        }
    }
    rep(i, 1, n){
        rep(j, 1, m){
            rep(k, 1, n){
                rep(l, 1, m){
                    if (map[i][j] && map[k][l] && (i != k || j != l) && dist(i, j, k, l) <= d * d){
                        add(getPoint(m, i, j) * 2, getPoint(m, k, l) * 2 - 1, INF);
                    }
                }
            }
        }
    }
    rep(i, 1, n){
        rep(j, 1, m){
            if (map[i][j] && (i <= d || n - i + 1 <=d || j <= d || m - j + 1 <= d)){
                add(getPoint(m, i, j) * 2, ed, INF);
            }
        }
    }
    int ans = dinic(n, m, st, ed);
    printf("%d\n",tot - ans);
    return 0;
}
评论
添加红包

请填写红包祝福语或标题

红包个数最小为10个

红包金额最低5元

当前余额3.43前往充值 >
需支付:10.00
成就一亿技术人!
领取后你会自动成为博主和红包主的粉丝 规则
hope_wisdom
发出的红包
实付
使用余额支付
点击重新获取
扫码支付
钱包余额 0

抵扣说明:

1.余额是钱包充值的虚拟货币,按照1:1的比例进行支付金额的抵扣。
2.余额无法直接购买下载,可以购买VIP、付费专栏及课程。

余额充值