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 1：

120 51840

Sample Output 1：

7

Sample Input 2：

42 42

Sample Output 2：

0

Sample Input 3：

48 72

Sample Output 3:

1

题目分析：

该题要判断n是否可以通过乘以2或者乘以3的方式得到m，那么就可以反过来对m进行因式分解，看m是否可以分解成若干个2或者3和n的乘积组合的形式，如果可以那么输出因数的个数，否则说明不能得到输出-1。

• 首先判断m是否能整除n，如果不能，无法分解，直接输出-1。
• 如果能整除，那么令m=m/n,对新的m进行2或者3的因式分解，直到最小单元1或者2或者3，分解完成，输出分解需要次数，如果m既不能整除2又不能整除3，说明无法按题意分解，输出-1。
  （注意m=n的情况，这个时候不需要ans++)

代码：

#include<iostream>
#include<stdio.h>
using namespace std;
bool solve(int m)//判断是不是质数 
{
	int i;
	for(i=2;i<=m;i++)
	{
		if(m%i==0)
		{
			break;
		}
	}
	if(i==m&&m!=1)
	{
		return true;
	}
	return false; 
}
int main()
{
	int n,m;
	cin>>n>>m;
	int ans=0;
	if(m%n!=0)
	{
		cout<<-1<<endl;	
	}
	else
	{
		m=m/n;
		if(m!=1)
		{
			ans++;
		}
		while(m!=2&&m!=1&&m!=3)
		{
			if(m%3==0)
			{
				m=m/3;
				ans++;
			}
			else if(m%2==0)
			{
				m=m/2;
				ans++;
			}
			else
			{
				ans=-1;
				break;
			}
		}
		cout<<ans<<endl;
	}
	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 2 3 4 5
3 2 1 4 5

Sample Output:

5 3

题目分析：

• 对于LIS，要求一个序列中的最长上升子序列，利用动态规划的思想维护一个dp_LIS数组，dp_LIS[i]表示以第i个元素为结尾的最长上升子序列的长度，初始化dp_LIS[i]=1（即本身）。然后对于对于每一个 i ，枚举在i 之前的每一个元素 j ，然后对于每一个 dp_LIS[j],如果a[i] 大于a[j]，那么就可以考虑继承，更新dp_LIS[i]=dp_LIS[j]+1。最后对所有的dp_LIS[i]，求出max最大值即为最长上升子序列长度（注意max(dp_LIS[i])并不一定是dp_LIS[n]哦!!!）
• 对于LCS，要求最长的公共子序列，利用动态规划的思想维护一个二维数组dp_LCS[i][j]，表示的含义是第一个序列的前i位与第二个序列的前j位的公共序列长度，遍历两个序列，如果a[i]=b[j]，那么说明有新的公共元素，需要更新dp_LCS[i][j]，dp_LCS[i][j]=max(dp_LCS[i][j],dp_LCS[i-1][j-1]+1);否则的话，说明没有新的公共元素，令dp_LCS[i][j]=max(dp_LCS[i-1][j],dp_LCS[i][j-1])，最后输出dp_LCS[n][m]即为最长的公共子序列长度。

关于LIS&LCS的详解以及算法的优化见动态规划-各种子序列问题

代码：

#include<iostream>
#include<stdio.h>
#include<string.h>
#include<algorithm>
using namespace std;
const int maxn=5e3+5;
int dp_LIS[maxn];
int dp_LCS[maxn][maxn];
int a[maxn],b[maxn];
int main()
{
	int n,m;
	cin>>n>>m;
	for(int i=1;i<=n;i++)
	{
		cin>>a[i];
	}
	for(int i=1;i<=m;i++)
	{
		cin>>b[i];
	}
	for(int i=1;i<=n;i++)//初始化 
	{
		dp_LIS[i]=1;
		for(int j=1;j<i;j++)//枚举i之前的每一个j 
		{
			if(a[j]<a[i]&&dp_LIS[i]<dp_LIS[j]+1)
			{
				dp_LIS[i]=dp_LIS[j]+1;
			}		
		}
	} 
	int ans1=0;
	for(int i=1;i<=n;i++)//初始化 
	{
		ans1=max(ans1,dp_LIS[i]);
	} 
	memset(dp_LCS,0,sizeof(dp_LCS));
	for(int i=1;i<=n;i++)
	{
		for(int j=1;j<=m;j++)
		{
			dp_LCS[i][j]=max(dp_LCS[i-1][j],dp_LCS[i][j-1]);
			if(a[i]==b[j])
			{
				dp_LCS[i][j]=max(dp_LCS[i][j],dp_LCS[i-1][j-1]+1);
			}
		}
	} 
	cout<<ans1<<" "<<dp_LCS[n][m]<<endl;
	return 0;
}

C-拿数问题||

题目：

YJQ 上完第10周的程序设计思维与实践后，想到一个绝妙的主意，他对拿数问题做了一点小修改，使得这道题变成了 拿数问题 II。

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

本题和课上讲的有些许不一样，但是核心是一样，需要你自己思考。

Input：

第一行包含一个整数 n (1 ≤ n ≤ 105)，表示数字里的元素的个数。
第二行包含n个整数a1, a2, …, an (1 ≤ ai ≤ 105)。

Output：

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

Input:

2 
1 2 

Output:

2 

Input:

3 
1 2 3 

Output:

4 

Input:

9 
1 2 1 3 2 2 2 2 3 

Output:

10 

Hint

对于第三个样例：先选任何一个值为2的元素，最后数组内剩下4个2。然后4次选择2，最终得到10分。

题目分析：

该题在原来的拿数问题上进行了修改，关键在于理解x+1 和 x-1 (如果有 Aj = x+1 或 Aj = x-1 存在) 就会变得不可拿（但是有 Aj = x 的话可以继续拿这个 x）这一句话，这句话说明不再是按照位置进行选或者不选，而是按照数的大小进行决策，所以我们的dp[i]数组针对的不再是数的下标而是数的大小。

• 首先对于输入的n个数，进行每个数的个数的统计并且用cnt数组记录，同时记录下数的大小的左右边界minn,maxx。
• 对于minn-maxx区间内的数进行遍历，根据题意，每个数有两种情况（1）i-1已经被选，那么i将不能再被选；（2）i-1未被选择，那么在i-2的基础上，再选择i（注意选择了i代表所有的数的大小等于i的数都将被选择，所以需要乘上它的出现次数）；所以综合二者可得dp[i]=max(dp[i-1],dp[i-2]+i*cnt[i]);。
• 最后输出dp[i]的最大值即为最大分数值。

(注意数据范围，数的个数e5,数的大小e5,那么最大分值可能会超过int取值范围，所以需要long long int)

代码：

#include<iostream>
#include<stdio.h>
#include<string.h>
using namespace std;
const int maxn=1e5+5;
const int INF=1e8+8;
long long int dp[maxn];
long long int cnt[maxn];
long long int minn=INF;long long int maxx=0; 
int main()
{
	int n;
	cin>>n;
	long long int *a=new long long int[n];
	for(int i=1;i<=n;i++)
	{
		cin>>a[i];
		minn=min(minn,a[i]);
		maxx=max(maxx,a[i]);
	}
	memset(cnt,0,sizeof(cnt));
	for(int i=1;i<=n;i++)//记录每个数出现个数 
	{
		cnt[a[i]]++;
	}
	for(int i=minn;i<=maxx;i++)
	{
		dp[i]=max(dp[i-1],dp[i-2]+i*cnt[i]);
	}
	long long int ans=0;
	for(int i=minn;i<=maxx;i++)
	{
		ans=max(ans,dp[i]);
	}
	cout<<ans<<endl;
	return 0;
}