Codeforces Round #573 (Div. 2) F. Tokitsukaze and Strange Rectangle

F. Tokitsukaze and Strange Rectangle

time limit per test:3 seconds
memory limit per test:256 megabytes
input:standard input
output:standard output

There are $n$ points on the plane, the $i$-th of which is at ($x$_i,$y$_i)

. Tokitsukaze wants to draw a strange rectangular area and pick all the points in the area.

The strange area is enclosed by three lines, $x$=$l$, $y$=$a$ and $x$=$r$, as its left side, its bottom side and its right side respectively, where $l$,$r$ and $a$ can be any real numbers satisfying that $l<r$. The upper side of the area is boundless, which you can regard as a line parallel to the $x$-axis at infinity. The following figure shows a strange rectangular area.

A point $(x$_i,$y$_i$)$ is in the strange rectangular area if and only if $l<x$_i$<r$ and $y$_i$>a$. For example, in the above figure, $p$₁ is in the area while $p$₂ is not.
Tokitsukaze wants to know how many different non-empty sets she can obtain by picking all the points in a strange rectangular area, where we think two sets are different if there exists at least one point in one set of them but not in the other.

Input

The first line contains a single integer $n$(1≤$n$≤2×10⁵) — the number of points on the plane.
The $i$-th of the next $n$ lines contains two integers $x$_i, $y$_i (1≤$x$_i,$y$_i≤10⁹) — the coordinates of the $i$-th point.
All points are distinct.
Output
Print a single integer — the number of different non-empty sets of points she can obtain.
题意
二维平面上有$n$个点，第i个点的坐标以$（x$_i，$y$_i$）$的形式给出，保证每个点的坐标都不同且均在区间$[1,10$⁹$]$内。画一个由$x=l,x=r,y=a$，上不封顶的矩形（其中至少包含一个点），求能够画出的包含的点的集合不同的矩形个数（若至少有一个点存在于集合A而不存在于集合B则认为集合A与集合B不同）
思路
显然，由于这个矩形上不封顶的性质，若有两个点A，B的x坐标相同，那么是无法做到画出的矩形中包含y坐标较小的点且不包含y坐标较大的点的。故先将点按x坐标排序，用map离散化得出x轴坐标不同的点的个数(设其为$cnt$)，再将点按y坐标排序，y坐标相同时按x坐标排序。每次处理当前所有y坐标最小的点。
对于仅有一个点的情况，设该点左边有$l$个点，右边有$r$个点，故对于左边r个点，可以选取$0$~$r$共$（r+1）$种方案，同理，右边有$（l+1）$种方案，共$（l+1）\ast （r+1）$种方案。
对于多个点的情况，处理第二个点时就不能直接计算为$（l+1）\ast （r+1）$了：假设第一个点前有$a$个$x$坐标不同的点，第二个点与第一个点间有$b$个$x$坐标不同的点，第二个点后有$c$个$x$坐标不同的点，若此时选取前$a$个点之一作为起点，后$c$个点之一作为终点，那么就与处理第一个点时的结果重复了（因为处理第一个点时也能分别以这两个点为起点和终点）。
故令$m[k]$表示$x=k$是第$m[k]$小的、至少包含一个点满足$x=k$的$x$坐标。令$sum$_ $t(i)$表示有多少个满足$x<=i$，且至少包含一个未被处理/正在被处理的点的$x$坐标。则处理任意一个点$m_i$时我们可以用$(sum$_ $t(m_i)$-$sum$_ $t(m$_i-1$))\ast ($$sum$_ $t(cnt)$-$sum$_ $t(m_i))$来表示该点且不选取任何一个与该点$y$坐标相同且x坐标小于该点的点的情况下能得到的矩形数。
将这些点处理完后，我们需要从$sum$_$t(i)$中删去这些点（因为每次被处理的点都是当前剩余点中$y$坐标最小的所有点）。
重复以上步骤直到没有剩余的点。
AC代码

#include<stdio.h>
#include<map>
#include<queue>
#include<iostream>
#include<algorithm>
#include<vector>
#define MOD 10000000007
#define N 200005
using namespace std;
map<int,int> m;
typedef struct node{
	int x,y;
}node;
node data[N];
int ct1[N],ct2[N];//ct1用于记录每个x坐标剩余的点的数量，ct2用于记录每个x坐标是否存在点 
int q[N];
int n,x,y;
long long temp,cot,ans,cnt;
inline int lowbit(int x){return (x&(-x));}//树状数组
inline long long sum_t(int n)
{
	long long result=0;
	while(n)
	{
		result+=ct2[n];n-=lowbit(n);
	}
	return result;
}//树状数组
inline void update(int pos,int val)
{
	do
	{
		ct2[pos]+=val;
		pos+=lowbit(pos);
	}while (pos<=n);
	return;
}//树状数组 
int cmp1(node a,node b)
{
	return a.x<b.x;
}
int cmp2(node a,node b)
{
	return a.y<b.y||a.y==b.y&&a.x<b.x;
}
int main()
{
	scanf("%d",&n);
	for(int i=1;i<=n;i++)
	{
		scanf("%d%d",&x,&y);
		data[i].x=x;data[i].y=y;
	}
	sort(data+1,data+n+1,cmp1);//按x坐标排序 
	for(int i=1;i<=n;i++)
	{
		if(m.find(data[i].x)==m.end()) 
		{
			m[data[i].x]=++cnt;//使用map实现离散化 
		}
		ct1[m[data[i].x]]++;
	}
	for(int i=1;i<=cnt;i++) update(i,1);
	sort(data+1,data+n+1,cmp2);//按y坐标优先，x坐标其次排序 
	int i=0,j=1;
	while(i<n&&j<=n)
	{
		i=j;int te=1;
		for(j=i;data[j].y==data[i].y;j++,te++)
		{
			q[te]=data[j].x;//寻找当前所有可行点 
		}
		for(int k=1;k<te;k++)
		ans+=(sum_t(m[q[k]])-sum_t(m[q[k-1]]))*(sum_t(cnt)-sum_t(m[q[k]])+1);//计算每个点的可行方案数 
		for(int k=1;k<te;k++)
		{
			ct1[m[q[k]]]--;//删去这个点 
			if(ct1[m[q[k]]]==0) update(m[q[k]],-1);//删去不存在点的x坐标 
		}
	}
	printf("%I64d",ans);
	return 0;
}