BZOJ3726: PA2014Final Wykladzina
来源:互联网 发布:淘宝注册需要打电话吗 编辑:程序博客网 时间:2024/05/21 13:54
先不考虑可以包含坏格
如果求不包含坏格的最大子矩形,有一种办法是枚举矩形的下界,用一个单调队列扫过去,队列里维护左端点和对应的高度,显然随左端点递增,高度递增,预处理u[i][j]代表(i,j)上方第一个坏格的位置,枚举右端点j扫过去时维护队列里的高度,如果当前队尾不合法,那他在j-1的时候合法且面积最大,更新一下答案,更新一下队列
那么包含一个坏格的最大子矩形就多维护一个队列就行了,每次从不包含坏格队列里出来的子矩形可以把他扔进第二个队列,用u[u[i][j]][j]更新一下它的高度
细节有点多qwq,建议拍一下再交
code:
#include<set>#include<map>#include<deque>#include<queue>#include<stack>#include<cmath>#include<ctime>#include<bitset>#include<string>#include<vector>#include<cstdio>#include<cstdlib>#include<cstring>#include<climits>#include<complex>#include<iostream>#include<algorithm>#define ll long longusing namespace std;const int maxn = 2100;int n,m;char str[maxn];bool v[maxn][maxn];int L[maxn][maxn],U[maxn][maxn];struct node{int p,u;}q[maxn],q2[maxn],t[maxn]; int tp;int head,tail,head2,tail2;int re;int cal(int l,int r,int d,int u) { return (r-l+1)*(u-d+1); }int main(){ scanf("%d%d",&n,&m); for(int i=1;i<=n;i++) { scanf("%s",str+1); for(int j=1;j<=m;j++) if(str[j]=='#') v[i][j]=true; } for(int i=1;i<=n;i++)for(int j=1;j<=m;j++) { L[i][j]=L[i][j-1]; if(v[i][j]) L[i][j]=j; } for(int j=1;j<=m;j++) for(int i=1;i<=n;i++) { U[i][j]=U[i-1][j]; if(v[i][j]) U[i][j]=i; } for(int i=1;i<=n;i++) { head=1,tail=0; head2=1,tail2=0; for(int j=1;j<=m;j++) { int u=U[i][j],uu=0; if(u) uu=U[u-1][j]; int la=0; while(head2<=tail2&&q2[tail2].u<=u) { re=max(re,cal(q2[tail2].p,j-1,q2[tail2].u,i)); la=q2[tail2--].p; } if(!(head2<=tail2&&q[tail2].u==u+1)&&la) q2[++tail2]=(node){la,u+1}; la=0; int la2=0; while(head<=tail&&q[tail].u<=u) { re=max(re,cal(q[tail].p,j-1,q[tail].u,i)); if(q[tail].u>uu) t[++tp]=q[tail]; else la2=q[tail].p; la=q[tail--].p; } if(!(head<=tail&&q[tail].u==u+1)&&la) q[++tail]=(node){la,u+1}; else q[++tail]=(node){j,u+1}; while(tp) { if(!la2||t[tp].p<la2) { while(head2<=tail2&&q2[tail2].p>=t[tp].p) tail2--; q2[++tail2]=t[tp]; } tp--; } if(!(head2<=tail2&&q2[tail2].u==uu+1)&&la2) { while(head2<=tail2&&q2[tail2].p>=la2) tail2--; q2[++tail2]=(node){la2,uu+1}; } else q2[++tail2]=(node){j,uu+1}; } while(head<=tail) re=max(re,cal(q[head].p,m,q[head].u,i)),head++; while(head2<=tail2) re=max(re,cal(q2[head2].p,m,q2[head2].u,i)),head2++; } printf("%d\n",re); return 0;}
阅读全文
0 0
- 【PA2014Final】【BZOJ3726】Wykladzina
- BZOJ3726: PA2014Final Wykladzina
- 【PA2014Final】【BZOJ3721】Bazarek
- 【bzoj3721】 PA2014Final Bazarek 贪心
- 【bzoj3728】 PA2014Final Zarowki
- 【bzoj3728】PA2014Final Zarowki
- BZOJ3728: PA2014Final Zarowki
- BZOJ3724: PA2014Final Krolestwo
- 3728: PA2014Final Zarowki set+priority_queue
- BZOJ 3728: PA2014Final Zarowki|贪心
- BZOJ 3728 PA2014Final Zarowki 堆+贪心
- 【PA】PA2014Final Wykładzina【单调栈+双指针+均摊分析】
- Unity3D C#对字符串的操作
- Apache-Kafka集群安装
- 省略perl命令来运行你常用的perl脚本吧
- Eclipse优化和快捷键大全
- SnippetsLab for Mac(代码片段管理工具)破解版 V1.7激活版SnippetsLab
- BZOJ3726: PA2014Final Wykladzina
- 学习笔记-dubbo介绍和使用
- java.sql.SQLException: ORA-01810: 格式代码出现两次
- leetcode 567. Permutation in String 滑动窗口
- Navicat for MySQL 快捷键
- 目标检测之SSD(single shot multibox detector)的pytorch代码阅读总结
- BioSolveIT Leadit.2.3.2.Win32 1CD
- Unsupported major.minor version 52.0
- codeforces 337A Puzzles