腐女的生日

题目描述

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

扫描线

除了第一步会左走
其余情况均不会左走。
用扫描线保存当前x坐标下到每点的最短路。
遇到矩形边界时讨论一下。

#include<cstdio>
#include<algorithm>
#define fo(i,a,b) for(i=a;i<=b;i++)
using namespace std;
const int maxn=100000+10,mx=1000002,maxm=2000003;
int sx[maxm*4+100],gc[maxm*4+100];
bool bz[maxm*4+100],pd[maxm*4+100];
struct dong{
    int x,ca,y0,y1;
} ask[maxn*2];
int i,j,k,l,r,t,n,m,tot,top,x,x0,y0,x1,y1,x2,y2,ans;
int read(){
    int x=0,f=1;
    char ch=getchar();
    while (ch<'0'||ch>'9'){
        if (ch=='-') f=-1;
        ch=getchar();
    }
    while (ch>='0'&&ch<='9'){
        x=x*10+ch-'0';
        ch=getchar();
    }
    return x*f;
}
bool cmp(dong a,dong b){
    return a.x<b.x||a.x==b.x&&a.ca>b.ca;
}
void markad(int p,int v,int d){
    if (pd[p]) return;
    sx[p]+=v;
    gc[p]+=d;
}
void markcl(int p){
    bz[p]=1;
    sx[p]=gc[p]=0;
}
void down(int p,int l,int r){
    int mid=(l+r)/2;
    if (bz[p]){
        markcl(p*2);
        markcl(p*2+1);
        bz[p]=0;
    }
    if (sx[p]||gc[p]){
        markad(p*2,sx[p],gc[p]);
        markad(p*2+1,sx[p]+(mid-l+1)*gc[p],gc[p]);
        sx[p]=gc[p]=0;
    }
}
void change(int p,int l,int r,int a,int b,int v,int d){
    if (pd[p]||a>b) return;
    if (l==a&&r==b){
        markad(p,v,d);
        return;
    }
    down(p,l,r);
    int mid=(l+r)/2;
    if (b<=mid) change(p*2,l,mid,a,b,v,d);
    else if (a>mid) change(p*2+1,mid+1,r,a,b,v,d);
    else{
        change(p*2,l,mid,a,mid,v,d);
        change(p*2+1,mid+1,r,mid+1,b,v+d*(mid-a+1),d);
    }
}
void clear(int p,int l,int r,int a,int b,int v){
    if (l==a&&r==b){
        pd[p]=v;
        if (v) markcl(p);
        return;
    }
    down(p,l,r);
    int mid=(l+r)/2;
    if (b<=mid) clear(p*2,l,mid,a,b,v);
    else if (a>mid) clear(p*2+1,mid+1,r,a,b,v);
    else{
        clear(p*2,l,mid,a,mid,v);
        clear(p*2+1,mid+1,r,mid+1,b,v);
    }
}
int query(int p,int l,int r,int a){
    if (l==r) return sx[p];
    down(p,l,r);
    int mid=(l+r)/2;
    if (a<=mid) return query(p*2,l,mid,a);else return query(p*2+1,mid+1,r,a);
}
int main(){
    freopen("bl.in","r",stdin);freopen("bl.out","w",stdout);
    x0=read();y0=read();
    ask[top=1].x=x0;ask[1].y0=y0;
    n=read();
    fo(i,1,n){
        x1=read();y1=read();x2=read();y2=read();
        if (x1>x2){
            swap(x1,x2);
            swap(y1,y2);
        }
        if (y1>y2) swap(y1,y2);
        ask[++top].x=x1;
        ask[top].ca=1;
        ask[top].y0=y1;ask[top].y1=y2;
        ask[++top].x=x2+1;
        ask[top].ca=2;
        ask[top].y0=y1;ask[top].y1=y2;
    }
    sort(ask+1,ask+top+1,cmp);
    x=0;
    change(1,1,maxm,mx+1,maxm,1,1);
    change(1,1,maxm,1,mx-1,mx-1,-1);
    fo(i,1,top){
        change(1,1,maxm,1,maxm,ask[i].x-x,0);
        x=ask[i].x;
        if (ask[i].ca==0){
            ans=query(1,1,maxm,ask[i].y0+mx);
            break;
        }
        else if (ask[i].ca==1){
            clear(1,1,maxm,ask[i].y0+mx,ask[i].y1+mx,1);
        }
        else if (ask[i].ca==2){
            clear(1,1,maxm,ask[i].y0+mx,ask[i].y1+mx,0);
            j=query(1,1,maxm,ask[i].y0-1+mx);
            k=query(1,1,maxm,ask[i].y1+1+mx);
            l=ask[i].y1-ask[i].y0+1;
            t=(k+l+1-j)/2;
            if (t<1) change(1,1,maxm,ask[i].y0+mx,ask[i].y1+mx,k+l,-1);
            else if (t>l) change(1,1,maxm,ask[i].y0+mx,ask[i].y1+mx,j+1,1);
            else{
                change(1,1,maxm,ask[i].y0+mx,ask[i].y0+t-1+mx,j+1,1);
                change(1,1,maxm,ask[i].y0+t+mx,ask[i].y1+mx,k+l-t,-1);
            }
        }
    }
    printf("%d\n",ans);
}