2017年10月5提高组T1 独立集

Description

---

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

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

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最大独立集的大小。但是事情不止于此。顺旺基同学有时候心情会不爽，这个时候他就会要求你再回答多一个问题：最大独立集可能不是唯一的，但有些点是一定要选的，问哪些点一定会在最大独立集里。今天恰好他不爽，被他问到的同学就求助于你了。

Input

---

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

Output

---

两行。第一行输出最大独立集的大小，第二行从小到大输出一定在最大独立集的点的编号(输入时的序号)。

Hint

---

30%的数据满足 N<=16

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

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

Source

---

YPX

Solution

---

可以发现实际上是在找一个最长下降子序列。唯一性我们正反做两次就可以了。

Code

---

#include <stdio.h>
#include <string.h>
#include <math.h>
#include <iostream>
#include <algorithm>
#define rep(i, st, ed) for (int i = st; i <= ed; i += 1)
#define drp(i, st, ed) for (int i = st; i >= ed; i -= 1)
#define erg(i, st) for (int i = ls[st]; i; i = e[i].next)
#define fill(x, t) memset(x, t, sizeof(x))
#define max(x, y) ((x)>(y)?(x):(y))
#define min(x, y) ((x)<(y)?(x):(y))
#define abs(x) (x)<(0)?(-(x)):(x)
#define INF 0x3f3f3f
#define N 100001
#define E 1001
#define L 1001
int t[N], f[N], s[N], d[N], g[N], n;
inline int read() {
    int x = 0; char ch = getchar();
    for(; ch<'0'||ch>'9'; ch=getchar());
    for(; ch<='9'&&ch>='0'; (x*=10)+=ch-'0',ch=getchar());
    return x;
}
int main(void) {
    n = read();
    rep(i, 1, n) {
        t[i] = read();
    }
    fill(s, 31);
    int ans = 0;
    s[0] = -1;
    s[1] = t[1];
    f[1] = 1;
    rep(i, 2, n) {
        int l = 0;
        int r = n;
        while (l <= r) {
            int mid = (l + r) >> 1;
            if (s[mid] <= t[i]) {
                l = mid + 1;
            }else {
                r = mid - 1;
            }
        }
        f[i] = r + 1;
        s[f[i]] = min(s[f[i]], t[i]);
        ans = max(ans, f[i]);
    }
    printf("%d\n", ans);
    fill(s, -31);
    s[0] = INF;
    s[1] = t[n];
    g[n] = 1;
    drp(i, n - 1, 1) {
        int l = 0;
        int r = n;
        while (l <= r) {
            int mid = (l + r) >> 1;
            if (s[mid] <= t[i]) {
                r = mid - 1;
            }else {
                l = mid + 1;
            }
        }
        g[i] = r + 1;
        s[g[i]] = max(s[g[i]], t[i]);
    }
    rep(i, 1, n) {
        if (f[i] + g[i] == ans + 1) {
            d[f[i]] += 1;
        }
    }
    rep(i, 1, n) {
        if (d[f[i]] == 1 && f[i] + g[i] == ans + 1) {
            printf("%d ", i);
        }
    }
    return 0;
}