动态规划--最长上升子序列

问题描述      

一个数的序列bi，当b1 < b2 < ... < bS 的时候，我们称这个序列是上升的。对于给定的一个序列(a1, a2, ..., aN)，我们可以得到一些上升的子序列(ai1, ai2, ..., aiK)，这里1 <= i1 < i2 < ... <iK <= N。

       比如，对于序列(1, 7, 3, 5, 9, 4, 8)，有它的一些上升子序列，如(1, 7), (3, 4, 8)等等。这些子序列中最长的长度是4，比如子序列(1, 3, 5, 8).

你的任务，就是对于给定的序列，求出最长上升子序列的长度。

输入要求

输入的第一行是序列的长度N(1 <= N <= 1000)。第二行给出序列中的N个整数，这些整数的取值范围都在0到10000。

输出要求

最长上升子序列的长度。

输入样例

7

1 7 3 5 9 4 8

输出样例

4

如何把这个问题分解成子问题呢？经过分析，发现“求以ak（k=1,2, 3…N）为终点的最长上升子序列的长度”是个好的子问题。

 由上所述的子问题只和一个变量相关，就是数字的位置。因此序列中数的位置 k 就是“状态”，而状态 k 对应的“值”，就是以 ak 做为“终点”的最长上升子序列的长度。这个问题的状态一共有 N 个。状态定义出来后，转移方程就不难想了。

假定MaxLen (k)表示以ak做为“终点”的最长上升子序列的长度，那么：

MaxLen (1) = 1

MaxLen (k) = Max { MaxLen (i)：1<i< k 且 ai < ak且 k≠1 } + 1

实际实现:

假设数列5 1 2 4 3在MaxLen[5]的时候，找到b[5]=3上升的数，即b2=1,和b3=2,然后找这些对应的MaxLen最大的一个加1存在MaxLen[5]中，以此类推。

#include <stdio.h>  
#include<string.h>  
const int MAXN=1000;  
  
int b[MAXN],aMaxlen[MAXN+10];  
int main(int argc, char *argv[])  
{  
    int N;  
    scanf("%d",&N);  
    for (int i=1;i<=N;i++)  
       scanf("%d",&b[i]);  
    aMaxlen[1]=1;  
    for (int i=2;i<=N;i++)  
    {  
        int nTmp=0;  
        for (int j=1;j<i;j++)  
          if(b[i]>b[j])  
         {  
           if(nTmp<aMaxlen[j])  
              nTmp=aMaxlen[j];  
          }  
          aMaxlen[i]=nTmp+1;  
    }     
    int nMax=-1;  
    for (int i=1;i<=N;i++)  
      if(nMax<aMaxlen[i])  
         nMax=aMaxlen[i];  
    printf("%d\n",nMax);       
    return 0;  
}