dp之lis//SDNU1040+1221+1292

第一行——咕咕咕。
第二行——我现在才发现我竟然欠了这么多算法不会(果然讲完了不做题就会忘得一干二净(自学ing

第三行:22.3.28我自己写的博客我自己都看不懂

子串:连续且顺序与文本串中的字符顺序一致

子序列:可不连续但顺序需一致

最长上升子序列:序列不一定唯一,但长度是唯一的

方法1:

1. dp就是从上一阶段的最优解推到先阶段最优解。定义d[i] (i∈[1,n])来表示前 i 个数以a[i]结尾的最长上升子序列长度,即以第i个数结尾的最长上升子序列是前i-1个数的最长子序列依次跟a[i]比较判断,若a[i]大于a[k],dp[i]=dp[k]+1,若a[i]小于所有a[k],则dp[i]=1,最终记录dp的最大值即为最长上升子序列的长度。

ex:数列为:2 7 1 5 6

        前1个数 d(1)=1 子序列为2;

  前2个数 7前面有2小于7 d(2)=d(1)+1=2 子序列为2 7

  前3个数 在1前面没有比1更小的,1自身组成长度为1的子序列 d(3)=1 子序列为1

  前4个数 5前面有2小于5 d(4)=d(1)+1=2 子序列为2 5

  前5个数 6前面有2 5小于6 d(5)=d(4)+1=3 子序列为2 5 6

2.其实如果只让求最长长度是很简单的,分暴力+模拟栈(二分两种,时间复杂度是O(n^2)和O(nlogn)(应该是的(不过模拟栈没法记录路径
3.然而重点是如何记录路径(好的我不会

 for(int i = 0; i < tot; ++i)
    {
        dp[i] = 1;//截止到第i个数,最长子序列的长度(初始为1,只有a[i]
        for(int j = 0; j < i; ++j)
        {
            if(a[j]>=a[i]&&dp[i]<dp[j]+1)
            {
                dp[i] = dp[j]+1;
            }
        }
        res = max(res,dp[i]);
    }

方法2:

贪心:对于一个上升子序列,其结尾元素越小,越有利于在后面接其他的元素,也就越可能变得更长。因为涉及到对子序列的替换导致最终得到的子序列的长度一致但元素非应该求得的最长上升子序列的元素。
初始时dp[0] = a[0], 从i = 1时遍历原数列, 将每个数与dp数列的末尾的数进行比较, 如果大于末尾的数, 就把a[i]放在dp数列的最后, 如果小于末尾的数,那么就用a[i]替换掉dp中第一个比a[i]大的数。

ex:数列为:2 7 1 5 6

        前1个数 d(1)=1 子序列为2;

  前2个数 7前面有2小于7 d(2)=d(1)+1=2 子序列为2 7

  前3个数 1比2小,用2替换1,最长子序列长度为2,但序列为1 7

  前4个数 同上1 5

  前5个数 6比5小直接加到末尾变为1 5 6

dp[0] = a[0];
len = 0;
for(i = 1; i < n; i++)
{
    if(a[i] > dp[len])
        dp[++len] = a[i];//加在数列的最后
    else 
        *lower_bound(dp, dp+len, a[i])=a[i];//用二分法找dp中第一个大于a[i]的数
}
len++;

记录路径:

在方法1的基础上,第i个数是基于前i-1个数求出来的,所以可以定义一个结构体,把上一位的下标储存在这一位里,类似于构建一个指针,然后从末尾开始保存到数组中再逆向输出即可。

本代码以SDNU1040为例

别人思路:
拦截系统与上升子序列有关,每有一个上升序列的元素,就需要有额外的拦截系统进行拦截,若全为上升序列,则序列有几个即需要几个拦截系统,到这里需要先求出最长上升子序列。
在序列里,最长上升子序列之外,分别有前面的大于上升序列首元素,后面的小于上升序列的末尾元素,中间的逆序的元素,这些元素,一定可以与上升序列的某个元素组成下降序列,被同一系统拦截。一定是上升子序列,因为非下降子序列(即包含相同元素),因为相同的可以算一个元素被一个系统拦截。
综上,只需找出最长子序列的长度-1即为额外所需的导弹拦截系统的个数

我的思路:

将所有的元素存到一个容器里,一次次找最长非上升子序列,将子序列中的元素从容器中删除,剩余元素再次找,直至容器中没有元素。

本代码思路:

#include <iostream>
#include <cstdio>
#include <cmath>

using namespace std;

int main()
{
    int a[25];
    int tot = 0;
    int ans = 1;
    int res = 1;
    int dp[25];
    while(~scanf("%d",&a[tot]))
    {
        tot++;
        if(getchar()=='\n')
            break;
    }
    for(int i = 0; i < tot; ++i)
    {
        dp[i] = 1;//截止到第i个数,最长子序列的长度(初始为1,只有a[i]
        for(int j = 0; j < i; ++j)
        {
            if(a[j]>=a[i]&&dp[i]<dp[j]+1)
            {
                dp[i] = dp[j]+1;
            }
        }
        res = max(res,dp[i]);
    }
    for(int i = 0; i < tot; ++i)
    {
        dp[i] = 1;
        for(int j = 0; j < i; ++j)
        {
            if(a[j]<a[i]&&dp[i]<dp[j]+1)
            {
                dp[i] = dp[j]+1;
            }
        }
        ans = max(ans,dp[i]);
    }
    printf("%d,%d\n",res,ans-1);
    return 0;
}

附上SDNU其他lis

SDNU1292

SDNU1221

  • 0
    点赞
  • 1
    收藏
    觉得还不错? 一键收藏
  • 0
    评论

“相关推荐”对你有帮助么?

  • 非常没帮助
  • 没帮助
  • 一般
  • 有帮助
  • 非常有帮助
提交
评论
添加红包

请填写红包祝福语或标题

红包个数最小为10个

红包金额最低5元

当前余额3.43前往充值 >
需支付:10.00
成就一亿技术人!
领取后你会自动成为博主和红包主的粉丝 规则
hope_wisdom
发出的红包
实付
使用余额支付
点击重新获取
扫码支付
钱包余额 0

抵扣说明:

1.余额是钱包充值的虚拟货币,按照1:1的比例进行支付金额的抵扣。
2.余额无法直接购买下载,可以购买VIP、付费专栏及课程。

余额充值