CF480E Parking Lot(单调队列+dp然鹅并不是优化)
(全英文题面所以直接放化简题意)
题意:在一个二维平面内,初始有一些点,然后每个时间点加入一些点,对每个时间点求平面内最大的无障碍正方形
(这次的题目是真的神仙啊。。。)
- 首先,考虑暴力,如果对每一个加点进行一遍扫描,那么,可以跑到天荒地老了。。。
- 其次,如果像以前的dp一样跑呢?因为是动态的,所以不行。。。
- 很容易想到,这个面积一定是单调不增的,于是,就可以倒序来,变成加点,离线做。
- 那么,有了加点就可以跑暴力了吗?很显然不行。。。
- 那怎么办呢?
solution:
(以下初始思路来自老师)
这个正方形可能出现在哪里?
- 一个点左上
- 左下
- 右上
- 右下
一个正方形由长和宽构成,要找最大边长,其实就是要维护两个东西:长和宽。
分别考虑长宽如何维护(好吧差不多都是边长)
宽:它可能经过一个点,也就是这个点在宽这条边上。
于是不太容易想到:维护两个dp数组,存每个点向左右能拓展的长度。
但是,如果依旧暴力更新,还是会炸成渣滓。。。
瞪大眼看题,当一个点更新时,只有它所在的一行会变,所以就可以O(n)更新dp数组了。
于是,只要在每个点加上后更新dp数组,然后查出最大正方形,就完事了
下面考虑列:
能不能像刚刚那样继续dp嘞?(像cdq那样dp套dp呢???)如果这样做,那可能会得到一个错位的正方形。。。很麻烦.....
况且我们要去最大的不是吗?这好像有什么提示.....
想想我学过的维护最值的数据结构(qsort rmq 线段树 单调队列 主席树 平衡树 树套树)
再结合一下问题实际:找长度,不太容易想到,对于每一个答案,进行一次check。这个check,有点像二分,对于每一个可行答案进行判断然后更新最优解。而这个可行答案,就是看是否能找到更大的正方形。
对此,只需要枚举是否能找到比当前更大的边长就行了。因为它是正方形,所以长宽相等,这里列问题就转化为了:
对于每一个定长区间,找最大的一个之前处理出的dp数组里是否有一个长度大于等于它的值,也就是最大值。
这就很明朗了,线段树,单调队列。
两者差距在哪里呢?线段树O(nlogn)单调队列O(n);
我用了单调队列。
这份代码给我的感觉就是:代码还可以继续化简,逻辑运算符可以有效增加效率和增加美观度
代码:
#include<bits/stdc++.h>
using namespace std;
const int maxn=;
char s[maxn];
int n,m,k;
int a[maxn][maxn];
struct node
{
int x,y,res;
}ans[maxn];
int f[maxn][maxn];
int l[maxn][maxn];
int r[maxn][maxn]; void init()
{
scanf("%d%d%d",&n,&m,&k);
for(int i=;i<=n;i++)
{
scanf("%s",s+);
for(int j=;j<=m;j++)
{
a[i][j]=(s[j]=='.');
}
}
for(int i=;i<=k;i++)
{
scanf("%d%d",&ans[i].x,&ans[i].y);
a[ans[i].x][ans[i].y]=;
}
} void dp()
{
for(int i=;i<=n;i++)
{
for(int j=m;j>=;j--)
{
if(!a[i][j])
continue;
if(i==||j==m)
{
f[i][j]=;
ans[k].res=max(ans[k].res,);
continue;
}
f[i][j]=min(f[i-][j],f[i][j+]);
if(a[i-f[i][j]][j+f[i][j]])f[i][j]++;
ans[k].res=max(ans[k].res,f[i][j]);
}
}
for(int i=;i<=n;i++)
{
for(int j=;j<=m;j++)
{
l[i][j]=a[i][j]*(l[i][j-]+);
}
}
for(int i=;i<=n;i++)
{
for(int j=m;j>=;j--)
{
r[i][j]=a[i][j]*(r[i][j+]+);
}
}
} int tmp[maxn];
int q[maxn]; bool check(int x,int y)
{
int head=,tail=;
memset(tmp,,sizeof(tmp));
for(int i=;i<=n;i++)
{
while(head<=tail&&l[q[tail]][y]>=l[i][y])
tail--;
q[++tail]=i;
while(q[head]<=i-x)
head++;
tmp[i]+=l[q[head]][y];
}
head=;tail=;
for(int i=;i<=n;i++)
{
while(head<=tail&&r[q[tail]][y]>=r[i][y])
tail--;
q[++tail]=i;
while(q[head]<=i-x)
head++;
tmp[i]+=r[q[head]][y]-;
}
for(int i=x;i<=n;i++)
{
if(tmp[i]>=x)
return ;
}
return ;
} void solve()
{
int tem=ans[k].res;
for(int i=k;i>=;i--)
{
ans[i].res=tem;
int x=ans[i].x;
int y=ans[i].y;
a[x][y]=;
for(int j=;j<=m;j++)
{
l[x][j]=a[x][j]*(l[x][j-]+);
}
for(int j=m;j>=;j--)
{
r[x][j]=a[x][j]*(r[x][j+]+);
}
while(check(tem+,y))
tem++;
}
ans[].res=tem;
for(int i=;i<=k;i++)
{
printf("%d\n",ans[i].res);
}
} int main()
{
init();
dp();
solve();
return ;
}
(完)
CF480E Parking Lot(单调队列+dp然鹅并不是优化)的更多相关文章
- POJ 3017 单调队列dp
Cut the Sequence Time Limit: 2000MS Memory Limit: 131072K Total Submissions: 8764 Accepted: 2576 ...
- [TyvjP1313] [NOIP2010初赛]烽火传递(单调队列 + DP)
传送门 就是个单调队列+DP嘛. ——代码 #include <cstdio> ; , t = , ans = ~( << ); int q[MAXN], a[MAXN], f ...
- zstu 4237 马里奥的求救——(单调队列DP)
题目链接:http://oj.acm.zstu.edu.cn/JudgeOnline/problem.php?id=4237 这题可以转化为每次可以走g~d+x步,求最大分数,且最大分数的步数最少. ...
- 1304F2 - Animal Observation (hard version) 线段树or单调队列 +DP
1304F2 - Animal Observation (hard version) 线段树or单调队列 +DP 题意 用摄像机观察动物,有两个摄像机,一个可以放在奇数天,一个可以放在偶数天.摄像机在 ...
- vijos P1243 生产产品(单调队列+DP)
P1243生产产品 描述 在经过一段时间的经营后,dd_engi的OI商店不满足于从别的供货商那里购买产 品放上货架,而要开始自己生产产品了!产品的生产需要M个步骤,每一个步骤都可以在N台机器 ...
- POJ 1821 单调队列+dp
题目大意:有K个工人,有n个墙,现在要给墙涂色.然后每个工人坐在Si上,他能刷的最大范围是Li,且必须是一个连续子区间,而且必须过Si,他刷完后能获得Pi钱 思路:定义dp[i][j]表示前i个人,涂 ...
- codeforces 1077F2. Pictures with Kittens (hard version)单调队列+dp
被队友催着上(xun)分(lian),div3挑战一场蓝,大号给基佬紫了,结果从D开始他开始疯狂教我做人??表演如何AKdiv3???? 比赛场上:A 2 分钟,B题蜜汁乱计数,结果想得绕进去了20多 ...
- 【LOJ#10180】烽火传递 单调队列+dp
题目大意:给定一个 N 个非负整数数组成的序列,每个点有一个贡献值,现选出其中若干数,使得每连续的 K 个数中至少有一个数被选,要求选出的数贡献值最小. 题解:设 \(dp[i]\) 表示考虑了序列前 ...
- BZOJ 5281--[Usaco2018 Open]Talent Show(分数规划&单调队列&DP)
5281: [Usaco2018 Open]Talent Show Time Limit: 10 Sec Memory Limit: 128 MBSubmit: 79 Solved: 58[Sub ...
随机推荐
- 设置VMWare CentOS7虚拟机上网(配置静态地址)
针对CentOS安装后设置虚拟机上网,参考网上相关资料和实际操作经验总结如下.本人亲测有效,进入主题. 设置虚拟机上网步骤 1.虚拟机设置->网络适配器 2.编辑->虚拟机网络编辑器-&g ...
- python入门之jieba库的使用
对于一段英文,如果希望提取其中的的单词,只需要使用字符串处理的split()方法即可,例如“China is a great country”. 然而对于中文文本,中文单词之间缺少分隔符,这是中文 ...
- 手把手教你吧Python应用到实际开发 不再空谈悟法☝☝☝
手把手教你吧Python应用到实际开发 不再空谈悟法☝☝☝ 想用python做机器学习吗,是不是在为从哪开始挠头?这里我假定你是新手,这篇文章里咱们一起用Python完成第一个机器学习项目.我会手把手 ...
- Integer类型与int的==比较
前言 Java中有两种类型 基本类型 基本数据类类型存的是数值本身 引用类型 引用类型变量在内存放的是数据的引用 基本类型通过==比较的是他们的值大小,而引用类型比较的是他们的引用地址 正文 在一些特 ...
- 攻防世界(XCTF)WEB(进阶区)write up(一)
cat ics-05 ics-06 lottery Cat XCTF 4th-WHCTF-2017 输入域名 输入普通域名无果 输入127.0.0.1返回了ping码的结果 有可能是命令执行 ...
- win10系统plsql卡顿、菜单闪烁解决办法
右键快捷方式--属性--兼容性: 设置为以win7模式运行,以管理员模式运行.如图:
- cinatra::http框架编译使用
cinatra 一个高效易用的c++ http框架 1.下载源码 https://github.com/qicosmos/cinatra 2. vs2017 编译boost库 [msvc-14.1] ...
- python属性的默认值
python类的构造函数中属性可以设置默认值,实例化出来的对象如果属性使用默认值,默认值的地址是相同的. class A: def __init__(self, name = []): self.__ ...
- Linux Centos7 基于Docker 搭建 Nexus私服搭建
创建Blob Stores[本地文件存储目录,统一管理] 1.设置名称和工作路径: ps[注意事项]: 1.storage name:自定义名称 2.storage path:存储路径,默认[/nex ...
- 三维动画形变算法(Mixed Finite Elements)
混合有限元方法通入引入辅助变量后可以将高阶微分问题变成一系列低阶微分问题的组合.在三维网格形变问题中,我们考虑如下泛函极值问题: 其中u: Ω0 → R3是变形体的空间坐标,上述泛函极值问题对应的欧拉 ...