1. 程式人生 > >BZOJ1047 [HAOI2007]理想的正方形

BZOJ1047 [HAOI2007]理想的正方形

com mat strong light ade main 延伸 tro end

本文版權歸ljh2000和博客園共有,歡迎轉載,但須保留此聲明,並給出原文鏈接,謝謝合作。

本文作者:ljh2000
作者博客:http://www.cnblogs.com/ljh2000-jump/
轉載請註明出處,侵權必究,保留最終解釋權!

題目鏈接:BZOJ1047

正解:單調隊列

解題報告:

  考慮以每個點為左下角的一個正方形,那麽我們需要支持快速查詢一個矩形內的最大值、最小值,暴力做似乎必須帶$log$。

  不難想到算出每個點往右延伸$n$個位置的最大值、最小值很好算,那麽對於同一列維護一個單調隊列,這樣就能快速得到一個矩陣的最大值、最小值了。

  貌似這玩意可以叫“二維單調隊列”?好吧,其實就是把二維拍扁了變成了一維...

//It is made by ljh2000
//有誌者,事竟成,破釜沈舟,百二秦關終屬楚;苦心人,天不負,臥薪嘗膽,三千越甲可吞吳。
#include <algorithm>
#include <iostream>
#include <cstring>
#include <vector>
#include <cstdio>
#include <string>
#include <queue>
#include <cmath>
using namespace std;
typedef long long LL;
const int MAXN = 1011;
int n,m,L,Max[MAXN][MAXN],Min[MAXN][MAXN];
int a[MAXN][MAXN],ans;
int l1,r1,l2,r2,d1[MAXN],d2[MAXN];
int L1[MAXN],R1[MAXN],L2[MAXN],R2[MAXN],D1[MAXN][MAXN],D2[MAXN][MAXN];

inline int getint(){
    int w=0,q=0; char c=getchar(); while((c<‘0‘||c>‘9‘) && c!=‘-‘) c=getchar();
    if(c==‘-‘) q=1,c=getchar(); while (c>=‘0‘&&c<=‘9‘) w=w*10+c-‘0‘,c=getchar(); return q?-w:w;
}

inline void work(){
	n=getint(); m=getint(); L=getint(); for(int i=1;i<=n;i++) for(int j=1;j<=m;j++) a[i][j]=getint();
	ans=((1<<30)-1)<<1|1;
	for(int i=1;i<=m;i++) L1[i]=L2[i]=1,R1[i]=R2[i]=0;//!!!
	for(int i=1;i<=n;i++) {
		l1=l2=1; r1=r2=0;
		for(int j=m;j>=1;j--) {
			Max[i][j]=Min[i][j]=a[i][j];
			while(l1<=r1 && d1[l1]-j>=L) l1++;
			while(l2<=r2 && d2[l2]-j>=L) l2++;

			while(r1>=l1 && a[i][j]>a[i][d1[r1]]) r1--;
			d1[++r1]=j;
			if(l1<=r1) Max[i][j]=max(Max[i][j],a[i][d1[l1]]);

			while(r2>=l2 && a[i][j]<a[i][d2[r2]]) r2--;
			d2[++r2]=j;
			if(l2<=r2) Min[i][j]=min(Min[i][j],a[i][d2[l2]]);
		}

		for(int j=m-L+1;j>=1;j--) {
			while(L1[j]<=R1[j] && i-D1[j][L1[j]]>=L) L1[j]++;
			while(L2[j]<=R2[j] && i-D2[j][L2[j]]>=L) L2[j]++;
			while(L1[j]<=R1[j] && Max[i][j]>Max[ D1[j][R1[j]] ][j]) R1[j]--;
			while(L2[j]<=R2[j] && Min[i][j]<Min[ D2[j][R2[j]] ][j]) R2[j]--;
			D1[j][++R1[j]]=i; D2[j][++R2[j]]=i;
			if(i>=L) {//!!!
				if(L1[j]<=R1[j] && L2[j]<=R2[j])
					ans=min(ans,Max[ D1[j][L1[j]] ][j]-Min[ D2[j][L2[j]] ][j]);
			}
		}
	}
	printf("%d",ans);
}

int main()
{
#ifndef ONLINE_JUDGE
	freopen("1047.in","r",stdin);
	freopen("1047.out","w",stdout);
#endif
    work();
    return 0;
}
//有誌者,事竟成,破釜沈舟,百二秦關終屬楚;苦心人,天不負,臥薪嘗膽,三千越甲可吞吳。

  

BZOJ1047 [HAOI2007]理想的正方形