洛谷-1169 [ZJOI2007]棋盘制作

题目描述
国际象棋是世界上最古老的博弈游戏之一，和中国的围棋、象棋以及日本的将棋同享盛名。据说国际象棋起源于易经的思想，棋盘是一个8×8大小的黑白相间的方阵，对应八八六十四卦，黑白对应阴阳。
而我们的主人公小Q，正是国际象棋的狂热爱好者。作为一个顶尖高手，他已不满足于普通的棋盘与规则，于是他跟他的好朋友小W决定将棋盘扩大以适应他们的新规则。
小Q找到了一张由N×M个正方形的格子组成的矩形纸片，每个格子被涂有黑白两种颜色之一。小Q想在这种纸中裁减一部分作为新棋盘，当然，他希望这个棋盘尽可能的大。
不过小Q还没有决定是找一个正方形的棋盘还是一个矩形的棋盘（当然，不管哪种，棋盘必须都黑白相间，即相邻的格子不同色），所以他希望可以找到最大的正方形棋盘面积和最大的矩形棋盘面积，从而决定哪个更好一些。
于是小Q找到了即将参加全国信息学竞赛的你，你能帮助他么？
输入输出格式
输入格式：
包含两个整数N和M，分别表示矩形纸片的长和宽。接下来的N行包含一个N ×M的01矩阵，表示这张矩形纸片的颜色（0表示白色，1表示黑色）。
输出格式：
包含两行，每行包含一个整数。第一行为可以找到的最大正方形棋盘的面积，第二行为可以找到的最大矩形棋盘的面积（注意正方形和矩形是可以相交或者包含的）。

输入输出样例
输入样例#1：
3 3
1 0 1
0 1 0
1 0 0

输出样例#1：
4
6

说明
对于20%的数据，N, M ≤ 80
对于40%的数据，N, M ≤ 400
对于100%的数据，N, M ≤ 2000

解释：我们首先预处理出来每点向上，向左，向右可以延伸的最大长度，然后我们枚举行，这样抽象成最大矩形面积的经典问题，最后特殊处理一下边界就好

#include<iostream>
#include<cstring>
#include<stack> 
using namespace std;
int G[2003][2003]={0};
int sum[2003][2003]={0};
int L[2003]={0},R[2003]={0};
int LL[2003][2003]={0},RR[2003][2003]={0};
int n=0,m=0;
int ans1=0,ans2=0;
stack<int> s;
int main(){
	ios::sync_with_stdio(false);
	cin>>n>>m;
	for(int i=1;i<=n;i++) for(int j=1;j<=m;j++) cin>>G[i][j];
	for(int i=1;i<=n;i++){
		for(int j=1;j<=m;j++){
			if(G[i][j]!=G[i-1][j]) sum[i][j]=sum[i-1][j]+1;
			else sum[i][j]=1;
		}
	}
	for(int i=1;i<=m;i++){
		for(int j=1;j<=n;j++){
			if(G[j][i]!=G[j][i-1]) LL[j][i]=LL[j][i-1]+1;
			else LL[j][i]=1;
		}
	}
	for(int i=m;i>=1;i--){
		for(int j=1;j<=n;j++){
			if(G[j][i]!=G[j][i+1]) RR[j][i]=RR[j][i+1]+1;
			else RR[j][i]=1;
		}
	}
	for(int i=1;i<=n;i++){
		memset(L,0,sizeof(L));
		memset(R,0,sizeof(R));
		for(int j=1;j<=m;j++){
			if(s.empty()||sum[i][s.top()]<=sum[i][j]) s.push(j);
			else{
				while(s.size()&&sum[i][s.top()]>sum[i][j]){
					R[s.top()]=j-1;s.pop();
				}
				s.push(j);
			}
		}
		while(s.size()){R[s.top()]=m;s.pop();}
		for(int j=m;j>=1;j--){
			if(s.empty()||sum[i][s.top()]<=sum[i][j]) s.push(j);
			else{
				while(s.size()&&sum[i][s.top()]>sum[i][j]){
					L[s.top()]=j+1;s.pop();
				}
				s.push(j);
			}
		}
		while(s.size()){L[s.top()]=1;s.pop();}
		for(int j=1;j<=m;j++){
			ans1=max(ans1,(min(RR[i][j]+j-1,R[j])-max(1+j-LL[i][j],L[j])+1)*sum[i][j]);
			ans2=max(ans2,min(sum[i][j],min(RR[i][j]+j-1,R[j])-max(1+j-LL[i][j],L[j])+1)*min(sum[i][j],min(RR[i][j]+j-1,R[j])-max(1+j-LL[i][j],L[j])+1));
		}
	}
	cout<<ans2<<endl;
	cout<<ans1<<endl;
	return 0;
}