洛谷.2219.[HAOI2007]修筑绿化带(单调队列)

题目链接 洛谷

COGS.24

对于大的矩阵可以枚举；对于小的矩阵，需要在满足条件的区域求一个矩形和的最小值

预处理S2[i][j]表示以(i,j)为右下角的C\(*\)D的矩阵和，

然后对于求矩形区域的最小值，可以先将每行看做一个数列，对于每个点y，得到一个[y-(B-3),y]的最小值

处理完行后得到Minr[][]，对每列的进行同样的操作，就可以得到Min[x][y]表示([x-A+3,x],[y-B+3,y])的最小矩形和

但是注意单调队列处理的是S2，S2表示的是C\(*\)D的和,not a single point!所以端点应该是([x-A+C+1],[y-B+D+1])

原先做过这样的套路题，还是费了近一下午==不行效率太低了

#include <cstdio>
#include <cctype>
//#define gc() getchar()
#define gc() (SS==TT&&(TT=(SS=IN)+fread(IN,1,MAXIN,stdin),SS==TT)?EOF:*SS++)
const int N=1005,MAXIN=6e5;

int n,m,A,B,C,D,S

,S1

,S2

,q
,Minr

;
char IN[MAXIN],*SS=IN,*TT=IN;

inline int read()
{
int now=0,f=1;register char c=gc();
for(;!isdigit(c);c=gc()) if(c=='-') f=-1;
for(;isdigit(c);now=now*10+c-'0',c=gc());
return now*f;
}
//void Debug(int a

,int s,int s2)
//{
//  puts("Debug");
//  for(int i=s; i<=n; ++i,putchar('\n'))
//      for(int j=s2; j<=m; ++j) printf("%d ",a[i][j]);
//  putchar('\n');
//}
void Pre()
{
for(int i=C+1; i<n; ++i)//第1~C,n行用不到
{
int h=1,t=0;
for(int j=D+1; j<m; ++j)
{
while(h<=t && S2[i][q[t]]>=S2[i][j]) --t;
q[++t]=j;
if(q[h]<=j-B+D+1) ++h;
Minr[i][j]=S2[i][q[h]];
}
}
//  Debug(S,1,1);
//  Debug(S1,A,B);
//  Debug(S2,C,D);
//  Debug(Minr,2,2);
}
void Solve()
{
int res=0;
for(int j=D+1; j<m; ++j)//枚举小矩形可放范围，行+1列+1即大矩形的位置
{
int h=1,t=0;
for(int i=C+1; i<n; ++i)
{
while(h<=t && Minr[q[t]][j]>=Minr[i][j]) --t;
q[++t]=i;
if(q[h]<=i-A+C+1) ++h;
if(res<S1[i+1][j+1]-Minr[q[h]][j]) res=S1[i+1][j+1]-Minr[q[h]][j];
}
}
printf("%d",res);
}

int main()
{
n=read(),m=read(),A=read(),B=read(),C=read(),D=read();
if(n<=2||m<=2) {putchar('0'); return 0;}
for(int i=1; i<=n; ++i)
for(int j=1; j<=m; ++j)
S[i][j]=read()+S[i][j-1]+S[i-1][j]-S[i-1][j-1];
for(int i=A; i<=n; ++i)
for(int j=B; j<=m; ++j)
S1[i][j]=S[i][j]-S[i][j-B]-S[i-A][j]+S[i-A][j-B];
for(int i=C+1; i<=n; ++i)
for(int j=D+1; j<=m; ++j)
S2[i][j]=S[i][j]-S[i][j-D]-S[i-C][j]+S[i-C][j-D];
Pre(), Solve();

return 0;
}