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<cstdio>
#include<iostream>
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<<lg[c])+1][lg[c]],f1[x2-(1<<lg[c])+1][y1][lg[c]],f1[x2-(1<<lg[c])+1][y2-(1<<lg[c])+1][lg[c]]);
}

inline int get2(int x1,int y1,int x2,int y2)
{
	return Min(f2[x1][y1][lg[c]],f2[x1][y2-(1<<lg[c])+1][lg[c]],f2[x2-(1<<lg[c])+1][y1][lg[c]],f2[x2-(1<<lg[c])+1][y2-(1<<lg[c])+1][lg[c]]);
}
int main()
{
	n=read(),m=read(),c=read();
	lg[0]=-1;
	for(int i=1;i<=c;i++)
		lg[i]=lg[i>>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<<k)+1;i++)   //大意在这儿没减（1<<k）
			for(int j=1;j<=m-(1<<k)+1;j++)
				f1[i][j][k]=Max(f1[i][j][k-1],f1[i][j+(1<<(k-1))][k-1],f1[i+(1<<(k-1))][j][k-1],f1[i+(1<<(k-1))][j+(1<<(k-1))][k-1]),
				f2[i][j][k]=Min(f2[i][j][k-1],f2[i][j+(1<<(k-1))][k-1],f2[i+(1<<(k-1))][j][k-1],f2[i+(1<<(k-1))][j+(1<<(k-1))][k-1]);
	ans=1e9;
	for(int i=1;i<=n-c+1;i++)
		for(int j=1;j<=m-c+1;j++)
			ans=min(ans,get1(i,j,i+c-1,j+c-1)-get2(i,j,i+c-1,j+c-1));
	printf("%d\n",ans);
	return 0;
}

解二：单调队列

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

#include<cstdio>
#include<iostream>
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]<a[i][q[r]]) r--;
            q[++r]=j;
            if(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]<m2[q2[r2]][j]) r2--;
            q2[++r2]=i;
            if(i>=c) ans=min(ans,m1[q1[l1]][j]-m2[q2[l2]][j]);
        }
    }
    
    printf("%d\n",ans);
    return 0;
}