bzoj3144 [HNOI2013]切糕(最小割)

bzoj Luogu

题面描述见上

题解时间

一开始我真就把这玩意所说的切面当成了平面来做的

事实上只是说相邻的切点高度差都不超过 $ d $

对于一条 $ z $ 轴方向的线,把原题的点看成边,每个原题的点两端看成两个点就好(就是说一条线上有 $ r+1 $ 个点 $ r $ 条边),底端每一个点有一条由 $ S $ 连向它的不能断开( $ inf $ )的边,顶端每个点同理连向 $ T $

之后考虑处理相邻两点之间高度差不超过 $ d $

假设我们已经选了线 $ l_1 $ 上的一个原图点 $ p_1 $,那么与之相邻的一条线 $ l_2 $ 上选择的原图点 $ p_2 $ 纵坐标必须在一个区间内

从 $ p_1 $ 两端向 $ p_2 $ 可选范围两端连不能断开( $ inf $ )的边就行

然后直接跑最小割。

#include<bits/stdc++.h>
using namespace std;
template<typename TP>inline void read(TP &tar)
{
TP ret=0,f=1;char ch=getchar();
while(ch<'0'||ch>'9'){if(ch=='-')f=-1;ch=getchar();}
while(ch>='0'&&ch<='9'){ret=ret*10+ch-'0';ch=getchar();}
tar=ret*f;
}
namespace LarjaIX
{
const int N=70011,M=1000011,inf=0x3f3f3f3f;
struct sumireko{int to,ne,w;}e[M<<1];int he[N],ecnt=1;
void addline(int f,int t,int w)
{
e[++ecnt].to=t,e[ecnt].w=w;
e[ecnt].ne=he[f],he[f]=ecnt;
e[++ecnt].to=f,e[ecnt].w=0;
e[ecnt].ne=he[t],he[t]=ecnt;
}
int head[N],dep[N];bool ins[N];
queue<int> qu;
bool bfs(int sp,int ep)
{
memcpy(head,he,sizeof(head));
memset(dep,0x3f,sizeof(dep));
dep[sp]=1,ins[sp]=1,qu.push(sp);
while(!qu.empty())
{
int x=qu.front();qu.pop();ins[x]=0;
for(int ei=he[x],t=e[ei].to;ei;ei=e[ei].ne,t=e[ei].to)
if(e[ei].w&&dep[t]>dep[x]+1)
{
dep[t]=dep[x]+1;
if(!ins[t]) qu.push(t),ins[t]=1;
}
}
return dep[ep]!=inf;
}
int dfs(int x,int lim,int ep)
{
if(x==ep||!lim) return lim;
int ret=0,tmp=0;
for(int ei=head[x],t=e[ei].to;ei;ei=e[ei].ne,t=e[ei].to)
{
head[x]=ei;
if(dep[t]==dep[x]+1) if(tmp=dfs(t,min(lim,e[ei].w),ep))
{
lim-=tmp,ret+=tmp;
e[ei].w-=tmp,e[ei^1].w+=tmp;
}
}
return ret;
}
int dinic(int sp,int ep)
{
int ret=0;
while(bfs(sp,ep))
ret+=dfs(sp,inf,ep);
return ret;
}
int p,q,r,d,a[44][44][44],id[44][44][44],ic,S,E;
int maid()
{
read(p),read(q),read(r),read(d);
for(int z=1;z<=r;z++)for(int x=1;x<=p;x++)for(int y=1;y<=q;y++) read(a[x][y][z]);
r++;
for(int x=1;x<=p;x++)for(int y=1;y<=q;y++)for(int z=1;z<=r;z++) id[x][y][z]=++ic;S=++ic,E=++ic;
for(int x=1;x<=p;x++)for(int y=1;y<=q;y++)
{
addline(S,id[x][y][1],inf),addline(id[x][y][r],E,inf);
for(int z=1;z<r;z++)
{
addline(id[x][y][z],id[x][y][z+1],a[x][y][z]);
if(x+1<=p) addline(id[x][y][z],id[x+1][y][max(1,z-d)],inf),addline(id[x+1][y][min(r,z+1+d)],id[x][y][z+1],inf);
if(x-1) addline(id[x][y][z],id[x-1][y][max(1,z-d)],inf),addline(id[x-1][y][min(r,z+1+d)],id[x][y][z+1],inf);
if(y+1<=q) addline(id[x][y][z],id[x][y+1][max(1,z-d)],inf),addline(id[x][y+1][min(r,z+1+d)],id[x][y][z+1],inf);
if(y-1) addline(id[x][y][z],id[x][y-1][max(1,z-d)],inf),addline(id[x][y-1][min(r,z+1+d)],id[x][y][z+1],inf);
}
}
printf("%d\n",dinic(S,E));
return 0;
}
}
int main(){return LarjaIX::maid();}

