bzoj 1949 [Ceoi2006]Walk dp 线段树

woc这题做得蛋疼。。。
一开始自己yy了一下，发现一个矩形到另一个矩形的转移一定是从一个的左上/左下角走到另一个的左上/左下角。如图：

然后设$f[i][0/1]$ 表示原点到第i个矩形左上/左下角的最短路。
然后按左边界排序，从右往左扫一遍维护每个矩形的左上/左下角向右到达的第一个矩形。
然后再更新一下可以直接从原点走到的矩形的f值。
但问题是无法确定最后一个矩形的左上/左下角是怎么走到终点的。
自己搞了几个结论然后写得整个人都不好了。
然后看了一下题解发现只需要从终点往起点走就好了。
因为所有矩形都在原点右侧因此就没有那么多特判了。。。。。。。

#include <bits/stdc++.h>
using namespace std;
#define N 110000
#define M 2100000
#define A 1000001
#define ll long long
#define ls l,mid,now<<1
#define rs mid+1,r,now<<1|1
int X,Y,n;
int bj[M<<2],tar[N][2];
ll f[N][2],ans;
void upd(ll &x,ll y){x=min(x,y);}
struct node
{
    int x1,x2,y1,y2;
    void read()
    {
        scanf("%d%d%d%d",&x1,&y1,&x2,&y2);
        if(x1>x2)swap(x1,x2);
        if(y1>y2)swap(y1,y2);
    }
    friend bool operator < (const node &r1,const node &r2)
        {return r1.x2>r2.x2;}
}a[N];
void pushdown(int x)
{
    if(bj[x])
    {
        bj[x<<1]=bj[x<<1|1]=bj[x];
        bj[x]=0;
    }
}
int query(int l,int r,int now,int pos)
{
    if(l==r)return bj[now];
    int mid=(l+r)>>1;
    pushdown(now);
    if(mid>=pos)return query(ls,pos);
    else return query(rs,pos);
}
void update(int l,int r,int now,int lq,int rq,int v)
{
    if(lq<=l&&r<=rq)
        {bj[now]=v;return;}
    int mid=(l+r)>>1;
    pushdown(now);
    if(mid>=lq)update(ls,lq,rq,v);
    if(mid<rq) update(rs,lq,rq,v);
}
void down(int l,int r,int now)
{
    if(bj[now])
    {
        int t=bj[now];
        upd(f[t][0],abs(l-A-Y)+abs(l-A-(a[t].y1-1))+abs(X-(a[t].x2+1)));
        upd(f[t][1],abs(r-A-Y)+abs(r-A-(a[t].y2+1))+abs(X-(a[t].x2+1)));
        return;
    }
    if(l==r)
    {
        if(l==Y+A){printf("%d\n",X+Y);exit(0);}
        return;
    }
    int mid=(l+r)>>1;
    down(ls);down(rs);
}
int main()
{
    //freopen("tt.in","r",stdin);
    scanf("%d%d",&X,&Y);
    scanf("%d",&n);
    for(int i=1;i<=n;i++)
        a[i].read();
    sort(a+1,a+1+n);
    for(int i=n;i>=1;i--)
    {
        if(a[i].x2>=X)continue;
        tar[i][0]=query(0,A*2,1,a[i].y1+A-1);
        tar[i][1]=query(0,A*2,1,a[i].y2+A+1);
        update(0,A*2,1,a[i].y1+A,a[i].y2+A,i);
    }
    memset(f,0x3f,sizeof(f));ans=f[0][0];
    down(0,A*2,1);
    for(int i=1;i<=n;i++)
    {
        if(tar[i][0])
        {
            int t=tar[i][0];
            upd(f[t][0],f[i][0]+abs(a[i].y1-a[t].y1)+abs(a[i].x2-a[t].x2));
            upd(f[t][1],f[i][0]+abs(a[i].y1-a[t].y2-2)+abs(a[i].x2-a[t].x2));
        }
        else upd(ans,f[i][0]+abs(a[i].y1-1)+a[i].x2+1);
        if(tar[i][1])
        {
            int t=tar[i][1];
            upd(f[t][0],f[i][1]+abs(a[i].y2-a[t].y1+2)+abs(a[i].x2-a[t].x2));
            upd(f[t][1],f[i][1]+abs(a[i].y2-a[t].y2)+abs(a[i].x2-a[t].x2));
        }
        else upd(ans,f[i][1]+abs(a[i].y2+1)+a[i].x2+1);
    }
    printf("%lld\n",ans);
    return 0;
}