动态规划题目一：最长单调递增子序列
习题3-1：最长单调递增子序列
[算法设计与实验题解 page 63]
问题描述：
所谓子序列，就是在原序列里删掉若干个元素后剩下的序列，以字符串"abcdefg"为例子，去掉bde得到子序列"acfg"
现在的问题是，给你一个数字序列，你要求出它最长的单调递增子序列。
输入：
多组测试数据，每组测试数据第一行是n(1<=n<=10000)，下一行是n个比1e9小的非负整数
输出：
对于每组测试数据输出一行，每行内容是最长的单调递增子序列的长度
样例输入：
5
1 2 4 8 16
5
1 10 4 9 7
9
0 0 0 1 1 1 5 5 5
样例输出：
5
3
3
我的Java 代码一：
/**
* 算法时间复杂度是O(n^2)
* 用数组b[0:n-1]记录以a[i],0<=i<n为结尾的元素的最长的递增子序列长度
* 序列a的最长的递增子序列的长度是max{b[i]}(0<=i<n)，b[i]满足最优子结构，可以递归的定义：
* b[0]=1;b[i]=max{b[k]}+1 (0<=k<i,a[k]<a[i])
*/
import java.util.Scanner;
public class LIS {
public static int l;
public static int maxl;
public static int[] a;
public static int[] b;
public static void main(String[] args) {
Scanner sin=new Scanner(System.in);
while(sin.hasNext()){
l=sin.nextInt();
a=new int[l];
b=new int[l];//记住数组一定要初始化给定一个大小，不然后面赋值时会报空指针错误
for(int i=0;i<l;i++){
a[i]=sin.nextInt();
}
maxl=dp();
System.out.println(maxl);
}
}
private static int dp() {
int i,j,k,t;
t=0;
b[0]=1;
for(i=0;i<l;i++){
k=0;
for(j=0;j<i;j++){
if(a[i]>a[j] && k<b[j]){
k=b[j];
}
}
b[i]=k+1;
}
for(i=0;i<l;i++){
if(t<b[i]){
t=b[i];
}
}
return t;
}
}
我的Java 代码二：
/**
* 算法时间复杂度是O(nlog(n))
* 通过归纳假设得到问题变为：已知计算序列a的最长的递增子序列的长度k以及序列
* a中所有长度为k的递增子序列中的最小结尾元素b[k]
*/
import java.util.Scanner;
public class LIS2 {
public static int l;
public static int maxl;
public static int[] a;
public static int[] b;
public static void main(String[] args) {
Scanner sin=new Scanner(System.in);
while(sin.hasNext()){
l=sin.nextInt();
a=new int[l];
b=new int[l+1];
for(int i=0;i<l;i++){
a[i]=sin.nextInt();
}
dp();
System.out.println(maxl);
}
}
private static void dp() {
int i,low,up,mid;
b[1]=a[0];
maxl=1;
for(i=1;i<l;i++){
if(a[i]>b[maxl]){//此时最大长度要加1，并且要改变最小元素值
maxl++;
b[maxl]=a[i];
}else if(a[i]<b[1]){//如果比第一个元素还要小，那么就换一下最小的元素值
b[1]=a[i];
}else{//这种情况下要找到合适的位置将a[i]插入到b中的某个位置
low=1;
up=maxl;
// while(low!=up-1){//原书的写法
// mid=(low+up)/2;
// //System.out.println(low+" "+up+" "+mid);
// //System.out.println(b[mid]+" "+a[i]);
// if(b[mid]<=a[i]){//找到满足b[j-1]<=a[i]<b[j]的j(=up)
// low=mid;
// }else{
// up=mid;
// }
// //System.out.println(low+" "+up+" "+mid);
// }
// b[up]=a[i];
while(low<up){//另一种写法，更加符合二分查找
mid=(low+up)/2;//注意条件：low等于up时就要退出了
//System.out.println(low+" "+up+" "+mid);
//System.out.println(b[mid]+" "+a[i]);
if(b[mid]<=a[i]){//找到满足b[j-1]<=a[i]<b[j]的j(=up)
low=mid+1;
}else{
up=mid-1;
}
//System.out.println(low+" "+up+" "+mid);
}
b[up]=a[i];
}
}
}
}
其他的资源：
可供参考的资料：
1.这道题目的两种解题算法总结博客：http://old.blog.edu.cn/user2/43845/archives/2006/1115309.shtml
2.国家集训队论文 李睿的论文：〈二分法和统计研究〉：http://wenku.baidu.com/view/313fbcd126fff705cc170a18.html
3.论坛中的其他解题方法：http://yzfy.org/dis/listpost.php?tid=101