bzoj3144 [HNOI2013]切糕(最小割)的更多相关文章

  1. BZOJ3144[Hnoi2013]切糕——最小割

    题目描述 输入 第一行是三个正整数P,Q,R,表示切糕的长P. 宽Q.高R.第二行有一个非负整数D,表示光滑性要求.接下来是R个P行Q列的矩阵,第z个 矩阵的第x行第y列是v(x,y,z) (1≤x≤ ...

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

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

  3. 【BZOJ-3144】切糕 最小割-最大流

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

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

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

  5. Luogu P3227 [HNOI2013]切糕 最小割

    首先推荐一个写的很好的题解,个人水平有限只能写流水账,还请见谅. 经典的最小割模型,很多人都说这个题是水题,但我还是被卡了=_= 技巧:加边表示限制 在没有距离\(<=d\)的限制时候,我们对每 ...

  6. bzoj 3144 [Hnoi2013]切糕——最小割

    题目:https://www.lydsy.com/JudgeOnline/problem.php?id=3144 一根纵轴上切一个点,可以把一根纵轴上的点连成一串来体现.自己的写法是每个点连向前一个点 ...

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

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

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

    题目 输入格式 第一行是三个正整数P,Q,R,表示切糕的长P. 宽Q.高R.第二行有一个非负整数D,表示光滑性要求.接下来是R个P行Q列的矩阵,第z个 矩阵的第x行第y列是v(x,y,z) (1≤x≤ ...

  9. BZOJ3144/LG3227 「HNOI2013」切糕 最小割离散变量模型

    问题描述 BZOJ3144 LG3227 还想粘下样例 输入: 2 2 2 1 6 1 6 1 2 6 2 6 输出: 6 题解 关于离散变量模型,我不想再抄一遍,所以: 对于样例,可以建立出这样的图 ...

随机推荐

  1. PHP获取用户IP地址

    PHP获取访问者IP地址 这是一段 PHP 代码,演示了如何获得来访者的IP address. <?php//打印出IP地址:echo (GetIP());function GetIP()  / ...

  2. python-xlutils模块-修改excel

    前言 xlrd库用于读取excel文件中的数据,xlwt库用于将数据写入excel文件,修改用xlutils模块: xlutils库也仅仅是通过复制一个副本进行操作后保存一个新文件,像是xlrd库和x ...

  3. PentestBox在win10里打不开工具 显示无系统命令的解决方法

    PentestBox详细安装过程:http://www.cnblogs.com/ESHLkangi/p/8336398.html 在使用PentestBox的时候出现了打不开工具的问题,最后看到一个大 ...

  4. mapreduce 中 groupingComparator 用法

    groupingComparator是对reduce输入的数据进行分组,比如 public int compare(WritableComparable a, WritableComparable b ...

  5. Qt:QNetworkReply

    0.说明 QNetworkReply对象包含了Manager发送的请求头和返回的数据. 它继承自QIODevice,所以可以用各种read获取其中返回的数据: QByteArray data = re ...

  6. c++刷leetcode记录

    #include<iostream> #include<sstream> #include<vector> std::vector<int> split ...

  7. 探索TiDB Lightning的源码来解决发现的bug

    背景 上一篇<记一次简单的Oracle离线数据迁移至TiDB过程>说到在使用Lightning导入csv文件到TiDB的时候发现了一个bug,是这样一个过程. Oracle源库中表名都是大 ...

  8. 1.java连接pulsar服务

    目录 是什么 安装 Java客户端 1.引入GAV 2.创建配置项 3.验证测试 其他及代码下载 是什么 Pulsar 是一个用于服务器到服务器的消息系统,具有多租户.高性能等优势.详见 安装 本文主 ...

  9. linux访问控制列表 ACL实现文件权限设置

    ACL:Access Control List,实现灵活的文件权限管理 除了文件的所有者,所属组和其它人,可以对更多的用户设置权限 CentOS7 默认创建的xfs和ext4文件系统具有ACL功能 A ...

  10. tp5 (自写) 实现redis消息队列 + 排行榜

    1:小皮开启redis, 控制器按Ctrl 点击new Redis 进入 redis.php 进行封装 //向队列添加数据 // LPUSH key value1 [value2] //将一个或多个值 ...