题目链接:BZOJ - 1066

题目分析

题目限制了高度为 x 的石柱最多可以有 x 只蜥蜴从上面跳起,那么就可以用网络流中的边的容量来限制。我们把每个石柱看作一个点,每个点拆成 i1, i2,从 i1 到 i2 连一条边,容量为这个石柱 i 的高度,即跳跃次数限制。来到这个石柱就是向 i1 连边,从这个石柱跳起就是从 i2 向外连边,这样只要从石柱 i 跳起就一定会消耗 i1 到 i2 的边的容量。如果 i 有蜥蜴,就从 S 到 i1 连一条容量为 1 的边,如果从石柱 i 能跳出边界,就从 i2 到 T 连 INF 的边。如果从石柱 i 能跳到石柱 j ,就从 i2 到 j1 连 INF 的边。这样求出的最大流就是能够跳出的蜥蜴的数量。

题目还限制了“任何时刻不能有两只蜥蜴在同一个石柱上”,但是想一想就会发现,这个限制其实并没有影响,只要我们按照合适的顺序让蜥蜴们先后跳跃,就一样可以让那些应该跳出的蜥蜴跳出。

另外,题目要求输出不能跳出的蜥蜴的数量,结果我直接输出了能跳出的数量...于是对着数据调了1h...最后发现是输出的东西不对...QAQ

代码

#include <iostream>
#include <cstdio>
#include <cstdlib>
#include <algorithm>
#include <cmath>
#include <cstring>
 
using namespace std;
 
const int MaxR = 20 + 5, MaxN = 400 * 2 + 5, MaxM = 500000 * 2 + 5, INF = 0x3fffffff;
 
int n, m, Dis, Top, S, T, Tot, Ans, Sum;
int Num[MaxN], d[MaxN];
 
bool Be[MaxR][MaxR];
 
char Str[MaxR];
 
struct Edge
{
    int v, w;
    Edge *Next, *Other;
} E[MaxM], *P = E, *Point[MaxN], *Last[MaxN];
 
inline void AddEdge(int x, int y, int z) {
    Edge *Q = ++P; ++P;
    P -> v = y; P -> w = z;
    P -> Next = Point[x]; Point[x] = P; P -> Other = Q;
    Q -> v = x; Q -> w = 0;
    Q -> Next = Point[y]; Point[y] = Q; Q -> Other = P;
}
 
struct ES
{
    int x, y, z;
    ES() {}
    ES(int a, int b, int c) {
        x = a; y = b; z = c;
    }
} JF[MaxN];
 
typedef double DB;
 
bool Check(int x, int y) {
    if ((x - 1) < Dis || (y - 1) < Dis) return true;
    if ((n - x) < Dis || (m - y) < Dis) return true;
    return false;
}
 
inline DB Sqr(DB x) {return x * x;}
 
DB Calc(DB x, DB y, DB xx, DB yy) {
    return sqrt(Sqr(x - xx) + Sqr(y - yy));
}
 
inline int gmin(int a, int b) {return a < b ? a : b;}
 
int DFS(int Now, int Flow) {
    if (Now == T) return Flow;
    int ret = 0;
    for (Edge *j = Last[Now]; j; j = j -> Next) {
        if (j -> w && d[Now] == d[j -> v] + 1) {
            Last[Now] = j;
            int p = DFS(j -> v, gmin(j -> w, Flow - ret));
            j -> w -= p; j -> Other -> w += p; ret += p;
            if (ret == Flow) return ret;
        }
    }
    if (d[S] >= Tot) return ret;
    if (--Num[d[Now]] == 0) d[S] = Tot;
    ++Num[++d[Now]];
    Last[Now] = Point[Now];
    return ret;
}
 
int main()
{
    scanf("%d%d%d", &n, &m, &Dis);
    Top = 0;
    for (int i = 1; i <= n; ++i) {
        scanf("%s", Str);
        for (int j = 1; j <= m; ++j) {
            if (Str[j - 1] != '0') {
                JF[++Top] = ES(i, j, Str[j - 1] - '0');
            }
        }  
    }
    Sum = 0;
    for (int i = 1; i <= n; ++i) {
        scanf("%s", Str);
        for (int j = 1; j <= m; ++j) {
            if (Str[j - 1] == 'L') {
                Be[i][j] = true;
                ++Sum;
            }
            else Be[i][j] = false;
        }
    }
    S = Top * 2 + 1; T = S + 1; Tot = T;
    for (int i = 1; i <= Top; ++i) {
        int xx = JF[i].x, yy = JF[i].y;
        AddEdge(i, Top + i, JF[i].z);
        if (Check(xx, yy)) AddEdge(Top + i, T, INF);
        if (Be[xx][yy]) AddEdge(S, i, 1);
        for (int j = 1; j <= Top; ++j) {
            if (j == i) continue;
            if (Calc(xx, yy, JF[j].x, JF[j].y) <= (DB)Dis)
                AddEdge(Top + i, j, INF);
        }
    }
    for (int i = 1; i <= Tot; ++i) Last[i] = Point[i];
    memset(Num, 0, sizeof(Num)); Num[0] = Tot;
    memset(d, 0, sizeof(d));
    Ans = 0;
    while (d[S] < Tot) Ans += DFS(S, INF);
    printf("%d\n", Sum - Ans);
    return 0;
}

  

