#729287 (C++14) No.1804 Intersection of LIS

• 問題: No.1804 Intersection of LIS
• ユーザー: 🍮かんプリン
• 提出日時: 2022-01-08 18:29:54
• 言語: C++14 (gcc 12.3.0 + boost 1.83.0)
• 結果: AC
• 実行時間: 435 ms / 2,000 ms
• コード長: 4,230 bytes
• コンパイル時間: 2,239 ms
• コンパイル使用メモリ: 183,488 KB
• 実行使用メモリ: 25,408 KB
• 最終ジャッジ日時: 2024-11-14 09:42:46
• 合計ジャッジ時間: 11,518 ms
• ジャッジサーバーID （参考情報）: judge4 / judge5

入力	結果	実行時間; 実行使用メモリ
testcase_00	AC	2 ms; 6,820 KB
testcase_01	AC	2 ms; 6,816 KB
testcase_02	AC	2 ms; 6,816 KB
testcase_03	AC	435 ms; 23,176 KB
testcase_04	AC	433 ms; 23,204 KB
testcase_05	AC	425 ms; 23,180 KB
testcase_06	AC	409 ms; 23,172 KB
testcase_07	AC	416 ms; 23,224 KB
testcase_08	AC	414 ms; 23,172 KB
testcase_09	AC	414 ms; 23,196 KB
testcase_10	AC	412 ms; 23,160 KB
testcase_11	AC	408 ms; 23,148 KB
testcase_12	AC	427 ms; 23,164 KB
testcase_13	AC	413 ms; 23,120 KB
testcase_14	AC	419 ms; 23,256 KB
testcase_15	AC	418 ms; 23,344 KB
testcase_16	AC	417 ms; 23,180 KB
testcase_17	AC	418 ms; 23,256 KB
testcase_18	AC	5 ms; 6,820 KB
testcase_19	AC	4 ms; 6,816 KB
testcase_20	AC	4 ms; 6,816 KB
testcase_21	AC	4 ms; 6,816 KB
testcase_22	AC	4 ms; 6,820 KB
testcase_23	AC	5 ms; 6,816 KB
testcase_24	AC	5 ms; 6,820 KB
testcase_25	AC	5 ms; 6,816 KB
testcase_26	AC	5 ms; 6,820 KB
testcase_27	AC	5 ms; 6,816 KB
testcase_28	AC	2 ms; 6,816 KB
testcase_29	AC	2 ms; 6,820 KB
testcase_30	AC	2 ms; 6,816 KB
testcase_31	AC	2 ms; 6,816 KB
testcase_32	AC	2 ms; 6,816 KB
testcase_33	AC	1 ms; 6,816 KB
testcase_34	AC	2 ms; 6,816 KB
testcase_35	AC	347 ms; 25,404 KB
testcase_36	AC	354 ms; 25,408 KB
testcase_37	AC	316 ms; 23,112 KB
testcase_38	AC	324 ms; 23,104 KB
testcase_39	AC	2 ms; 6,820 KB

ソースコード

/**
 *   @FileName	a.cpp
 *   @Author	kanpurin
 *   @Created	2022.01.08 18:29:52
**/

#include "bits/stdc++.h" 
using namespace std; 
typedef long long ll;


template < class Monoid >
struct SegmentTree {
private:
    using Func = std::function< Monoid(Monoid, Monoid) >;
    Func F;
    Monoid UNITY;
    int n;
    std::vector< Monoid > node;

    int _binary_search(int a, int b, const std::function<bool(Monoid)> &f, Monoid &s, int k = 0, int l = 0, int r = -1) {
        if (r < 0) r = n;
        if (r <= a || b <= l) return n;
        if (a <= l && r <= b && !f(F(s,node[k]))) {
            s = F(s,node[k]);
            return n;
        }
        if (l == r - 1) {s = F(s,node[k]); return l;}
        int ret = _binary_search(a, b, f, s, 2 * k + 1, l, (r - l) / 2 + l);
        return ret != n ? ret : _binary_search(a, b, f, s, 2 * k + 2, (r - l) / 2 + l, r);
    }

public:
    SegmentTree() {}

