首页 > 其他分享 >洛谷 P1336 最佳课题选择 题解

洛谷 P1336 最佳课题选择 题解

时间:2023-08-07 16:45:20浏览次数:40  
标签:背包 洛谷 int 题解 LL array cost P1336 转移

P1336 最佳课题选择 题解

状态:考虑\(f_{i,j}\)表示前\(i\)种论文里面,一共写了\(j\)篇,的最少花费时间。

转移策略:我们一次考虑每一种论文写多少篇。假设写\(k\)篇,\(k \in [0,j] \cap \mathbb{Z}\) ,有转移方程:

\[f_{i,j} = min(f_{i-1,j-k} + cost(i,k)), k \in [0,j] \cap \mathbb{Z} \]

其中

\[cost(i,k) = a_i \times k ^ {b_i} \]

可以做记忆化优化,因为多次访问同一个\(cost(i,k)\)。

初始化:我们考虑选0篇无论怎样花费为0,选0种论文但多于0篇一定是不可能的,花费为Inf,再加上转移中要取min的缘故,我们初始化为:

\[\begin{equation*} f_{i,j}=\left\{ \begin{array}{} 0 & ,j=0\\ + \inf &, j \neq 0\\ \end{array} \right . \end{equation*} \]

相关变量范围见代码。

最后,我们就可以开滚动数组,或者倒序\(j\)优化掉一维。

代码在文末。


这里和01背包及无限背包做个对比。

01背包的转移选择策略,是一个一个选(同时也是一种一种选,因为一种里面只有一个)。转移方程:

\[f_{i,j} = max(f_{i-1,j}, f_{i-1,j-c_i} + v_i) \]

无限背包的转移选择策略,是在一种里面一个个选。转移方程:

\[f_{i,j} = max(f_{i-1,j}, f_{i,j-c_i} + v_i) \]

可以看到两者细微的差别在,选物品的时候,01背包访问\(i-1\),无限背包访问\(i\),所以可以把其中一个维度压缩掉之后,01反着做,无限正着做。

本题目有两点不同:

第一,本题相当于“选择超过某价值的物品,使得容量需要量最小”,把状态的内外对换了一对。

第二,本题目中,如果选了前面一个,那么后面再选的时候,贡献会不一样,因为时间贡献是幂级的而非线性的。类比到无限背包中就是“选的越多,价值越少”,已经选了多少会影响到后面选的东西。选同一个东西的贡献,在不同情况下不一样。

于是,按照无限背包“一个一个”选的策略,就会涉及到之前的选择情况,我们尝试把第\(i\)种选择数目\(k\)加入到状态中,有

\[\begin{equation*} f_{i,j,k(k \leq j)}=\left\{ \begin{array}{} f_{i,j-1,k-1} - cost(i,k-1) + cost(i,k) & ,k \neq 0\\ min(f_{i-1,j,l}), l \in [0,j] \cap \mathbb{Z} &, k = 0\\ \end{array} \right . \end{equation*} \]

这样是正确的,但是我们发现大多数的状态被浪费了,因为\(k \neq 0\) 的时候只有一种选择。

还有一种方式是把\(f_{i,j}\)中选择了多少个第i种论文,另存为\(g_{i,j} = k\),于是我们有转移方程(错误的):

\[\begin{equation*} f_{i,j}=min\left( \begin{array}{} f_{i,j-1} - cost(i,g_{i,j-1}) + cost(i,g_{i,j}), g_{i,j} = g_{i,j-1} + 1 ,\\ f_{i-1,j}, g(i,j) = 0 \end{array} \right) \end{equation*} \]

这个转移方程看上去好像是对的,因为就是按照无限背包的思路改的,但实际上是错误的,提交后只能有10分。因为理论上这种转移方程只能处理线性的\(cost(i,k)\),即\(b_{i} = 1\)或\(b_{i} = 0\)

我们和无限背包进行对比:假设有\(f_{1,1} = 12\),\(f_{1,2} = 24\),\(f_{2,1} = 2\),那么对于无限背包问题,有:

\[f_{2,2} = max(f_{1,2},f_{2,1} + v_{2})\\ \because f_{2,1} = max(f_{1,1},f_{2,0} + v_{2}) \neq f_{1,1}\\ \therefore f_{2,1} \geq f_{1,1}\\ \therefore f_{2,1} + v_{2} \geq f_{1,1} + v_{2} \]

因此\(f_{1,1} + v_{2}\)不会影响\(f_{2,1} + v_{2}\)对\(f_{2,2}\)可能的贡献,也就是不可能贡献(被max掉了)。

但是在这个问题中,我们令\(w_{i,k} = cost(i,k) - cost(i,k-1)\),有:

\[f_{2,2} = min(f_{1,2},f_{2,1} + w_{2,1})\\ \because f_{2,1} = min(f_{1,1},f_{2,0} + w_{2,2}) \neq f_{1,1}\\ \therefore f_{2,1} \leq f_{1,1}\\ \therefore f_{2,1} + w_{2,2} \quad? \quad f_{1,1} + w_{2,1} \]

可以发现我们无法判断\(f_{1,1} + w_{2,1}\)对于\(f_{2,2}\)是否有贡献,其根源在于对于同一个\(i\),\(w_{i,k}\)是可能不同的,因此不等式失效,即在这种状态和转移方程下,不满足最优子结构。因此DP失效。

(这也说明,一个可DP的问题是存在“最优子结构”“无后效性”,只有特定的方式去“成立”这些性质才能DP,不是随便搞一搞就一定可以DP的)

也就是说,在无限背包问题中,我们把一个物品拿走,从规模更小的状态中转移,这个物品的贡献是确定的,决定“拿走”这一选择中的最优解,就只能是承接之前的最优解。

然而在这个问题中,如果我们把一个物品拿走,从规模更小的状态中转移时,这个物品的贡献是不确定的,也就是说除了最优解以外,次有解也有可能贡献出更好的结果,因此“最优子结构失效了”。

因此,我们在考虑状态和状态转移方程时,可以灵活更改状态、转移方式和策略、记录信息等等,同时也要按照普适状况考虑转移的合法性和具体转移方法。


标准AC代码:

#include <bits/stdc++.h>

#define N (int)(205)
#define M (int)(25)

using namespace std;

typedef long long LL;

int n,m;
LL f[N];
LL a[M],b[M];
LL p[M][N];

LL cost(int i, int k)
{
	if(p[i][k] == -1)
		return p[i][k] = a[i] * pow(k,b[i]);
	else
		return p[i][k];
}

int main()
{
	memset(p,-1,sizeof(p));
	memset(f,0x7f,sizeof(f));
	ios::sync_with_stdio(false);
	cin.tie(0);cout.tie(0);
	cin >> n >> m;
	for(int i = 1;i <= m;i++) cin >> a[i] >> b[i];
	f[0] = 0;
	for(int i = 1;i <= m;i++)
	{
		for(int j = n;j >= 1;j--)
		{
			for(int k = 0;k <= j;k++)
			{
				f[j] = min(f[j],f[j-k] + cost(i,k));
			}
		}
	}
	cout << f[n];
	return 0;
}

另一种AC的代码:

#include <bits/stdc++.h>

#define N (int)(205)
#define M (int)(25)

using namespace std;

typedef long long LL;

int n,m;
LL f[M][N][N];
LL a[M],b[M];
LL p[M][N];

LL cost(int i, int k)
{
	if(p[i][k] == -1)
		return p[i][k] = a[i] * pow(k,b[i]);
	else
		return p[i][k];
}

int main()
{
	memset(p,-1,sizeof(p));
	memset(f,0x7f,sizeof(f));
	ios::sync_with_stdio(false);
	cin.tie(0);cout.tie(0);
	cin >> n >> m;
	for(int i = 1;i <= m;i++) cin >> a[i] >> b[i];
	for(int i = 0;i <= m;i++) f[i][0][0] = 0;
	for(int i = 1;i <= m;i++)
	{
		for(int j = 1;j <= n;j++)
		{
			for(int l = 0;l <= j;l++)
			{
				f[i][j][0] = min(f[i][j][0],f[i-1][j][l]);
			}
			for(int k = 1;k <= j;k++)
			{
				f[i][j][k] = f[i][j-1][k-1] - cost(i,k-1) + cost(i,k);
			}
		}
	}

	LL ans = 1e18;
	for(int k = 0;k <= n;k++)
		ans = min(ans,f[m][n][k]);
	cout << ans;
	return 0;
}

标签:背包,洛谷,int,题解,LL,array,cost,P1336,转移
From: https://www.cnblogs.com/l-cacherr/p/17611828.html

相关文章

  • 洛谷 P6010 - [USACO20JAN] Falling Portals P
    先考虑怎么对一组询问求解答案。容易想到一种贪心策略:如果\(a_{q_i}<a_i\),那么我们肯定希望自己能够尽量快地下落,因此如果遇到一个下落速度大于当前世界下落速度的传送门我们肯定就会进那个世界,如此走下去知道能够传送到世界\(q_i\)为止。对于\(a_{q_i}>a_i\)的情况也类似,只......
  • RTSP/Onvif视频服务器LntonNVR(源码版)视频平台无法通过Onvif控制摄像头云台的问题解决
    LntonNVR视频边缘计算网关平台是我们推出的软硬一体的视频平台,既有软件版本,又有硬件版本。LntonNVR与摄像头连接时,可以通过平台自带的Onvif探测进行设备探测、连接,还能实现对摄像头的PTZ云台控制,包括镜头转向、变焦等操作。通过Onvif控制云台是非常实用的功能,在很多用户实际项目中......
  • 国标GB28181视频平台LntonGBS(源码版)国标平台出现录像无法播放并存在RTMP重复推流现象
    LntonGBS国标视频云服务通过支持国标GB28181协议,实现了设备接入、实时监控直播、录像、语音对讲、云存储、告警、级联等功能。同时,它还支持将接入的视频流以多种格式(包括RTSP、RTMP、FLV、HLS、WebRTC)进行全终端、全平台分发,实现了无插件播放在Web浏览器、手机浏览器、微信端、PC客......
  • 洛谷 P8500 - [NOI2022] 冒泡排序
    显然将权值离散化是没有问题的,因为必然存在一组最优解,满足每个\(a_i\)都取自于某个\(V_i\),于是不管三七二十一先将\(V_i\)离散化了再说。考虑从部分分入手逐步分析这道题:特殊性质A:\(V_i=1\)相当于这个区间中的数必须是\(1\),先将这些数去掉不管,紧接着考虑\(V_i=0\)的......
  • zookeeper常见问题解决
     注意:自zk3.5.5版本以后,已编译的jar包,尾部有bin标识,应该使用的是apache-zookeeper-3.x.x-bin.tar.gz 错误一:Startingzookeeper…FAILEDTOSTART版本问题,自3.5以上的版本,随着版本的更新,3.5版本以后的压缩包分成了两种我们需要使用文件名带有bin的那个压缩包,例如:ap......
  • [ABC313] C~E 题解
    [ABC313]C~E题解C-ApproximateEqualization2让所有的数字都尽量接近平均数,先算出平均数,然后把所有数字分成两份,一份要加,一份要减,因为平均数有余数,余数肯定给最大的几个,所以这样计算总共需要加减多少个,然后在加减里面取\(\max\)即可。时间复杂度:\(O(n\logn)\)#include......
  • 题解 P6831 - [IOI2020] 嘉年华奖券
    小清新IOI题。首先考虑怎么求出答案。等价于我选择\(\dfrac{nk}{2}\)个数令它们系数为\(1\),再选\(\dfrac{nk}{2}\)个数令它们系数为\(-1\),最大化每个数的值乘以系数之和,并且要求每个奖券选择的数的个数恰好是\(k\)个。考虑先令每个奖券的前\(k\)个数系数为\(-1\),然......
  • HHKB2020 D 题解
    problem&blog。特判一下\(a+b>n\)时为\(0\)。正难则反,计算重叠的方案数。重叠即\(x,y\)轴均有交集,于是转化为两条线段相交的方案数(两条线段认为是不同的)。正难则反,计算不相交的方案数。第一条线段有\((n-b+1)-a+1=n-a-b+2\)种可能,第二条有\((n-a-b+1)\)种可能,故方案......
  • 题解 [POI2012] OKR-A Horrible Poem
    题目链接询问循环节的“模板题”?首先,有一个经典结论:若存在一长度为\(len\)的循环节,则\(s[l\simr-len]=s[l+len\simr]\),简单来说就是利用移位,说明是否是循环节。有了这个结论,再使用字符串哈希,我们就可以做到\(O(1)\)判断。因为:字符串长度=循环节长度\(\times\)循环......
  • [国家集训队] Tree II 题解报告
    [国家集训队]TreeII一道·真·板子·题就是练习LCT懒标记的题目除了翻转标记以外还要维护乘法标记和加法标记注意加法标记和乘法标记的维护!!!加法标记因为splay的区间大小不是固定的,所以我们要维护size,并且子树的sum要加上size乘上标记其他的就只用直接加上即可voidpusha......