BZOJ_1818_[Cqoi2010]内部白点 _扫描线+树状数组

BZOJ_1818_[Cqoi2010]内部白点 _扫描线+树状数组

Description

无限大正方形网格里有n个黑色的顶点，所有其他顶点都是白色的（网格的顶点即坐标为整数的点，又称整点）。每秒钟，所有内部白点同时变黑，直到不存在内部白点为止。你的任务是统计最后网格中的黑点个数。 内部白点的定义：一个白色的整点P(x,y)是内部白点当且仅当P在水平线的左边和右边各至少有一个黑点（即存在x1 < x < x2使得(x1,y)和(x2,y)都是黑点），且在竖直线的上边和下边各至少有一个黑点（即存在y1 < y < y2使得(x,y1)和(x,y2)都是黑点）。

Input

输入第一行包含一个整数n，即初始黑点个数。以下n行每行包含两个整数(x,y)，即一个黑点的坐标。没有两个黑点的坐标相同，坐标的绝对值均不超过109。

Output

输出仅一行，包含黑点的最终数目。如果变色过程永不终止，输出-1。

Sample Input

4
0 2
2 0
-2 0
0 -2

Sample Output

5

数据范围
36%的数据满足：n < = 500
64%的数据满足：n < = 30000
100%的数据满足：n < = 100000

---

先排两边序，取出所有两个端点是黑点且极长的横/竖线。

可以发现新增的黑点一定是这些线的交点。

于是把横线拆成两个。

全部按横坐标排序，遇到横线决定+1还是-1，遇到竖线就查一次区间和，这个可以用树状数组搞一搞。

 

代码：

#include <cstdio>
#include <cstring>
#include <algorithm>
using namespace std;
#define N 300050
int n,V[N],c[N];
struct P {
	int x,y;
}a[N];
inline bool cmp1(const P &x,const P &y) {return x.x==y.x?x.y<y.y:x.x<y.x;}
inline bool cmp2(const P &x,const P &y) {return x.y==y.y?x.x<y.x:x.y<y.y;}
struct A {
	int x,y,z,opt,pri;
	bool operator < (const A &u) const {
		return x==u.x?pri<u.pri:x<u.x;
	}
}q[N];
void fix(int x,int v) {for(;x<=n;x+=x&(-x)) c[x]+=v;}
int inq(int x) {int re=0;for(;x;x-=x&(-x)) re+=c[x]; return re;}
int main() {
	scanf("%d",&n);
	int i,cnt=0;
	for(i=1;i<=n;i++) scanf("%d%d",&a[i].x,&a[i].y),V[i]=a[i].y;
	sort(V+1,V+n+1);
	int ww=unique(V+1,V+n+1)-V-1;
	for(i=1;i<=n;i++) a[i].y=lower_bound(V+1,V+ww+1,a[i].y)-V;
	sort(a+1,a+n+1,cmp1);
	for(i=1;i<=n;i++) {
		int j=i;
		for(;i<n&&a[i+1].x==a[j].x;i++);
		q[++cnt]=(A){a[j].x,a[j].y,a[i].y,1,2};
	}
	sort(a+1,a+n+1,cmp2);
	for(i=1;i<=n;i++) {
		int j=i;
		for(;i<n&&a[i+1].y==a[j].y;i++);
		q[++cnt]=(A){a[j].x,a[j].y,1,0,1};
		q[++cnt]=(A){a[i].x,a[j].y,-1,0,3};
	}
	sort(q+1,q+cnt+1);
	int ans=0;
	for(i=1;i<=cnt;i++) {
		if(q[i].opt==0) fix(q[i].y,q[i].z);
		else {
			ans+=inq(q[i].z)-inq(q[i].y-1);
		}
	}
	printf("%d\n",ans);
}