最长不下降序列（C++）

设有由n个不相同的整数组成的数列，记为:b(1)、b(2)、……、b(n)且b(i)<>b(j)  (i<>j)，若存在i1<i2<i3< … < ie 且有b(i1)<b(i2)< … <b(ie)则称为长度为e的不下降序列。程序要求，当原数列出之后，求出最长的不下降序列。

例如13，7，9，16，38，24，37，18，44，19，21，22，63，15。例中13，16，18，19，21，22，63就是一个长度为7的不下降序列，同时也有7 ，9，16，18，19，21，22，63长度为8的不下降序列。

算法分析：

根据动态规划的原理，由后往前进行搜索(当然从前往后也一样)。

1、对b(n)来说，由于它是最后一个数，所以当从b(n)开始查找时，只存在长度为1的不下降序列；

2、若从b(n-1)开始查找，则存在下面的两种可能性：

①若b(n-1)<b(n)则存在长度为2的不下降序列b(n-1)，b(n)。

②若b(n-1)>b(n)则存在长度为1的不下降序列b(n-1)或b(n)。

3、一般若从b(i)开始，此时最长不下降序列应该按下列方法求出:

在b(i+1),b(i+2),…,b(n)中，找出一个比b(i)大的且最长的不下降序列，作为它的后继。

数据结构】

为算法上的需要，定义一个整数类型二维数组b(N,3)

1、b(I,1)表示第I个数的数值本身；

2、b(I,2)表示从I位置到达N的最长不下降序列长度

3、b(I,3)表示从I位置开始最长不下降序列的下一个位置，若b[I,3]=0则表示后面没有连接项。

求解过程：

①从倒数第二项开始计算，后面仅有1项，比较一次，因63>15，不符合要求，长度仍为1。

②从倒数第三项开始其后有2项，需做两次比较，得到目前最长的不下降序列为2，如下表：

	11	12	13	14	……		11	12	13	14
		22	63	15	……		21	22	63	15
		2	1	1	……		3	2	1	1
		13	0	0	……		12	13	0	0

一般处理过程是：

①在i+1,i+2,…,n项中，找出比b[I,1]大的最长长度L以及位置K;

②若L>0，则b[I,2]:=L+1;b[I,3]:=k;

最后本题经过计算，其数据存储表如下：

1	2	3	4	5	6	7	8	9	10	11	12	13	14
13	7	9	16	38	24	37	18	44	19	21	22	63	15
7	8	7	6	3	4	3	5	2	4	3	2	1	1
4	3	4	8	9	7	9	10	13	11	12	13	0	0

初始化:
for (i=1;i<=n;i++)
{
       cin>>b[i][1];
       b[i][2]=1;b[i][3]=0;
}
下面给出求最长不下降序列的算法:
for (i=n-1;i>=1;i--)
{
       l=0;k=0;
       for (j=i+1;j<=n;j++)
             if((b[j][1]>b[i][1])&&(b[j][2]>l))
             {
                   l=b[j][2];
                   k=j;
             }
       if (l>0) 
       {
             b[i][2]=l+1;
             b[i][3]=k;
       }
}
下面找出最长不下降序列：
k=1;
for (j=1;j<=n;j++)
      if (b[j][2]>b[k][2])  k=j;
最长不下降序列长度为b[k][2]序列
while (k!=0)  
{
      cout<<’’<<b[k][1];
      k=b[k][3];
}
程序代码：

#include<iostream>
using namespace std;
int main()
{
   int n,i,j,l,k,b[200][10]; 
   cout<<"input n:"<<endl;
   cin>>n;
   for (i=1;i<=n;i++)                          //输入序列的初始值
    {
      cin>>b[i][1];
      b[i][2]=1;b[i][3]=0;
    }
     for (i=n-1;i>=1;i--)                      //求最长不下降序列
     {
       l=0;k=0;
        for (j=i+1;j<=n;j++)
         if ((b[j][1]>b[i][1])&&(b[j][2]>l))
          {
             l=b[j][2];
             k=j;
          }
        if (l>0)
         {
           b[i][2]=l+1;b[i][3]=k;
         }  
     }
     k=1;
     for (j=1;j<=n;j++)                        //求最长不下降序列的起始位置
      if (b[j][2]>b[k][2]) k=j;
     cout<<“max=”<<b[k][2]<<endl;  //输出结果
     while (k!=0)                                //输出最长不下降序列
     {
        cout<<' '<<b[k][1];
        k=b[k][3];   
      } 
}

程序运行结果：

输入：14

13 7 9 16 38 24 37 18 44 19 21 22 63 15

输出：max=8

7  9  16  18  19  21  22  63