2019多校第二场：Beauty Of Unimodal Sequence(LIS变形)

题目链接：http://acm.hdu.edu.cn/showproblem.php?pid=6592

解题心得：

• 刚开始看到的时候没什么思路，其实是忽略了一个东西，就是$n\ast logn$复杂度的$LIS$中二分查某个数找出来的位置就是以这个数为终点能够得到的最长不降子序列的长度。
• 这样正跑一遍$LIS$，反跑一遍$LIS$，记录下每个位置的答案，这样就能得到一个一个的峰值。第二个需要处理的问题就是字典序最大和最小的问题，仔细想想，如果要字典序最小，当有多个长度相同都是最长长度的满足要求的序列，那么选择的肯定是峰值最靠近左边的，因为如果如果选择右边的，那么字典序肯定会更大，或者用右边的选择替换更新就会更新成最左边那个，所以选择的一定是左边的那个峰值。然后峰值左边就只要满足长度要求并且值不大于峰值一直更新就行了，右边则直接替换。字典序最大则找最右边的那个峰值，然后操作策略和最小左右替换一下就行了。

---

#include <bits/stdc++.h>
using namespace std;
const int maxn = 5e5+100;
const int INF = 0x3f3f3f3f;

int dp[maxn], len[2][maxn], n, num[maxn], pos[maxn];
//len[0]是从左到右跑的LIS，len[1]是从右往左跑的LIS，
//pos[i]记录i这个位置放的是哪个数

vector <int> ve_min, ve_max;

void init() {
    for(int i=1;i<=n;i++) scanf("%d", &num[i]);
    memset(dp, 0x3f, sizeof(dp));
    ve_min.clear();
    ve_max.clear();
    for(int i=1;i<=n;i++) {
        len[0][i] = lower_bound(dp + 1, dp + 1 + n, num[i]) - dp;
        dp[len[0][i]] = num[i];
    }

    memset(dp, 0x3f, sizeof(dp));
    for(int i=n;i>=1;i--) {
        len[1][i] = lower_bound(dp+1, dp+1+n, num[i]) - dp;
        dp[len[1][i]] = num[i];
    }
}

void get_min() {
    memset(pos , 0x3f, sizeof pos);
    int index = 0, max_len = 0;
    for(int i=1;i<=n;i++) {
        int sum = len[0][i] + len[1][i];
        if(sum > max_len) {
            max_len = sum;
            index = i;
        }
    }

    stack <int> st;
    pos[index] = max_len;
    for(int i=index-1;i>=1;i--) {
        if(len[0][i] > len[0][index]) continue;
        if(num[i] >= pos[len[0][i] + 1]) continue;

        while(!st.empty() && len[0][i] >= len[0][st.top()]) {
            pos[len[0][st.top()]] = INF;
            st.pop();
        }
        st.push(i);
        pos[len[0][i]] = num[i];
    }
    while(!st.empty()) {
        ve_min.push_back(st.top());
        st.pop();
    }

    ve_min.push_back(index);
    int last = index;
    for(int i=index+1;i<=n;i++) {
        if(len[1][i] == len[1][last] - 1 && num[last] > num[i]) {
            last = i;
            ve_min.push_back(i);
        }
    }
}

void get_max() {
    memset(pos, 0x3f, sizeof(pos));
    int index = 0, max_len = 0;
    for(int i=1;i<=n;i++) {
        int temp = len[0][i] + len[1][i];
        if(temp >= max_len) {
            max_len = temp;
            index = i;
        }
    }

    int last = index;
    ve_max.push_back(index);
    for(int i=index-1;i>=1;i--) {
        if(len[0][i] == len[0][last] - 1 && num[i] < num[last]) {
            last = i;
            ve_max.push_back(i);
        }
    }
    reverse(ve_max.begin(), ve_max.end());

    stack <int> st;
    for(int i=index+1;i<=n;i++) {
        if(len[1][i] >= len[1][index]) continue;
        if(num[i] >= pos[len[1][i]+1]) continue;

        while(!st.empty() && len[1][i] >= len[1][st.top()]) {
            pos[len[1][st.top()]] = INF;
            st.pop();
        }
        st.push(i);
        pos[len[1][i]] = num[i];
    }
    int now = ve_max.size();
    while(!st.empty()) {
        ve_max.push_back(st.top());
        st.pop();
    }
    reverse(ve_max.begin()+now, ve_max.end());
}

void Print() {
    for(int i=0;i<ve_min.size();i++) {
        printf("%d", ve_min[i]);
        if(i != ve_min.size()-1) printf(" ");
        else printf("\n");
    }
    for(int i=0;i<ve_max.size();i++) {
        printf("%d", ve_max[i]);
        if(i != ve_max.size()-1) printf(" ");
        else printf("\n");
    }
}

int main() {
//    freopen("1.in.txt", "r", stdin);
    while(scanf("%d", &n)!= EOF) {
        init();
        get_min();//得到字典序最小
        get_max();//得到字典序最大
        Print();
    }
    return 0;
}