NOIP2012开车旅行

开车旅行

【题目分析】
看了几遍题目，发现这是一道模拟题，可普通的去模拟，超时是妥妥的。首先，我们一定要先预处理出在当前这个点，A和B会去向那个城市（否则每次在模拟的时候，还要先知道要去哪里，这也太鬼畜了）。我们可以用一个对来维护高度有序，A和B要去的点一定是大于等于h[i]的两个和小于h[i]的两个点。
其次，我们也可以知道需要花费的时间。我们基于以往做LCA等在树上操作的题目的想法，知道了可以用倍增来”跳”（花费的时间也同理可以处理出来）。这一题也就解决了。

【代码】

#include<set>
#include<cmath>
#include<cstdio>
#include<cstring>
#include<cstdlib>
#define M 1000005
#include<iostream>
#include<algorithm>
#define ll long long
using namespace std;
void Rd(int &res){
    char c;res=0;int k=1;
    while(c=getchar(),!isdigit(c)&&c!='-');
    if(c=='-')k=-1,c=getchar();
    do{res=(res<<3)+(res<<1)+(c^48);}while(c=getchar(),isdigit(c));
    res*=k;
}
template<class T>void Pt(T x){
    if(x==0)return;
    Pt(x/10);
    putchar(x%10^48);
}
template<class T>void Ps(T x){
    if(x==0)putchar('0');
    else Pt(x);
}
struct City{
    int h,p;
    bool operator<(const City rhs)const{return h<rhs.h;}
}h[M];
set<City>S;
set<City>::iterator it,it1;
int n,x0,m,next[2][M],dist[2][M],g[M][21];
ll f[M][21][2];
void update(City x,City y){
    ll cha=abs(x.h-y.h);
    if(!next[0][x.p]){
        next[0][x.p]=y.p;
        dist[0][x.p]=cha;
    }
    else if(cha<dist[0][x.p]||(cha==dist[0][x.p]&&y.h<h[next[0][x.p]].h)){
        next[1][x.p]=next[0][x.p];
        dist[1][x.p]=dist[0][x.p];
        next[0][x.p]=y.p;
        dist[0][x.p]=cha;
    }
    else if(cha<dist[1][x.p]||(cha==dist[1][x.p]&&y.h<h[next[1][x.p]].h)){
        next[1][x.p]=y.p;
        dist[1][x.p]=cha;
    }
    else if(!next[1][x.p]){
        next[1][x.p]=y.p;
        dist[1][x.p]=cha;
    }
}
void ask(int s,int x,ll &dista,ll &distb){
    for(int i=18;i>=0;i--){
        if(f[s][i][0]+f[s][i][1]<=x&&g[s][i]){
            dista+=f[s][i][0];
            distb+=f[s][i][1];
            x-=f[s][i][0]+f[s][i][1];
            s=g[s][i];
        }
    }if(next[1][s]&&dist[1][s]<=x)dista+=dist[1][s];
}
void Init(){
    Rd(n);
    for(int i=1;i<=n;i++)Rd(h[i].h),h[i].p=i;
    for(int i=n;i;i--){
        S.insert(h[i]);
        it=S.find(h[i]);it1=it;
        if(it!=S.begin()){
            it--;//前一个 
            update(h[i],*it);
            if(it!=S.begin()){
                it--;//再前一个 
                update(h[i],*it);
            }
        }
        if((++it1)!=S.end()){//后一个，等等 
            update(h[i],*it1);
            if((++it1)!=S.end())update(h[i],*it1);
        }
    }
    for(int i=1;i<=n;i++){
        f[i][0][0]=dist[1][i];
        g[i][0]=next[0][next[1][i]];
        f[i][0][1]=dist[0][next[1][i]];
    }
    for(int j=1;j<=18;j++){
        for(int i=1;i<=n;i++){
            g[i][j]=g[g[i][j-1]][j-1];
            f[i][j][0]=f[i][j-1][0]+f[g[i][j-1]][j-1][0];
            f[i][j][1]=f[i][j-1][1]+f[g[i][j-1]][j-1][1];
        }
    }
}
int main(){
    Init();
    int s0=1;
    ll a=1e15,b=0;
    Rd(x0);
    for(int i=1;i<=n;i++){
        ll dista=0,distb=0;
        ask(i,x0,dista,distb);
        if(!distb){
            if(!b&&h[i].h>h[s0].h)s0=i;
        }//else if(!b)
        else if(distb&&(!s0||1.0*a*distb>1.0*b*dista))s0=i,a=dista,b=distb;//防止精度不够，故用乘 
        else if(distb&&(!s0||1.0*a*distb==1.0*b*dista)){
            if(h[i].h>h[s0].h)s0=i;
        }   
    }Ps(s0);putchar('\n');

    Rd(m);
    while(m--){
        int s,x;
        Rd(s),Rd(x);
        ll dista=0,distb=0;
        ask(s,x,dista,distb);
        Ps(dista),putchar(' '),Ps(distb),putchar('\n');
    }return 0;
}