博客
关于我
强烈建议你试试无所不能的chatGPT,快点击我
[ZJOI2007]棋盘制作 【最大同色矩形】
阅读量:5060 次
发布时间:2019-06-12

本文共 2876 字,大约阅读时间需要 9 分钟。

题目描述

国际象棋是世界上最古老的博弈游戏之一,和中国的围棋、象棋以及日本的将棋同享盛名。据说国际象棋起源于易经的思想,棋盘是一个8*8大小的黑白相间的方阵,对应八八六十四卦,黑白对应阴阳。

而我们的主人公小Q,正是国际象棋的狂热爱好者。作为一个顶尖高手,他已不满足于普通的棋盘与规则,于是他跟他的好朋友小W决定将棋盘扩大以适应他们的新规则。

小Q找到了一张由N*M个正方形的格子组成的矩形纸片,每个格子被涂有黑白两种颜色之一。小Q想在这种纸中裁减一部分作为新棋盘,当然,他希望这个棋盘尽可能的大。

不过小Q还没有决定是找一个正方形的棋盘还是一个矩形的棋盘(当然,不管哪种,棋盘必须都黑白相间,即相邻的格子不同色),所以他希望可以找到最大的正方形棋盘面积和最大的矩形棋盘面积,从而决定哪个更好一些。

于是小Q找到了即将参加全国信息学竞赛的你,你能帮助他么?

输入输出格式

输入格式:

包含两个整数N和M,分别表示矩形纸片的长和宽。接下来的N行包含一个N * M的01矩阵,表示这张矩形纸片的颜色(0表示白色,1表示黑色)。

输出格式:

包含两行,每行包含一个整数。第一行为可以找到的最大正方形棋盘的面积,第二行为可以找到的最大矩形棋盘的面积(注意正方形和矩形是可以相交或者包含的)。

输入输出样例

输入样例#1:

3 31 0 10 1 01 0 0
输出样例#1:

46

说明

对于20%的数据,N, M ≤ 80

对于40%的数据,N, M ≤ 400

对于100%的数据,N, M ≤ 2000

题解

首先,我们选出所有01联通块,只需要一个很简单的做法,就是隔一个反转一个数字,具体来讲,就是将i+j为偶数的(i,j)反转
反转之后呢,就是  最大同色矩形问题

最大同色矩形问题

最大同色矩形问题,顾名思义,就是求一个矩阵中的最大同色矩形
用到 动归 + 并查集 + 排序
具体而言,我们先求出h[i][j]为从(i,j)点向右扩展最远的距离
利用dp可以在O(n^2)内求出
for (int i = 1; i <= N; i++){		h[i][M] = 1;		for (int j = M - 1; j >= 1; j--){			h[i][j] = (A[i][j] == A[i][j + 1] ? h[i][j + 1] + 1 : 1);		}	}
这样,我们就可以很快求出横向边的长度,如何求出纵向边呢?
由于是纵向边,我们就在每一列讨论:
对于j这一列,我们按h的大小从大到小枚举,枚举到当前行i时,包含i的矩形的横向边长度一定是i,而纵向边的长度就是当前从i出发向上向下能经过已访问点的最长长度【因为只有已访问的点的h比i大,才能使h[i][j]作为横向长度】,具体用两个并查集l[i],r[i]实现
#include
#include
#include
#include
#define LL long long intusing namespace std;const int maxn = 2005,maxm = 100005,INF = 2000000000;inline int read(){ int out = 0,flag = 1;char c = getchar(); while (c < 48 || c > 57) {if (c == '-') flag = -1;c = getchar();} while (c >= 48 &&c <= 57) {out = out * 10 + c - 48;c = getchar();} return out * flag;}int N,M,A[maxn][maxn],l[maxn],r[maxn],h[maxn][maxn];bool vis[maxn];struct node{ int len,id;}tmp[maxn];inline bool operator < (const node& a,const node& b){ return a.len > b.len;}inline int findl(int u) {return u == l[u] ? u : l[u] = findl(l[u]);}inline int findr(int u) {return u == r[u] ? u : r[u] = findr(r[u]);}void init(){ N = read(); M = read(); for (int i = 1; i <= N; i++) for (int j = 1; j <= M; j++) A[i][j] = read()^((i^j) & 1);}void solve(){ int ans1 = 0,ans2 = 0,a,b; for (int i = 1; i <= N; i++){ h[i][M] = 1; for (int j = M - 1; j >= 1; j--){ h[i][j] = (A[i][j] == A[i][j + 1] ? h[i][j + 1] + 1 : 1); } } for (int j = 1; j <= M; j++){ for (int i = 1; i <= N; i++){ tmp[i].len = h[i][j]; tmp[i].id = i; vis[i] = false; l[i] = r[i] = i; } sort(tmp + 1,tmp + 1 + N); for (int i = 1; i <= N; i++){ int k = tmp[i].id; vis[k] = true; if (k > 1 && vis[k - 1] && A[k - 1][j] == A[k][j]){ l[k] = k - 1; r[k - 1] = k; } if (k < N && vis[k + 1] && A[k + 1][j] == A[k][j]){ r[k] = k + 1; l[k + 1] = k; } a = findr(k) - findl(k) + 1; b = tmp[i].len; ans1 = max(ans1,min(a,b) * min(a,b)); ans2 = max(ans2,a * b); } } printf("%d\n%d\n",ans1,ans2);}int main(){ init(); solve(); return 0;}

转载于:https://www.cnblogs.com/Mychael/p/8282863.html

你可能感兴趣的文章
eggs
查看>>
一步步学习微软InfoPath2010和SP2010--第七章节--从SP列表和业务数据连接接收数据(4)--外部项目选取器和业务数据连接...
查看>>
如何增强你的SharePoint 团队网站首页
查看>>
FZU 1914 Funny Positive Sequence(线性算法)
查看>>
oracle 报错ORA-12514: TNS:listener does not currently know of service requested in connec
查看>>
基于grunt构建的前端集成开发环境
查看>>
MySQL服务读取参数文件my.cnf的规律研究探索
查看>>
java string(转)
查看>>
__all__有趣的属性
查看>>
写博客
查看>>
利用循环播放dataurl的视频来防止锁屏:NoSleep.js
查看>>
python3 生成器与迭代器
查看>>
java编写提升性能的代码
查看>>
ios封装静态库技巧两则
查看>>
Educational Codeforces Round 46 (Rated for Div. 2)
查看>>
Abstract Factory Pattern
查看>>
C# 实现Bresenham算法(vs2010)
查看>>
基于iSCSI的SQL Server 2012群集测试(一)--SQL群集安装
查看>>
list 容器 排序函数.xml
查看>>
存储开头结尾使用begin tran,rollback tran作用?
查看>>