[dp]答疑

题目描述

在一个数轴上,有 n n n 个同学等待老师的答疑。老师最先站在 k k k 号同学的位置。老师需要解决所有同学的问题,解决一个问题都只需要 0 0 0 秒。学生很没有耐心,每等一秒钟就会增加 w i w_i wi 的不开心值。老师以 1 m / s 1m/s 1m/s 的速度开始行动,他想知道解决完所有同学的问题最少需要多少的不开心值。

输入格式

第一行输入两个整数 n n n 和 k k k,表示等待答疑的学生数量和老师最先站在几号同学的旁边。

接下来 n n n 行,每行两个整数 x x x 和 y y y,表示第 i i i 个同学的位置和每秒增加的不开心值。

输出格式

输出一行一个整数表示最小的不开心值。

数据范围

对于 40 % 40\% 40% 的数据, 2 ≤ n ≤ 7 2 \le n \le 7 2≤n≤7。

对于 100 % 100\% 100% 的数据, 2 ≤ n ≤ 1000 2 \le n \le 1000 2≤n≤1000, 1 ≤ k ≤ 1000 1 \le k \le 1000 1≤k≤1000 0 ≤ x , y ≤ 1000 0 \le x,y \le 1000 0≤x,y≤1000。

样例

样例输入1

4 3
2 2
5 8
6 1
8 7

样例输出1

56

题解

1

对于 40 % 40\% 40% 的数据,我们可以直接暴力枚举,枚举到每个学生的顺序,计算出不开心值。

进而,我们可以将枚举进化,可以得到 60 60 60 分。

我们还是 dfs 出每个位置,边枚举边计算不开心值,进行剪枝。

cpp 复制代码
在 dfs 前,可以先排序,在实践中发现排序后一般比不排序快
x 表示解决的数量,z 表示现在的不开心值(即每秒增加的不开心),s 表示不开心值,o 记录上一次走到的位置
inline void dfs(int x, int z, int s, int o){
	//剪枝
	if(s >= ans){
		return;
	}
	if(x > n){
		ans = min(ans, s);
		return;
	}
	for(int i = o - 1; i >= 1; -- i){
		if(!f[i]){
			f[i] = 1;
			dfs(x + 1, z - u[i].b, s + abs(u[o].a - u[i].a) * z, i);
			f[i] = 0;
		}
	}
	for(int i = o + 1; i <= n; ++ i){
		if(!f[i]){
			f[i] = 1;
			dfs(x + 1, z - u[i].b, s + abs(u[o].a - u[i].a) * z, i);
			f[i] = 0;
		}
	}
}

这里 dfs 也可以优化到 100 100 100 分,跟下面做法类似。

2

考虑使用 区间dp,记录从 i i i 到 j j j 的最小不开心值( i ≤ j i \le j i≤j)。

从 i i i 到 j j j 的区间,有两种情况:

  1. 从 i i i 开始到 j j j,到 j j j 停止。
  2. 从 j j j 开始到 i i i,到 i i i 停止。

定义 f i , j , 0 / 1 f_{i,j,0/1} fi,j,0/1 表示从 i i i 到 j j j 和从 j j j 到 i i i 的最小不开心值。

最开始 f k , k , 0 f_{k,k,0} fk,k,0 和 f k , k , 1 f_{k,k,1} fk,k,1 都为 0 0 0(因为老师最开始就在 k k k 同学位置),其他位置都为无穷大。

考虑转移 f i , j , 0 f_{i,j,0} fi,j,0,它可以从 f i , j − 1 , 0 f_{i,j-1,0} fi,j−1,0 转移过来,即从 i i i 到 j j j,可以先从 i i i 到 j − 1 j - 1 j−1,再从 j − 1 j - 1 j−1 到 j j j;也可以从 f i , j − 1 , 1 f_{i,j-1,1} fi,j−1,1 转移过来,即从 j − 1 j - 1 j−1 到 i i i,再从 i i i 到 j j j。尽管第二个方法看起来不是很优,但不开心值也可能比第一种小。

