【模板】动态规划之最长上升子序列LIS问题、概率DP问题（烟花题解）

最新推荐文章于 2020-11-01 18:59:43 发布

coderwait

最新推荐文章于 2020-11-01 18:59:43 发布

阅读量627

点赞数

分类专栏：算法理论和模板

本文链接：https://blog.csdn.net/coderwait/article/details/89679684

版权

算法理论和模板专栏收录该内容

31 篇文章 0 订阅

订阅专栏

最长上升子序列

Description:

老师有n个题目，他希望出一张考试试卷，从中选取一定数量的题目，在不改变给定题目顺序的情况下，要求选取的题目难度严格递增，为了防止有人AK，老师希望在考试中出尽可能多的题目，求最大题目数量。

Input:

每个测试点只有一组测试数据。

第一行一个整数n表示题目数量，第二行n个整数ai表示题目难度。

Output:

输出一个整数ans，一场考试中题目数量的最大值

Sample Input:

1 2 3 1 5

Sample Output:

Hint:

可以选取第1、2、3、5题

//最大上升子序列
#include<iostream>
#include<cstring>
const int MAX = 10e5+7;
using namespace std;
int dp[MAX];//表示最长上升子序列以a[i]结束的情况下，其元素个数 
int a[MAX];
int main(void)
{
	int n;
	cin>>n;
	for(int i=1;i<=n;i++) cin>>a[i];
	memset(dp,0,sizeof(dp));
	for(int i=1;i<=n;i++){
		dp[i]=1;//初始值为1,至少为1 
		for(int j=1;j<i;j++){
			if(a[j]<a[i]) dp[i]=max(dp[i],dp[j]+1);
		}
	} 
	int ans =1;//默认为1 
	for(int i=1;i<=n;i++) ans=max(dp[i],ans); 
	cout<<ans<<endl;
	return 0;
 }

LIS和最大字段和的比较：都是定住一侧范围，然后求最大值。

如LIS，dp[]表示最长上升子序列以a[i]结束的情况下其元素个数，可见时间复杂度为 o(n) 。而最大字段和，是用C[i]表示以必须第i个元素为结束的最大子段和，而C[i]又具有最优子结构性质，可由C[i-1]+A[i]或者A[i]转移而来，时间复杂度为 o(n) 。

概率DP问题

烟花

题目来源

注意点：

1. 边界条件。前0种烟花有i种升天的概率dp[0][i]=0,前i种烟花有0种升天的概率循环求解，其中dp[0][0]=1。

2. 输出形式。

#include<iostream>
#include<cstring>
using namespace std;
double dp[201][201];//前i种烟花有j种升天的概率
double p[201];
int main(void)
{
    memset(dp,0,sizeof(dp));
    memset(p,0,sizeof(p));
    int n,k;//n种烟花 
    cin>>n>>k;
    for(int i=1;i<=n;i++) cin>>p[i];
    dp[0][0]=1; //出错点 
    for(int i=1;i<=n;i++){
        dp[i][0]=dp[i-1][0]*(1-p[i]);
    }
    for(int i=1;i<=n;i++){
        for(int j=1;j<=k;j++){
            dp[i][j]=dp[i-1][j]*(1-p[i])+dp[i-1][j-1]*p[i];
        }
    }
    printf("%.4f",dp[n][k]);
    return 0;
    
 }

优化代码（一维数组）：

#include<iostream>
#include<cstring>
using namespace std;
double dp[201];//前i种烟花有j种升天的概率
double p[201];
int main(void)
{
    memset(dp,0,sizeof(dp));
    memset(p,0,sizeof(p));
    int n,k;//n种烟花 
    cin>>n>>k;
    for(int i=1;i<=n;i++) cin>>p[i];
    dp[0]=1;
    for(int i=1;i<=n;i++){
        for(int j=k;j>=0;j--){
            if(j==0) dp[j]=dp[j]*(1-p[i]);
            else dp[j] = dp[j-1]*p[i] + dp[j]*(1-p[i]);
        }
    }
    printf("%.4f\n",dp[k]);
    return 0;
}