GDOI2016模拟8.18的士

题目
Bessie为农场上的其他奶牛提供的士服务。奶牛们在一条长为M（1<=M<=1,000,000,000）的栅栏的不同位置上。不幸的是，它们厌倦了它们现在所在的位置而想要去栅栏上其他的位置。Bessie必须把她每一个朋友从它们各自的起始地接上车然后送它们到目的地。但Bessie的车太小了，所以她每次只能运送一只奶牛。奶牛们上车下车是瞬间的事情。

为了省油钱，Bessie想要使她的驾驶量最小。给出N只奶牛（1<=N<=100,000）每一只的起始地和目的地，计算Bessie最少需要的驾驶量。Bessie意识到她需要偶尔把牛放在某一个地方而不是把它送到目的地才能省油钱。

Bessie在栅栏最左端（位置0）开始工作，而且必须在最右端（位置M）结束她的工作。

这题由于标解太抽象了（虽然效率快得惊人），我们这里讲一种简单易懂的方法。

好了，我们首先看看题目让我们做些什么：
确定最少驾驶量

这个东西由什么组成呢？很多来来去去的路线，想着很乱（一般我们会全局观察，然后想到这里）

到这里，我们会发现一起考虑会很乱，因而会想着去简化它：变成很多段路线，每一段可以往返多次，答案就是每一段路线的距离乘上他们的最小通过次数。

那么限制着通过次数的量有：向右要通过这段路的牛的个数k和向左要通过这段路的牛的个数p。
这时候分情况考虑：
1：（假如什么都没有）k=0、p=0，这时我们知道必须向右走（因为最终要到达最右端），所以通过次数为1；
2：k=0，p＞0，由于我们知道一开始是向右走，走到右边p个牛那里，我们会折返才能送一头牛，这样我们得走2*p次，然而由于我们最终得走到右端，所以还要加上1，走2*p+1次，（p=0,k＞0同理可得为2*k-1）
3：k＞0，p＞0，这时，若k＞p我们要走2*k-1次，若k≤p，我们要走2*p+1次，这个我们可以视作2情况中k、p取max值（大的那边肯定可以捎上小的那边）

分析完后，我们可以发现，段与段之间可以相互独立（因为我们考虑的是必须走的，为下界），同时由于都是下界，那么合起来就一定是答案的下界，也就是最终的答案。

具体过程可以离线坐标，然后用+1表示从第i个位置开始多了一头牛向右走，-1表示有一头牛走完了（向左向右要分开处理），然后前缀和就可得到到第i个点有多少头牛向左或向右走，然后按上述情况考虑即可。

贴代码

#include<iostream>
#include<algorithm>
#include<cstdio>
#define N 200101
using namespace std;
long long a[N][2],b[N][2],c[N*2],d[N*2],e[N*2];
int n,m,topa,topb;
long long ans;
void init(){
    static int x,y;
    scanf("%d %d",&n,&m);
    for (int i=1;i<=n;i++){
        scanf("%d %d",&x,&y);
        if (x<y)
            a[++topa][0]=x,a[topa][1]=y;
        else
            b[++topb][0]=x,b[topb][1]=y;
        c[i]=x,c[i+n]=y;
    }
    c[0]=n+n+2;
    c[n+n+1]=0,c[n+n+2]=m;
}
int er(int x){
    static int l,r,mid;
    l=1,r=c[0];
    while (l<=r)
        if (c[mid=(l+r)/2]>x)r=mid-1;else l=mid+1;
    if (!l)l=1;
    if (l>c[0])l=c[0];
    while (c[l]>x)l--;
    while (c[l]<x)l++;
    return l;
}
void pre(){
    static int x;
    sort(c+1,c+c[0]+1);
    x=c[0];
    c[0]=1;
    for (int i=2;i<=x;i++)
        if (c[i]!=c[c[0]])c[++c[0]]=c[i];
    for (int i=1;i<=topa;i++)
        d[er(a[i][0])]++,d[er(a[i][1])]--;
    for (int i=1;i<=topb;i++)
        e[er(b[i][0])]++,e[er(b[i][1])]--;
}
void work(){
    static int suma,sumb;
    suma=0,sumb=0;
    for (int i=c[0];i;i--)
        sumb+=e[i];
    for (int i=1;i<c[0];i++){
        suma+=d[i],sumb-=e[i];
        if (!suma){
            if (!sumb)ans+=(c[i+1]-c[i]);
            else
                ans+=(c[i+1]-c[i])*(sumb*2+1);
        }else
        if (!sumb)ans+=(c[i+1]-c[i])*(suma*2-1);
        else
            if (suma>sumb)
                ans+=(c[i+1]-c[i])*(suma*2-1);
            else
                ans+=(c[i+1]-c[i])*(sumb*2+1);
    }
}
void write(){
    printf("%lld",ans);
}
int main(){
    init();
    pre();
    work();
    write();
    return 0;
}