POJ2374 Fence Obstacle Course DP+线段树优化

　　题目链接：http://poj.org/problem?id=2374

　　这个题目坑了好久啊，题意看错了，bs自己。

　　其实不难的，状态方程容易相处来，每个fence只有两种状态，即左边下去和右边下去，如果要知道当前状态的最优情况，就必须知道前面所有状态的最优情况，最坏情况下复杂度O(n^2)。但显然可以发现前面有很多的状态都不能到达当前状态于，因为前面的有些状态被覆盖了，于是我们可以用线段树来优化，在log(n)的时间里找到可以到达当前状态的fence。虽然在log(n)的时间里找出了状态，但是还是有比较多的状态，还可以继续优化。注意到编号为 i 的fence和所有编号为 1,2,,,,i-1 的fence，编号为 i 的fence只能有前面一个到达，因此我们可以倒叙来DP。

//STATUS:C++_AC_329MS_4044KB
#include<stdio.h>
#include<stdlib.h>
#include<string.h>
#include<math.h>
#include<iostream>
#include<string>
#include<algorithm>
#include<vector>
#include<queue>
#include<stack>
#include<map>
using namespace std;
#define LL __int64
#define pii pair<int,int>
#define Max(a,b) ((a)>(b)?(a):(b))
#define Min(a,b) ((a)<(b)?(a):(b))
#define mem(a,b) memset(a,b,sizeof(a))
#define lson l,mid,rt<<1
#define rson mid+1,r,rt<<1|1
const int N=50010,M=100000,INF=0x3f3f3f3f,MOD=100000000;
const double DNF=100000000000;

int l[N][2],f[N][2],sta[200010<<2];
int n,s,w,a,b;

void update(int l,int r,int rt)
{
    if(a<=l && r<=b){
        sta[rt]=w;
        return;
    }
    if(sta[rt]>=0){
        sta[rt<<1]=sta[rt<<1|1]=sta[rt];
        sta[rt]=-1;
    }
    int mid=(l+r)>>1;
    if(a<=mid)update(lson);
    if(b>mid)update(rson);
}

void query(int l,int r,int rt)
{
    if(sta[rt]!=-1){
        w=sta[rt];
        return;
    }
    int mid=(l+r)>>1;
    if(a<=mid)query(lson);
    else query(rson);
}

int main()
{
 //   freopen("in.txt","r",stdin);
    int i,j,ans;
    while(~scanf("%d%d",&n,&s))
    {
        mem(f,INF);
        mem(sta,0);
        for(i=1;i<=n;i++){
            scanf("%d%d",&l[i][0],&l[i][1]);
            l[i][0]+=M;l[i][1]+=M;
        }
        f[0][0]=f[0][1]=0;
        l[0][0]=l[0][1]=M;
        for(i=1;i<=n;i++){
            w=0;a=l[i][0];
            query(0,M<<1,1);
            f[i][0]=Min(f[w][0]+abs(l[i][0]-l[w][0]),f[w][1]+abs(l[i][0]-l[w][1]));
            w=0;a=l[i][1];
            query(0,M<<1,1);
            f[i][1]=Min(f[w][0]+abs(l[i][1]-l[w][0]),f[w][1]+abs(l[i][1]-l[w][1]));
            a=l[i][0],b=l[i][1];w=i;
            update(0,M<<1,1);
        }

        ans=Min(f[n][0]+abs(s+M-l[n][0]),f[n][1]+abs(s+M-l[n][1]));
        printf("%d\n",ans);
    }
    return 0;
}