P1081 [NOIP2012 提高组] 开车旅行

思路：

首先令 \(nxt1_i\) 表示右侧最近的城市距离（\(id1_i\) 为编号），令 \(nxt2_i\) 表示右侧第二近的城市编号（\(id2_i\) 为编号）；可以使用 set 找出离这个城市最近的 \(4\) 个城市（前面两个，后面两个）。

定义：

• \(f_{i,j}\) 表示从 \(i\) 点出发走 \(2^j\) 轮最后到达的位置。

• \(dp1_{i,j}\) 表示从 \(i\) 点出发走 \(2^j\) 轮最后 A 走过的距离。

• \(dp2_{i,j}\) 表示从 \(i\) 点出发走 \(2^j\) 轮最后 B 走过的距离。

初始化：

\[f_{i,0} = id1_{id2_i}
\]

\[dp1_{i,0} = nxt2_{i}
\]

\[dp2_{i,0} = nxt1_{id2_i}
\]

状态转移方程为：

\[f_{i,j} = f_{f_{i,j-1},j-1}
\]

\[dp1_{i,j} = dp1_{i,j-1} + dp1_{f_{i,j-1},j-1}
\]

\[dp2_{i,j} = dp2_{i,j-1} + dp2_{f_{i,j-1},j-1}
\]

此时对于询问 \(1\) 和询问 \(2\)：

• 本质上是求出从每个城市出发后 \(A\) 走的距离与 \(B\) 走的距离。

• 那么考虑从高位到低位贪心，即设当前跳到了 \(s\) 点，若 \(dp1_{s,i} + dp2_{s,i} \le x\)，可以从 \(s\) 跳到 \(f_{s,i}\)，需要令 \(x \gets x - (dp1_{s,i} + dp2_{s,i})\)，然后继续遍历 \(i-1\) 位。

• 因为是 A 先开车，所以 A 可能会在最后一轮结束后还能再开上一次，需要特判。

时间复杂度为 \(O((N+Q) \log N)\)。

完整代码：

#include<bits/stdc++.h>
#define Add(x,y) (x+y>=mod)?(x+y-mod):(x+y)
#define lowbit(x) x&(-x)
#define pi pair<ll,ll>
#define pii pair<ll,pair<ll,ll>>
#define iip pair<pair<ll,ll>,ll>
#define ppii pair<pair<ll,ll>,pair<ll,ll>>
#define fi first
#define se second
#define full(l,r,x) for(auto it=l;it!=r;it++) (*it)=x
#define Full(a) memset(a,0,sizeof(a))
#define open(s1,s2) freopen(s1,"r",stdin),freopen(s2,"w",stdout);
using namespace std;
typedef double db;
typedef unsigned long long ull;
typedef long long ll;
bool Begin;
const ll N=1e5+10,M=18,INF=1e18;
inline ll read(){
    ll x=0,f=1;
    char c=getchar();
    while(c<'0'||c>'9'){
        if(c=='-')
          f=-1;
        c=getchar();
    }
    while(c>='0'&&c<='9'){
        x=(x<<1)+(x<<3)+(c^48);
        c=getchar();
    }
    return x*f;
}
inline void write(ll x){
	if(x<0){
		putchar('-');
		x=-x;
	}
	if(x>9)
	  write(x/10);
	putchar(x%10+'0');
}
struct Node{
	ll a,b;
	ll id;
	bool operator<(const Node&rhs)const{
		if(a!=rhs.a)
		  return a<rhs.a;
		return b<rhs.b;
	}
};
ll n,m,s,x,x0,s0,cnt,dis1,dis2,disa,disb;
ll a[N],nxt1[N],nxt2[N],id1[N],id2[N];
ll f[N][M],dp1[N][M],dp2[N][M];
Node h[N];
vector<Node> V;
multiset<Node> S;
void solve(ll s,ll x){
	dis1=dis2=0;
	for(int i=M-1;i>=0;i--){
		if(dp1[s][i]+dp2[s][i]<=x&&f[s][i]){
			dis1+=dp1[s][i],dis2+=dp2[s][i];
			x-=dp1[s][i]+dp2[s][i];
			s=f[s][i];
		}
	}
	if(nxt2[s]<=x)
	  dis1+=nxt2[s];
}
bool End;
int main(){
	n=read();
	for(int i=1;i<=n;i++){
		a[i]=read();
		h[i]={a[i],a[i],i};
	}
	S.insert({INF,INF,0});
	S.insert({INF-1,INF-1,0});
	S.insert({-INF,-INF,0});
	S.insert({-INF+1,-INF+1,0});
	for(int i=n;i>=1;i--){
		V.clear();
		V.push_back(*--S.lower_bound(h[i]));
		V.push_back(*--S.lower_bound(V[0]));
		V.push_back(*S.upper_bound(h[i]));
		V.push_back(*S.upper_bound(V[2]));
		for(auto &v:V)
		  v.a=abs(h[i].a-v.a);
		sort(V.begin(),V.end());
		nxt1[i]=V[0].a,nxt2[i]=V[1].a;
		id1[i]=V[0].id,id2[i]=V[1].id;
		cerr<<id1[i]<<' '<<id2[i]<<'\n';
		S.insert(h[i]);
	}
	for(int i=1;i<=n;i++){
		f[i][0]=id1[id2[i]];
		dp1[i][0]=nxt2[i];
		dp2[i][0]=nxt1[id2[i]];
	}
	for(int j=1;j<M;j++){
		for(int i=n;i>=1;i--){
			f[i][j]=f[f[i][j-1]][j-1];
			dp1[i][j]=dp1[i][j-1]+dp1[f[i][j-1]][j-1];
			dp2[i][j]=dp2[i][j-1]+dp2[f[i][j-1]][j-1];
		}
	}
	x0=read();
	for(int i=1;i<=n;i++){
		solve(i,x0);
		if(dis2&&(!s0||disa*dis2>disb*dis1)){
			s0=i;
			disa=dis1,disb=dis2;
		}
	}
	write(s0);
	putchar('\n');
	m=read();
	while(m--){
		s=read(),x=read();
		solve(s,x);
		write(dis1);
		putchar(' ');
		write(dis2);
		putchar('\n');
	}
	cerr<<'\n'<<abs(&Begin-&End)/1048576<<"MB";
	return 0;
}