[BZOJ2658][ZJOI2012]小蓝的好友(mrx)-Treap

小蓝的好友

Description

终于到达了这次选拔赛的最后一题，想必你已经厌倦了小蓝和小白的故事，为了回馈各位比赛选手，此题的主角是贯穿这次比赛的关键人物——小蓝的好友。
在帮小蓝确定了旅游路线后，小蓝的好友也不会浪费这个难得的暑假。与小蓝不同，小蓝的好友并不想将时间花在旅游上，而是盯上了最近发行的即时战略游戏——SangoCraft。但在前往通关之路的道路上，一个小游戏挡住了小蓝的好友的步伐。
“国家的战争其本质是抢夺资源的战争”是整款游戏的核心理念，这个小游戏也不例外。简单来说，用户需要在给定的长方形土地上选出一块子矩形，而系统随机生成了N个资源点，位于用户所选的长方形土地上的资源点越多，给予用户的奖励也越多。悲剧的是，小蓝的好友虽然拥有着极其优秀的能力，但同时也有着极差的RP，小蓝的好友所选的区域总是没有一个资源点。
终于有一天，小蓝的好友决定投诉这款游戏的制造厂商，为了搜集证据，小蓝的好友想算出至少包含一个资源点的区域的数量。作为小蓝的好友，这自然是你分内之事。

Input

每个输入文件中仅包含一个测试数据。
第一行包含两个由空格隔开的正整数R,C,N，表示游戏在一块[1,R]X[1,C]的地图上生成了N个资源点。
接下来有N行，每行包含两个整数 x,y，表示这个资源点的坐标(1<=x<=R,1<=Y<=c)。

Output

输出文件应仅包含一个整数，表示至少包含一个资源点的区域的数量。具体的说，设N个资源点的坐标为(i=1..n),你需要计算有多少个四元组(LB,DB,RB,UB)满足1<=LB<=RB<=R,1<=DB<=UB<=C，且存在一个i使得LB<=Xi<=RB,DB<=Yi<=UB均成立。

Sample Input

5 5 4
1 2
2 3
3 5
4 1

Sample Output

139

HINT

【数据范围】

对于100%的数据，R,C<=40000,N<=100000，资源点的位置两两不同，且位置为随机生成。

---

一道被搬进校内互测的神题，现场无人AC……
当时搬题人给了个70pts的$O(n^2)$部分分，其他dalao都用了单调栈，就咱这种数据结构学傻选手强行用主席树加了一个$log$……

---

思路:
考虑补集转化，也就是计算不合法的矩形数量。

考虑逐行扫描，每次计算以这一行为底边的空矩形个数。
可以发现这跟每一个位置向上第一个遇到的资源点位置有关。

那么可以发现每一时刻咱们要维护的都是这样一个东西:

于是把这个奇怪的东西找个数据结构维护一下就好了~
咱校内互测的时候就用了主席树去维护~

然而主席树会T。
考虑换成平衡树，以高度为关键字建出出一棵笛卡尔树。
此时记每个节点高度为$h[i]$，子树大小为$siz[i]$，那么单个节点的贡献为:

$${siz\times (siz+1) \over 2}\times(h[i]-h[fa])$$

这样只需要支持整棵树$h[i]+1$，和将某个节点的$h[i]$改为$0$两个操作辣~

这样做貌似复杂度依旧不对。
但是，由于数据随机，这其实也是棵Treap……
于是复杂度得到了保证。$O(nlogn)$。

代码比咱的主席树还好写qwq

#include<bits/stdc++.h>
using namespace std;

typedef long long ll;
const int N=40009;
const int Q=100009;

inline int read()
{
    int x=0;char ch=getchar();
    while(ch<'0' || '9'<ch)ch=getchar();
    while('0'<=ch && ch<='9')x=x*10+(ch^48),ch=getchar();
    return x;
}

int rt;
namespace treap
{
    const int M=N;

    int ch[M][2],siz[M],fa[M];
    ll val[M],sum[M],tag[M];

    inline void calc(int x)
    {
        sum[x]=sum[ch[x][0]]+sum[ch[x][1]]+(ll)(siz[x]+1ll)*siz[x]/2*(val[x]-val[fa[x]]);
    }

    inline void mark(int x,int v){tag[x]+=v;val[x]+=v;calc(x);}

    inline void update(int x)
    {
        if(ch[x][0])calc(ch[x][0]);
        if(ch[x][1])calc(ch[x][1]);
        siz[x]=siz[ch[x][0]]+siz[ch[x][1]]+1;
        calc(x);
    }

    inline void push(int x)
    {
        if(tag[x])
        {
            if(ch[x][0])mark(ch[x][0],tag[x]);
            if(ch[x][1])mark(ch[x][1],tag[x]);
            tag[x]=0;
        }
    }

    inline void rotate(int x)
    {
        int y=fa[x],z=fa[y];
        int l=(ch[y][1]==x);
        if(z)ch[z][ch[z][1]==y]=x;
        fa[x]=z;fa[y]=x;fa[ch[x][l^1]]=y;
        ch[y][l]=ch[x][l^1];ch[x][l^1]=y;
        update(y);update(x);
        if(y==rt)rt=x;
    }

    inline int build(int l,int r)
    {
        if(l>r)return 0;
        int mid=l+r>>1,lc,rc;
        if(ch[mid][0]=build(l,mid-1))
            fa[ch[mid][0]]=mid;
        if(ch[mid][1]=build(mid+1,r))
            fa[ch[mid][1]]=mid;
        siz[mid]=r-l+1;
        return mid;
    }

    inline void modify(int x)
    {
        static int stk[M],top;
        stk[top=1]=x;
        for(int i=x;i!=rt;i=fa[i])stk[++top]=fa[i];
        while(top)push(stk[top--]);

        val[x]=0;
        while(x!=rt)rotate(x);
        update(x);
    }
}

int n,m,q;
vector<int> g[N];

int main()
{
    n=read();m=read();q=read();
    for(int i=1;i<=q;i++)
    {
        int x=read(),y=read();
        g[x].push_back(y);
    }
    for(int i=1;i<=n;i++)
        if(g[i].size())
            sort(g[i].begin(),g[i].end());
    rt=treap::build(1,m);

    ll ans=0;
    for(int i=1;i<=n;i++)
    {
        treap::mark(rt,1);
        for(int j=0,ej=g[i].size();j<ej;j++)
            treap::modify(g[i][j]);
        ans+=treap::sum[rt];
    }

    printf("%lld\n",(n*(n+1ll)/2ll)*(m*(m+1ll)/2ll)-ans);
    return 0;
}