未知：独立集——题解

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

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。
//结束。

那么我们要算出这个无向图G最大独立集的大小。但是事情不止于此。顺旺基同学有时候心情会不爽，这个时候他就会要求你再回答多一个问题：最大独立集可能不是唯一的，但有些点是一定要选的，问哪些点一定会在最大独立集里。今天恰好他不爽，被他问到的同学就求助于你了。
输入
两行。第一行为N，第二行为1到N的一个全排列。
输出
两行。第一行输出最大独立集的大小，第二行从小到大输出一定在最大独立集的点的编号(输入时的序号)。
样例输入
3
3 1 2
样例输出
2
2 3
数据范围

30%的数据满足 N<=16

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

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

这题难想是难想，但是……
标程更难理解好不？！
所以做法直接扔标程里。
所以可能和路由器的代码风格不一致。

#include <cstdio>
#include <iostream>
#include <algorithm>
#include <map>
#include <cstring>
using namespace std;
const int MAXN = 100000;
const int INF = 1e9+7;
int B[MAXN+10], LenB;
int a[MAXN+10], b[MAXN+10];
int A[MAXN+10], LenA;
int L[MAXN+10], R[MAXN+10];
map<pair<int, int>, int> check;
int main(){
    int n;
    scanf("%d", &n);
    for(int i=1;i<=n;i++){
        scanf("%d", &a[i]);
        b[i] = -a[i];//变求下降为求上升 
    }
    //下面求出最长上升子序列 
    for(int i=1;i<=n;i++){
        int pos = lower_bound(A+1, A+1+LenA, a[i]) - A;//找到a[i]放在A的位置 
        L[i] = pos;
        LenA += int(A[pos] == 0);//计算A的长度 
        A[pos] = a[i];
    }
    for(int i=1;i<=n;i++) B[i] = -INF;
    //下面求出最长下降子序列 
    for(int i=n;i>=1;i--){
        int pos = lower_bound(B+1, B+1+LenB, b[i]) - B;
        R[i] = pos;
        LenB += int(B[pos] == -INF);
        B[pos] = b[i];
    }
    printf("%d\n", LenA);//最长上升子序列长
    for(int i=1;i<=n;i++)
        if(L[i] + R[i] - 1 == LenA)
            check[make_pair(L[i], R[i])]++;//创建二元组 
    LenB=0;
    for(int i=1;i<=n;i++){
        //printf("%d\n",check[make_pair(L[i], R[i])]);
        if(L[i] + R[i] - 1 == LenA && check[make_pair(L[i], R[i])] == 1){
            //有一个奇怪的定理
            //当一个数的向左最长下降长度与向右最长上升长度
            //对于所有数的左最长下降长度与向右最长上升长度中是独一无二的
            //那么这个点就符合条件
            //简单的证明：
            //如果这个点并非唯一
            //那么这个点就可以被另一个和它相同的点所替代 
            B[++LenB] = i;
        }
    }
    for(int i=1;i<=LenB;i++)
        printf("%d ", B[i]);
    return 0;
}