二维rmq

一维的推广，本来想一个小时之内ac，结果先是开四维，re，再是重编三维，少考虑一种情况，wa，终于ac还是2小时之后的事。

题目数据很弱，暴力也可以过，题目中b是确定的，所以四维中表示矩阵大小的两维完全可以用一维表示。

预处理：

f[i,j,k]:=min{f[i,j,k-1],f[i+d[k-1],j,k-1],f[i,j+d[k-1],k-1],f[i+d[k-1],j+d[k-1],k-1]}

询问：

si:=i+b-1;sj:=j+b-1;

ans:=min{f[i,j,k],f[si-d[k]+1,j,k],f[i,sj-d[k]+1,k],f[si-d[k]+1,sj-d[k]+1,k]}

var d:array[0..10]of longint; n,b,p,l:longint; f,g:array[1..400,1..400,0..8]of longint;function min(x,y:longint):longint;begin if x<y then exit(x) else exit(y)end;function max(x,y:longint):longint;begin if x>y then exit(x) else exit(y)end;procedure init;var i,j,k,sx,sy,x,y,ans1,ans2:longint;begin readln(n,b,p); fillchar(f,sizeof(f),127); for i:=1 to n do begin for j:=1 to n do begin read(f[i,j,0]);g[i,j,0]:=f[i,j,0] end; readln end; d[0]:=1;l:=trunc(ln(b)/ln(2)); for i:=1 to 8 do d[i]:=d[i-1]<<1; for k:=1 to l do for i:=1 to n do for j:=1 to n do begin f[i,j,k]:=min(f[i,j,k-1],f[i+d[k-1],j,k-1]); f[i,j,k]:=min(f[i,j+d[k-1],k-1],f[i,j,k]); f[i,j,k]:=min(f[i+d[k-1],j+d[k-1],k-1],f[i,j,k]); g[i,j,k]:=max(g[i,j,k-1],g[i+d[k-1],j,k-1]); g[i,j,k]:=max(g[i,j+d[k-1],k-1],g[i,j,k]); g[i,j,k]:=max(g[i+d[k-1],j+d[k-1],k-1],g[i,j,k]); end; for i:=1 to p do begin readln(x,y);sx:=x+b-1;sy:=y+b-1; ans1:=max(g[x,y,k],g[x,sy-d[k]+1,k]);ans1:=max(ans1,g[sx-d[k]+1,y,k]); ans1:=max(ans1,g[sx-d[k]+1,sy-d[k]+1,k]); ans2:=min(f[x,y,k],f[x,sy-d[k]+1,k]);ans2:=min(ans2,f[sx-d[k]+1,y,k]); ans2:=min(ans2,f[sx-d[k]+1,sy-d[k]+1,k]); writeln(ans1-ans2) endend;beginassign(input,'2019.in');reset(input); init;close(input)end.