记录编号 |
239444 |
评测结果 |
AAAAAAAAAA |
题目名称 |
[HAOI 2007]理想的正方形 |
最终得分 |
100 |
用户昵称 |
liu_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;
}