关闭

最长不降子序列

标签: 动态规划C++
125人阅读 评论(0) 收藏 举报
分类:

比如arr={1,5,8,2,3,4}的最长不降子序列是1,2,3,4

  • 动态规划的思路:用P[i]去记录以i结尾的最长的不降子序列的长度,比如arr中,P[0]=1;P[1]=2;P[2]=3;P[3]=2;P[4]=3;P[5]=4;

  • 那P是如何生成的:首先P数组初始化每个值都为1,因为每个数自身也算做长度为1的不降子序列;其次在算P[i+1]时要去找从0-i中哪些数比a[i+1]小,从这些元素的P中找到一个最大的P[j],则P[i+1]=P[j]+1就是当前以i+1结尾的最大的不降子序列长度了;最终在P中找到最大的值就是当前数组的最大的不降子序列的长度,同时可以输出这个子序列。

  • 代码如下:

int longest_increasing_sequence1(int a[], int n)
//最长不降子序列 O(N^2),能输出子序列的值
{
    int D[10] = { 1,1,1,1,1,1,1,1,1,1 };
    //记录到第i个数位置,最长子序列的长度
    for (int i = 1; i < n; i++)
    {
        for (int j = 0; j < i; j++)
        {
            if (a[j] <= a[i])//从头找比当前数小的数进行比较
            //<=考虑到了相同的数字也组成了不降序列
            {
                D[i] = ((D[j] + 1)>D[i]) ? (D[j] + 1) : D[i];
                //保证能得到当前最长的子序列
            }
        }
    }
    int max = 1;
    for (int i = 0; i < 10; i++)
        max = ((D[i]>max) ? D[i] : max);
    return max;
}
  • 另一种方法也是从网上得知的,比较巧妙,复杂度为O(nlogn),但是无法输出具体的子序列。思路如下:
    – 用数组D[size]去存长度为size的不降子序列,最后一位数最小的值是D[size];例如a[i]元素:

        if a[i]>=D[size]:size++;D[size]=a[i];
        else:find(D[j-1]<=a[i]<D[j]);D[j]=a[i];
    

    – 也就是说,如果要比较的数比当前最长不降子序列最后一位数的最小值要大,说明这个数可以被加入到这个子序列中去,那么直接将子序列的长度size+1且更新size+1时子序列最后一位的值;如果要比较的数小于子序列最后一位数的最小值,说明这个数只能去找较小的size的子序列,然后插入到这个较小的子序列中去,寻找较小的size=j的方法是二分查找的思想:从0-size中去二分,直到找到了这样一个j,使得D[j-1]<=a[i]< D[j],这样a[i]就可以放在D[j]上,更新D[j]=a[i],但此时size并不会更新;之所以等号放在D[j-1]这里是为了保证当a[i]与D[j-1]相等时,仍构成一个不降的序列。

  • 代码如下:

int b_search(int b[], int m, int k)//二分查找,找长度为m的b中,最小的大于l的数的坐标
//记住二分查找的两个写法:递归和非递归,递归情况形参有四个包括left和right指针,用if(l<=r)判断条件;非递归情况形参只有三个包括size即可,用while(l<=r)判断条件
{
    int l = 0;
    int r = m - 1;
    while (1 <= r)
    {
        int mid = (l + r) / 2;
        if ((b[mid - 1] <= k) && (b[mid] > k))
        //不能写成a<b<c的形式,这样即使出现了重复,也会将重复的数计入到总长度里
            return mid;
        else if (b[mid - 1]>k)
            r = mid - 1;
//即使l=b[mid],也可以让l=mid+1,因为最终还会比较l与l-1两者的范围
            l = mid + 1;
    }
    return 0;//如果比第一个数还要小,则返回-1
}

int longest_increasing_sequence2(int a[],int n)
//最长不降子序列O(nlogn) 无法输出子序列的值
{
    int D[10];
    D[0] = a[0];
    //D[[i]记录最长子序列的长度为i+1,且此序列最后一位最小的数为D[i]
    int size = 1;//记录最长子序列的长度为size=i+1
    for (int i = 1; i < n; i++)
    {
//如果比当前最大长度的最小值大,则长度加一
        {
            D[size] = a[i];
            size++;
        }
        if (a[i] < D[size - 1])
        //如果比当前最大长度的最小值小,去长度较小的里面找能否符合
        {
            int j = b_search(D, size, a[i]);
            D[j] = a[i];
        }
    }
    return size;
}
0
0

查看评论
* 以上用户言论只代表其个人观点,不代表CSDN网站的观点或立场
    个人资料
    • 访问:942次
    • 积分:75
    • 等级:
    • 排名:千里之外
    • 原创:7篇
    • 转载:0篇
    • 译文:0篇
    • 评论:0条
    文章存档