记录编号 239444 评测结果 AAAAAAAAAA
题目名称 [HAOI 2007]理想的正方形 最终得分 100
用户昵称 Gravatarliu_runda 是否通过 通过
代码语言 C++ 运行时间 1.173 s
提交时间 2016-03-20 11:19:42 内存使用 17.60 MiB
显示代码纯文本
#include<cstdio>
#include<cstring>
const int maxn=1005;
int m1[maxn][maxn];//原始矩阵 
int m2a[maxn][maxn];//对于每一行,单调队列求每个长为n区间的最小值 以区间右端点代表区间 
int m2b[maxn][maxn];//对于每一行,单调队列求每个长为n区间的最大值 
int m3a[maxn][maxn],m3b[maxn][maxn];
int a,b,n;
bool high(int a,int b){
	return a>=b;
}
bool low(int a,int b){
	return a<=b;
}
struct queue{
	int q[1005],index[1005];
	int head,tail;
	void clear(){
		memset(q,0,sizeof(q));
		head=tail=0;
	}
	void add(int x,int i,bool high(int,int)){
		int lim=i+1-n;
		while(head<tail&&index[head]<lim)head++;
		while(head<tail&&high(q[tail-1],x))tail--;
		index[tail]=i;
		q[tail++]=x;
	}
	int min(){
		return q[head];
	}
}q1;
int main(){
	freopen("square.in","r",stdin);
	freopen("square.out","w",stdout);
	scanf("%d%d%d",&a,&b,&n);
	for(int i=1;i<=a;++i)
		for(int j=1;j<=b;++j)
			scanf("%d",&m1[i][j]);
	for(int i=1;i<=a;++i){//获取最小值 
		q1.clear();
		for(int j=1;j<n;++j)q1.add(m1[i][j],j,high);
		for(int j=n;j<=b;++j){
			q1.add(m1[i][j],j,high);
			m2a[i][j]=q1.min();
		}
	}
	for(int i=1;i<=a;++i){//获取最大值 
		q1.clear();
		for(int j=1;j<n;++j)q1.add(m1[i][j],j,low);
		for(int j=n;j<=b;++j){
			q1.add(m1[i][j],j,low);
			m2b[i][j]=q1.min();
		}
	}
	for(int i=n;i<=b;++i){
		q1.clear();
		for(int j=1;j<n;++j)q1.add(m2a[j][i],j,high);
		for(int j=n;j<=a;++j){
			q1.add(m2a[j][i],j,high);
			m3a[j][i]=q1.min();
		}
	}
	for(int i=n;i<=b;++i){
		q1.clear();
		for(int j=1;j<n;++j)q1.add(m2b[j][i],j,low);
		for(int j=n;j<=a;++j){
			q1.add(m2b[j][i],j,low);
			m3b[j][i]=q1.min();
		}
	}
	int ans=m3b[a][b]-m3a[a][b];
	for(int i=n;i<=a;++i){
		for(int j=n;j<=b;++j){
			if(m3b[i][j]-m3a[i][j]<ans)ans=m3b[i][j]-m3a[i][j];
		}
	}
	printf("%d",ans);
	fclose(stdin);fclose(stdout);
	return 0;
}