【单调队列】COGS 308. [HAOI2007] 理想的正方形

orz大神，附上链接：http://blog.csdn.net/orpinex/article/details/7179078

题目链接：http://cojs.tk/cogs/problem/problem.php?pid=308

分析：

对每一列用单调队列维护从这个点开始接下来n个数里面最大的和最小的
那么对于一列上的n个数，我们把它压成了一个数

lmax[i][j],lmin[i][j]分别表示以这个点为起点往下n个数里面最大和最小

那么我们对这个两个数组横向的也用单调队列维护

那么对于每一行，我们又把它压成了一个数
hmax[i][j],hmin[i][j]就表示这个点为左上角的正方形里面最大和最小的了
然后暴力O(ab)的找答案就可以了
代码：

#include<cstdio>
#include<iostream>
#include<queue>
using namespace std;

int lmax[1008][1008],lmin[1008][1008],
hmax[1008][1008],hmin[1008][1008],
n,m,r,c,map[1008][1008],ans=2000000000;
typedef pair<int,int> pii;
deque<pii> qmax,qmin;
void calc(int col)
{
qmax.clear();qmin.clear();
for (int i=1;i<=r;i++)
{
while (!qmax.empty()&&i-qmax.front().second>=n) qmax.pop_front();
while (!qmin.empty()&&i-qmin.front().second>=n) qmin.pop_front();
while (!qmax.empty()&&qmax.back().first<=map[i][col]) qmax.pop_back();
while (!qmin.empty()&&qmin.back().first>=map[i][col]) qmin.pop_back();
qmax.push_back(pii(map[i][col],i));
qmin.push_back(pii(map[i][col],i));
if (i>=n)
lmax[i-n+1][col]=qmax.front().first,
lmin[i-n+1][col]=qmin.front().first;
}
}
void calc2(int col)
{
qmax.clear();qmin.clear();
for (int i=1;i<=c;i++)
{
while (!qmax.empty()&&i-qmax.front().second>=n) qmax.pop_front();
while (!qmin.empty()&&i-qmin.front().second>=n) qmin.pop_front();
while (!qmax.empty()&&qmax.back().first<=lmax[col][i]) qmax.pop_back();
while (!qmin.empty()&&qmin.back().first>=lmin[col][i]) qmin.pop_back();
qmax.push_back(pii(lmax[col][i],i));
qmin.push_back(pii(lmin[col][i],i));
if (i>=n)
hmax[col][i-n+1]=qmax.front().first,
hmin[col][i-n+1]=qmin.front().first;
}
}

int main()
{
freopen("square.in","r",stdin);
freopen("square.out","w",stdout);
cin>>r>>c>>n;
for(int i=1;i<=r;i++)
for (int j=1;j<=c;j++)cin>>map[i][j];
for (int i=1;i<=c;i++)  calc(i);
for (int j=1;j<=r;j++)  calc2(j);
for (int i=1;i<=r-n+1;i++)
for (int j=1;j<=c-n+1;j++)
ans=min(ans,hmax[i][j]-hmin[i][j]);
cout<<ans<<endl;
return 0;
}