RQNOJ 合唱队形（LIS最大上升序列解析）

题目链接 http://www.rqnoj.cn/Problem_26.html
题目要求最小的出列人数，换个想法就是求满足题目要求队形的人数最多是多少。题目中要求他们的身高满足T1<...<Ti>Ti+1>…>TK(1<=i<=K)，那么我们对于Ti可以求左边的最大上升序列，右边就是求最大下降序列，算出他们的总和，然后用总人数减去这个总和就能得出最后结果。当然我们求最大下降序列时可以完全转化为从数组尾部开始的最大上升序列。
关于LIS介绍两种算法，一数组a[N](a[1]开始), D[i]代表以a[i]结尾的最长上升子序列的个数
D[i]=MAX(D[j])+1 (1=<j<i且a[i]>a[j]);显然这种算法会有两重循环，时间复杂度为O(n^2)。如果采用了这种思路因为外层还有一个数组遍历的循环，时间复杂度会达到O(n^3),估计是题目测试数据太少也能通过，但不是最好的方法。
采用二分的策略
开辟一个栈b，每次取栈顶元素s和读到的元素a做比较，如果a>s，则置为栈顶；如果a<s，则二分查找栈中的大于等于a的第1个数，并替换。最终栈的大小即为最长递增子序列为长度。时间复杂度为O(nlogn)。题目的时间复杂度就为O(n^2logn)。
下面举个例子
原序列为2，4，6，3，4，5

原栈为2，4，6
1 读到3，则用3替换4，得到栈中元素为2，3，6
2 读到4，用4替换6，变为2，3，4
3 读到5，栈顶加一个元素变为2，3，4，5
最终得到最大上升序列的长度为 4，最大元素为 5。为什么一定要替换呢，我们看例子假如没有发生替换，栈中元素一直为2，4，6。那么读入最后一个元素5时，最大序列长度便不会发生变化，发生错误。替换意义就在于看后边的元素能不能组成更长的上升序列。
假如序列为2,4,6,3
只发生了第一步替换栈中元素为 2，3，6 。2，3，6在这个序列中是最大上升序列吗？并不是，应该是2，4，6。所以由于发生了替换不能求出最大上升序列的每个元素。只能求出最后一个元素和序列的长度。
 O(n^3)

#include <stdio.h>
#define MAX  102

int N, a[MAX], a1[MAX], b[MAX], Max_Lis, MAX_sum;

int Lis(int c[], int n)
{
	for(int i = 1; i <= n; i++)
	{
		Max_Lis = 0;
		for(int j = 1; j <= i; j++)
		{
			if(c[i] > c[j] && Max_Lis < b[j]) Max_Lis = b[j];
		}
		b[i] = Max_Lis + 1; 
	}
	return b[n];
}

int main()
{
	scanf("%d", &N);
	for(int i = 1; i <= N; i++) 
	{
		scanf("%d", &a[i]);
		a1[N - i + 1] = a[i];
		b[i] = 1;
	}
	MAX_sum = 0;
	for(int i = 1; i <= N; i++)
	{
		if(MAX_sum < Lis(a, i) + Lis(a1, N - i + 1) - 1)
		{
			MAX_sum = Lis(a, i) + Lis(a1, N - i + 1) - 1;
		}
	}
	printf("%d\n",N - MAX_sum);
	return 0;
}

O(n^2logn)

#include <stdio.h>
#include <string.h>
#define MAX  102

int N, a[MAX], a1[MAX], b[MAX], s[MAX];
int Max_Lis, MAX_sum;

int Lis(int c[], int n)
{
	int ans = 0;
	memset(s, 0, sizeof(s));
	for(int i = 1; i <= n; i++)
	{
		int low = 1, high = ans;
		while(low <= high)
		{
			int mid = (low + high) >> 1;
			if(s[mid] < c[i]) low = mid + 1;
			else high = mid - 1;
		}
		s[low] = c[i];
		if(low > ans) ans++;
	}
	return ans;	
}

int main()
{
	scanf("%d", &N);
	for(int i = 1; i <= N; i++) 
	{
		scanf("%d", &a[i]);
		a1[N - i + 1] = a[i];
		b[i] = 1;
	}
	MAX_sum = 0;
	for(int i = 1; i <= N; i++)
	{
		if(MAX_sum < Lis(a, i) + Lis(a1, N - i + 1) - 1)
		{
			MAX_sum = Lis(a, i) + Lis(a1, N - i + 1) - 1;
		}
	}
	printf("%d\n",N - MAX_sum);
	return 0;
}