[luoguP4142]洞穴遇险

https://www.zybuluo.com/ysner/note/1240792

题面

戳我

解析

这种用来拼接的奇形怪状的东西,要不就是轮廓线\(DP\),要不就是网络流。

为了表示奇数点(即\((x+y)\%2=1\))的危险值,把该点拆为两个点,连一条边长为该点危险值相反数的边(两点分别称为起点和终点)。

鉴于一根柱子跨越\(3\)个格子,其中一点为奇数点,另外两个点都是偶数点,不能区分。

于是也要把偶数点分为两类(不用拆点)。一类连源点,一类连汇点。连汇点的一类连奇数点的起点,另一类连终点。(让源点出发能到汇点就成)

然后思考如何表示柱子。

如果强行给柱子规定方向,则有\(8\)个方向。
表示出来,有两种情况,一是\(x\)轴方向出,\(y\)轴方向进;另一种是\(y\)轴方向进,\(x\)轴方向出。
于是两个奇数点分别反映一种情况,同时注意相邻的偶数点连奇数点中的起点、还是终点即可。

还要注意的是,最小费用最大流模板求出的是在最大流前提下的最小流,在后期,可能为了得到最大流而付出更多费用(在本题中就是为了放更多柱子而增加不稳定度)。在费用开始非负时(开始退流时)记得\(break\)

唯一一种能让网络流TLE的方式就是cnt=0

