bzoj 1818 Cqoi2010 内部白点 扫描线

[Cqoi2010]内部白点

Time Limit: 10 Sec  Memory Limit: 64 MB
Submit: 1126  Solved: 530
[Submit][Status][Discuss]

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

然后就是扫描线求交点了。

 #include<cstring>
 #include<cmath>
 #include<algorithm>
 #include<iostream>
 #include<cstdio>

 #define N 100001
 #define M 1000001
 #define ll long long
 using namespace std;
 inline int read()
 {
     int x=,f=;char ch=getchar();
     while(ch>''||ch<''){if (ch=='-') f=-;ch=getchar();}
     while(ch<=''&&ch>=''){x=(x<<)+(x<<)+ch-'';ch=getchar();}
     return x*f;
 }

 int n,cnt,ans,hash[N],tr[N];
 struct point{int x,y;}a[N];
 struct seg{int k,x,y,r;}s[M];

 inline bool cmp1(point a,point b){if(a.x==b.x){return a.y<b.y;}return a.x<b.x;}
 inline bool cmp2(point a,point b){if(a.y==b.y){return a.x<b.x;}return a.y<b.y;}
 inline bool cmp3(seg a,seg b)
 {
     if(a.y==b.y)return a.k<b.k;
     return a.y<b.y;
 }
 int find(int x)
 {
     int l=,r=n,mid;
     while(l<=r)
     {
         int mid=(l+r)>>;
         if(hash[mid]<x)l=mid+;
         else if(hash[mid]>x)r=mid-;
         else return mid;
     }
 }
 void insert(int k,int l,int r,int t)//0横线,1竖线
 {
     if(!k){s[++cnt].x=find(l);s[cnt].r=find(r);s[cnt].y=t;}
     else{
         s[++cnt].x=find(t);s[cnt].y=l;s[cnt].k=;
         s[++cnt].x=find(t);s[cnt].y=r;s[cnt].k=-;
     }
 }
 void build()
 {
     sort(a+,a+n+,cmp1);
     for(int i=;i<=n;i++)
         if(a[i].x==a[i-].x)
             insert(,a[i-].y,a[i].y,a[i].x);
     sort(a+,a+n+,cmp2);
     for(int i=;i<=n;i++)
         if(a[i].y==a[i-].y)
             insert(,a[i-].x,a[i].x,a[i].y);
 }
 int lowbit(int x){return x&(-x);}
 void update(int x,int y)
 {
     while(x<=n)
     {
         tr[x]+=y;
         x+=lowbit(x);
     }
 }
 int ask(int x)
 {
     int s=;
     while(x)
     {
         s+=tr[x];
         x-=lowbit(x);
     }
     return s;
 }
 void work()
 {
     for(int i=;i<=cnt;i++)
     {
         if(!s[i].k)ans+=ask(s[i].r-)-ask(s[i].x);
         else update(s[i].x,s[i].k);
     }
 }
 int main()
 {
     scanf("%d",&n);
     for(int i=;i<=n;i++)
     {
         scanf("%d%d",&a[i].x,&a[i].y);
         hash[i]=a[i].x;
     }
     sort(hash+,hash+n+);
     build();
     sort(s+,s+cnt+,cmp3);
     work();
     printf("%d",ans+n);
 }