[HNOI2013]切糕
题目描述
网址: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]切糕的更多相关文章
- BZOJ 3144: [Hnoi2013]切糕
3144: [Hnoi2013]切糕 Time Limit: 10 Sec Memory Limit: 128 MBSubmit: 1495 Solved: 819[Submit][Status] ...
- bzoj 3144: [Hnoi2013]切糕 最小割
3144: [Hnoi2013]切糕 Time Limit: 10 Sec Memory Limit: 128 MBSubmit: 681 Solved: 375[Submit][Status] ...
- BZOJ_3144_[Hnoi2013]切糕_最小割
BZOJ_3144_[Hnoi2013]切糕_最小割 Description Input 第一行是三个正整数P,Q,R,表示切糕的长P. 宽Q.高R.第二行有一个非负整数D,表示光滑性要求.接下来是R ...
- bzoj千题计划142:bzoj3144: [Hnoi2013]切糕
http://www.lydsy.com/JudgeOnline/problem.php?id=3144 如果D=2 ,两个点,高度为4,建图如下 #include<queue> #inc ...
- 【BZOJ3144】[HNOI2013]切糕
[BZOJ3144][HNOI2013]切糕 题面 题目描述 经过千辛万苦小 A 得到了一块切糕,切糕的形状是长方体,小 A 打算拦腰将切糕切成两半分给小 B.出于美观考虑,小 A 希望切面能尽量光滑 ...
- 【BZOJ 3144】 3144: [Hnoi2013]切糕 (最小割模型)
3144: [Hnoi2013]切糕 Time Limit: 10 Sec Memory Limit: 128 MBSubmit: 1764 Solved: 965 Description Inp ...
- BZOJ3144 Hnoi2013 切糕 【网络流】*
BZOJ3144 Hnoi2013 切糕 Description Input 第一行是三个正整数P,Q,R,表示切糕的长P. 宽Q.高R.第二行有一个非负整数D,表示光滑性要求.接下来是R个P行Q列的 ...
- 【BZOJ3144】[Hnoi2013]切糕 最小割
[BZOJ3144][Hnoi2013]切糕 Description Input 第一行是三个正整数P,Q,R,表示切糕的长P. 宽Q.高R.第二行有一个非负整数D,表示光滑性要求.接下来是R个P行Q ...
- 3144: [Hnoi2013]切糕
3144: [Hnoi2013]切糕 Time Limit: 10 Sec Memory Limit: 128 MBSubmit: 1526 Solved: 827[Submit][Status] ...
- bzoj3144 [HNOI2013]切糕(最小割)
bzoj3144 [HNOI2013]切糕(最小割) bzoj Luogu 题面描述见上 题解时间 一开始我真就把这玩意所说的切面当成了平面来做的 事实上只是说相邻的切点高度差都不超过 $ d $ 对 ...
随机推荐
- java程序设计原则
前言: 前言:java这种面向对象的的语言其实是很符合人的生活思维模式的,因为它以对象作为自己的研究对象,所谓"万物皆对象".一个人.一颗苹果.一只蚂蚁都是对象.所以它的设计原则和 ...
- XAMPP简介、安转、使用
虽然没有写Mac安装方法及使用, 但方法也都大相径庭, 殊途同归而已. XAMPP简介 XAMPP是一款开源.免费的网络服务器软件,经过简单安装后,就可以在个人电脑上搭建服务器环境.本文为大家介绍Wi ...
- php分布式redis实现session共享
方法一:找到配置文件php.ini,修改为下面内容,保存并重启服务 session.save_handler = redis session.save_path = "tcp://127.0 ...
- 启动MySql提示:The server quit without updating PID file(…)失败
1.可能是/usr/local/mysql/data/rekfan.pid文件没有写的权限解决方法 :给予权限,执行 "chown -R mysql:mysql /var/data" ...
- MysqL错误之_ERROR! MySQL server PID file could not be found!
在配置Mysql主从GTID模式下,启动Mysql服务时出现报错,搜索了一番,找到了一个简单可靠的方法,直接成功.如果遇到相同问题没有解决的童鞋,那就去试一下很多其他方案,如,强制杀掉进程重启,修改其 ...
- 使用C#开发C/S框架高级版添加新项目实例
操作步骤: 1.新建一个项目CSFramework3.test,在CSFramework3.test项目下新建一个FormMan窗口(此时不需要继承). 2.按F7打开类,替换一下引用 using S ...
- MacOS下SVN迁移Git踩坑记
1. First Blood 之前在Windows环境下进行svn到git的迁移是很简单的,(参考官方文档:https://git-scm.com/book/zh/v1/Git-%E4%B8%8E%E ...
- JDBC各种数据库连接方式
1)连接Oracle 8/8i/9i/10g/11g(thin模式) Class.forName("oracle.JDBC.driver.OracleDriver").newIns ...
- ORACLE 中NUMBER类型默认的精度和Scale问题
在ORACLE数据库中,NUMBER(P,S)是最常见的数字类型,可以存放数据范围为10^-130~10^126(不包含此值),需要1~22字节(BYTE)不等的存储空间.P 是Precison的英文 ...
- 提升R代码运算效率的11个实用方法——并行、效率
转载于36大数据,原文作者:Selva Prabhakaran 译者:fibears 众所周知,当我们利用R语言处理大型数据集时,for循环语句的运算效率非常低.有许多种方法可以提升你的代码运算效率 ...