4245. 【五校联考6day2】er

Description

小明在业余时间喜欢打电子游戏，不是星际和魔兽这些，是赛尔号一类的游戏。最近小明在玩一款新出的游戏，叫做■■■■■■■■。小明觉得游戏里自己的装备太垃圾了，每次都被大神虐，一怒之下充了■■元准备强化装备。
这个游戏中用于强化装备的道具叫做强化符文。有以下3 种：
1. 赋值强化符文，对某个装备使用这个符文以后，装备威力值会变为一个常数。因为这个功能很IMBA，可以让一个垃圾装备变得非常牛■，所以它在游戏里很稀有，市场上最多能见到一个。
2. 加法强化符文，对某个装备使用后，威力值加上一个常数。
3. 乘法强化符文，对某个装备使用后，威力值乘上一个常数。
市场上有M 个不同强化符文可以购买，小明有N 件装备准备强化，他只能购买K 个强化符文，然后以任意顺序使用这些符文，强化他的任意装备（一个装备可以不强化也可以强化多次）。根据游戏的设定，所有装备威力值乘积为总效果。请为他设计一个购买和强化方案，使强化后的所有装备总效果值最大。
由于小明RP 不太好，打BOSS 都不掉神装，所以他的装备不超过两件。

Input

第一行3 个正整数N;M;K, 含义见题面。
第二行N 个正整数Ai，表示他的每个装备的初始威力值。
第三行开始共M 行，每行两个正整数Type_i;Ci，描述一个强化符文。Type_i表示符文类型，1 表示赋值，2 表示加法，3 表示乘法。Ci 是对应的常数值。

Output

一个数，表示最大的总效果值。由于这个数可能很大，请输出它的自然对数，保留3 位小数。

Sample Input

2 5 3
0 1
2 3
2 1
2 4
3 4
3 2

Sample Output

4.159

Data Constraint

对于20% 的数据，N = 1；
对于全部数据M,K ≤ 100;N ≤ 2，最多一个Type_i = 1。
输入数据中所有数不超过2000。

Solution

贪心+DP。

首先我们根据贪心策略可以得出肯定先用赋值强化符文，再用加法强化符文，最后用乘法强化符文是最优的。

我们考虑n=1的情况。那么赋值强化符文可以看做加法，若要a1变成b实际上就可以看成将a1加上了(b-a1)，所以可以跟加法一起做。对于加法符文，我们肯定优先选大的，所以可以将加法符文按从大到小排序后做一个前缀和。对于乘法符文则排序后做一个前缀积。然后我们只要枚举一个i(可为0)，使得（a1+i个最大的加法符文）*（k-i个乘积最大的乘法符文）。

当n=2时，就要用dp做了。先不考虑赋值强化符文的情况。可以发现其实只要使a1,a2以及用了i个加法符文就可以算出答案了(答案为a1*a2*使用的乘法符文)。所以我们可以设f[ i ][ j ]表示选了前i个加法符文，a1的值为 j 时，a2的最大值。那么f [ i ][ j ]转移时就会有两种情况。

• 将第i个符文给a1用。f [ i ][ j ]=max{ f [ i - 1 ][ j ] - 加法[ i ] ] }
• 将第i个符文给a2用。f [ i ][ j ]=max{ f [ i - 1 ][ j ]+加法[ i ] }

最后答案就是max{ j * f [ i ][ j ] * k-i个乘积最大的乘法符文 }。

特别特别要注意边界条件！！！因为有关精度尤其还是这神马辣鸡自然对数问题时就会有各种神奇问题影响答案。

那么如果有赋值强化符文的时候呢？其实可以将a1和a2中更小的那个直接变成赋值强化符文，表示直接使用了一个符文，然后再用新的a1,a2做一遍上述dp就可以了。

最后，因为答案会很大，所以出题人良心不让我们大高精度，改成了求答案的自然对数，然而你看了超链接也没懂，简单来说就是 2.718...^a=b,我们就说b的自然对数为a(可能为小数)，实际上就是C++的 log(x) 运算，同时，因为不能将答案算出来再求自然对数，我们有一个性质，log(a*b)=log(a)+log(b)。证明可以自行查资料，这里就不再证明了。通过这个性质我们就可以在算出加法符文的贡献之后直接将答案加起来而不用乘了。

Code

#include<cstdio>
#include<cmath>
#include<cstring>
#include<algorithm>
#define ll long long
#define ld long double
using namespace std;
int n,m,ma,k,a1,a2,x,y,f[110][202010],t1,s1,s2,s3,t2[110],t3[110];
ld ans,s[110];
bool cmp(ll x,ll y){return x>y;}
void dp(int a1,int a2){
	memset(f,129,sizeof(f));
	f[0][a1]=a2;
	ma=a1;
	for(int i=1;i<=s2;i++){
		if(i>k) break;
		for(int j=0;j<=ma+t2[i];j++){
			f[i][j]=max(f[i][j],f[i-1][j]+t2[i]);
			if(j>=t2[i]) f[i][j]=max(f[i][j],f[i-1][j-t2[i]]);
			if(f[i][j]>0) ans=max(ans,(ld)log(j)+(ld)log(f[i][j])+s[k-i]);
		}
		ma+=t2[i];
	}
	ans=max(ans,log(a1)+log(a2)+s[k]);
}
int main(){
	scanf("%d%d%d",&n,&m,&k);
	if(n==1) scanf("%d",&a1);
	else scanf("%d%d",&a1,&a2);
	for(int i=1;i<=m;i++){
		scanf("%d%d",&x,&y);
		if(x==1){t1=y;s1=1;}
		if(x==2){t2[++s2]=y;}
		if(x==3){t3[++s3]=y;}
	}
	sort(t2+1,t2+1+s2,cmp);
	sort(t3+1,t3+1+s3,cmp);
	if(n==1){
		for(int i=1;i<=s2;i++) t2[i]=t2[i-1]+t2[i];
		for(int i=1;i<=s3;i++) s[i]=s[i-1]+log(t3[i]);
		for(int i=0;i<=min(k,s2);i++) ans=max(ans,log(a1+t2[i])+s[k-i]);
	}
	else{
		for(int i=1;i<=s3;i++) s[i]=s[i-1]+log(t3[i]);
		dp(a1,a2);
		if(s1){
			k--;
			if(a1<a2) dp(t1,a2);
			else dp(a1,t1);
		}
	}
	printf("%.3Lf",ans);
	return 0;
}

---

---

作者：zsjzliziyang 
QQ:1634151125 
转载及修改请注明 
本文地址:https://blog.csdn.net/zsjzliziyang/article/details/86695662