JZOJ4605。ソート（セグメントツリーのマージおよび分割）

効果の件名：

---

たびに昇順または降順にゾーンは、最終的にはそう何点を尋ねました。

ソリューション：

---

あなただけ確かに古典的な問題である一点を、せれば、それは半分のセグメントツリーのメンテナンス01順番に答えます。

PTYは、リアルタイムでのツリーラインのシーケンス全体のメンテナンスと組み合わせることができるスプリットにそこから取得します。

ルートにマークされ、この番号間隔ツリーラインは、正または逆で入れてしまったの並べ替えを考えてみましょう。

あなたがツリーに取得したい場合はポテンシャルエネルギー、複雑さやの分析によると、スプリットを使用してマージする必要があります\（O（N〜N-ログ）\） 。

コード：

#include<bits/stdc++.h>
#define fo(i, x, y) for(int i = x, B = y; i <= B; i ++)
#define ff(i, x, y) for(int i = x, B = y; i <  B; i ++)
#define fd(i, x, y) for(int i = x, B = y; i >= B; i --)
#define ll long long
#define pp printf
#define hh pp("\n")
using namespace std;

const int N = 1e5 + 5;

int n, m, a[N];
int op, l, r;

struct tree {
    int l, r, x;
} t[N * 120];
#define i0 t[i].l
#define i1 t[i].r
int tot, pl, pr, px;
int rt[N * 120], rev[N * 120];

void upd(int i) {
    t[i].x = t[i0].x + t[i1].x;
}

void add(int &i, int x, int y) {
    i = ++ tot;
    t[i].x = 1;
    if(x == y) return;
    int m = x + y >> 1;
    if(pl <= m) add(i0, x, m); else add(i1, m + 1, y);
}

void merge(int &i, int j) {
    if(!i || !j) { i = i + j; return;}
    merge(t[i].l, t[j].l);
    merge(t[i].r, t[j].r);
    upd(i);
}

void ft(int i, int x, int y) {
    if(x == y) { if(t[i].x == pl) px = x; return;}
    int m = x + y >> 1;
    if(t[i0].x >= pl) ft(i0, x, m); else
        pl -= t[i0].x, ft(i1, m + 1, y);
}

void sp(int i, int &j1, int &j2, int x, int y, int k) {
    if(!i) return;
    if(y <= k) { j1 = i, j2 = 0; return;}
    if(x > k)  { j1 = 0, j2 = i; return;} 
    int m = x + y >> 1;
    if(!j1) j1 = ++ tot;
    if(!j2) j2 = ++ tot;
    sp(i0, t[j1].l, t[j2].l, x, m, k);
    sp(i1, t[j1].r, t[j2].r, m + 1, y, k);
    upd(j1); upd(j2);
}

int bz[N];
void dg(int i, int x, int y) {
    if(!t[i].x) return;
    if(x == y) { bz[x] ++; return;}
    int m = x + y >> 1;
    dg(i0, x, m); dg(i1, m + 1, y);
}

void split(int &i, int &j, int k) {
    pl = rev[i] ? t[i].x - k : k;
    px = 0; ft(i, 0, n);
    if(rev[i]) sp(i, j, i, 0, n, px); else sp(i, i, j, 0, n, px);
}

set<int> s;
set<int> :: iterator st, en;
int d[N], d0;

int main() {
    scanf("%d %d", &n, &m);
    fo(i, 1, n) scanf("%d", &a[i]);
    fo(i, 0, n + 1) s.insert(i);
    fo(i, 0, n) pl = pr = a[i], add(rt[i], 0, n);
    fo(i, 1, m) {
        scanf("%d %d %d", &op, &l, &r);
        en = s.upper_bound(r);
        st = -- s.lower_bound(l);
        d0 = 0;
        for(; st != en; st ++) d[++ d0] = *st;
        d[++ d0] = *en;
        int rt1 = 0, rt2 = 0;
        split(rt[d[1]], rt1, l - d[1]);
        if(t[rt1].x) {
            rt[l] = rt1;
            rev[rt[l]] = rev[rt[d[1]]];
            d[1] = l;
        }
        if(d[d0] != r + 1) {
            split(rt[d[d0 - 1]], rt2, r + 1 - d[d0 - 1]);
            rt[r + 1] = rt2;
            rev[rt[r + 1]] = rev[rt[d[d0 - 1]]];
        }
        fo(j, 2, d0 - 1) {
            if(l != d[j]) merge(rt[l], rt[d[j]]), rt[d[j]] = 0;
            s.erase(d[j]);
        }
        rev[rt[l]] = op;
        s.insert(l); s.insert(r + 1);
    }
    scanf("%d", &l);
    l ++;
    for(st = s.begin(); st != s.end(); st ++) {
        int i = *st;
        if(t[rt[i]].x >= l) {
            pl = rev[rt[i]] ? t[rt[i]].x - l + 1 : l;
            px = 0;
            ft(rt[i], 0, n);
            pp("%d\n", px);
            return 0;
        } else l -= t[rt[i]].x;
    }
}