POJ1156:A STRIP OF LAND

题面
大意：在矩阵中求一个最大的子矩阵，使其中最大值与最小值的差不超过给定的C，且宽不超过100

乍一看：这不是暴力枚举可以搞定的吗？
稍微算算：$O(100n^{2}m)$的枚举，加上$O(100n)$的check，也就$O(10^4{n}^{3}m)$
其实我觉得应该算$O(n^3{m}^3)$
那我们肯定不用这么笨的方法

正文

我不知道原作者是如何想到这么妙的方法的，所以我就不用以往的优化法，直接解析了。
我们用两个循环$l,r$枚举左右边界，使用$maxx$和$minn$两个数组储存每一行中在$l,r$中间的最大、最小值。
然后从上往下扫一个下边界，然后用两个单调队列实现$O(1)$求最值。
具体看代码中的注释

#include<cstdio>
/*/
#include<deque>
/*/
int max(int a,int b)
{return a>b?a:b;}
int min(int a,int b)
{return a<b?a:b;}
template<typename tp>
struct deque
{
	int list[11000],head,tail;
	deque()
	{
		head=tail=1;
	}
	int size(){return tail-head;}
	
	int front(){return list[head];}
	void push_front(int x){list[--head]=x;}
	void pop_front(){head++;}
	
	int back(){return list[tail-1];}
	void push_back(int x){list[tail++]=x;}
	void pop_back(){tail--;}
};//手动模拟deque可加速，并且易于debug
//*/
using namespace std;
int maxx[510],minn[510],a[510][510];
int main()
{
	int n,m,c,ans=0;
	scanf("%d%d%d",&m,&n,&c);
	for(int i=1;i<=n;i++)
		for(int j=1;j<=m;j++)
			scanf("%d",&a[i][j]);
	for(int l=1;l<=m;l++)//枚举左边界
	{
		for(int i=1;i<=n;i++)
			maxx[i]=minn[i]=a[i][l];//初始化每一行的目前最值
		int rmax=min(m,l+99);
		for(int r=l+1;r<=rmax;r++)//枚举右边界
		{
			for(int i=1;i<=n;i++)//更新最值
				maxx[i]=max(maxx[i],a[i][r]),
				minn[i]=min(minn[i],a[i][r]);
			deque<int>maxq,minq;//单调队列：记录当前最值的行号
			int w=r-l+1,high=1,low=1;//w为宽度，high为上边界，low为下边界
			while(low<=n&&(n-high+1)*w>ans)//可以拓展且有拓展价值
			{
				while(maxq.size()&&maxx[maxq.back()]<=maxx[low])maxq.pop_back();
				maxq.push_back(low);
				while(minq.size()&&minn[minq.back()]>=minn[low])minq.pop_back();
				minq.push_back(low);
				//加入当前一行的值
				while(high<=low&&maxx[maxq.front()]-minn[minq.front()]>c)
				//当我们还有一片区域且我们的区域中最值之差不符合要求时，向下移动上边界，
				//尝试把导致不符合要求的元素剔除
				{
					++high;
					if(maxq.front()<high)maxq.pop_front();
					if(minq.front()<high)minq.pop_front();
				}
				ans=max(ans,w*(low-high+1));//记录该矩阵
				++low;//下一个下边界
			}
		}
	}
	printf("%d",ans);
	return 0;
}