    SegmentTree(const std::vector< Monoid > &v, const Func f, const Monoid &unity) {
        F = f;
        UNITY = unity;
        int sz = v.size();
        n = 1;
        while (n < sz) n <<= 1;
        node.resize(n * 2 - 1, UNITY);
        for (int i = 0; i < sz; i++) node[i + n - 1] = v[i];
        build();
    }

    SegmentTree(int m, const Monoid &val, const Func f, const Monoid &unity) {
        F = f;
        UNITY = unity;
        n = 1;
        while (n < m) n <<= 1;
        node.resize(n * 2 - 1, UNITY);
        if (val != UNITY) {
            for (int i = 0; i < m; i++) node[i + n - 1] = val;
            build();
        }
    }

    void set(int k, const Monoid &x) {
        node[n + k - 1] = x;
    }

    void build() {
        for (int i = n - 2; i >= 0; i--) node[i] = F(node[2 * i + 1], node[2 * i + 2]);
    }

    void update_query(int x, const Monoid &val) {
        if (x >= n || x < 0) return;
        x += n - 1;
        node[x] = val;
        while (x > 0) {
            x = (x - 1) >> 1;
            node[x] = F(node[2 * x + 1], node[2 * x + 2]);
        }
    }

    Monoid get_query(int l, int r) {
        Monoid L = UNITY, R = UNITY;
        if (l < 0) l = 0;
        if (r < 0) return UNITY;
        if (l >= n) return UNITY;
        if (r >= n) r = n;
        for (l += n, r += n; l < r; l >>= 1, r >>= 1) {
            if (l & 1) L = F(L, node[l++ - 1]);
            if (r & 1) R = F(node[--r - 1], R);
        }
        return F(L, R);
    }

    int binary_search(int l, int r, const std::function<bool(Monoid)> &f) {
        Monoid s = UNITY;
        int ret = _binary_search(l,r,f,s);
        return ret != n ? ret : -1;
    }

    Monoid operator[](int x) const {
        return node[n + x - 1];
    }

    int size() {
        return n;
    }

    void print() {
        for (int i = 0; i < n; i++) {
            std::cout << i << "\t: " << node[n + i - 1] << std::endl;
        }
    }
};

int main() {
    int n;cin >> n;
    vector<int> p(n);
    for (int i = 0; i < n; i++) {
        cin >> p[i];p[i]--;
    }
    vector<int> idx(n);
    iota(idx.begin(), idx.end(),0);
    sort(idx.begin(), idx.end(),[&](int a,int b){return p[a] < p[b];});
    SegmentTree<int> seg1(n,0,[](int a,int b){return max(a,b);},0);
    SegmentTree<int> seg2(n,0,[](int a,int b){return max(a,b);},0);
    vector<int> val(n,-1);
    for (int i = 0; i < n; i++) {
        int v = seg1.get_query(0,idx[i])+1;
        val[i] += v;
        seg1.update_query(idx[i],v);
    }
    for (int i = n-1; i >= 0; i--) {
        int v = seg2.get_query(idx[i]+1,n)+1;
        val[i] += v;
        seg2.update_query(idx[i],v);
    }
    SegmentTree<int> seg3(n,0,[](int a,int b){return max(a,b);},0);
    SegmentTree<int> seg4(n,0,[](int a,int b){return max(a,b);},0);
    vector<bool> isOK(n,true);
    vector<int> ans;
    for (int i = 0; i < n; i++) {
        if (val[i] <= seg3.get_query(idx[i]+1,n)) isOK[i] = false;
        seg3.update_query(idx[i],val[i]);
    }
    for (int i = n-1; i >= 0; i--) {
        if (val[i] <= seg4.get_query(0,idx[i])) isOK[i] = false;
        seg4.update_query(idx[i],val[i]);
    }
    for (int i = 0; i < n; i++) {
        if (isOK[i]) ans.push_back(i+1);
    }
    cout << ans.size() << endl;
    for (int &a : ans) cout << a << " ";
    cout << endl;
    return 0;
}