洛谷P1081——开车旅行

传送门：QAQQAQ

题意注意点：

1.是从前往后走，不能回头

2.小A小B轮流开，先小A开，而小A是到第二近的点（这点调试的时候查了好久）

3.若绝对值差相同海拔低的更近，而第一个询问若比值相同是海拔高的更优

思路：我们先预处理出离i点最近和次近点的编号，没有就是-1（链表可以达到O(n)，但是用STL O(nlog(n)) 绰绰有余，代码量还小）

然后用倍增预处理出从i号点走2^j轮A走的路程，B走的路程，和走完这些轮最终停下的点（定义AB各走一次为一轮）

注意可能走半轮，即最后一下A走B不走。

写题过程：

大佬MYY告诉我不要心急，先打70分暴力，然后打了2个小时暴力（STL真的好难查啊QAQ） 打完暴力以后告诉我主程序要全部删掉重写。。。。。

看到直接旁边的HCY没打暴力直接正解AC，心态崩了。。。

然后就在没有特判掉走不到的情况和AB顺序一半写正，一般写反上调了好久，不过最终做出来还是蛮有成就感的（感觉自己很弱啊。。）

代码：

#include<bits/stdc++.h>
#define m_k make_pair
#define ll long long
using namespace std;
const ll N=;
const ll inf=(ll)5e18;

ll f[N][],n,m,a[N],x0;
//0:closest 1:second closest
ll dp[N][][],Bits[],p[N][];
//0:disB 1:disA p:finalpos AB为一轮，走2^t轮 

ll _abs(ll x)
{
    if(x<) return -x;
    else return x;
}

set<pair<ll,ll> > st;
vector<pair<ll,ll> > tmp;
void init()
{
    memset(f,-,sizeof(f));
    st.insert(m_k(a[n],n));
    st.insert(m_k(a[n-],n-));
    f[n-][]=n;
    for(ll i=n-;i>=;i--)
    {
        ll bl=;
        tmp.clear();
        set<pair<ll,ll> > :: iterator it;
        it=st.lower_bound(m_k(a[i],i));
        if(it!=st.end()) tmp.push_back(*it),it++,bl++;
        if(it!=st.end()) tmp.push_back(*it);
        if(bl) it--;
        it--;
        set<pair<ll,ll> > :: iterator be;
        be=st.begin(); be--;
        if(it!=be) tmp.push_back(*it),it--;
        if(it!=be) tmp.push_back(*it);
        ll ans1=inf,ans2=inf;
        for(ll j=;j<(ll)tmp.size();j++)
        {
            if(ans1>_abs(a[i]-tmp[j].first)||(ans1==_abs(a[i]-tmp[j].first)&&a[f[i][]]>tmp[j].first))
            {
                f[i][]=f[i][]; ans2=ans1;
                f[i][]=tmp[j].second; ans1=_abs(a[i]-tmp[j].first);
            }
            else if(ans2>_abs(a[i]-tmp[j].first)||(ans2==_abs(a[i]-tmp[j].first)&&a[f[i][]]>tmp[j].first))
            {
                f[i][]=tmp[j].second; ans2=_abs(a[i]-tmp[j].first);
            }
        }
        st.insert(m_k(a[i],i));
    }
}

void ready()
{
    memset(p,-,sizeof(p));
    memset(dp,-,sizeof(dp));
    for(ll i=;i<=n;i++)
    {
        ll pos=f[i][]; if(pos==-) continue;
        dp[i][][]=_abs(a[pos]-a[i]);
        if(f[pos][]==-) continue;
        dp[i][][]=_abs(a[pos]-a[f[pos][]]);
        p[i][]=f[pos][];
    }
    for(ll j=;j<;j++)
    {
        for(ll i=;i<=n;i++)
        {
            for(ll k=;k<=;k++)
                if(dp[i][j-][k]!=-&&dp[p[i][j-]][j-][k]!=-) dp[i][j][k]=dp[i][j-][k]+dp[p[i][j-]][j-][k];
            if(p[i][j-]!=-) p[i][j]=p[p[i][j-]][j-];

        }
    }
}

ll suma,sumb;
double calc(ll x,ll y)
{
    if(y==) return inf*1.0;
    return(x*1.0)/(y*1.0);
}

int main()
{
    scanf("%lld",&n);
    for(ll i=;i<=n;i++) scanf("%lld",&a[i]);
    init();
    ready();
    scanf("%lld",&x0);
    double ans=inf*1.0;
    ll now=;
    a[]=-inf;
    for(int i=;i<=n;i++)
    {
        ll suma=,sumb=,pos=i;
        for(int j=;j>=;j--)
        {
            if(dp[pos][j][]==-||dp[pos][j][]==-||p[pos][j]==-) continue;
            if(dp[pos][j][]+dp[pos][j][]+suma+sumb>x0) continue;
            suma+=dp[pos][j][];
            sumb+=dp[pos][j][];
            pos=p[pos][j];
        }
        if(dp[pos][][]+suma+sumb<=x0&&dp[pos][][]!=-) suma+=dp[pos][][];
        if(ans>calc(suma,sumb)||(ans==calc(suma,sumb)&&a[i]>a[now]))
        {
            now=i;
            ans=calc(suma,sumb);
        }
    }
    cout<<now<<endl; ll pos;
    scanf("%lld",&m);
    for(int i=;i<=m;i++)
    {
        scanf("%lld%lld",&pos,&x0);
        suma=; sumb=;
        for(int j=;j>=;j--)
        {
            if(dp[pos][j][]==-||dp[pos][j][]==-||p[pos][j]==-) continue;
            if(dp[pos][j][]+dp[pos][j][]+suma+sumb>x0) continue;
            suma+=dp[pos][j][];
            sumb+=dp[pos][j][];
            pos=p[pos][j];
        }
        if(dp[pos][][]+suma+sumb<=x0&&dp[pos][][]!=-) suma+=dp[pos][][];
        printf("%lld %lld\n",suma,sumb);
    }
    return ;
}