题目描述

网址:https://daniu.luogu.org/problemnew/show/3227
大意:
平面上有一长方体,目标为将其切割为上下两半。
切割点为\((x,y,z)\)的点,每个点有一个不和谐值v,现在有两个要求:

  • [1]任何两个平面相邻的切割点之间的高度差不能超过D
  • [2]要使得最终的不和谐值最小。
    求解最小的$ ∑vi $

题目解法

直接建图,化点为边。将最下面的点与S相连,最上面的与T相连。
然后跑最小割即为答案。
现在关键在于如何限制高度差不超过D。
以\((x,y,z)\)与\((x+1,y,z)\)为例。
建立\((x,y,z) —>(x+1,y,z-D)\) 与 \((x+1,y,z+D+1)—>(x,y,z+1)\)两条边
那么在割断\((x,y,z)->(x,y,z+1)\)时,就可保证\((x+1,y)\)上割的为题目所限制范围。
动手画一下就明白了。 如果割超出范围的边,S、T依旧联通。
所以用上述方法建图,然后跑最小割即可得到答案。

实现代码

#include<cstdio>
#include<cstdlib>
#include<iostream>
#include<cstring>
#include<cmath>
#include<queue>
#include<algorithm>
#define ll long long
#define RG register
#define IL inline
#define maxn 125000
#define INF 1e16+7
using namespace std;

IL ll gi(){
    RG ll date = 0, m = 1;  RG char ch = 0;
    while(ch!='-'&&(ch<'0'||ch>'9'))ch = getchar();
    if(ch == '-'){m = -1; ch = getchar();}
    while(ch>='0' && ch<='9')
       {date=date*10+ch-'0'; ch = getchar();}
    return date*m;
}

ll cnt,Q,P,R,S,T,D;
ll vis[maxn],dep[maxn],head[maxn],cur[maxn];
ll mx[4] = {0,0,-1,1},my[4] = {-1,1,0,0};
struct Road { ll to,next,cap,flow; }t[20*maxn];

IL bool Bfs(){
    for(RG ll i = 0; i <= T; i ++)vis[i] = 0;
    dep[S] = 0; vis[S] = 1; queue<ll>Que; Que.push(S);
    while(!Que.empty()){
        RG ll u = Que.front(); Que.pop();
        for(RG ll i = head[u]; i!= -1; i = t[i].next){
            RG ll v = t[i].to;
            if(!vis[v] && t[i].cap > t[i].flow){
                dep[v] = dep[u] + 1;
                vis[v] = true; Que.push(v);
            }
        }
    }return vis[T];
}

int Dfs(RG ll u,RG ll Lim){
    if(u == T || Lim == 0)return Lim;
    RG ll flow =0,f = 0;
    for(RG ll &i = cur[u]; i != -1; i = t[i].next){
        RG ll v = t[i].to;
        if(dep[v] == dep[u] + 1 && (f = Dfs(v,min(Lim,t[i].cap-t[i].flow))) > 0){
            flow += f; Lim -= f;
            t[i].flow += f; t[(i^1)].flow -= f;
            if(Lim == 0)break;
        }
    }return flow;
}

IL int Dinic(){
    RG ll Flow = 0;
    while(Bfs()){
        for(RG ll i = 0; i <= T; i ++)cur[i] = head[i];
        Flow += Dfs(S,INF);
    }return Flow;
}

IL ll Cod(RG ll z,RG ll x,RG ll y){
    return z*(P*Q) + (x-1)*Q + y;
}

IL void Add(ll uu,ll vv,ll cc){
    t[cnt++] = (Road){vv,head[uu],cc,0}; head[uu] = cnt - 1;
    t[cnt++] = (Road){uu,head[vv],0,0};  head[vv] = cnt - 1;
}

int main(){
    freopen("testdate.in","r",stdin);
    P = gi(); Q = gi(); R = gi(); D = gi();
    S = 0; T = R*P*Q+P*Q+1;
    for(RG ll i = 0; i <= T; i ++)head[i] = -1;
    for(RG ll h = 1; h <= R; h ++)
        for(RG ll i = 1; i <= P; i ++)
            for(RG ll j = 1; j <= Q; j ++)
                Add( Cod(h,i,j) , Cod(h-1,i,j) , gi());
    for(RG ll i = 1; i <= P; i ++)
        for(RG ll j = 1; j <= Q; j ++)
            Add(Cod(0,i,j),T,INF);
    for(RG ll i = 1; i <= P; i ++)
        for(RG ll j = 1; j <= Q; j ++)
            Add(S,Cod(R,i,j),INF);
    for(RG ll h = 1; h <= R; h ++)
        for(RG ll i = 1; i <= P; i ++)
            for(RG ll j = 1; j <= Q; j ++)
                for(RG ll f = 0; f < 4; f ++){
                    RG int x = mx[f]+i,y = my[f]+j;
                    if(x<1||y<1||x>P||y>Q)continue;
                    if(h-D>=0)Add(Cod(h-D,x,y),Cod(h,i,j),INF);
                    if(h+D<=R)Add(Cod(h,i,j),Cod(h+D,x,y),INF);
                }

    RG ll Ans = Dinic(); printf("%lld",Ans);
    return 0;
}

