Codeforces1190 D. Tokitsukaze and Strange Rectangle（扫描线+树状数组）

题意：

数据范围：n<=2e5，1<=xi,yi<=1e9

解法：

考虑枚举矩形下边界(横线),之后再考虑左右端点.
对点按照y从大到小排序,y相等时按照x从小到大排序.
因为点只有n个,显然横线也最多只有n种取法,
从上到下枚举横线y时,开树状数组丢入纵坐标<=当前y的所有点的横坐标(相同的只统计一次).

1.假设纵坐标=当前y的点只有一个,设当前点为(x,y),因为当前点(x,y)必须选,
那么满足条件的左端点数量为T.ask(x-1)+1,满足条件的右端点数量为T.ask(n)-T.ask(x)+1.
相乘就是当前横线的贡献.
2.假设纵坐标=当前y的点不止一个,那么稍微麻烦一点,不能直接按上面那样计算了,
首先将纵坐标=当前y的点(xi,y)的所有横坐标xi加入树状数组.
接下来看图示:

假设当前树状数组的桶图是这样的：

对于x1，满足条件的左端点数量为2，满足条件的右端点数量为6，贡献为2*6=12。
对于x2，x2右端点集合是之前x1右端点集合的子集，x1左端点集合是x2左端点集合的子集，
如果还是按照上面的方法计算的话，会出现贡献重复计算的情况（左右端点相同），
因此x2左端点的取值数量只有：x2左端点集合对x1左端点集合的差，即[x1+1,x2-1]范围内点的数量+1.

code：

#include <bits/stdc++.h>
using namespace std;
#define int long long
#define PI pair<int,int>
const int maxm=2e5+5;
struct Node{
    int x,y;
}e[maxm];
struct BIT{
    int c[maxm];
    inline int lowbit(int i){
        return i&-i;
    }
    void add(int i,int t){
        while(i<maxm){
            c[i]+=t;
            i+=lowbit(i);
        }
    }
    int ask(int i){
        int ans=0;
        while(i){
            ans+=c[i];
            i-=lowbit(i);
        }
        return ans;
    }
}T;
int mark[maxm];
int b[maxm];
int n;
bool cmp(Node a,Node b){//按y从大到小排序
    if(a.y!=b.y)return a.y>b.y;
    return a.x<b.x;
}
signed main(){
    ios::sync_with_stdio(0);cin.tie(0);
    cin>>n;
    for(int i=1;i<=n;i++){
        cin>>e[i].x>>e[i].y;
        b[i]=e[i].x;
    }
    sort(e+1,e+1+n,cmp);
    sort(b+1,b+1+n);
    int num=unique(b+1,b+1+n)-b-1;
    for(int i=1;i<=n;i++){
        e[i].x=lower_bound(b+1,b+1+num,e[i].x)-b;
    }
    int ans=0;
    for(int i=1;i<=n;i++){//枚举横线y
        //存在相同的y
        int j=i;
        while(j<=n&&e[j].y==e[i].y)j++;
        j--;
        for(int k=i;k<=j;k++){
            if(!mark[e[k].x]){//因为相同的只算一个,所以要打标记
                T.add(e[k].x,1);
                mark[e[k].x]=1;
            }
        }
        //
        int l=T.ask(e[i].x-1)+1;
        int r=T.ask(num)-T.ask(e[i].x)+1;
        ans+=l*r;
        for(int k=i+1;k<=j;k++){
            l=max(0LL,T.ask(e[k].x-1)-T.ask(e[k-1].x))+1;
            r=T.ask(num)-T.ask(e[k].x)+1;
            ans+=l*r;
        }
        //
        i=j;
    }
    cout<<ans<<endl;
    return 0;
}

---