「NOIP 2012 D1T3」开车旅行

传送门


problem

在这里插入图片描述
数据范围: 1 ≤ n ≤ 1 0 5 1≤n≤10^5 1n105 1 ≤ m ≤ 1 0 4 1≤m≤10^4 1m104


solution

感觉这道题挺妙的。

首先要解决的一个问题是,怎么找到题目中定义的最近城市次近城市

一个思路就是从 n ∼ 1 n\sim 1 n1 倒着扫,处理到 i i i 时,我们想知道 [ i + 1 , n ] [i+1, n] [i+1,n] 中, h i h_i hi前驱,后继,前驱的前驱,后继的后继,因为只最近和次近肯定在这 4 4 4 个点中,排序判一下具体是哪两个即可。

由于我们要支持插入查询前驱后继这些操作,可以很方便地用 set 维护。

接下来就是询问了。两种询问实质上是一样的,因为第一个询问可以通过枚举每一个点转换成第二个询问。

暴力的做法就是直接模拟,先走一步 A,然后走 B,再走 A ⋯ \cdots ,直到不能走为止。

怎么优化这个暴力呢,不妨考虑倍增A,B 各走一次为一轮,暴力跳看能走几轮,最后特判一下 A 是否可以再走一次即可。

时间复杂度 O ( n log ⁡ n ) O(n\log n) O(nlogn)

PS:代码中倍增是以 1 1 1 为下标的, 0 0 0 记的是 A 跳一次的信息。


code

#include<bits/stdc++.h>
#define ll long long
using namespace std;
namespace IO{
	const int Rlen=1<<22|1;
	char buf[Rlen],*p1,*p2;
	inline char gc(){
		return (p1==p2)&&(p2=(p1=buf)+fread(buf,1,Rlen,stdin),p1==p2)?EOF:*p1++;
	}
	template<typename T>
	inline T Read(){
		char c=gc();T x=0,f=1;
		while(!isdigit(c))  f^=(c=='-'),c=gc();
		while( isdigit(c))  x=(x+(x<<2)<<1)+(c^48),c=gc();
		return f?x:-x;
	}
	inline int in()  {return Read<int>();}
}
using IO::in;
const int N=1e5+5;
int n,m,pos,X0,h[N];
ll disA[N],disB[N],posA[N],posB[N],inf=1e15;
ll p[N][21],A[N][21],B[N][21];
set<ll>S;
unordered_map<ll,int>Map;
struct node{ll h,dis;}tmp[5];
bool operator<(const node &p,const node &q)  {return (p.dis==q.dis)?p.h<q.h:p.dis<q.dis;}
ll Abs(ll x)  {return x>0?x:-x;}
void prework(){
	S.insert(inf),S.insert(inf-1),S.insert(-inf),S.insert(-inf+1);
	for(int i=n;i>=1;--i){
		S.insert(h[i]);
		tmp[1].h=*--S.lower_bound(h[i]);
		tmp[2].h=*--S.lower_bound(tmp[1].h);
		tmp[3].h=*S.upper_bound(h[i]);
		tmp[4].h=*S.upper_bound(tmp[3].h);
		for(int j=1;j<=4;++j)  tmp[j].dis=Abs(tmp[j].h-h[i]);
		sort(tmp+1,tmp+5);
		disA[i]=tmp[2].dis,posA[i]=Map[tmp[2].h];
		disB[i]=tmp[1].dis,posB[i]=Map[tmp[1].h];
		if(posA[i])  A[i][0]=disA[i],p[i][0]=posA[i];
		if(posB[posA[i]])  A[i][1]=disA[i],B[i][1]=disB[posA[i]],p[i][1]=posB[posA[i]];
		for(int j=2;j<=20;++j){
			p[i][j]=p[p[i][j-1]][j-1];
			A[i][j]=A[i][j-1]+A[p[i][j-1]][j-1];
			B[i][j]=B[i][j-1]+B[p[i][j-1]][j-1];
		}
	}
}
void calc(int S,int X,ll &ansA,ll &ansB){
	ansA=ansB=0;
	for(int i=20;i>=1;--i){
		if(p[S][i]&&ansA+ansB+A[S][i]+B[S][i]<=X){
			ansA+=A[S][i];
			ansB+=B[S][i];
			S=p[S][i];
		}
	}
	if(p[S][0]&&ansA+ansB+A[S][0]<=X)  ansA+=A[S][0];
}
ll ansA,ansB;
int main(){
	n=in();
	for(int i=1;i<=n;++i)  h[i]=in(),Map[h[i]]=i;
	prework(),X0=in(),m=in();
	double Min=1e18;
	for(int i=1;i<=n;++i){
		calc(i,X0,ansA,ansB);
		double now=ansB?1.0*ansA/ansB:1e18;
		if(Min>now)  Min=now,pos=i;
	}
	printf("%d\n",pos);
	while(m--){
		int S=in(),X=in();
		calc(S,X,ansA,ansB),printf("%lld %lld\n",ansA,ansB);
	}
	return 0;
}
  • 0
    点赞
  • 0
    收藏
    觉得还不错? 一键收藏
  • 0
    评论
你好!感谢你的提问。根据你的要求,我将为你解答关于NOIP2017普及t3问题的内容。 在这道题目中,我们需要解决一个关于棋盘的问题。题目描述了一个n×m的棋盘,每个格子上有一个非负整数。开始时,小明站在左上角的格子上,他可以向右或向下移动一步,每次移动的代价为目标格子上的数值。我们需要找到一条从左上角到右下角的路径,使得移动的总代价最小。 解决这个问题的一个常见的方法是使用动态规划(Dynamic Programming)。我们可以创建一个二维数组dp,其中dp[i][j]表示从起点到达坐标为(i, j)的格子时的最小代价。然后,我们可以按照从左上角到右下角的顺序依次计算dp数组的值。 具体的计算方法如下: 1. 首先,我们可以初始化dp数组的第一行和第一列,即dp[0][j]和dp[i][0],它们表示从起点到达第一行和第一列的格子时的最小代价。初始化的方法是累加前面的格子的代价。 2. 接下来,我们可以使用一个双重循环,从(1, 1)开始遍历整个棋盘。对于每个格子(i, j),我们可以选择从上方格子(i-1, j)或左方格子(i, j-1)中选择一个代价较小的路径,并加上当前格子的代价。即dp[i][j] = min(dp[i-1][j], dp[i][j-1]) + grid[i][j]。 3. 最后,当我们计算完dp数组的所有值后,dp[n-1][m-1]即为从起点到达右下角的格子时的最小代价。 这样,我们就可以得到从左上角到右下角的最小代价。希望能对你的问题有所帮助!如果你还有其他问题,请随时提问。

“相关推荐”对你有帮助么?

  • 非常没帮助
  • 没帮助
  • 一般
  • 有帮助
  • 非常有帮助
提交
评论
添加红包

请填写红包祝福语或标题

红包个数最小为10个

红包金额最低5元

当前余额3.43前往充值 >
需支付:10.00
成就一亿技术人!
领取后你会自动成为博主和红包主的粉丝 规则
hope_wisdom
发出的红包
实付
使用余额支付
点击重新获取
扫码支付
钱包余额 0

抵扣说明:

1.余额是钱包充值的虚拟货币,按照1:1的比例进行支付金额的抵扣。
2.余额无法直接购买下载,可以购买VIP、付费专栏及课程。

余额充值