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

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

题面

bzoj

洛谷

题解

二维\(st\)表，代码是以前的

#include<iostream>
#include<cstdio>
#include<cstdlib>
#include<cstring>
#include<string>
#include<cmath>
#include<vector>
#include<queue>
#include<map>
#include<set>
#include<algorithm>
#include<ctime>
#define RG register
#define file(x) freopen(x".in","r",stdin);
using namespace std;

#define ll long long
inline int gi(){
    int data=0,w=1;
    char ch=0;
    while(ch!='-'&&(ch>'9'||ch<'0')) ch=getchar();
    if(ch=='-') w=-1,ch=getchar();
    while(ch>='0'&&ch<='9') data=data*10+ch-'0',ch=getchar();
    return w*data;
}

#define N 1010
int a,b,n;
int val[N][N],dp1[N][N][11],dp2[N][N][11],lg[N];
int MX(int e,int f,int g,int h){
	return max(max(e,f),max(g,h));
}
int MI(int e,int f,int g,int h){
	return min(min(e,f),min(g,h));
}
void Pre(){
	for(RG int i=1;i<=a;i++)
	  for(RG int j=1;j<=b;j++)
	    dp1[i][j][0]=dp2[i][j][0]=val[i][j];

	for(RG int k=1;k<=10;k++)
	  for(RG int i=1;i<=a-(1<<k-1);i++)
	    for(RG int j=1;j<=b-(1<<k-1);j++){
	    dp1[i][j][k]=MX (dp1[i][j][k-1],
	                     dp1[i+(1<<k-1)][j][k-1],
						 dp1[i][j+(1<<k-1)][k-1],
						 dp1[i+(1<<k-1)][j+(1<<k-1)][k-1]);

		dp2[i][j][k]=MI (dp2[i][j][k-1],
	                     dp2[i+(1<<k-1)][j][k-1],
						 dp2[i][j+(1<<k-1)][k-1],
						 dp2[i+(1<<k-1)][j+(1<<k-1)][k-1]);
	}
}
int query(int i,int j,int dis){
	int t=log(n)/log(2);
	int sum;
	sum=(MX(dp1[i][j][t],dp1[n+i-(1<<t)][j][t],dp1[i][j+n-(1<<t)][t],dp1[i+n-(1<<t)][j+n-(1<<t)][t])-
	     MI(dp2[i][j][t],dp2[i+n-(1<<t)][j][t],dp2[i][j+n-(1<<t)][t],dp2[i+n-(1<<t)][j+n-(1<<t)][t]));
	return sum;
}
int main(){
	a=gi();b=gi();n=gi();
	for(RG int i=1;i<=a;i++)
	  for(RG int j=1;j<=b;j++)
	    val[i][j]=gi();
	Pre();

	int MIN=0x7fffffff;;
	for(RG int i=1;i<=a-n+1;i++)
	   for(RG int j=1;j<=b-n+1;j++)
	       MIN=min(MIN,query(i,j,n));
    printf("%d\n",MIN);
	return 0;
}