[BZOJ1047][HAOI2007]理想的正方形

#include <cstdio>
#include <algorithm>
#include <cmath>
using namespace std;
const int N=1010;
const int LOGN=11;
int a,b,n,ans=2147483647;
int lg2
,two[LOGN]={1};
int M

;
int fmi

[LOGN],fmx

[LOGN];

int main(){
scanf("%d%d%d",&a,&b,&n);
for(int i=1;i<=n;i++) lg2[i]=(int)log2(i);
for(int i=1;i<LOGN;i++) two[i]=two[i-1]<<1;
for(int i=1;i<=a;i++){
for(int j=1;j<=b;j++){
scanf("%d",&M[i][j]);
fmi[i][j][0]=fmx[i][j][0]=M[i][j];
}
}
for(int i=1;i<LOGN;i++){
for(int j=1;j+two[i-1]<=a;j++){
for(int k=1;k+two[i-1]<=b;k++){
fmi[j][k][i]=fmi[j][k][i-1];
fmi[j][k][i]=min(fmi[j][k][i],fmi[j+two[i-1]][k][i-1]);
fmi[j][k][i]=min(fmi[j][k][i],fmi[j][k+two[i-1]][i-1]);
fmi[j][k][i]=min(fmi[j][k][i],fmi[j+two[i-1]][k+two[i-1]][i-1]);
fmx[j][k][i]=fmx[j][k][i-1];
fmx[j][k][i]=max(fmx[j][k][i],fmx[j+two[i-1]][k][i-1]);
fmx[j][k][i]=max(fmx[j][k][i],fmx[j][k+two[i-1]][i-1]);
fmx[j][k][i]=max(fmx[j][k][i],fmx[j+two[i-1]][k+two[i-1]][i-1]);
}
}
}
int d=two[lg2
];
for(int i=1;i+n-1<=a;i++){
for(int j=1;j+n-1<=b;j++){
int mi=fmi[i][j][lg2
],mx=fmx[i][j][lg2
];
mi=min(mi,min(fmi[i+n-d][j][lg2
],min(fmi[i][j+n-d][lg2
],fmi[i+n-d][j+n-d][lg2
])));
mx=max(mx,max(fmx[i+n-d][j][lg2
],max(fmx[i][j+n-d][lg2
],fmx[i+n-d][j+n-d][lg2
])));
ans=min(ans,mx-mi);
}
}
printf("%d\n",ans);

return 0;
}