【开车旅行】题解(NOIP2012提高组)
分析
首先我们可以发现，两个询问都可以通过一个子程序来求。
接着，如果每到一个城市再找下一个城市，显然是行不通的。所以首先先预处理从每一个城市开始，小A和小B要去的城市。预处理的方法很多，我用的是双向链表：按高度把城市排序，用双向链表把每个城市的相邻的城市里连起来，再从每个城市按双向链表往两边搜，如果搜到的城市在这个城市的西边就删掉，否则记录。左右分别记录两个城市，排个序就可以的出小A和小B要去的城市了。
然后我们就可以发现这是一棵树，所以倍增。
我们设g[i][j]指从i城市走j2轮走到的点（一轮指小A和小B各走一次，小A先走），fa[i][j]指指从i城市走j2轮的小A走过路程，fb[i][j]指指从i城市走j^2轮的小B走过路程。转移显然。
#include <cmath>
#include <iostream>
#include <cstdio>
#include <cstdlib>
#include <cstring>
#include <algorithm>
#include <queue>
const long long maxlongint=2147483647;
using namespace std;
struct lyd
{
long long last;
long long next;
};
lyd a[200000];
long long h[200000][2],c[200000][5][3],qu[200000][3],yh[200000],lf[3],n,m,tot,sum,num,g[150000][17],fa[150000][17],fb[150000][17];
double ans;
void q(long long l,long long r)
{
long long i=l,j=r,mid=h[(l+r)/2][1],e;
while(i<j)
{
while(h[i][1]<mid) i++;
while(h[j][1]>mid) j--;
if(i<=j)
{
e=h[i][0];
h[i][0]=h[j][0];
h[j][0]=e;
e=h[i][1];
h[i][1]=h[j][1];
h[j][1]=e;
i++;
j--;
}
}
if(i<r) q(i,r);
if(l<j) q(l,j);
}
int work(long long x,long long value)
{
long long i,j,k=0,l;
for(i=log2(n);i>=0;i--)
{
if(fa[x][i]+fb[x][i]+k<=value)
{
k+=fa[x][i]+fb[x][i];
lf[0]+=fa[x][i];
lf[1]+=fb[x][i];
x=g[x][i];
}
}
if(c[x][2][1]+k<=value)
{
k+=c[x][2][1];
lf[0]+=c[x][2][1];
}
}
int main()
{
long long i,j,k,l,x,y;
scanf("%lld",&n);
for(i=1;i<=n;i++)
{
scanf("%lld",&h[i][1]);
yh[i]=h[i][1];
h[i][0]=i;
}
scanf("%lld",&qu[0][2]);
scanf("%lld",&m);
for(i=1;i<=m;i++)
{
scanf("%lld%lld",&qu[i][1],&qu[i][2]);
}
q(1,n);
for(i=1;i<=n;i++)
{
a[h[i][0]].last=h[i-1][0];
a[h[i][0]].next=h[i+1][0];
}
for(i=1;i<=n;i++)
{
j=i;
while(a[j].last!=0)
{
j=a[j].last;
if(j<i)
{
a[a[j].next].last=a[j].last;
a[a[j].last].next=a[j].next;
j=a[j].last;
}
if(j!=0)
{
c[i][0][0]++;
c[i][c[i][0][0]][0]=j;
c[i][c[i][0][0]][1]=abs(yh[i]-yh[j]);
c[i][c[i][0][0]][2]=yh[j];
if(c[i][0][0]==2) break;
}
}
j=i;
while(a[j].next!=0)
{
j=a[j].next;
if(j<i)
{
a[a[j].next].last=a[j].last;
a[a[j].last].next=a[j].next;
j=a[j].next;
}
if(j!=0)
{
c[i][0][0]++;
c[i][c[i][0][0]][0]=j;
c[i][c[i][0][0]][1]=abs(yh[i]-yh[j]);
c[i][c[i][0][0]][2]=yh[j];
if(c[i][0][0]==4) break;
}
}
for(j=1;j<=c[i][0][0];j++)
{
for(k=1;k<=c[i][0][0];k++)
if(c[i][j][1]<c[i][k][1] || c[i][j][1]==c[i][k][1] && c[i][j][2]<c[i][k][2])
{
l=c[i][k][1];
c[i][k][1]=c[i][j][1];
c[i][j][1]=l;
l=c[i][k][0];
c[i][k][0]=c[i][j][0];
c[i][j][0]=l;
l=c[i][k][2];
c[i][k][2]=c[i][j][2];
c[i][j][2]=l;
}
}
}
for(i=1;i<=n;i++)
{
g[i][0]=c[c[i][2][0]][1][0];
fa[i][0]=c[i][2][1];
fb[i][0]=c[c[i][2][0]][1][1];
}
for(j=1;j<=log2(n);j++)
{
for(i=1;i<=n;i++)
{
g[i][j]=g[g[i][j-1]][j-1];
fa[i][j]=fa[i][j-1]+fa[g[i][j-1]][j-1];
fb[i][j]=fb[i][j-1]+fb[g[i][j-1]][j-1];
}
}
ans=double(maxlongint);
for(i=1;i<=n;i++)
{
lf[0]=lf[1]=0;
work(i,qu[0][2]);
if(lf[1]!=0)
if(double(lf[0])/double(lf[1])<ans)
{
ans=double(lf[0])/double(lf[1]);
num=i;
}
}
cout<<num<<endl;
for(i=1;i<=m;i++)
{
lf[0]=lf[1]=0;
work(qu[i][1],qu[i][2]);
printf("%lld %lld\n",lf[0],lf[1]);
}
}
