[第一波模拟\day3\T2]{独立集}(bubble.cpp)

【问题描述】

有一天，一个名叫顺旺基的程序员从石头里诞生了。又有一天，他学会了冒泡排序和独立集。在一个图里，独立集就是一个点集，满足任意两个点之间没有边。于是他就想把这两个东西结合在一起。众所周知，独立集是需要一个图的。那么顺旺基同学创造了一个算法，从冒泡排序中产生一个无向图。

这个算法不标准的伪代码如下：

Pascal版本	C/C++版本
procedure bubblesortgraph(n, a[]) ： /*输入：点数n，1到n的全排列a。 输出：一个点数为n的无向图G。*/ 创建一个有n个点，0条边的无向图G。 repeat swapped = false for i 从 1 到 n-1 ： if a[i] > a[i + 1] ： 在G中连接点a[i]和点a[i + 1] 交换a[i]和a[i + 1] swapped = true until not swapped 输出图G。 //结束。	void bubblesortgraph(n,a[]) //输入：点数n，1到n的全排列a //输出：一个点数为n的无向图G {  创建一个有n个点，0条边的无向图G。 do{  swapped=false for i 从1 到n-1 if(a[i]>a[i+1]) {  在G中连接点a[i]和点a[i+1] 交换a[i]和a[i+1] swapped =true } }while(swapped); 输出图G。 } //结束。

那么我们要算出这个无向图G最大独立集的大小。但是事情不止于此。顺旺基同学有时候心情会不爽，这个时候他就会要求你再回答多一个问题：最大独立集可能不是唯一的，但有些点是一定要选的，问哪些点一定会在最大独立集里。今天恰好他不爽，被他问到的同学就求助于你了。

【输入】

输入包含两行，第一行为N，第二行为1到N的一个全排列。

【输出】

输出包含两行，第一行输出最大独立集的大小，第二行从小到大输出一定在最大独立集的点的编号。

【输入输出样例】

bubble.in	bubble.out
3 3 1 2	2 2 3

【样例说明】

如上图，顶点1和2一定在最大独立集中，其对应的编号为2和3。

【数据范围】

30%的数据满足 N<=16

60%的数据满足 N<=1,000

100%的数据满足 N<=100,000

【题解】

很经典的操作，求最长上升子序列的不动点和唯一点的交

#include <stdio.h>
#include <string.h>
#define RG register
#define oo 0x7fffffff
inline void read(RG int &x)
	{
	RG int c = getchar();
	while(c < '0' || c > '9') c = getchar();
	for(x = 0; c >= '0' && c <= '9'; c = getchar()) x = (x << 1) + (x << 3) + c - '0';
	}
int n, a[100003], d[100003], len, f[100003], g[100003], hash[100003];
int main()
	{
	freopen("bubble.in", "r", stdin), freopen("bubble.out", "w", stdout);
	read(n);
	for(RG int i = 1; i <= n; ++i) read(a[i]);
	for(RG int i = 1; i <= n; ++i)
		{
		if(a[i] > d[len]) d[++len] = a[i], f[i] = len;
		else
			{
			RG int L = 1, R = len, mid;
			while(L <= R)
				{
				mid = L + R >> 1;
				if(d[mid] >= a[i]) R = mid - 1;
				else L = mid + 1;
				}
			d[L] = a[i];
			f[i] = L;
			}
		}
	len = 0, d[0] = oo;
	for(RG int i = n; i; --i)
		{
		if(a[i] < d[len]) d[++len] = a[i], g[i] = len;
		else
			{
			RG int L = 1, R = len, mid;
			while(L <= R)
				{
				mid = L + R >> 1;
				if(d[mid] <= a[i]) R = mid - 1;
				else L = mid + 1;
				}
			d[L] = a[i];
			g[i] = L;
			}
		}
	for(RG int i = 1; i <= n; ++i)
		if(f[i] + g[i] == len + 1) ++hash[f[i]];
	printf("%d\n", len);
	for(RG int i = 1; i <= n; ++i)
		if(f[i] + g[i] == len + 1 && hash[f[i]] == 1)
			printf("%d ", i);
	fclose(stdin), fclose(stdout);
	return 0;
	}