洛谷1081 (NOIp2012) 开车旅行——倍增预处理

题目：https://www.luogu.org/problemnew/show/P1081

预处理从每个点开始a能走多少、b能走多少。可以像dp一样从后往前推。

但有X的限制。所以该数组可以变成倍增的样子。

预处理第一步的找下一个点可以从后往前弄，在树状数组上二分。然后正常转移即可。

树状数组上的二分有细节：

　　找第一个 f 值为需要的值k的，可以if( query( mid )>=k )ans=mid,r=mid-1；也可以if( query(mid)<=k-1 )ans=mid,l=mid+1; ans++；

倍增数组也有细节，大约就是当目的地为0的时候就不要给距离赋值了。

#include<iostream>
#include<cstdio>
#include<cstring>
#include<algorithm>
#include<cmath>
#define ll long long
using namespace std;
const int N=1e5+;
const double eps=1e-,INF=2e9+;
int n,m,h[N],tp[N],f[N],ps[N],S;
ll X,dp[N][][],mp[N],s[N],ca,cb;
double tc=INF;
int rdn()
{
    int ret=,fx=;char ch=getchar();
    while(ch>''||ch<''){if(ch=='-')fx=-;ch=getchar();}
    while(ch>=''&&ch<='')(ret*=)+=ch-'',ch=getchar();
    return ret*fx;
}
ll rdl()
{
    ll ret=,fx=;char ch=getchar();
    while(ch>''||ch<''){if(ch=='-')fx=-;ch=getchar();}
    while(ch>=''&&ch<='')(ret*=)+=ch-'',ch=getchar();
    return ret*fx;
}
void add(int x,int k){for(;x<=m;x+=(x&-x))f[x]+=k;}
int query(int x){int ret=;for(;x;x-=(x&-x))ret+=f[x];return ret;}
void find(int i,int &ps0,int &ps1)
{
    ps0=;ps1=;
    int k=query(h[i]),l=,r=m,tp0=,tp1=,antp=;
//    printf("i=%d k=%d\n",i,k);
    if(n-i+>k)
    {
        while(l<=r)
        {
            int mid=l+r>>;
            if(query(mid)<=k)tp0=mid,l=mid+;
            else r=mid-;
        }
        tp0++;
    }
    if(k>)
    {
        l=;r=m;
        while(l<=r)
        {
            int mid=l+r>>;
            if(query(mid)>=k-)tp1=mid,r=mid-;
            else l=mid+;
        }
    }
    if(!tp0&&tp1)ps0=ps[tp1],antp=tp1;else if(tp0&&!tp1)ps0=ps[tp0],antp=tp0;
    else if(tp0&&tp1)
        {if(tp[h[i]]-tp[tp1]>tp[tp0]-tp[h[i]])ps0=ps[tp0],antp=tp0;else ps0=ps[tp1],antp=tp1;}
//    printf("tp0=%d tp1=%d ps0=%d\n",tp0,tp1,ps0);
    if(antp)add(antp,-);
    k=query(h[i]);tp0=tp1=;
//    printf("i=%d k=%d(antp=%d)\n",i,k,antp);
    if(n-i>k)
    {
        l=;r=m;
        while(l<=r)
        {
            int mid=l+r>>;
            if(query(mid)<=k)tp0=mid,l=mid+;
            else r=mid-;
        }
        tp0++;
    }
    if(k>)
    {
        l=;r=m;
        while(l<=r)
        {
            int mid=l+r>>;
            if(query(mid)>=k-)tp1=mid,r=mid-;
            else l=mid+;
        }
    }
    if(!tp0&&tp1)ps1=ps[tp1];else if(tp0&&!tp1)ps1=ps[tp0];
    else if(tp0&&tp1){if(tp[h[i]]-tp[tp1]>tp[tp0]-tp[h[i]])ps1=ps[tp0];else ps1=ps[tp1];}
//    printf("tp0=%d tp1=%d ps1=%d\n",tp0,tp1,ps1);
    if(antp)add(antp,);
}
void cz(int cr)
{
    int a,b;find(cr,b,a);
    dp[cr][][]=a;mp[cr]=b;
    if(a)dp[cr][][]=abs(tp[h[a]]-tp[h[cr]]);

    int ta=mp[a];
    if(ta)
    {
        dp[cr][][]=ta;
        dp[cr][][]=dp[cr][][];  if(ta)dp[cr][][]=abs(tp[h[ta]]-tp[h[a]]);
    }
//    printf("cr=%d a=%lld b=%lld(a=%d ta=%d)\n",cr,dp[cr][1][0],dp[cr][1][1],a,ta);

    for(int i=;i<=;i++)
    {
        a=dp[cr][i-][];if(!dp[a][i-][])break;//
        dp[cr][i][]=dp[cr][i-][]+dp[a][i-][];
        dp[cr][i][]=dp[cr][i-][]+dp[a][i-][];
        dp[cr][i][]=dp[a][i-][];
    }
}
int main()
{
    n=rdn();
    for(int i=;i<=n;i++)tp[i]=h[i]=rdn();
    sort(tp+,tp+n+);m=unique(tp+,tp+n+)-tp-;
    for(int i=n;i;i--)
    {
        h[i]=lower_bound(tp+,tp+n+,h[i])-tp;
        add(h[i],);ps[h[i]]=i;cz(i);
    }
    X=rdl();ll yX=X;
    for(int i=;i<=n;i++)
    {
//        printf("i=%d\n",i);
        int now=i;ca=;cb=;X=yX;
        for(int j=;j;j--)
//        {
//            if(j<=3)printf("X=%lld now=%d j=%d a=%lld b=%lld\n"
//                ,X,now,j,dp[now][j][0],dp[now][j][1]);
            if(dp[now][j][]+dp[now][j][]<=X)
            {
                X-=dp[now][j][]+dp[now][j][];
                ca+=dp[now][j][];cb+=dp[now][j][];
                if(dp[now][j][])now=dp[now][j][];
            }
//        }
//        printf("X=%lld now=%d j=0 a=%lld\n",X,now,dp[now][0][0]);
        if(now&&dp[now][][]<=X)ca+=dp[now][][];
//        printf("i=%d ca=%lld cb=%lld\n",i,ca,cb);
        if(cb)
        {
            double w=(double)ca/cb;
            if(w<tc||(fabs(w-tc)<=eps&&h[i]>h[S]))tc=w,S=i;
        }
//        printf("i=%d S=%d\n",i,S);
    }
    printf("%d\n",S);
    m=rdn();
    while(m--)
    {
        S=rdn();X=rdn();
        int now=S;ca=;cb=;
        for(int j=;j;j--)
//        {
//            if(j<=3)printf("X=%lld now=%d j=%d a=%lld b=%lld\n"
//                ,X,now,j,dp[now][j][0],dp[now][j][1]);
            if(dp[now][j][]+dp[now][j][]<=X)
            {
                X-=dp[now][j][]+dp[now][j][];
                ca+=dp[now][j][];cb+=dp[now][j][];
                if(dp[now][j][])now=dp[now][j][];
            }
//        }
//        printf("X=%lld now=%d j=0 a=%lld\n",X,now,dp[now][0][0]);
        if(now&&dp[now][][]<=X)ca+=dp[now][][];
        printf("%lld %lld\n",ca,cb);
    }
    return ;
}