实际上这应当是贪心算法,d[k]表示[0..i]中长度为k的子串中末尾最小的子串的末尾,因为选择末尾最小的子串,x[i+1]连接到他后面更容易一些,那么,为了让d[k]尽可能小,就要让d[k-1]尽可能小。。。。依次类推,于是需要将小一点的数取代前面的数,以保证后面的可以最小。d[k]是不减的,因此可以用二分查找。
具体请看:http://blog.csdn.net/zsc09_leaf/article/details/6536802
即使看了博客,写出代码还是用了很久,在递归终止条件上花了很多时间
#include<stdio.h>
#define max 1000
int biSearch(int d[],int x,int l,int r)/*返回小于等于x的d[m]中,最大的m*/
{
int m=(l+r)/2;
if (x>=d[r])/*合并了r=l的情况,因为d[0]=负无穷,保证一定有解*/
return r;
if ((r-l)==1)
return l;
if (x>=d[m])
return biSearch(d,x,m,r);
else
return biSearch(d,x,l,m-1);
}
int LS(int x[],int d[],int n)
{
int i,k=1,tmp;
d[0]=0x80000000;/*为了简化递归终止条件,将d[0]设为负无穷*/
d[1]=x[0];
for (i=1;i<n;i++)
if (x[i]>=d[k])
d[++k]=x[i];
else
{
tmp=biSearch(d,x[i],0,k-1);
d[tmp+1]=x[i];
}
return k;
}
int main(void)
{
int x[max]={33,123,55,3,78,54,32,67,4,29},n=10,d[max+1],ans;
ans=LS(x,d,n);
printf("%d\n",ans);
int i;
for (i=1;i<=n;i++)
printf("%d ",d[i]);
printf("\n");
return 0;
}