BZOJ 1047: [HAOI2007]理想的正方形

Description

　　有一个a*b的整数组成的矩阵，现请你从中找出一个n*n的正方形区域，使得该区域所有数中的最大值和最小值
的差最小。

Input

　　第一行为3个整数，分别表示a,b,n的值第二行至第a+1行每行为b个非负整数，表示矩阵中相应位置上的数。每
行相邻两数之间用一空格分隔。
100%的数据2<=a,b<=1000,n<=a,n<=b,n<=1000

Output

　　仅一个整数，为a*b矩阵中所有“n*n正方形区域中的最大整数和最小整数的差值”的最小值。

Sample Input

5 4 2
1 2 5 6
0 17 16 0
16 17 2 1
2 10 2 1
1 2 2 2

Sample Output

1

 

根据伟大的大佬zyh的话

一题要多解

这题可以多解

解一：二维ST表

n*m*c，卡卡能过

（代码巨丑，横向巨长，不爽）

#include
#include
using namespace std;
int read()
{
	int ret=0;
	char ch=getchar();
	while(ch<'0'||ch>'9') ch=getchar();
	while(ch>='0'&&ch<='9')
		ret=(ret<<1)+(ret<<3)+ch-'0',ch=getchar();
	return ret;
}

int n,m,c,ans;
const int N=1005,C=105;
int lg[N],f1[N][N][20],f2[N][N][20];

inline int Max(int x1,int x2,int x3,int x4)
{
	return max(x1,max(x2,max(x3,x4)));
}

inline int Min(int x1,int x2,int x3,int x4)
{
	return min(x1,min(x2,min(x3,x4)));
}

inline int get1(int x1,int y1,int x2,int y2)
{
	return Max(f1[x1][y1][lg[c]],f1[x1][y2-(1<>1]+1;
	for(int i=1;i<=n;i++)
		for(int j=1;j<=m;j++)
			f1[i][j][0]=f2[i][j][0]=read();
	for(int k=1;k<=lg[c];k++)
		for(int i=1;i<=n-(1<

解二：单调队列

代码更长，不过可以复制，打一部分就好

#include
#include
using namespace std;

int n,m,c,l,r,l1,r1,l2,r2,ans;
const int N=1005;
int a[N][N],q[N],q1[N],q2[N],m1[N][N],m2[N][N];

inline int read()
{
    int ret=0;
    char ch=getchar();
    while(ch<'0'||ch>'9') ch=getchar();
    while(ch>='0'&&ch<='9')
        ret=(ret<<1)+(ret<<3)+ch-'0',ch=getchar();
    return ret;
}

int main()
{
    n=read(),m=read(),c=read();
    for(int i=1;i<=n;i++)
        for(int j=1;j<=m;j++) 
            a[i][j]=read(); 
    
    for(int i=1;i<=n;i++)
    {
        l=1;r=0;
        for(int j=1;j<=m;j++)
        {
            while(l<=r&&j-q[l]+1>c) l++;
            while(l<=r&&a[i][j]>a[i][q[r]]) r--;
            q[++r]=j;
            if(j>=c)m1[i][j]=a[i][q[l]];
        }
        
        l=1;r=0;
        for(int j=1;j<=m;j++)
        {
            while(l<=r&&j-q[l]+1>c) l++;
            while(l<=r&&a[i][j]=c)m2[i][j]=a[i][q[l]];
        }
    }
    ans=1e9;
    for(int j=c;j<=m;j++)
    {
        l1=1;r1=0; l2=1,r2=0; 
        for(int i=1;i<=n;i++)
        {
            while(l1<=r1&&i-q1[l1]+1>c) l1++;
            while(l1<=r1&&m1[i][j]>m1[q1[r1]][j]) r1--;
            q1[++r1]=i;
            while(l2<=r2&&i-q2[l2]+1>c) l2++;
            while(l2<=r2&&m2[i][j]=c) ans=min(ans,m1[q1[l1]][j]-m2[q2[l2]][j]);
        }
    }
    
    printf("%d\n",ans);
    return 0;
}