倍增+双向链表——洛谷P1081 开车旅行

https://daniu.luogu.org/problem/show?pid=1081#sub
提高组的第三题；
哇；

首先题目很恶心；
然而发现也没什么歧义，也不难看懂；
70分的暴力做法；
n*n预处理出A[i],B[i]两个数组；
标表示在i点时小A和小B开车分别会走到哪里；
这个只要暴力高就好了；
指的一提的是负最大初始值1e9不够；
要3e9;
然后我们就可以发现也已搞一下优化；
我们把读入排序；
拍好后搞成一个双向链表；
然后按照点的先后顺序去枚举；
对于k点；
显然A[i],B[i]；
出自双向链表的
k-1,k-2,k+1,k+2里面；
我们统计号好k之后直接把k删掉；
这样对于后面的k，链表里面的非k点的顺序一定比k后；
这样预处理时间就变成了nlogn;
然后我们把小A一天+小B一天合为一步；
用这个去搞倍增表；
然后直接像lca一样去倍增就好啦；

#include<iostream>#include<cstdio>#include<cstdlib>#include<algorithm>using namespace std;const int N=1e5+100;struct L2{    int v,l,r,num;//双向链表 }L[N];struct zz{//一个答案结构体     int a,b;    zz(){a=b=0;}};int a[N],A[N],B[N],c[N];int f[N][20],AA[N][20],BB[N][20];//倍增表，AA BB是这段倍增区间A B的花费 int x,y,z,n,m;bool cmp(L2 a,L2 b){return a.v<b.v;}void make(int k){    int mst=3e9,med=3e9,aa=0,bb=0,v,m;    m=L[k].l;    if(m){        v=abs(L[k].v-L[m].v);        if(v<mst||(v==mst&&L[k].v>L[m].v)){            med=mst; bb=aa; mst=v; aa=L[m].num;        }else        if(v<med||(v==med&&L[k].v>L[m].v)){            med=v; bb=L[m].num;        }        m=L[m].l;        if(m){            v=abs(L[k].v-L[m].v);            if(v<mst||(v==mst&&L[k].v>L[m].v)){                med=mst; bb=aa; mst=v; aa=L[m].num;            }else            if(v<med||(v==med&&L[k].v>L[m].v)){                med=v; bb=L[m].num;            }        }    }    m=L[k].r;    if(m){        v=abs(L[k].v-L[m].v);        if(v<mst||(v==mst&&L[k].v>L[m].v)){            med=mst; bb=aa; mst=v; aa=L[m].num;        }else        if(v<med||(v==med&&L[k].v>L[m].v)){            med=v; bb=L[m].num;        }        m=L[m].r;        if(m){            v=abs(L[k].v-L[m].v);            if(v<mst||(v==mst&&L[k].v>L[m].v)){                med=mst; bb=aa; mst=v; aa=L[m].num;            }else            if(v<med||(v==med&&L[k].v>L[m].v)){                med=v; bb=L[m].num;            }        }    }    A[L[k].num]=bb;//注意，我的写法比较怪异，A和aa翻一下的     B[L[k].num]=aa;    L[L[k].r].l=L[k].l;//删除     L[L[k].l].r=L[k].r;}void makebz(){    for(int i=1;i<=n;i++){        f[i][0]=B[A[i]];        AA[i][0]=abs(a[i]-a[A[i]]);        BB[i][0]=abs(a[A[i]]-a[B[A[i]]]);    }    for(int j=1;(1<<j)<=n;j++)        for(int i=1;i<=n;i++)        if(f[i][j-1]){            f[i][j]=f[f[i][j-1]][j-1];            AA[i][j]=AA[i][j-1]+AA[f[i][j-1]][j-1];            BB[i][j]=BB[i][j-1]+BB[f[i][j-1]][j-1];        }}zz work(int now,int x){    int aa=0,bb=0,j;    while(1){        j=0;        for(;f[now][j]&&(aa+AA[now][j]+bb+BB[now][j]<=x);j++);        if(!j)break;        j--;        aa+=AA[now][j];        bb+=BB[now][j];        now=f[now][j];            }    if(A[now]&&(aa+bb+abs(a[A[now]]-a[now])<=x))aa+=abs(a[A[now]]-a[now]);    zz ans;                 //最后A可能还会走一下     ans.a=aa;    ans.b=bb;    return ans;}void outit1(){    int x,ans;    double Ans=1e9+5,temp;    zz c;    scanf("%d",&x);    for(int i=1;i<=n;i++){        c=work(i,x);        if(!c.b)temp=1e9;else temp=(double)c.a/c.b;        if(Ans>temp||(Ans==temp&&a[ans]<a[i]))Ans=temp,ans=i;    }    printf("%d\n",ans);}void outit2(){    int m,x,y;    zz c;    scanf("%d",&m);    for(int i=1;i<=m;i++){        scanf("%d%d",&x,&y);        c=work(x,y);        printf("%d %d\n",c.a,c.b);    }}int main(){    scanf("%d",&n);    for(int i=1;i<=n;i++)scanf("%d",&a[i]),L[i].v=a[i],L[i].num=i;    sort(L+1,L+n+1,cmp);    for(int i=1;i<=n;i++){//创建一个链表         L[i].l=i-1;        L[i].r=i+1;        c[L[i].num]=i;    }    L[n].r=0;    for(int i=1;i<=n;i++)make(c[i]);//按顺序来处理点     makebz();//做倍增表     outit1();    outit2();}