BZOJ 1047 理想的正方形

可以先合一维再合一维

两个方向单调队列分别搞一下即可

蒟蒻我居然想不到单调队列

#include <cstdio>

const int MAXN=1011;
const int MAXM=1011;
const int INF=1034567890;

int N, M, K;
int Map[MAXN][MAXM], Temp[MAXN][MAXM], Max[MAXN][MAXM], Min[MAXN][MAXM];

int Q[MAXN+MAXM], Head, Tail;

int main(){
	
	scanf("%d%d%d", &N, &M, &K);
	for(int i=1;i<=N;++i)
		for(int j=1;j<=M;++j)
			scanf("%d", &Map[i][j]);
	
	for(int i=1;i<=N;++i){
		Head=1;Tail=1;
		for(int j=1;j<=M;++j){
			while(Head<Tail && Map[i][Q[Tail-1]]<=Map[i][j])	--Tail;
			Q[Tail++]=j;
			while(Q[Head]<=j-K)	++Head;
			Temp[i][j]=Map[i][Q[Head]];
		}
	}
	
	for(int j=1;j<=M;++j){
		Head=1;Tail=1;
		for(int i=1;i<=N;++i){
			while(Head<Tail && Temp[Q[Tail-1]][j]<=Temp[i][j])	--Tail;
			Q[Tail++]=i;
			while(Q[Head]<=i-K)	++Head;
			Max[i][j]=Temp[Q[Head]][j];
		}
	}
	
	for(int i=1;i<=N;++i){
		Head=1;Tail=1;
		for(int j=1;j<=M;++j){
			while(Head<Tail && Map[i][Q[Tail-1]]>=Map[i][j])	--Tail;
			Q[Tail++]=j;
			while(Q[Head]<=j-K)	++Head;
			Temp[i][j]=Map[i][Q[Head]];
		}
	}
	
	for(int j=1;j<=M;++j){
		Head=1;Tail=1;
		for(int i=1;i<=N;++i){
			while(Head<Tail && Temp[Q[Tail-1]][j]>=Temp[i][j])	--Tail;
			Q[Tail++]=i;
			while(Q[Head]<=i-K)	++Head;
			Min[i][j]=Temp[Q[Head]][j];
		}
	}
	
	for(int i=1;i<=N;++i)
		for(int j=1;j<=M;++j)
			Temp[i][j]=Max[i][j]-Min[i][j];
	
	int ANS=INF;
	
	for(int i=K;i<=N;++i)
		for(int j=K;j<=M;++j)
			if(Temp[i][j]<ANS)	ANS=Temp[i][j];
	
	printf("%d\n", ANS);
	
	return 0;
}

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

1

*/