Weekly10-动态规划part1-拿数问题/LIS&LCS

A - 签到题

题目

东东在玩游戏“Game23”。

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

Input

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

Output

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

Simple Input 1
120 51840

Simple Output 1
7

Simple Input 2
42 42

Simple Output 2
0

Simple Input 3
48 72

Simple Output 3
-1

思路

• 首先判断m能否整除n，若不能，则m无法转化成n
• 然后让m不断地去整除2和3，记录次数，直到余数不为0为止
• 若结果为1，则证明可以转化，否则，m无法转化成n

总结

注意数据范围

实现代码

#include<iostream>
using namespace std;

int main()
{
	long long n,m,res=0;
	cin>>n>>m;
	if(m%n!=0)
	{
		cout<<"-1"<<endl;
		return 0;
	}
	else
	{
		long long cur=m/n;
		while(cur%2==0)
		{
			cur=cur/2;
			res++;
		}
		while(cur%3==0)
		{
			cur=cur/3;
			res++;
		}
		if(cur!=1)
		{
			cout<<"-1"<<endl;
			return 0;
		}
		else cout<<res<<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的长度

Simple Input
5 5
1 3 2 5 4
2 4 3 1 5

Simple Output
3 2

思路

LIS的求解

• 定义fa[i]为以a[i]为结尾的LIS的长度
• 初始化fa[i]=1
• fa[i]=max(fa[i],fa[j]+1) (a[j]<a[i])
• 答案为max(fa[i])

LCS的求解

• 定义fb[i][j]为$A_1$,……，$A_i$和$B_1$,……，$B_i$的LCS长度
• 初始化:fb[0][0]=fb[1][0]=fb[0][1]=0
• 当$A_i$==$B_j$时 fb[i][j]=fb[i-1][j-1]+1 否则 fb[i][j]=max(fb[i-1][j],fb[i][j-1])
• 答案为fb[n][m]

总结

比较基础的动态规划例题

实现代码

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

int m,n;
long long a[5005],fa[5005];
long long b[5005],fb[5005][5005];
long long ans1,ans2;

int main()
{
	memset(fa,0,sizeof(fa));
	memset(fb,0,sizeof(fb));
	cin>>n>>m;
	for(int i=0;i<n;i++)
		cin>>a[i],fa[i]=1;
	for(int i=0;i<m;i++)
		cin>>b[i];
	for(int i=0;i<n;i++)
	{
		for(int j=0;j<i;j++)
		{
			if(a[j]<a[i]) fa[i]=max(fa[i],fa[j]+1);
		}
	}
	for(int i=0;i<n;i++)
	{
		ans1=max(ans1,fa[i]);
	}
    for(int i=1;i<=n;i++)
    {
        for(int j=1;j<=m;j++)
        {
            if(a[i-1]==b[j-1])
            {
                fb[i][j]=fb[i-1][j-1]+1;
            }
            else
            {
                fb[i][j]=max(fb[i-1][j],fb[i][j-1]);
            }
        }
    }
	ans2=fb[n][m];
	cout<<ans1<<" "<<ans2<<endl;
	return 0;
}

C - 拿数问题 II

题目

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你能得到最大分值。

Example

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分。

思路

• 和上课讲的方法类似
• 定义dp[i]为仅考虑a[1…i]能拿到的最大分数，b[i]为得分为i的个数
• dp[i]=max（dp[i-2]+i*b[i]，dp[i-1]）
• 答案为dp[max(a[i])]

总结

得分可能比较大，开long long类型

实现代码

#include<iostream>
#include<algorithm>
using namespace std;

int n,a[100005];
long long int b[100005],dp[100005];

int main()
{
	cin>>n;
	for(int i=0;i<n;i++)
	{
		cin>>a[i];
		b[a[i]]++;
	}
	sort(a,a+n);
	int Max=a[n-1];
	dp[1]=b[1];
	dp[2]=max(b[1],2*b[2]);
	for(int i=3;i<=Max;i++)
		dp[i]=max(dp[i-1],dp[i-2]+b[i]*i);
	cout<<dp[Max]<<endl;
	/*for(int i=0;i<=Max;i++)
		cout<<dp[i]<<endl;*/
	return 0;
}