2015 UESTC Training for Dynamic Programming N - 导弹拦截 LIS nlog(n)+打印字典序最小的路径

N - 导弹拦截

Time Limit:1000MS     Memory Limit:65535KB     64bit IO Format:%lld & %llu

Submit  Status  Practice  UESTC 251

Appoint description:  System Crawler  (2016-03-21)

Description

某国为了防御敌国的导弹袭击，开发出一种导弹拦截系统。但是这种导弹拦截系统有一个缺陷：虽然它的第一发炮弹能够到达任意的高度，但是以后每一发炮弹都要高于前一发的高度。某天，雷达捕捉到敌国的导弹来袭，并观测到导弹依次飞来的高度，请计算这套系统最多能拦截多少导弹，同时，司令部想知道拦截下来的导弹的高度。拦截来袭导弹时，必须按来袭导弹袭击的时间顺序，不允许先拦截后面的导弹，再拦截前面的导弹。

Input

第一行是一个整数 t,代表case数。 
对于每一个case，第一行是一个整数 n(1≤n≤100000)； 
第二行是 n个非负整数，表示第 n枚导弹的高度，按来袭导弹的袭击时间顺序给出，以空格分隔。数据保证高度不会超过 100000.

Output

对于每一个case,第一行输出最多能拦截的导弹数，第二行按来袭顺序输出拦截下来的导弹的高度构成的序列，以一个空格隔开。若有不止一种方法可以拦截最多的导弹，输出字典序最小的。

Sample Input

1 
5 
1 6 3 5 7

Sample Output

4 
1 3 5 7

Hint

思路：

b[i]存储了对应长度为i的LIS所对应的最小末尾; dp[i]存储了结尾为a[i]的LIS

求LIS的值挺简单的，难就难在打印路径，而且该路径还是字典序最小的。可以通过证明从后往前所获得的路径字典序最小；

证明：当i<j时，当a[i]<=a[j]时，b[i]<=b[j],这里无论从从后往前或者从前往后都一样; 而当 a[i]>a[j]时,b[i]=b[j],然而你从前往后遍历的时候，记录了a[i]而没有记录a[j]了,很明显比从后往前记录a[j]所得到的字典序要大，因为a[i]>a[j]。

代码：

#include <iostream>
#include <cstdio>
#include <cstring>
#include <algorithm>
#include <vector>
#include <string>
#include <cmath>
#define ll long long
using namespace std;
#define maxn 110000
#define inf 0x3f3f3f3f
int a[maxn],dp[maxn],b[maxn],path[maxn],t,n;
//b[i]存储了对应长度为i的LIS所对应的最小末尾
//dp[i]存储了结尾为a[i]的LIS
//path[i]存储了LIS字典序最小的路径
int main()
{
#ifdef CDZSC_June
	freopen("t.txt","r",stdin);
#endif
	int max1,L,R,mid;
	scanf("%d",&t);
	while(t--)
	{
		memset(b,inf,sizeof(b));//将b全部初始化为正无穷
		scanf("%d",&n);
		max1 =1;
		for(int i = 1; i<=n; i++)
		{
			scanf("%d",&a[i]);
			L = (int)(lower_bound(b+1,b+n+1,a[i]) - b);//L为b数组中第一个>=a[i]的下标(即dp[i]的值)，用来更新b[L]的值
			b[L] = a[i];
			dp[i] = L;
			max1 = max(max1,L);//更新最大值从而获得最长的LIS
		}
		int last = inf,sum = max1;
		for(int i = n; i>=1; i--)//从后前所得到的序列边上字典序最小的LIS路径
		{
			if(!sum)break;
			if(dp[i] == sum && last > a[i])
			{
				last = a[i];
				path[sum--] = a[i];
			}
		}
		printf("%d\n%d",max1,path[1]);
		for(int i = 2; i<=max1; i++)
		{
			printf(" %d",path[i]);
		}
		puts("");
	}
	return 0;
}