[HAOI2007]理想的正方形

这道题本来想能不能用STL平衡树水过……

一看值的范围在1e9范围内，hash不能……

好吧

MS这道题很老的样子……

正解是对每一列用单调队列维护从这个点开始接下来n个数里面最大的和最小的

那么对于一列上的n个数，我们把它压成了一个数

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

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

那么对于每一行，我们又把它压成了一个数

hmax[i][j],hmin[i][j]就表示这个点为左上角的正方形里面最大和最小的了

然后暴力O(ab)的找答案就可以了

单调队列用的deque，超级慢……其实没什么必要……现在养成了基本数据结构都用stl的习惯……

//Lib #include<cstdio> #include<cstring> #include<cstdlib> #include<cmath> #include<ctime> #include<iostream> #include<algorithm> #include<vector> #include<string> #include<queue> using namespace std; //Macro #define rep(i,a,b) for(int i=a,tt=b;i<=tt;++i) #define drep(i,a,b) for(int i=a,tt=b;i>=tt;--i) #define erep(i,e,x) for(int i=x;i;i=e[i].next) #define irep(i,x) for(__typedef(x.begin()) i=x.begin();i!=x.end();i++) #define read() (strtol(ipos,&ipos,10)) #define sqr(x) ((x)*(x)) #define pb push_back #define PS system("pause"); typedef long long ll; typedef pair<int,int> pii; const int oo=~0U>>1; const double inf=1e100; const double eps=1e-6; string name="", in=".in", out=".out"; //Var int n,m,R,C,ans=oo; int map[1008][1008]; int lmax[1008][1008],lmin[1008][1008],hmax[1008][1008],hmin[1008][1008]; deque<pii> qmax,qmin; void Calc(int col) { qmax.clear();qmin.clear(); rep(i,1,R) { 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.pb(pii(map[i][col],i)); qmin.pb(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 row) { qmax.clear();qmin.clear(); rep(i,1,C) { 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[row][i])qmax.pop_back(); while(!qmin.empty()&&qmin.back().first>=lmin[row][i])qmin.pop_back(); qmax.pb(pii(lmax[row][i],i)); qmin.pb(pii(lmin[row][i],i)); if(i>=n) hmax[row][i-n+1]=qmax.front().first, hmin[row][i-n+1]=qmin.front().first; } } void Init() { scanf("%d%d%d",&R,&C,&n); rep(i,1,R)rep(j,1,C){scanf("%d",&map[i][j]);} } void Work() { rep(i,1,C) Calc(i); rep(i,1,R)Calc2(i); rep(i,1,R-n+1)rep(j,1,C-n+1)ans=min(ans,hmax[i][j]-hmin[i][j]); cout<<ans<<endl; } int main() { // freopen((name+in).c_str(),"r",stdin); // freopen((name+out).c_str(),"w",stdout); Init(); Work(); return 0; }