算法学习笔记 最长不下降子序列LIS

仍是一道经典的动态规划问题，设原数组为a[]，用数组dp[i]记录数组a中以点i为结尾的最长不下降子序列的长度，用pre[]数组来记录每一个节点在以它为结尾的最长不下降子列中的前趋，易知初始条件：dp[i]=1(每一个节点自身就是一个lis)，pre[i]=-1（每一个节点自身的lis没有前趋）。

可知若点i之前存在a[j]<=a[i]，则dp[i]=max(dp[j]+1,dp[i])(j<i)，那么pre[i]={j:max(dp[j]+1)(j<i)}容易求出以i为结尾的序列的最长不下降子列的长度dp[i]，再找到dp[i]的最大值即为所求的最大长度。同时根据最大值点的前趋不断回溯可求得路径，代码如下：

#include<iostream>
#include<stack>
const int maxn=110;
using namespace std;
int main(){
    int a[maxn];
    int dp[maxn];
    int pre[maxn];//用于回溯子串
    int i=0;
    while(cin>>a[i]){
            dp[i]=1;
            pre[i]=-1;//边界条件：每个元素自身都是一个长度为1的非递减序列,不存在前趋元素
            i++;
    }
    int subsize=1,n=i-1,j=0;
    for(i=1;i<=n;++i){//核心语句！！
        for(j=0;j<i;++j){
            if(a[i]>=a[j]&&dp[i]<dp[j]+1)
            {
                dp[i]=dp[j]+1;
                pre[i]=j;
            }
        }
    }
    for(i=1;i<=n;++i){
        if(dp[i]>subsize){
            subsize=dp[i];
            j=i;
        }
    }
    cout<<subsize<<endl;
    stack<int> s;
    while(pre[j]!=-1){//回溯求路径
        s.push(a[j]);
        j=pre[j];
    }
    s.push(a[j]);
    while(!s.empty()){
        cout<<s.top()<<" ";
        s.pop();
    }
   }