【noip2012】开车旅行 [倍增]
开车旅行
倍增
可以发现每个起点出发后面选择的城市都是一定的 所以预处理出\(to[i][j],da[i][j],db[i][j]\)表示从\(i\)出发小\(A\)和小\(B\)经过\(2^j\)轮后到达的地点、小\(A\)走的路程、小\(B\)走的路程
预处理时询问第一近和第二近的地点用双向链表
双向链表
不仅记录当前元素的下一个元素\((next)\)还记录当前元素的上一个元素\((pre)\)
数组版：
int head,tail,tot;
struct Node{
int value;
int pre,nxt;
}node[SIZE];
int initialize(){//新建链表
tot=2;
head=1,tail=2;
node[head].nxt=tail,node[tail].pre=head;
}
int insert(int p,int val){//在p后插入包含数据val的新节点
q=++tot;
node[q].value=val,node[node[q].nxt].pre=q,
node[q].nxt=node[p].nxt,node[p].nxt=q,node[q].pre=p;
}
void remove(int p){
node[node[p].pre].nxt=node[p].nxt,
node[node[p].nxt].pre=node[p].pre;
}
void clear(){
memset(node,0,sizeof(node));
head=tail=tot=0;
}
struct Node{
int value;
Node *prev,next;
};
Node *head,*tail;
void initialize(){
head=new Node(),tail=new Node();
head->next=tail,tail->prev=head;
}
void insert(Node *p,int val){
Node q=new Node(),q->value=val;
p->next->prev=q,q->next=p->next,
p->next=q,q->prev=p;
}
void remove(Node *p){
p->prev->next=p->next,p->next->prev=p->prev;
}
void recycle(){//链表内存回收
while(head!=tail) head=head->next,delete head->prev;
delete tail;
}
邻值查找：给定一个长度为\(n\)的序列\(A\)，\(A\)中的数各不相同。对于\(A\)中的每一个数 \(A_i\)，求：\(min_{1≤j<i}|Ai−Aj|\)
以及令上式取到最小值的\(j\)（记为\(P_i\)）。若最小值点不唯一，则选择使 AjAj 较小的那个。
int n,id[N],pre[N],nxt[N];
pair<int,int>ans[N];
template<class t>void rd(t &x){
x=0;int w=0;char ch=0;
while(!isdigit(ch)) w|=ch=='-',ch=getchar();
while(isdigit(ch)) x=(x<<1)+(x<<3)+(ch^48),ch=getchar();
x=w?-x:x;
}
struct node{int x,id;}a[N];
bool cmp(node x,node y){return x.x<y.x;}
int main(){
#ifndef ONLINE_JUDGE
freopen("in.txt","r",stdin);
#endif
rd(n);
for(int i=1;i<=n;++i) rd(a[i].x),a[i].id=i;
sort(a+1,a+n+1,cmp);a[0].x=-3000000000,a[n+1].x=3000000000;
for(int i=1;i<=n;++i) id[a[i].id]=i,pre[i]=i-1,nxt[i]=i+1;
for(int i=n,nw,l,r,lw,rw;i>=2;--i){
nw=id[i],l=pre[nw],r=nxt[nw];
lw=abs(a[nw].x-a[l].x),rw=abs(a[nw].x-a[r].x);
if(lw<=rw) ans[i]=make_pair(lw,a[l].id);
else if(lw>rw) ans[i]=make_pair(rw,a[r].id);
pre[nxt[nw]]=pre[nw],nxt[pre[nw]]=nxt[nw];
}
for(int i=2;i<=n;++i) printf("%d %d\n",ans[i].first,ans[i].second);
return 0;
}
code
struct node{int h,id;}a[N];
bool cmp(node x,node y){return x.h<y.h;}
int pos[N][2],to[N][20];
void upd(int nw,int x){
ll d=abs(h[nw]-a[x].h);
if(dis[nw][0]>d||(dis[nw][0]==d&&a[x].h<h[pos[nw][0]]))
dis[nw][1]=dis[nw][0],dis[nw][0]=d,pos[nw][1]=pos[nw][0],pos[nw][0]=a[x].id;
else if(dis[nw][1]>d||(dis[nw][1]==d&&a[x].h<h[pos[nw][1]]))
dis[nw][1]=d,pos[nw][1]=a[x].id;
}
void erase(int x){nxt[pre[x]]=nxt[x],pre[nxt[x]]=pre[x];}
void prepa(){
for(int i=1,nw,l1,l2,r1,r2;i<=n;++i){
nw=id[i],l1=pre[nw],l2=pre[l1],r1=nxt[nw],r2=nxt[r1];
pos[i][0]=pos[i][1]=0,dis[i][0]=dis[i][1]=inf;
if(l1&&l1!=n+1) upd(i,l1);
if(l2&&l2!=n+1) upd(i,l2);
if(r1&&r1!=n+1) upd(i,r1);
if(r2&&r2!=n+1) upd(i,r2);
erase(nw);
}
for(int i=1;i<=n;++i) to[i][0]=pos[pos[i][1]][0],da[i][0]=dis[i][1],db[i][0]=dis[pos[i][1]][0];
for(int j=1;j<=17;++j)
for(int i=1;i<=n;++i)
to[i][j]=to[to[i][j-1]][j-1],da[i][j]=da[i][j-1]+da[to[i][j-1]][j-1],db[i][j]=db[i][j-1]+db[to[i][j-1]][j-1];
}
void query(int s,int lim,ll &A,ll &B){
for(int i=17;i>=0;--i)
if(to[s][i]&&da[s][i]+db[s][i]<=lim)
lim-=da[s][i]+db[s][i],A+=da[s][i],B+=db[s][i],s=to[s][i];
if(pos[s][1]&&dis[s][1]<=lim) A+=dis[s][1];
}
int main(){
#ifndef ONLINE_JUDGE
freopen("in.txt","r",stdin);
#endif
rd(n);
for(int i=1;i<=n;++i) rd(h[i]),a[i]=(node){h[i],i};
sort(a+1,a+n+1,cmp);
for(int i=1;i<=n;++i) id[a[i].id]=i,pre[i]=i-1,nxt[i]=i+1;
rd(x0),rd(m),pre[0]=nxt[n+1]=0,pre[n+1]=n,nxt[0]=1;
prepa();
ll A,B;double mn=inf,ret;int ans;
for(int i=1;i<=n;++i){
query(i,x0,A=0,B=0);
if(!B) continue;
ret=(A*1.0)/(B*1.0);
if(mn>ret) mn=ret,ans=i;
else if(ret==mn) if(h[i]>h[ans]) ans=i;
}
printf("%d\n",ans);
for(int i=1,s,lim;i<=m;++i){
rd(s),rd(lim),A=0,B=0;
query(s,lim,A,B);
printf("%lld %lld\n",A,B);
}
return 0;
}
