【魔改】树状数组 牛客多校第五场I vcd 几何+阅读理解

https://www.nowcoder.com/acm/contest/143/I

vc-dimension

题解：分三种情况，组合数学算一下，其中一种要用树状数组维护

技巧（来自UESTC）：1.循环技巧i主j滑

2.树状数组：一个数列从左到右分别维护某个元素左边比它大num的与右边比他大的num时，从上往下扫， 对每个点的x坐标离散化累加1到X轴上，然后就会发现sum(p[i].x-1)就是左边比它大的，i-1-sum(p[i])就是右边比它大的。 注意y相同的点，需要一起更新（我已开始一个一个更新，根本写不出来）

orz太屌了

坑点：我循环写错了+公式里i,j写反了，wa了20发

比赛时x,y轴都看反了orz

#define _CRT_SECURE_NO_WARNINGS
#include<cmath>
#include<iostream>
#include<stdio.h>
#include<algorithm>
#include<cstring>
#include<stack>
#include<vector>
#include<string.h>
#include<queue>
#include<string>
#include<set>
using namespace std;
#define rep(i,t,n)  for(int i =(t);i<=(n);++i)
#define per(i,n,t)  for(int i =(n);i>=(t);--i)
#define mmm(a,b) memset(a,b,sizeof(a))
#define eps 1e-6
#define pb push_back
#define lowbit(x) ((x)&(-(x)))
const int maxn = 3e5 + ;
const int inf = 1e7 + ;//0x7fffffff;   //无限大
const int MOD = ;
typedef long long ll;
const int mod = ;
int power(int a, int b) {
    int c = ; a %= mod;
    while (b) {
        if (b & ) c = 1ll * c*a%mod;
        a = 1ll * a*a%mod; b >>= ;
    }
    return c;
}
 
struct node {
    int x, y;
    void sc() { scanf("%d%d", &x, &y); }
}p[maxn];
 
int n;
bool cmp(node a, node b) {
    return a.x > b.x;
}
 
bool cmpy(node a, node b) {
    return a.y < b.y;
}
 
int f[maxn];
int idx;
void add(int x, int y) { for (; x <= idx; x += lowbit(x)) f[x] += y; }
int sum(int x) { int ans = ; for (; x; x -= lowbit(x)) ans += f[x]; return ans; }
 
int main() {
    int div2 = power(, mod - );
    cin >> n;
    rep(i, , n) {
        p[i].sc();
    }
    int ans = n;
     ans +=  1ll*n * (n - )%mod*div2 % mod;
    if (ans >= mod)ans -= mod; else if (ans < )ans += mod;
    sort(p + , p +  + n, cmpy);
     idx = ;
    for (int i = , j; i <= n; i = j + ) {
        for (j = i ; j < n&&p[i].y == p[j+].y; j++);
        ++idx;
        rep(k, i, j)p[k].y = idx;
        ans -=1ll* (j - i)*(j - i + )%mod*div2%mod;
        if (ans >= mod)ans -= mod; else if (ans < )ans += mod;
    }
    mmm(f, );
    sort(p + , p +  + n, cmp);
    for (int i = , j; i <= n; i = j + ) {
        for (j = i ; j < n&&p[i].x == p[j+].x; j++);
        rep(k, i, j) {
            int t1 = sum(p[k].y - ); int t2 = i -  - sum(p[k].y);
            ans += 1ll * t1*t2%mod;
            if (ans >= mod)ans -= mod; else if (ans < )ans += mod;
        }
        rep(k, i, j)add(p[k].y, );
    }
    cout << ans << endl;
    cin >> n;
    return ;
}
 
/*
 
*/

一开始无脑分别前后维护

