【jzoj5049】【腐女的生日】【扫描线】【数据结构】

题目大意

腐女要过生日了，pty 想给腐女送礼物，但是腐女所在的教室离pty 的教室太远了，于是pty就拜托会动归和A星的djy帮忙送礼物。djy在学校建立了一个平面直角坐标系，他站在了（0，0）点，腐女在（x0,y0）点，djy每次只能往上下左右四个方向移动一步，中间有n栋矩形教学楼，每个教学楼给出两个对角的坐标，并且保证每栋教学楼的周围区域（如图所示）不会有别的教学楼，即djy可以绕一个教学楼走不会碰到任何障碍，现在djy 想知道从起点到终点不碰到任何教学楼，最短需要多少步。

解题思路

考虑扫描线，从左到右不会折返，用一颗线段树记录到达扫描线不同高度的纵坐标花费，出现矩阵把对应区间赋为正无穷，消失矩阵重新更新，发现只能从矩阵两端走过来，是两个一次函数，找出分界点用线段树维护即可。

code

#include<set>
#include<cstdio>
#include<cstring>
#include<algorithm>
#define LD double
#define LL long long
#define ULL unsigned long long
#define min(a,b) ((a<b)?a:b)
#define max(a,b) ((a>b)?a:b)
#define fo(i,j,k) for(int i=j;i<=k;i++)
#define fd(i,j,k) for(int i=j;i>=k;i--)
#define fr(i,j) for(int i=begin[j];i;i=next[i])
using namespace std;
int const mn=2*1e5+9,mp=8*1e6+9,con=1e5+1,inf=1e9;
int n,x0,y0,tag[mp][2];
struct rec{
    int x,l,r,op;
};
rec a[mn];
bool cmp(rec x,rec y){
    return (x.x<y.x)||((x.x==y.x)&&(x.op>y.op));
}
void retag(int now){
    int tmp=(tag[now][0]<tag[now][1])?1:-1;
    tag[now*2][0]=tag[now][0];
    tag[now*2][1]=(tag[now][0]+tag[now][1]+(tmp<0))/2;
    tag[now*2+1][0]=(tag[now][0]+tag[now][1]+(tmp<0))/2+tmp;
    tag[now*2+1][1]=tag[now][1];
    tag[now][0]=tag[now][1]=0;
}
void oper(int now,int l,int r,int u,int uu,int v,int vv){
    int mid=(l+r)>>1;
    if((l!=r)&&(tag[now][0]||tag[now][1]))retag(now);
    if((l==u)&&(r==v)){
        tag[now][0]=uu;tag[now][1]=vv;
        return;
    }
    if(v<=mid)oper(now*2,l,mid,u,uu,v,vv);
    else if(mid<u)oper(now*2+1,mid+1,r,u,uu,v,vv);
    else{
        oper(now*2,l,mid,u,uu,mid,((vv>uu)?1:-1)*(mid-u)+uu);
        oper(now*2+1,mid+1,r,mid+1,((vv>uu)?1:-1)*(mid-u+1)+uu,v,vv);
    }
}
int qury(int now,int l,int r,int v){
    int mid=(l+r)>>1;
    if(l==r)return tag[now][0];
    if((l!=r)&&(tag[now][0]||tag[now][1]))retag(now);
    if(v<=mid)return qury(now*2,l,mid,v);
    return qury(now*2+1,mid+1,r,v);
}
int main(){
    freopen("bl.in","r",stdin);
    freopen("bl.out","w",stdout);
    scanf("%d%d",&x0,&y0);
    scanf("%d",&n);int x1,y1,x2,y2,cnt=0;
    fo(i,1,n){
        scanf("%d%d%d%d",&x1,&y1,&x2,&y2);
        a[++cnt].x=x1,a[cnt].l=y1,a[cnt].r=y2,a[cnt].op=1;
        a[++cnt].x=x2,a[cnt].l=y1,a[cnt].r=y2,a[cnt].op=-1;
    }
    sort(a+1,a+cnt+1,cmp);
    oper(1,-con,con,-con,con,0,0);
    oper(1,-con,con,1,1,con,con);
    fo(i,1,cnt){
        if(a[i].x>x0)break;
        if(a[i].op==-1){
            int tmp1=qury(1,-con,con,a[i].l-1),tmp2=qury(1,-con,con,a[i].r+1),
                tmp=(a[i].l+a[i].r+tmp2-tmp1)/2;
            oper(1,-con,con,a[i].l-1,tmp1,tmp,tmp1+tmp-a[i].l+1);
            if(tmp+1<a[i].r+1)oper(1,-con,con,tmp+1,tmp2+a[i].r-tmp,a[i].r+1,tmp2);
        }
    }
    printf("%d",qury(1,-con,con,y0)+x0);
    return 0;
}