HAOI2007 理想的正方形 单调队列

单调队列

by   GeneralLiu

滑动窗口是比较裸的单调队列

理想的正方形   就拔高了一个层次（多了一维）

有一个a*b的整数组成的矩阵

现请你从中找出一个n*n的正方形区域

使得该区域所有数中的最大值和最小值的差最小

只写MAX了，MIN一个道理，懒 不写了

先横着跑单调队列

　　维护许多长度为 n 的 横着的MAX

　　存在数组 map1[][] 中

再对数组 map1[][] 竖着跑单调队列

　　就维护了 许多 n*n 的 矩阵的MAX

MIN同理

在竖着跑单调队列时

顺便MAX 与 MIN 作差

更新答案即可

代码

一个小地方debug了好久

代码上还有 

学 mzx dalao  “ 二分出错位置  输出中间答案 ” 的 “ 蛛丝马迹 ”

#include<bits/stdc++.h>
using namespace std;
#define A 1005
#define N 105
int an=1000000001,h2,t2,h1,t1,a,b,n;
int pos1[A],pos2[A],que1[A],que2[A];
int mip1[A][A],map1[A][A];
int mip2[A][A],map2[A][A];
int read(){  //  读入优化 
    int ans=0;
    char ch=getchar();
    for(;!isdigit(ch);ch=getchar());
    for(;isdigit(ch);ch=getchar())
        ans=ans*10+ch-'0';
    return ans;
}
int main(){
    a=read(),b=read(),n=read();
    for(int k,i=1;i<=a;i++){
        h1=h2=1;
        t1=t2=0;
        for(int j=1;j<=b;j++){
            k=read();
            
            for(;t1>=h1&&k>=que1[t1];t1--);
            que1[++t1]=k;
            pos1[t1]=j;
            
            for(;t2>=h2&&k<=que2[t2];t2--);
            que2[++t2]=k;
            pos2[t2]=j;
            
            if(j>=n){
                map1[i][j]=que1[h1];
                mip1[i][j]=que2[h2];
                if(pos1[h1]==j-n+1)h1++;
                if(pos2[h2]==j-n+1)h2++;
            }
        }
    }
    /*cout<<endl;
    for(int i=1;i<=a;i++){
      for(int j=1;j<=b;j++)
        cout<<map1[i][j]<<" ";
      cout<<endl;
    }
    cout<<endl;
    for(int i=1;i<=a;i++){
      for(int j=1;j<=b;j++)
        cout<<mip1[i][j]<<" ";
      cout<<endl;
    }
    cout<<endl;*/
    for(int k1,k2,i=n;i<=b;i++){
        h1=h2=1;
        t1=t2=0;
        for(int j=1;j<=a;j++){
            k1=map1[j][i];
            k2=mip1[j][i];
            
            for(;t1>=h1&&k1>=que1[t1];t1--);
            que1[++t1]=k1;
            pos1[t1]=j;
            
            for(;t2>=h2&&k2<=que2[t2];t2--);
            que2[++t2]=k2;
            pos2[t2]=j;
            
            
            if(j>=n){
              //map2[j][i]=que1[h1];
              //mip2[j][i]=que2[h2];
              an=min(an,que1[h1]-que2[h2]);
              if(pos1[h1]==j-n+1)h1++;
              if(pos2[h2]==j-n+1)h2++;
            }
        }
    }
    
    /*for(int i=1;i<=a;i++){
      for(int j=1;j<=b;j++)
        cout<<map2[i][j]<<" ";
      cout<<endl;
    }
    cout<<endl;
    for(int i=1;i<=a;i++){
      for(int j=1;j<=b;j++)
        cout<<mip2[i][j]<<" ";
      cout<<endl;
    }
    cout<<endl;*/
    
    cout<<an;
    
    return 0;
}/*5 4 2 
1 2 5 6 
0 17 16 0 
16 17 0 1 
2 10 2 1 
1 2 3 2

2
*/