NKOI 3500 独立集

输入格式

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

输出格式

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

样例输入

3 
3 1 2

样例输出

2
2 3

提示

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

仔细观察后发现，他说的冒泡排序并没有任何卵用(给你一个暴力的思路吧，大概),只要仔细观察就会发现答案其实就是最长上升子序列的长度

我们用f[i]表示以i结尾的最长上升子序列的长度，g[i]表示以i为开头的最长上升子序列的长度，那么对于一个点i,如果i在一个最长上升子序列上，则说明f[i]+g[i]-1==ans(ans为最长上升子序列长度)

并且如果这个点是独立集中所必须存在的点，说明cnt[f[i]]==1,这里的cnt[f[i]]记录的是值为f[i]的i的个数

显然这里要求两次LIs，必须用到nlogn的算法，这里用的是树状数组优化

#include<cstdio>
#include<iostream>
#include<cstring>
#define lowbit(xx) (xx&-xx)
using namespace std;
const int maxn=1e5+5,inf=0x3f3f3f3f;
inline void _read(int &x){
    char t=getchar();bool sign=true;
    while(t<'0'||t>'9')
    {if(t=='-')sign=false;t=getchar();}
    for(x=0;t>='0'&&t<='9';t=getchar())x=x*10+t-'0';
    if(!sign)x=-x;
}
int n,s[maxn],f[maxn],g[maxn],cnt[maxn],c[maxn],ans;
void insert(int x,int d){for(int i=x;i<=maxn;i+=lowbit(i))c[i]=max(d,c[i]);}
int getmax(int x){
	int maxx=0;
	for(int i=x;i;i-=lowbit(i))maxx=max(maxx,c[i]);
	return maxx; 
}
int main(){
	_read(n);
	int i,j;
	for(i=1;i<=n;i++)_read(s[i]);
	for(i=1;i<=n;i++){
		f[i]=getmax(s[i])+1;
		ans=max(ans,f[i]);
		insert(s[i],f[i]);
	}
	cout<<ans<<endl;
	memset(c,0,sizeof(c));
	for(i=n;i;i--){
		g[i]=getmax(maxn-s[i])+1;
		if(f[i]+g[i]-1==ans)cnt[f[i]]++;
		insert(maxn-s[i],g[i]);
	}
	for(i=1;i<=n;i++)
		if(f[i]+g[i]-1==ans&&cnt[f[i]]==1)printf("%d ",i);
}