Codeforces Round #619 (Div. 2)E思维+二维RMQ
题:https://codeforces.com/contest/1301/problem/E
题意:给个n*m的图形,q个询问,每次询问问询问区间最大的合法logo的面积是多少
分析:由于logo是由4个同色的小正方形组成的,所以我们先考虑一个数组val[i][j]表示[i][j]位置为‘R' , [i+1][j]位置为’Y‘,[i][j+1]位置为’G‘,[i+1][j+1]位置为’B'所能形成的最大图案边长是多少;
因为这样就只要预处理‘R’正方形沿着左上能形成的最大正方形,“Y”沿着左下能形成的最大正方形,’G‘沿着右上能形成的最大正方形,”B“沿着右下能形成的最大正方形,然后对于位置[i][j]就取min这四者;
对于询问,我们只要二分答案搜到最大符合条件的就行了,注意这里的二分的mid,假设询问[r1,c1,r2,c2],则搜索访问应该在[r1+mid-1,c1+mid-1,r2-mid,c2-mid]的范围内搜索才是合法的;
因为我们在二分时要求满足条件的mid就行,所以这里对val进行倍增处理;
合法性:因为小的合法图案一定在大的合法图案里面,所以对val取的图案也许是超过询问范围,但我们上面已经对询问范围进行如上缩减,所以是合法的
#include<bits/stdc++.h>
using namespace std;
const int M=;
const int N=;
int ex[][M][M];
int n,m,q;
int val[M][M],f[N][N][M][M],lg[M];
char s[M][M];
void RMQ(int n,int m){
for(int i=;i<n;i++)
for(int j=;j<m;j++)
f[][][i+][j+]=val[i][j];
for(int i=;i<M;i++)
lg[i]=+lg[i/];
for(int i=;i<=n;i++)
for(int k2=;(<<k2)<=m;k2++)
for(int j=;j<=m-(<<k2)+;j++)
f[][k2][i][j]=max(f[][k2-][i][j],f[][k2-][i][j+(<<(k2-))]);
for(int k1=;(<<k1)<=n;k1++)
for(int i=;i<=n-(<<k1)+;i++)
for(int k2=;(<<k2)<=m;k2++)
for(int j=;j<=m-(<<k2)+;j++)
f[k1][k2][i][j]=max(f[k1-][k2][i][j],f[k1-][k2][i+(<<(k1-))][j]);
} int query(int x1,int y1,int x2,int y2){
int k1=lg[x2-x1+];
int k2=lg[y2-y1+];
x2=x2-(<<k1)+;
y2=y2-(<<k2)+;
return max(f[k1][k2][x1][y1],max(f[k1][k2][x1][y2],max(f[k1][k2][x2][y1],f[k1][k2][x2][y2])));
}
bool check(int r1,int c1,int r2,int c2,int midd){
r1=r1+midd-;
c1=c1+midd-;
r2=r2-midd;
c2=c2-midd;
if(r1>r2||c1>c2)
return false;
return query(r1,c1,r2,c2)>=midd;
}
void color(int fx,int fy,char c,int k){///根据进来的参数不同,组成不同方向相同颜色的最大正方形
int stx=,sty=;
if(fx==-)stx=n-;
if(fy==-)sty=m-;
///while遍历二维数组
while(stx<n&&stx>=){
while(sty<m&&sty>=){
if(s[stx][sty]==c){
ex[k][stx][sty]=;
if(stx-fx>=&&stx-fx<n&&sty-fy>=&&sty-fy<m){///正方形区域取最小值
if(s[stx-fx][sty]==c&&s[stx][sty-fy]==c&&s[stx-fx][sty-fy]==c){
ex[k][stx][sty]=min(ex[k][stx-fx][sty],
min(ex[k][stx][sty-fy],ex[k][stx-fx][sty-fy]))+;
}
}
}
sty+=fy;
}
sty=;if(fy==-)sty=m-;
stx+=fx;
}
}
int main(){
scanf("%d%d%d",&n,&m,&q);
for(int i=;i<n;i++)
scanf("%s",s[i]);
//
color(,,'R',);
color(-,,'Y',);
color(,-,'G',);
color(-,-,'B',);
for(int i=;i<n-;i++){
for(int j=;j<m-;j++){
if(s[i][j]=='R'&&s[i+][j]=='Y'&&s[i][j+]=='G'&&s[i+][j+]=='B'){
val[i][j]=min(ex[][i][j],min(ex[][i+][j],min(ex[][i][j+],ex[][i+][j+])));
}
}
}
//
RMQ(n,m); //
for(int i=;i<q;i++){
int r1,c1,r2,c2,ans=;
scanf("%d%d%d%d",&r1,&c1,&r2,&c2);
int l=,r=min(r2-r1,c2-c1);
while(l<=r){
int midd=(l+r)>>;
if(check(r1,c1,r2,c2,midd)){
ans=midd;
l=midd+;
}
else
r=midd-;
}
printf("%d\n",*ans*ans);
}
return ;
}
Codeforces Round #619 (Div. 2)E思维+二维RMQ的更多相关文章
- Codeforces Round #368 (Div. 2) E. Garlands 二维树状数组 暴力
E. Garlands 题目连接: http://www.codeforces.com/contest/707/problem/E Description Like all children, Ale ...
- Codeforces Round #619 (Div. 2)
A. Three Strings 题意:给三个长度相同的非空字符串abc,依次将c中的每个字符和a或者b中对应位置的字符进行交换,交换必须进行,问能否使得ab相同. 思路:对于每一个位置,如果三个字符 ...
- codeforces 713D D. Animals and Puzzle 二分+二维rmq
题目链接 给一个01矩阵, 然后每个询问给出两个坐标(x1, y1), (x2, y2). 问你这个范围内的最大全1正方形的边长是多少. 我们dp算出以i, j为右下角的正方形边长最大值. 然后用二维 ...
- 【CF1256】Codeforces Round #598 (Div. 3) 【思维+贪心+DP】
https://codeforces.com/contest/1256 A:Payment Without Change[思维] 题意:给你a个价值n的物品和b个价值1的物品,问是否存在取物方案使得价 ...
- Codeforces Round #143 (Div. 2) (ABCD 思维场)
题目连链接:http://codeforces.com/contest/231 A. Team time limit per test:2 seconds memory limit per test: ...
- Codeforces Round #395 (Div. 2)(A.思维,B,水)
A. Taymyr is calling you time limit per test:1 second memory limit per test:256 megabytes input:stan ...
- Codeforces Round #416 (Div. 2)(A,思维题,暴力,B,思维题,暴力)
A. Vladik and Courtesy time limit per test:2 seconds memory limit per test:256 megabytes input:stand ...
- Codeforces Round #533 (Div. 2) C.思维dp D. 多源BFS
题目链接:https://codeforces.com/contest/1105 C. Ayoub and Lost Array 题目大意:一个长度为n的数组,数组的元素都在[L,R]之间,并且数组全 ...
- Codeforces Round #539 (Div. 2) D 思维
https://codeforces.com/contest/1113/problem/D 题意 将一个回文串切成一段一段,重新拼接,组成一个新的回文串,问最少切几刀 题解 首先无论奇偶串,最多只会切 ...
随机推荐
- 打开exe并传参
shellexecute(Application.Handle,'open',PWideChar('E:\控件\TMS.Scripter.Studio.Pro..6.0.2.0.Delphi.BCB. ...
- django ORM多对多正向查询时查询返回结果为None
表 class Books(models.Model): '''书籍''' id = models.AutoField(primary_key=True) name = models.CharFiel ...
- NOR Flash驱动
驱动程序 1 ] ] );81 ;83 }84 85 86 static void __exit nor_exit(void)87 {88 iounmap(nor_ ...
- 一百一十五、脱离SAP本体,通过ActiveX读取SAP表中数据
一.Sap自带有客户端,但是非常之臃肿卡顿,可以利用ActiveX的方式,脱离Sap本体,来读取Sap表中的内容进行插入等操作,非常之方便.代码如下: 二.界面如下,输入好相关内容,点击登录,提示登录 ...
- SQL分组后获取其中一个字段最大值的整条记录
SELECT * FROM( SELECT id,name,counts,createDate,row_number() OVER(partition BY name ORDER BY createD ...
- excel处理经纬度
=LEFT(A1,FIND("°",A1)-1)*1+MID(A1,FIND("°",A1)+1,2)/60+MID(A1,FIND("′" ...
- Spark SQL 笔记
Spark SQL 简介 SparkSQL 的前身是 Shark, SparkSQL 产生的根本原因是其完全脱离了 Hive 的限制.(Shark 底层依赖于 Hive 的解析器, 查询优化器) Sp ...
- Day 15:缓冲输入输出常用方法和小练习
以拷贝图片为例子,演示异常处理的代码: 拷贝一张图片 import java.io.File;import java.io.FileInputStream;import java.io.FileOut ...
- MFC 选择文件夹
WCHAR szPath[_MAX_PATH] = {}; BROWSEINFO bi; //指定父窗口,在对话框显示期间,父窗口将被禁用 bi.hwndOwner = this->GetSaf ...
- c++ 字符串转数字或数字转字符串
在C++中字符串转换为数字,或数字转换为字符串,用到如下函数: _itoa atoi.atof.itoa.itow _itoa_s 1.整形转换为字符串: wchar_t * _itot(int _V ...