Week10 - 程序设计思维与实践 - 动态规划（一）

A - 签到题

题意

东东在玩游戏“Game23”。

在一开始他有一个数字n，他的目标是把它转换成m，在每一步操作中，他可以将n乘以2或乘以3，他可以进行任意次操作。输出将n转换成m的操作次数，如果转换不了输出-1。

Input

输入的唯一一行包括两个整数n和m（1<=n<=m<=5*10^8).

Output

输出从n转换到m的操作次数，否则输出-1.

Sample Input

120 51840

Sample Output

7

代码实现

#include <iostream>
#include <cstdio>
#include <cstdlib>
using namespace std;

bool flag=false;
int n,m,ans=1e9;

void dfs(int x,int steps){
	if(x==m){
		ans=min(ans,steps);
		flag=true;
		return;
	}
	else if(x>m)
		return;
	else{
		dfs(x*2,steps+1);
		dfs(x*3,steps+1);
	}
}

int main()
{
	scanf("%d%d",&n,&m);
	dfs(n,0);
	if(flag)
		printf("%d\n",ans);
	else
		printf("-1\n");
	return 0;
}

B - LIS & LCS

题意

东东有两个序列A和B。

他想要知道序列A的LIS和序列AB的LCS的长度。

注意，LIS为严格递增的，即a1<a2<…<ak(ai<=1,000,000,000)。

Input

第一行两个数n，m（1<=n<=5,000,1<=m<=5,000）
第二行n个数，表示序列A
第三行m个数，表示序列B

Output

输出一行数据ans1和ans2，分别代表序列A的LIS和序列AB的LCS的长度

Sample Input

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

Sample Output

3 2

思路

最长上升子序列（LIS）

• 设 dp1[i] 是以第 i 个数为终点的最长上升子序列
• 初始状态 dp1[1]=1
• dp1[k] = max{ dp1[i] : 1≤i<k 且 ai<ak 且 k≠1 }+1
• 若找不到这样的 k 则 dp1[k]=1

最长公共子序列（LCS）

• 设 dp2[i][j] 表示 a 左边 i 个字符形成的子串，与 b 左边 j 个字符形成的子串的最长公共子序列的长度
• dp2[n][0] = 0 ( n = 0… len1)
• dp2[0][n] = 0 ( n = 0… len2)
• 若a[i]=b[j] 则 dp2[i][j]=dp2[i-1][j-1]+1;
• 否则 dp2[i][j] = max(dp2[i][j-1], dp2[i-1][j]);

代码实现

#include <iostream>
#include <cstdio>
#include <cstdlib>
using namespace std;
const int maxn=5010;
int n,m,dp1[maxn],dp2[maxn][maxn],a[maxn],b[maxn];

int main()
{
	scanf("%d%d",&n,&m);
	for(int i=1;i<=n;i++){
		scanf("%d",&a[i]);
		dp1[i]=1;
		dp2[i][0]=0;
	}	
	for(int i=1;i<=m;i++){
		scanf("%d",&b[i]);
		dp2[0][i]=0;
	}
	dp2[0][0]=0;	
	//LIS
	int ans1=1;
	for(int i=2;i<=n;i++){
		for(int j=1;j<i;j++){
			if(a[i]>a[j]){
				dp1[i]=max(dp1[i],dp1[j]+1);
				ans1=max(ans1,dp1[i]);
			}
		}
	}
	//LCS
	for(int i=1;i<=n;i++)
		for(int j=1;j<=m;j++)
			if(a[i]==b[j])
				dp2[i][j]=dp2[i-1][j-1]+1;
			else
				dp2[i][j]=max(dp2[i][j-1],dp2[i-1][j]);
	printf("%d %d\n",ans1,dp2[n][m]);
	return 0;
}

C - 拿数问题 II

题意

给一个序列，里边有 n 个数，每一步能拿走一个数，比如拿第 i 个数， Ai = x，得到相应的分数 x，但拿掉这个 Ai 后，x+1 和 x-1 (如果有 Aj = x+1 或 Aj = x-1 存在) 就会变得不可拿（但是有 Aj = x 的话可以继续拿这个 x）。求最大分数。

Input

第一行包含一个整数 n (1 ≤ n ≤ 10^5)，表示数字里的元素的个数

第二行包含n个整数a1, a2, …, an (1 ≤ ai ≤ 10^5)

Output

输出一个整数：n你能得到最大分值。

Sample Input

9
1 2 1 3 2 2 2 2 3

Sample Output

10

思路

和课上讲的 拿数问题Ⅰ 的状态转移方程类似，不过要先进行一步转换，即将相同项合并然后从小到大进行排序，这样就只需要考虑取了 x，不能取 x+1 的情况

• 设 dp(i) 表示以 b[i] 结尾的最大分值
• 若 b[i] = b[i-1]+1则 dp[i]=max(dp[i-1],dp[i-2]+b[i]*mp[b[i]]);
• 否则 dp[i]=dp[i-1]+b[i]*mp[b[i]];

代码实现

#include <iostream>
#include <cstdio>
#include <cstdlib>
#include <cstring>
#include <algorithm>
#include <map>
using namespace std;
typedef long long ll;
const int maxn=100010;
int n;
ll a[maxn],b[maxn],dp[maxn];
map<ll,ll> mp;

int main()
{
	scanf("%d",&n);
	for(int i=0;i<n;i++){
		scanf("%lld",&a[i]);
		mp[a[i]]++;
	}
	memset(dp,0,sizeof(dp));
	memset(b,0,sizeof(b)); 
	int idx=0;
	for(map<ll,ll>::iterator it=mp.begin();
			it!=mp.end();it++,idx++){
		b[idx]=it->first;
	}
	dp[0]=b[0]*mp[b[0]];
	dp[1]=(b[0]+1==b[1])?max(b[0]*mp[b[0]],b[1]*mp[b[1]]):(b[0]*mp[b[0]]+b[1]*mp[b[1]]);
	for(int i=2;i<idx;i++){
		if(b[i]==b[i-1]+1){
			dp[i]=max(dp[i-1],dp[i-2]+b[i]*mp[b[i]]);
		}
		else{
			dp[i]=dp[i-1]+b[i]*mp[b[i]];
		}
	}
	printf("%lld\n",dp[idx-1]);
	return 0;
}