[HNOI2013]切糕的更多相关文章

  1. BZOJ 3144: [Hnoi2013]切糕

    3144: [Hnoi2013]切糕 Time Limit: 10 Sec  Memory Limit: 128 MBSubmit: 1495  Solved: 819[Submit][Status] ...

  2. bzoj 3144: [Hnoi2013]切糕 最小割

    3144: [Hnoi2013]切糕 Time Limit: 10 Sec  Memory Limit: 128 MBSubmit: 681  Solved: 375[Submit][Status] ...

  3. BZOJ_3144_[Hnoi2013]切糕_最小割

    BZOJ_3144_[Hnoi2013]切糕_最小割 Description Input 第一行是三个正整数P,Q,R,表示切糕的长P. 宽Q.高R.第二行有一个非负整数D,表示光滑性要求.接下来是R ...

  4. bzoj千题计划142:bzoj3144: [Hnoi2013]切糕

    http://www.lydsy.com/JudgeOnline/problem.php?id=3144 如果D=2 ,两个点,高度为4,建图如下 #include<queue> #inc ...

  5. 【BZOJ3144】[HNOI2013]切糕

    [BZOJ3144][HNOI2013]切糕 题面 题目描述 经过千辛万苦小 A 得到了一块切糕,切糕的形状是长方体,小 A 打算拦腰将切糕切成两半分给小 B.出于美观考虑,小 A 希望切面能尽量光滑 ...

  6. 【BZOJ 3144】 3144: [Hnoi2013]切糕 (最小割模型)

    3144: [Hnoi2013]切糕 Time Limit: 10 Sec  Memory Limit: 128 MBSubmit: 1764  Solved: 965 Description Inp ...

  7. BZOJ3144 Hnoi2013 切糕 【网络流】*

    BZOJ3144 Hnoi2013 切糕 Description Input 第一行是三个正整数P,Q,R,表示切糕的长P. 宽Q.高R.第二行有一个非负整数D,表示光滑性要求.接下来是R个P行Q列的 ...

  8. 【BZOJ3144】[Hnoi2013]切糕 最小割

    [BZOJ3144][Hnoi2013]切糕 Description Input 第一行是三个正整数P,Q,R,表示切糕的长P. 宽Q.高R.第二行有一个非负整数D,表示光滑性要求.接下来是R个P行Q ...

  9. 3144: [Hnoi2013]切糕

    3144: [Hnoi2013]切糕 Time Limit: 10 Sec  Memory Limit: 128 MBSubmit: 1526  Solved: 827[Submit][Status] ...

  10. bzoj3144 [HNOI2013]切糕(最小割)

    bzoj3144 [HNOI2013]切糕(最小割) bzoj Luogu 题面描述见上 题解时间 一开始我真就把这玩意所说的切面当成了平面来做的 事实上只是说相邻的切点高度差都不超过 $ d $ 对 ...

随机推荐

  1. java程序设计原则

    前言: 前言:java这种面向对象的的语言其实是很符合人的生活思维模式的,因为它以对象作为自己的研究对象,所谓"万物皆对象".一个人.一颗苹果.一只蚂蚁都是对象.所以它的设计原则和 ...

  2. XAMPP简介、安转、使用

    虽然没有写Mac安装方法及使用, 但方法也都大相径庭, 殊途同归而已. XAMPP简介 XAMPP是一款开源.免费的网络服务器软件,经过简单安装后,就可以在个人电脑上搭建服务器环境.本文为大家介绍Wi ...

  3. php分布式redis实现session共享

    方法一:找到配置文件php.ini,修改为下面内容,保存并重启服务 session.save_handler = redis session.save_path = "tcp://127.0 ...

  4. 启动MySql提示:The server quit without updating PID file(…)失败

    1.可能是/usr/local/mysql/data/rekfan.pid文件没有写的权限解决方法 :给予权限,执行 "chown -R mysql:mysql /var/data" ...

  5. MysqL错误之_ERROR! MySQL server PID file could not be found!

    在配置Mysql主从GTID模式下,启动Mysql服务时出现报错,搜索了一番,找到了一个简单可靠的方法,直接成功.如果遇到相同问题没有解决的童鞋,那就去试一下很多其他方案,如,强制杀掉进程重启,修改其 ...

  6. 使用C#开发C/S框架高级版添加新项目实例

    操作步骤: 1.新建一个项目CSFramework3.test,在CSFramework3.test项目下新建一个FormMan窗口(此时不需要继承). 2.按F7打开类,替换一下引用 using S ...

  7. MacOS下SVN迁移Git踩坑记

    1. First Blood 之前在Windows环境下进行svn到git的迁移是很简单的,(参考官方文档:https://git-scm.com/book/zh/v1/Git-%E4%B8%8E%E ...

  8. JDBC各种数据库连接方式

    1)连接Oracle 8/8i/9i/10g/11g(thin模式) Class.forName("oracle.JDBC.driver.OracleDriver").newIns ...

  9. ORACLE 中NUMBER类型默认的精度和Scale问题

    在ORACLE数据库中,NUMBER(P,S)是最常见的数字类型,可以存放数据范围为10^-130~10^126(不包含此值),需要1~22字节(BYTE)不等的存储空间.P 是Precison的英文 ...

  10. 提升R代码运算效率的11个实用方法——并行、效率

    转载于36大数据,原文作者:Selva Prabhakaran  译者:fibears 众所周知,当我们利用R语言处理大型数据集时,for循环语句的运算效率非常低.有许多种方法可以提升你的代码运算效率 ...