#define _CRT_SECURE_NO_WARNINGS
#include<cmath>
#include<iostream>
#include<stdio.h>
#include<algorithm>
#include<cstring>
#include<stack>
#include<vector>
#include<string.h>
#include<queue>
#include<string>
#include<set>
using namespace std;
#define rep(i,t,n)  for(int i =(t);i<=(n);++i)
#define per(i,n,t)  for(int i =(n);i>=(t);--i)
#define mmm(a,b) memset(a,b,sizeof(a))
#define eps 1e-6
#define pb push_back
const int maxn = 1e5 + ;
const int inf = 1e7 + ;//0x7fffffff;   //无限大
const int MOD = ;
typedef long long ll;
const long long mod = ;
int power(int a,int b){
    int c=; a%=mod;
    while (b) {
        if (b&) c=1ll*c*a%mod;
        a=1ll*a*a%mod; b>>=;
    }
    return c;
}
 
ll div2 = power(, mod - );
struct node {
    int x, y;
    int id;
    void sc() { scanf("%d%d", &x, &y); }
}P[maxn];
 
ll n;
bool cmp(node a, node b) {
    if (a.x != b.x)
    return a.x < b.x;
    else return a.y < b.y;
}
bool cmpx(node a, node b) {
    if (a.x != b.x)
        return a.x > b.x;
    else return a.y < b.y;
}
 
bool cmpy(node a, node b) {
    return a.y < b.y;
}
int cnt[maxn];
ll l[maxn], r[maxn];
int d[maxn];
int lowbit(int x) { return x & (-x); }
int fr[maxn], bk[maxn];
void add(int x, int v) {//a[x]+=v;
    while (x <= maxn) {
        d[x] += v;
        x += lowbit(x);
    }
 
}
int query(int x) {
    int res = ;
    while (x) {
        res += d[x];
        x -= lowbit(x);
    }
    return res;
}
 
int main() {
    cin >> n;
    rep(i, , n) {
        P[i].sc();
        P[i].id = i;
    }
    sort(P + , P + n + , cmp);
 
    int idx = ;
    int now = ;
    rep(i, , n) {
        if (P[i].x != now) {
            cnt[++idx]++;
            now = P[i].x;
        }
        else {
            cnt[idx]++;
        }
    }
    ll ans = n;
    ll tmp = n * (n - ) %mod*div2;
    rep(i, , idx)if (cnt[i] > ) { tmp = (tmp - cnt[i] * (cnt[i] - ) % mod *div2) % mod; }
    ans = (ans + tmp) % mod;
    //维护每个点左右比它高的点
    sort(P + , P + n + , cmpy);
    idx = ;
    now = ;
    rep(i, , n) {
        if (P[i].y != now) {
            P[i].y = ++idx;
            now = P[i].y;
        }
        else {
            P[i].y = idx;
        }
    }
    sort(P + , P + n + , cmp);
     now = ;
    rep(i, , n)
    {
        if (P[i].x != P[i - ].x)fr[i] =now= i - ;
        else fr[i] = now;
    }
    now = ;
    per(i, n, )
    {
        if (P[i].x != P[i + ].x)bk[i] = now = i + ;
        else bk[i] = now;
    }
    rep(i, , n) {
        if(P[i].x!=P[i-].x)l[i] =ll( i--query(P[i].y));
        else {
            if (P[i].y >= P[fr[i]].y || P[i].y >= P[bk[i]].y)l[i] = ;
            else l[i] = l[i - ];
        }
        add(P[i].y, );
    }
    mmm(d, );
    //sort(P + 1, P + n + 1, cmpx);
    int first = ;
    per(i, n, ) {
        if (P[i].x != P[i + ].x)r[i] = ll(n-i - query(P[i].y)),first=;
        else {
            if (P[i].y >= P[fr[i]].y || P[i].y >= P[bk[i]].y)r[i] = ;
            else {
                if (first)
                {
                    first = ;
                    r[i] = ll(n - i -- query(P[i].y));
                }
                else r[i] = r[i + ];
            }
        }
        add(P[i].y, );
    }
    tmp = ;
    rep(i, , n) {
        tmp = (tmp + l[i] * r[i] % mod) % mod;
    }
    ans = (ans + tmp) % mod;
        cout << ans<<endl;
 
        //sort(P + 1, P + n + 1, cmp);
        rep(i, , n) {
            //cout << l[i] << ' ' << r[i] << endl;
        }
    cin >> n;
}
 
/*
 
*/