【bzoj1818】[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

题解
我们只需求出有多少个整点满足上下左右都存在黑点。
也就是计算线段的交点。
采用扫描线算法。大概就是把线段按照y端点排序，类似查分数组，一端+1，一端-1。然后查询x轴即可。

代码

#include<bits/stdc++.h>
#define N 100005
#define ll long long
#define mo 1000000007
using namespace std;
int t[N],n,cnt,ans,hash[100001];
struct node1{int x,y;}a[N];
struct node2{int x,y,r,k;}s[N*5];
inline bool cmp1(node1 a,node1 b)
{if(a.x==b.x){return a.y<b.y;}return a.x<b.x;}
inline bool cmp2(node1 a,node1 b)
{if(a.y==b.y){return a.x<b.x;}return a.y<b.y;}
inline bool cmp3(node2 a,node2 b)
{
    if (a.y==b.y) return a.k<b.k;
    return a.y<b.y;
}
inline int read()
{
    int x=0,f=1;char ch=getchar();
    while(ch<'0'||ch>'9'){if(ch=='-')f=-1;ch=getchar();}
    while(ch>='0'&&ch<='9'){x=x*10+ch-'0';ch=getchar();}
    return x*f;
}
int find(int x)
{
    int l=1,r=n,mid;
    while(l<=r)
    {
        int mid=(l+r)>>1;
        if(hash[mid]<x)l=mid+1;
        else if(hash[mid]>x)r=mid-1;
        else return mid;
    } 
}
void insert(int k,int l,int r,int p)
{
    if (!k)
    {
        s[++cnt].x=find(l);
        s[cnt].r=find(r);
        s[cnt].y=p;
    }
    else
    {
        s[++cnt].x=find(p);s[cnt].y=l;s[cnt].k=1;
        s[++cnt].x=find(p);s[cnt].y=r;s[cnt].k=-1;
    }
}
void build()
{
    sort(a+1,a+n+1,cmp1);
    for (int i=2;i<=n;i++)
        if (a[i].x==a[i-1].x) insert(1,a[i-1].y,a[i].y,a[i].x);
    sort(a+1,a+n+1,cmp2);
    for (int i=2;i<=n;i++)
        if (a[i].y==a[i-1].y) insert(0,a[i-1].x,a[i].x,a[i].y);
}
int lowbit(int x){return x&(-x);}
int getsum(int x)
{
    int ans=0;
    for (int i=x;i;i-=lowbit(i))ans+=t[i];
    return ans;
}
void update(int x,int k)
{
    for (int i=x;i<=n;i+=lowbit(i))t[i]+=k;
}
void work()
{
    for (int i=1;i<=cnt;i++)
    {
        if (!s[i].k)ans+=getsum(s[i].r-1)-getsum(s[i].x);
        else update(s[i].x,s[i].k);
    }
}
int main()
{
    n=read();
    for (int i=1;i<=n;i++)
    {
        a[i].x=read();a[i].y=read();
        hash[i]=a[i].x;
    }
    sort(hash+1,hash+n+1);
    build();
    sort(s+1,s+cnt+1,cmp3);
    work();
    cout<<ans+n<<endl;
    return 0;
}