[BZOJ 1066] [SCOI2007] 蜥蜴 【最大流】的更多相关文章

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

    结点容量..拆点然后随便写 --------------------------------------------------------------- #include<cstdio> ...

  2. poj 2711 Leapin' Lizards && BZOJ 1066: [SCOI2007]蜥蜴 最大流

    题目链接:http://poj.org/problem?id=2711 题目链接:http://www.lydsy.com/JudgeOnline/problem.php?id=1066 Your p ...

  3. POJ 2711 Leapin' Lizards / HDU 2732 Leapin' Lizards / BZOJ 1066 [SCOI2007]蜥蜴(网络流,最大流)

    POJ 2711 Leapin' Lizards / HDU 2732 Leapin' Lizards / BZOJ 1066 [SCOI2007]蜥蜴(网络流,最大流) Description Yo ...

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

    [题目链接] http://www.lydsy.com/JudgeOnline/problem.php?id=1066 [题目大意] 在一个r行c列的网格地图中有一些高度不同的石柱,一些石柱上站着一些 ...

  5. bzoj 1066: [SCOI2007] 蜥蜴

    这道题还是挺好想的,但我一开始还是想错了…… 把每个石柱拆成两个点,一个入度,一个出度,两个点连一条容量为高度的边,这样就可以限制从此石柱上经过的蜥蜴的数量.关于蜥蜴是否单独成点,我是单独当成了一个点 ...

  6. bzoj 1066 : [SCOI2007]蜥蜴 网络流

    题目链接 给一个n*m的图, 里面每一个点代表一个石柱, 石柱有一个高度. 初始时有些石柱上面有蜥蜴, 蜥蜴可以跳到距离他曼哈顿距离小于等于d的任意一个石柱上,跳完后, 他原来所在的石柱高度会减一, ...

  7. 1066: [SCOI2007]蜥蜴

    1066: [SCOI2007]蜥蜴 Time Limit: 1 Sec  Memory Limit: 162 MBSubmit: 3545  Solved: 1771[Submit][Status] ...

  8. [SCOI2007] 蜥蜴 (最大流)

    [SCOI2007] 蜥蜴 题目背景 07四川省选 题目描述 在一个r行c列的网格地图中有一些高度不同的石柱,一些石柱上站着一些蜥蜴,你的任务是让尽量多的蜥蜴逃到边界外. 每行每列中相邻石柱的距离为1 ...

  9. P2472 [SCOI2007]蜥蜴 (最大流)

    题目 P2472 [SCOI2007]蜥蜴 解析 这个题思路比较清晰,本(qi)来(shi)以(jiu)为(shi)无脑建图跑最大流,结果挂了,整了一个小时后重新建图才过的. 建立一个超级源点和一个超 ...

随机推荐

  1. Android color(颜色) 在XML文件和java代码中

    Android color(颜色) 在XML文件和java代码中,有需要的朋友可以参考下. 1.使用Color类的常量,如: int color = Color.BLUE;//创建一个蓝色 是使用An ...

  2. CentOS6.5解决中文乱码与设置字符集

    [ CleverCode发表在csdn博客中的原创作品,请勿转载,原创地址:http://blog.csdn.net/clevercode/article/details/46377577] 1)说明 ...

  3. Android 6.0 双卡拨号

    相关 api getCallCapablePhoneAccountsAdded in API level 23 Android 5.0 之前的版本 Call from second sim 获取 si ...

  4. eclipse4.3 kepler中安装maven

    1.软件准备 a:Eclipse 4.3 http://www.eclipse.org/downloads/ b:maven http://maven.apache.org/download.cgi ...

  5. MVC3.0 提交表单的方法

    在views中使用 @using (Html.BeginForm()) 来完成,其中放一个submit类型的提交按钮,如: @using (Html.BeginForm()) { @Html.Vali ...

  6. Oracle高版本导出dmp导入Oracle低版本报错:"不是有效的导出文件、头部验证失败"解决方法

    从Oracle高版本中导出dmp,然后导入到Oracle低版本时会报错:"不是有效的导出文件.头部验证失败",解决方法: 方法一:下载软件:AlxcTools,打开后选择要修改的文 ...

  7. C# 二叉查找树实现

    BuildTree 代码1次CODE完,没有BUG. 在画图地方debug了很多次.第一次画这种图. 一开始用treeview显示,但发现不是很好看出树结构,于是自己动手画了出来. using Sys ...

  8. asp.net对word文档进行修改 对于使用word文档做模板编辑比较适用

    最近做项目,需要多word文档进行编辑并导出一个新的word,在最初的word编辑中留下特定的字符串用来替换,然后在本地生成一个新的word文档,并且不修改服务器中的word文档,这样才能保证服务器中 ...

  9. 如何通过WiFi来进行Android的真机模拟

    我们知道,在使用模拟机模拟的时候会出现较多的问题,所以如果有一部Android手机的话进行真机模拟是极好的. 准备工作: 第一种方法:使用数据线,具体操作百度.略(非WIFI操作的真机模拟) 第二方法 ...

  10. windows 20003 扩展安装后不成功的原因

    windows扩展如果安装不成功(PHP扩展)很大的可能就是那个DLL的权限不够.需要分配:AdministratorAuthenticater UsersIIS_WPGSYSTEMUsers