【NOIP2012】旅行计划

题解

双向链表加倍增。。。

正写着不一定能写对2333

终于写对了。。。

然而我的双向链表和别人的都不一样。。。

瑟瑟发抖。。。

代码

//by 减维
#include<cstdio>
#include<iostream>
#include<cstring>
#include<queue>
#include<cstdlib>
#include<ctime>
#include<cmath>
#include<map>
#include<algorithm>
#define ll long long
#define inf (9223372036854775807LL)
using namespace std;

struct biao{
    int i,to,fr;
    ll v;
}e[];

int n,m,hest,ys[],zj[],cj[],f[][];
ll x,her,a[][],b[][],h[];

bool cmp(const biao&x,const biao&y){return x.v>y.v;}
bool cm2(const biao&x,const biao&y){return x.i<y.i;}

int main()
{
    scanf("%d",&n);
    for(int i=;i<=n;++i)
    {
        scanf("%lld",&e[i].v),e[i].i=i,h[i]=e[i].v;
        if(h[i]>her) her=h[i],hest=i;
    }
    sort(e+,e+n+,cmp);
    for(int i=;i<=n;++i)
        e[i].fr=e[i-].i,e[i].to=e[i+].i;
    sort(e+,e+n+,cm2);
    for(int i=;i<=n;++i)
    {
        ll m1=inf,m2=inf,hh[];int pos[],a1,a2;

        hh[]=abs(h[i]-h[e[i].fr]);                    pos[]=e[i].fr;
        hh[]=abs(h[i]-h[e[e[i].fr].fr]);              pos[]=e[e[i].fr].fr;
        hh[]=abs(h[i]-h[e[i].to]);                    pos[]=e[i].to;
        hh[]=abs(h[i]-h[e[e[i].to].to]);           pos[]=e[e[i].to].to;

        if(e[i].fr==)hh[]=inf;
        if(e[e[i].fr].fr==)hh[]=inf;
        if(e[i].to==)hh[]=inf;
        if(e[e[i].to].to==)hh[]=inf;

        for(int j=;j<=;++j)
            if(m1>hh[j]||(m1==hh[j]&&h[a1]>h[pos[j]]))m2=m1,a2=a1,m1=hh[j],a1=pos[j];
        else
            if(m2>hh[j]||(m2==hh[j]&&h[a2]>h[pos[j]]))m2=hh[j],a2=pos[j];

        if(m1!=inf)zj[i]=a1;
        if(m2!=inf)cj[i]=a2;
        if(e[i].fr!=)e[e[i].fr].to=e[i].to;
        if(e[i].to!=)e[e[i].to].fr=e[i].fr;

    }
    for(int i=;i<=n;++i)
    {
        f[i][]=zj[cj[i]];
        if(i<=n-)
            a[i][]=abs(h[i]-h[cj[i]]);
        if(i<=n-)
            b[i][]=abs(h[cj[i]]-h[zj[cj[i]]]);
    }
    for(int i=;i<=;++i)
        for(int j=;j<=n;++j)
        {
            f[j][i]=f[f[j][i-]][i-];
            a[j][i]=a[j][i-]+a[f[j][i-]][i-];
            b[j][i]=b[j][i-]+b[f[j][i-]][i-];
        }

    scanf("%lld",&x);
    int ans=;
    double m1=inf;
    ll mh=-;
    for(int i=;i<=n;++i)
    {
        int pos=i;
        ll l=x,disa=,disb=;
        for(int j=;j>=;--j)
        {
            if(f[pos][j]!=&&a[pos][j]+b[pos][j]<=l)
            {
                disa+=a[pos][j];
                disb+=b[pos][j];
                l-=a[pos][j]+b[pos][j];
                pos=f[pos][j];
            }
        }
        if(a[pos][]<=l){
            disa+=a[pos][];
            l-=a[pos][];
        }
        double c=inf;
        if(disb==) c=inf;
        else c=disa*1.0/disb;
        if(c<m1||(c==m1&&h[i]>mh))
        {
            m1=c;
            mh=h[i];
            ans=i;
        }
    }
    if(m1!=inf)printf("%d\n",ans);
    else printf("%d\n",hest);
    scanf("%d",&m);
    for(int i=;i<=m;++i)
    {
        int s;
        scanf("%d%lld",&s,&x);
        int pos=s;
        ll l=x,disa=,disb=;
        for(int j=;j>=;--j)
        {
            if(f[pos][j]!=&&a[pos][j]+b[pos][j]<=l)
            {
                disa+=a[pos][j];
                disb+=b[pos][j];
                l-=a[pos][j]+b[pos][j];
                pos=f[pos][j];
            }
        }
        if(a[pos][]<=l){
            disa+=a[pos][];
            l-=a[pos][];
        }
        printf("%lld %lld\n",disa,disb);
    }
}