// luogu-judger-enable-o2
#include<iostream>
#include<cstdio>
#include<cstdlib>
#include<cstring>
#include<cmath>
#include<algorithm>
#include<vector>
#include<queue>
#define re register
#define il inline
#define ll long long
#define max(a,b) ((a)>(b)?(a):(b))
#define min(a,b) ((a)<(b)?(a):(b))
#define fp(i,a,b) for(re int i=a;i<=b;i++)
#define fq(i,a,b) for(re int i=a;i>=b;i--)
using namespace std;
const int mod=1e9+7,N=5e3+100;
struct Edge{int to,nxt,w,c;}e[N*10];
int n,m,k,s,sum,h[N],d[55][55],ans=0,dis[N],S,T,pe[N],pv[N],tot,cnt=1,g;
bool vis[N],ban[55][55];
il void add(re int u,re int v,re int w,re int c)
{if(u==-1||v==-1) return ;e[++cnt]=(Edge){v,h[u],w,c};h[u]=cnt;e[++cnt]=(Edge){u,h[v],0,-c};h[v]=cnt;
}
queue<int>Q;
il ll gi()
{re ll x=0,t=1;re char ch=getchar();while(ch!='-'&&(ch<'0'||ch>'9')) ch=getchar();if(ch=='-') t=-1,ch=getchar();while(ch>='0'&&ch<='9') x=x*10+ch-48,ch=getchar();return x*t;
}
il int id(re int x,re int y){return (x-1)*n+y;}
il int SPFA()
{memset(dis,63,sizeof(dis));dis[S]=0;vis[S]=1;Q.push(S);while(!Q.empty()){re int u=Q.front();Q.pop();for(re int i=h[u];i+1;i=e[i].nxt){re int v=e[i].to;//printf("!!!%d %d\n",u,v);if(dis[v]>dis[u]+e[i].c&&e[i].w){dis[v]=dis[u]+e[i].c;pe[v]=i;pv[v]=u;if(!vis[v]) vis[v]=1,Q.push(v);}}vis[u]=0;}return dis[T]<dis[0];
}
int main()
{//freopen("marshland.in","r",stdin);//freopen("marshland.out","w",stdout);memset(h,-1,sizeof(h));n=gi();m=gi();k=gi();g=n*n;S=2*g+1;T=2*g+2;fp(i,1,n)fp(j,1,n)d[i][j]=gi(),s+=d[i][j];fp(i,1,k){re int u=gi(),v=gi();ban[u][v]=1;}fp(i,1,n)fp(j,1,n){if(ban[i][j]) continue;if((i+j)%2) add(id(i,j),id(i,j)+g,1,-d[i][j]);if((i+j)%2&&j%2==0){if(i>1&&!ban[i-1][j]) add(id(i,j)+g,id(i-1,j),1,0);if(i<n&&!ban[i+1][j]) add(id(i,j)+g,id(i+1,j),1,0);if(j>1&&!ban[i][j-1]) add(id(i,j-1),id(i,j),1,0);if(j<n&&!ban[i][j+1]) add(id(i,j+1),id(i,j),1,0);}if((i+j)%2&&j%2==1){if(i>1&&!ban[i-1][j]) add(id(i-1,j),id(i,j),1,0);if(i<n&&!ban[i+1][j]) add(id(i+1,j),id(i,j),1,0);if(j>1&&!ban[i][j-1]) add(id(i,j)+g,id(i,j-1),1,0);if(j<n&&!ban[i][j+1]) add(id(i,j)+g,id(i,j+1),1,0);}if((i+j)%2==0&&j%2==1) add(S,id(i,j),1,0);if((i+j)%2==0&&j%2==0) add(id(i,j),T,1,0);}while(SPFA()&&m){if(dis[T]>=0) break;sum=2e9;for(re int i=T;i!=S;i=pv[i])sum=min(sum,e[pe[i]].w);m-=sum;//printf("%d\n",sum);for(re int i=T;i!=S;i=pv[i])e[pe[i]].w-=sum,e[pe[i]^1].w+=sum,ans+=sum*e[pe[i]].c;}printf("%d\n",s+ans);fclose(stdin);fclose(stdout);return 0;
}

转载于:https://www.cnblogs.com/yanshannan/p/9433035.html

本文来自互联网用户投稿,该文观点仅代表作者本人,不代表本站立场。本站仅提供信息存储空间服务,不拥有所有权,不承担相关法律责任。如若转载,请注明出处:http://www.mzph.cn/news/388891.shtml

如若内容造成侵权/违法违规/事实不符,请联系多彩编程网进行投诉反馈email:809451989@qq.com,一经查实,立即删除!

相关文章

飞信虚拟机

做完了一个图片处理软件,突然想到上次上网看到C#程序脱离.NET FRAMEWORK运行的文章,于是决定自己动手试一下。 之前看到有用别的方法来实现的&#xff0c;但我还是选择了现在比较流行的软件飞信中带的VMDotNet&#xff0c;也就是所谓的.NET FRAMEWORK虚拟机吧。相信有很多人也已…

django的contenttype表

https://blog.csdn.net/aaronthon/article/details/81714496 这篇文章已经非常详细了,供自己以后忘了...回看...... 总结&#xff1a; 当一张表和多个表FK关联&#xff0c;并且多个FK中只能选择其中一个或其中n个时&#xff0c;可以利用contenttype&#xff0c;固定用三个字段…

视频播放问题和提高性能方案

1.Five symptoms of poor video performance 1.1 视频加载缓慢 ​Perceived Wait Time Time to first frame (TTFF): ​ 播放开始所需的adaptive bitrate(ABR)流媒体段的数量。(我们稍后将对此进行更详细的讨论。) ​ 视频请求发送到视频加载之间的时间(即接收到足够的数据…

rabbitmq 不同的消费者消费同一个队列_RabbitMQ 消费端限流、TTL、死信队列

消费端限流1. 为什么要对消费端限流假设一个场景&#xff0c;首先&#xff0c;我们 Rabbitmq 服务器积压了有上万条未处理的消息&#xff0c;我们随便打开一个消费者客户端&#xff0c;会出现这样情况: 巨量的消息瞬间全部推送过来&#xff0c;但是我们单个客户端无法同时处理这…

动量策略 python_在Python中使用动量通道进行交易

动量策略 pythonMost traders use Bollinger Bands. However, price is not normally distributed. That’s why only 42% of prices will close within one standard deviation. Please go ahead and read this article. However, I have some good news.大多数交易者使用布林…

css3 变换、过渡效果、动画

1 CSS3 选择器 1.1 基本选择器 1.2 层级 空格 > .itemli ~ .item~p 1.3 属性选择器 [attr] [attrvalue] [attr^value] [attr$value] [attr*value] [][][] 1.4 伪类选择器 :link :visited :hover :active :focus :first-child .list li:first-child :last-chi…

webservice 启用代理服务器

您会发现你写完了一个webservice在调用的时候发现怎也没办法调用&#xff0c;一个简单的webservice怎么不能使用&#xff0c;一肚子的怨恨&#xff0c;哈哈您可能没有为webservice设置代理。 下面就给您写个调用的用例和大家分享下。其实很简单&#xff0c;但是你没有想到的时…

mysql常用的存储引擎_Mysql存储引擎

什么是存储引擎&#xff1f;关系数据库表是用于存储和组织信息的数据结构&#xff0c;可以将表理解为由行和列组成的表格&#xff0c;类似于Excel的电子表格的形式。有的表简单&#xff0c;有的表复杂&#xff0c;有的表根本不用来存储任何长期的数据&#xff0c;有的表读取时非…

android studio设计模式和文本模式切换

转载于:https://www.cnblogs.com/judes/p/9437104.html

高斯模糊为什么叫高斯滤波_为什么高斯是所有发行之王?

高斯模糊为什么叫高斯滤波高斯分布及其主要特征&#xff1a; (Gaussian Distribution and its key characteristics:) Gaussian distribution is a continuous probability distribution with symmetrical sides around its center. 高斯分布是连续概率分布&#xff0c;其中心周…

C# webbrowser 代理

百度&#xff0c;google加自己理解后&#xff0c;将所得方法总结一下&#xff1a; 方法1&#xff1a;修改注册表Software//Microsoft//Windows//CurrentVersion//Internet Settings下 ProxyEnable和ProxyServer。这种方法适用于局域网用户&#xff0c;拨号用户无效。 1p…

C MySQL读写分离连接串_Mysql读写分离

一 什么是读写分离MySQL Proxy最强大的一项功能是实现“读写分离(Read/Write Splitting)”。基本的原理是让主数据库处理事务性查询&#xff0c;而从数据库处理SELECT查询。数据库复制被用来把事务性查询导致的变更同步到集群中的从数据库。当然&#xff0c;主服务器也可以提供…

golang 编写的在线redis 内存分析工具 rma4go

redis 内存分析工具 rma4go redis是一个很有名的内存型数据库&#xff0c;这里不做详细介绍。而rma4go (redis memory analyzer for golang) 是一个redis的内存分析工具&#xff0c;这个工具的主要作用是针对运行时期的redis进行内存的分析&#xff0c;统计redis中key的分布情…

从Jupyter Notebook到脚本

16 Aug: My second article: From Scripts To Prediction API8月16日&#xff1a;我的第二篇文章&#xff1a; 从脚本到预测API As advanced beginners, we know quite a lot: EDA, ML concepts, model architectures etc…… We can write a big Jupyter Notebook, click “Re…

【EasyNetQ】- 使用Future Publish调度事件

许多业务流程要求在将来某个日期安排事件。例如&#xff0c;在与客户进行初次销售联系后&#xff0c;我们可能希望在将来的某个时间安排跟进电话。EasyNetQ可以通过其Future Publish功能帮助您实现此功能。例如&#xff0c;这里我们使用FuturePublish扩展方法来安排未来一个月的…

Java这些多线程基础知识你会吗?

0、并发和并行、进程核线程、多进程和多线程的区别&#xff1a; &#xff08;这里的时间和时刻上的概念同物理上的一样&#xff09; 并发&#xff1a;在一段时间内多个任务同时执行&#xff0c;或者说是在一段很短的时间内可以执行多条程序指令&#xff0c;微观上看起来好像是可…

MySQL set names 命令_mysql set names 命令和 mysql 字符编码问题

先看下面的执行结果&#xff1a;(rootlocalhost)[(none)]mysql>show variables like character%;---------------------------------------------------------------------------------------| Variable_name | Value |---------------------------------------------------…

设置Proxy Server和SQL Server实现数据库安全

首先&#xff0c;我们需要了解一下SQL Server在WinSock上定义协议的步骤&#xff1a; 1. 在”启动”菜单上&#xff0c;指向”程序/Microsoft Proxy Server”&#xff0c;然后点击”Microsoft Management Console”。 2. 展开”Internet Information Service”,再展开运行Proxy…

Python django解决跨域请求的问题

解决方案 1.安装django-cors-headers pip3 install django-cors-headers 2.配置settings.py文件 INSTALLED_APPS [...corsheaders&#xff0c;...] MIDDLEWARE_CLASSES (...corsheaders.middleware.CorsMiddleware,django.middleware.common.CommonMiddleware, # 注意顺序...…

加勒比海兔_加勒比海海洋物种趋势

加勒比海兔Ok, here’s a million dollar question: is the Caribbean really dying? Or, more specifically, are marine species found on Caribbean reefs becoming less abundant?好吧&#xff0c;这是一个百万美元的问题&#xff1a;加勒比海真的死了吗&#xff1f; 或者…