LIS有N2的算法,还有NLOGN的算法,这里说一下NLOGN的算法。
我们设S[]为原始序列,F[i]为以S[i]结束的LIS的最大长度。在求LIS时,对于S[i],我们找一个最大的可能的F[j],(j<i)使S[i]尽可能大,在寻找一个F[i]时,我们可以用二分查找优化。
我们要通过枚举F[1...i-1]找到可作为F[i]前缀的LIS的最大长度,换一种思路,我们就可以枚举可能的长度——用二分——再检查它是否可以连接在S[i]的前面。
于是我们用C[i]记录长度为i的LIS的末位的最小值(这里有点贪心的想法),然后在处理S[i]时,通过二分,查找末尾数小于S[i]的LIS的最大长度,然后更新C[i]数组,最终通过查询C[i]数组得到结果。
编码的难点是二分查找,需要仔细体会。
对于二分:
A. 如果一个有序表中有元素X, 则通过二分查找必能找到X
B. 在二分中,我用v记录可行的一个LIS长度(即存在长度为v的可以连接在S[i]前面的LIS),然后查找是否有更大的v,于是要让mid=low+1,此时若存在更好的v一定可以找到。
/*
* TYVJ-1049 最长不下降子序列
* mike-w
* 2012-10-16
*/
#include<stdio.h>
#include<stdlib.h>
#include<string.h>
#define MAXSIZE 5555
#define INF 0x3ffffff
#define xLOL_DEBUG
int s[MAXSIZE];
int c[MAXSIZE];
int n;
int min(int e1, int e2)
{
return e1<e2?e1:e2;
}
int main(void)
{
int i, j, high, low, mid, v;
scanf("%d", &n);
for(i=1; i<=n; i++)
scanf("%d", s+i);
for(i=1; i<=n; i++)
c[i]=INF;
c[1]=s[1];
for(i=2; i<=n; i++)
{
high=i-1, low=1, v=-1;
while(low<=high)
{
mid=(low+high)/2;
if(c[mid]<=s[i])
low=mid+1, v=mid;
else
high=mid-1;
}
if(v==-1) /* s[i] is too small */
c[1]=s[i];
else /* update c[] */
c[v+1]=min(s[i], c[v+1]);
#ifdef LOL_DEBUG
printf("i=%d v=%d\n", i, v);
puts("c:");
for(j=1; j<=n; j++)
printf("%d ", c[j]);
putchar('\n');
#endif
}
for(i=n; i>=1; i--)
if(c[i]<INF)
break;
printf("%d\n", i);
return 0;
}
希望上述文字能解决大家的一点问题——貌似我有点言语混乱。