bzoj 1818 主席树

思路：主席树搞一搞。

#include<bits/stdc++.h>
#define LL long long
#define fi first
#define se second
#define mk make_pair
#define PII pair<int, int>
#define y1 skldjfskldjg
#define y2 skldfjsklejg
using namespace std;
 
const int N = 1e5 + ;
const int inf = 0x3f3f3f3f;
const LL INF = 0x3f3f3f3f3f3f3f3f;
const int mod = ;
 
int hs1[N], hs2[N], mnX[N], mxX[N], mnY[N], mxY[N], tot1, tot2, n;
int root[][N];
 
struct Point {
    int x, y;
    Point(int _x = , int _y = ) {
        x = _x; y = _y;
    }
} a[N];
 
struct Chairman_tree {
    int tot;
    struct node {
        int sum, l, r;
    } a[N * ];
    void update(int p, int l, int r, int &x, int y) {
        x = ++tot; a[x] = a[y]; a[x].sum++;
        if(l == r) return;
        int mid = l + r >> ;
        if(p <= mid) update(p, l, mid, a[x].l, a[y].l);
        else update(p, mid + , r, a[x].r, a[y].r);
    }
    int query(int L, int R, int l, int r, int x, int y) {
        if(L > R) return ;
        if(l >= L && r <= R) return a[x].sum - a[y].sum;
        int mid = l + r >> , ans = ;
        if(L <= mid) ans += query(L, R, l, mid, a[x].l, a[y].l);
        if(R > mid) ans += query(L, R, mid + , r, a[x].r, a[y].r);
        return ans;
    }
} ct[];
 
int main() {
    ct[].tot = ct[].tot = ;
    memset(mnX, inf, sizeof(mnX));
    memset(mnY, inf, sizeof(mnY));
    memset(mxX, , sizeof(mxX));
    memset(mxY, , sizeof(mxY));
    scanf("%d", &n);
    for(int i = ; i <= n; i++) {
        scanf("%d%d", &a[i].x, &a[i].y);
        hs1[++tot1] = a[i].x;
        hs2[++tot2] = a[i].y;
    }
    sort(hs1 + , hs1 +  + tot1);
    sort(hs2 + , hs2 +  + tot2);
    tot1 = unique(hs1 + , hs1 +  + tot1) - hs1 - ;
    tot2 = unique(hs2 + , hs2 +  + tot2) - hs2 - ;
 
    for(int i = ; i <= n; i++) {
        a[i].x = lower_bound(hs1 + , hs1 +  + tot1, a[i].x) - hs1;
        a[i].y = lower_bound(hs2 + , hs2 +  + tot2, a[i].y) - hs2;
        mnX[a[i].x] = min(mnX[a[i].x], a[i].y);
        mxX[a[i].x] = max(mxX[a[i].x], a[i].y);
        mnY[a[i].y] = min(mnY[a[i].y], a[i].x);
        mxY[a[i].y] = max(mxY[a[i].y], a[i].x);
    }
 
    for(int i = ; i <= tot2; i++) {
        ct[].update(mnY[i], , tot1, root[][i], root[][i - ]);
        ct[].update(mxY[i], , tot1, root[][i], root[][i - ]);
    }
 
    LL ans = ;
 
    for(int i = ; i <= tot1; i++) {
        ans += mxX[i] - mnX[i] + ;
        ans -= ct[].query(i + , tot1, , tot1, root[][mxX[i]], root[][mnX[i]-]);
        ans -= ct[].query(, i - , , tot1, root[][mxX[i]], root[][mnX[i]-]);
    }
 
    printf("%lld\n", ans);
    return ;
}