接下来就是计算 w w w 了。可以用一个前缀和数组 s u m sum sum 记录从 w 1 w_1 w1 到 w i w_i wi 的和。这下从 i i i 到 j j j 的区间,不开心值就是排除 i i i 到 j j j 的区间和,这里 j j j 是将要到的,并没有到,所以不计 j j j 的不开心值。

因此,我们就可以得出 f i , j , 0 f_{i,j,0} fi,j,0 的转移。
f i , j , 0 = min ⁡ ( f i , j , 0 , f i + 1 , j , 0 + ( ∣ x j + 1 . a − x j . a ∣ ) × ( s u m n − s u m k + s u m j ) , f j + 1 , k , 1 + ( ∣ x k . a − x j . a ∣ ) × ( s u m n − s u m k + s u m j ) ) f_{i,j,0} = \min(f_{i,j,0}, f_{i + 1, j, 0} + (\mid x_{j + 1}.a - x_{j}.a\mid) \times (sum_n - sum_k + sum_j), f_{j + 1, k, 1} + (\mid x_{k}.a - x_{j}.a \mid) \times (sum_{n} - sum_{k} + sum_{j})) fi,j,0=min(fi,j,0,fi+1,j,0+(∣xj+1.a−xj.a∣)×(sumn−sumk+sumj),fj+1,k,1+(∣xk.a−xj.a∣)×(sumn−sumk+sumj))

同理,我们可以的出 f i , j , 1 f_{i,j,1} fi,j,1 的转移。

答案就是 min ⁡ ( f 1 , n , 0 , f 1 , n , 1 ) \min(f_{1,n,0}, f_{1,n,1}) min(f1,n,0,f1,n,1) 的值。

cpp 复制代码
//前缀和
for(int i = 1; i <= n; ++ i){
	sum[i] = sum[i - 1] + x[i].b;
} 
//预处理 dp 数组
for(int i = 0; i <= n + 1; ++ i){
	for(int j = 0; j <= n + 1; ++ j){
		f[i][j][0] = f[i][j][1] = 1e16;
	}
}
f[m][m][0] = f[m][m][1] = 0;

//转移
for(int i = 1; i <= n; ++ i){
	for(int j = 1; j <= n - i; ++ j){
		//这里的 j 就是上文的 i, k 就是上文的 j
		int k = i + j;
		f[j][k][0] = min(f[j][k][0], min(f[j + 1][k][0] + abs(x[j + 1].a - x[j].a) * (sum[n] - sum[k] + sum[j]), f[j + 1][k][1] + abs(x[k].a - x[j].a) * (sum[n] - sum[k] + sum[j])));
		f[j][k][1] = min(f[j][k][1], min(f[j][k - 1][1] + abs(x[k].a - x[k - 1].a) * (sum[n] - sum[k - 1] + sum[j - 1]), f[j][k - 1][0] + abs(x[k].a - x[j].a) * (sum[n] - sum[k - 1] + sum[j - 1])));
	}
}

禁止抄袭!!!

相关推荐
咒法师无翅鱼20 分钟前
【定理证明工具调研】Coq, Isabelle and Lean.
算法
风清云淡_A1 小时前
【java基础系列】实现数字的首位交换算法
java·算法
涵涵子RUSH1 小时前
合并K个升序链表(最优解)
算法·leetcode
爱吃西瓜的小菜鸡1 小时前
【C语言】矩阵乘法
c语言·学习·算法
染指11102 小时前
50.第二阶段x86游戏实战2-lua获取本地寻路,跨地图寻路和获取当前地图id
c++·windows·lua·游戏安全·反游戏外挂·游戏逆向·luastudio
sjsjs112 小时前
【多维DP】力扣3122. 使矩阵满足条件的最少操作次数
算法·leetcode·矩阵
哲学之窗2 小时前
齐次矩阵包含平移和旋转
线性代数·算法·矩阵
Code out the future2 小时前
【C++——临时对象,const T&】
开发语言·c++
Sudo_Wang2 小时前
力扣150题
算法·leetcode·职场和发展
sam-zy2 小时前
MFC用List Control 和Picture控件实现界面切换效果
c++·mfc