1.最长单调递增子序列(从头部到尾部是严格递增的)
单调递增的辅助数组sort[],k表示数组sort[]目前的长度,算法完成后k的值即为LIS的长度。
初始化:sort[0] = A[0],k = 1
从前到后扫描数组A[],比较当前的数A[i]和sort[k-1]:
如果A[i]>sort[k-1],A[i]大于sort[]的最后一个元素,A[i]直接插至sort[]的末尾,sort[]的长度增加1,sort[k++]=A[i];
如果A[i]<sort[k-1],在sort[1]...sort[k]中二分查找第一个>A[i]的数sort[index],修改sort[index]=A[i]。(修改后序列可能不对,但是个数是正确的),将sort[index]变小,减小递增序列的值,这样使得递增序列sort可以容纳更多的元素。
二分法T(n) = O(nlgn)
#include <stdio.h>
#include <stdlib.h>
void print(int const *sort, int cnt)
{
int i;
for (i=0; i<cnt-1; ++i)
{
printf("%d ", sort[i]);
}
printf("%d\n", sort[i]);
}
/*在升序序列中二分查找第一个>A[i]的坐标*/
int biSearch(int const *sort, int low, int high, int key)
{
int mid;
while(low <= high)
{
mid = low + ((high - low) >> 1);
if (sort[mid] < key)
{
low = mid + 1;
}
else
{
high = mid - 1;
}
}
return low;
}
void LIS(int const *A, int n)
{
int i, index;
int cnt = 1;
int *sort = (int *)malloc(sizeof(int)*n);
int *dp = (int *)malloc(sizeof(int)*n); /*统计每个位置的个数*/
dp[0] = 0;
sort[0] = A[0];
for (i=1; i<n; ++i)
{
if (A[i] > sort[cnt-1]) /*严格递增,直接插至递增序列的尾部*/
{
dp[i] = cnt;
sort[cnt++] = A[i];
}
else
{
/*二分查找第一个>A[i]的index,修改sort[index]&#