题目大意: 有一个矩阵,每次操作前,矩阵每个格子内的权值都会 + 1 +1 +1,有三种操作:1、将某行权值归零;2、将某列权值归零;3、询问从 ( s x , s y ) (sx,sy) (sx,sy) 到 ( e x , e y ) (ex,ey) (ex,ey) 的最短路,要求不经过权值超过 k k k 的格子。
题解
显然,每次询问时,只有前 k k k 次的操作是有用的,于是只考虑前 k k k 次操作,然后分类讨论:
- ( s x , s y ) (sx,sy) (sx,sy) 和 ( e x , e y ) (ex,ey) (ex,ey),一个行被归零过,一个列被归零过,那么答案就是曼哈顿距离
- 行 s x sx sx ~ e x ex ex 都被归零过或列 s y sy sy ~ e y ey ey 都被归零过,答案也是曼哈顿距离
- ( s x , s y ) (sx,sy) (sx,sy) 和 ( e x , e y ) (ex,ey) (ex,ey) 两个行都被归零过,那么就是要找一个尽可能接近区间 [ s y , e y ] [sy,ey] [sy,ey] 的被归零过的列,假如恰好就在区间内,那么答案又是曼哈顿距离,否则就是曼哈顿距离加上多出来的长度
- 与情况 3 3 3 类似,考虑起点和终点的列都被归零过的情况,然后去找尽可能接近的行
于是维护两棵线段树,分别维护行和列,每一行或每一列的值表示它上一次被归零的时间,一开始都是 0 0 0。
然后需要维护的操作就是 区间求min
,区间求max
,找到第一个位置大于a且值大于b的位置
,找到第一个位置小于a且值大于b的位置
。
代码如下:
#include <cstdio>
#include <cstring>
#include <algorithm>
using namespace std;
#define inf 999999999
int n,m,q;
struct node{
int l,r,mid,ma,mi;node *zuo,*you;
node(int x,int y):l(x),r(y),mid(l+r>>1)
{if(x<y)zuo=new node(l,mid),you=new node(mid+1,r);}
void check(){ma=max(zuo->ma,you->ma);mi=min(zuo->mi,you->mi);}
void change(int x,int y){
if(l==r)return (void)(ma=mi=y);
if(x<=mid)zuo->change(x,y);else you->change(x,y);check();
}
int getmi(int x,int y)
{
if(l==x&&r==y)return mi;
if(y<=mid)return zuo->getmi(x,y);
else if(x>=mid+1)return you->getmi(x,y);
else return min(zuo->getmi(x,mid),you->getmi(mid+1,y));
}
int getma(int x,int y)
{
if(l==x&&r==y)return ma;
if(y<=mid)return zuo->getma(x,y);
else if(x>=mid+1)return you->getma(x,y);
else return max(zuo->getma(x,mid),you->getma(mid+1,y));
}
int find1(int x,int k)
{
if(ma<k)return -inf;if(l==r)return l;
if(x<=mid+1)return zuo->find1(x,k);
else if(getma(mid+1,x-1)>=k)return you->find1(x,k);
else return zuo->find1(mid+1,k);
}
int find2(int x,int k)
{
if(ma<k)return inf;if(l==r)return l;
if(x>=mid)return you->find2(x,k);
else if(getma(x+1,mid)>=k)return zuo->find2(x,k);
else return you->find2(mid,k);
}
}*line=NULL,*col=NULL;
inline char cn()
{
static char buf[1000010],*p1=buf,*p2=buf;
return p1==p2&&(p2=(p1=buf)+fread(buf,1,1000000,stdin),p1==p2)?EOF:*p1++;
}
template<class TY>void read(TY &x)
{
x=0;int f1=1;char ch=cn();
while(ch<'0'||ch>'9'){if(ch=='-')f1=-1;ch=cn();}
while(ch>='0'&&ch<='9')x=x*10+(ch-'0'),ch=cn(); x*=f1;
}
int main()
{
read(n);read(m);read(q);line=new node(1,n);col=new node(1,m);
for(int i=1,id,x,y,xx,yy,k;i<=q;i++)
{
read(id);
if(id==1)read(x),line->change(x,i);
if(id==2)read(x),col->change(x,i);
if(id==3)
{
read(x),read(y),read(xx),read(yy),read(k);k=i-k;
int X=line->getma(x,x),Y=col->getma(y,y),XX=line->getma(xx,xx),YY=col->getma(yy,yy);
#define manhattan {printf("%d\n",abs(yy-y)+abs(xx-x));continue;}
#define manhattan2 {printf("%d\n",abs(yy-y)+abs(xx-x)+2*ans);continue;}
if((X>=k&&YY>=k)||(Y>=k&&XX>=k))manhattan;//一个行,一个列
if(x>xx)swap(x,xx);if(y>yy)swap(y,yy);
if(line->getmi(x,xx)>=k||col->getmi(y,yy)>=k)manhattan;//一段行/一段列
if(X>=k&&XX>=k){//两个行
if(col->getma(y,yy)>=k)manhattan;
int ans=min(y-col->find1(y,k),col->find2(yy,k)-yy);
if(ans<n+m)manhattan2;
}
if(Y>=k&&YY>=k){//两个列
if(line->getma(x,xx)>=k)manhattan;
int ans=min(x-line->find1(x,k),line->find2(xx,k)-xx);
if(ans<n+m)manhattan2;
}
printf("-1\n");
}
}
}