【数据结构】 POJ3109 Inner Vertices 题解

题意

给你一个无限大的白色点网格，给出n个黑点的坐标，如果一个白点上下左右都有黑点（不一定相邻），它就会变成黑点，求图上最终一共有多少个黑点。（原题意有无限多输出-1，但是想想就知道不可能）
其中n<=10^5 坐标绝对值<=10^9

分析

首先，离散化是肯定的，不然10^9处理起来不方便。
然后就可以用到扫描线法，按坐标x从小到大，x相同y从小到大排序，然后对于x值相同的点，从上到下，比较这个点与上一个同样x的点，用树状数组记录它们之间的白点左右都有黑点的有多少个，加入答案，具体实现是一开始就把每一行上有多少个黑点sum记录下来，然后每一竖列上的点，经过上面的操作后，再挨个处理，通过一个记录第i行有多少个点now[i]，再和sum[i]比较，如果now[i]==1，就说明这是当前行的第一个点，加入树状数组中，如果now[i]==sum[i]，说明是最后一个点，从树状数组中删除，最后输出n+增加的点就可以了。

代码

#include<cstdio>
#include<algorithm>
#include<map>
#define MAXN 100006
using namespace std;
int bit[MAXN],cntx,cnty,n,ans,l,r,now[MAXN],sum[MAXN],last;
struct node
{
    int x,y;
}p[MAXN];
map<int,int> mapx,mapy;
int lowbit(int x)
{
    return x&(-x); 
} 
int addup(int x)
{
    int re=0;
    while(x)
    {
        re+=bit[x];
        x-=lowbit(x);
    }
    return re;
}
void change(int x,int d)
{
    while(x<=cnty)
    {
        bit[x]+=d;
        x+=lowbit(x);
    }
}
bool cmpx(node a,node b)
{
    if(a.x==b.x) return a.y<b.y;
    return a.x<b.x;
}
bool cmpy(node a,node b)
{
    if(a.y==b.y) return a.x<b.x;
    return a.y<b.y;
}
int main()
{
    scanf("%d",&n);
    for(int i=1;i<=n;i++)
        scanf("%d%d",&p[i].x,&p[i].y);
    sort(p+1,p+1+n,cmpy);
    for(int i=1;i<=n;i++)
    {
        if(i==1||p[i].y!=last) 
        {
            last=p[i].y;
            p[i].y=++cnty;
        }
        else p[i].y=p[i-1].y;
        sum[p[i].y]++;
    }
    sort(p+1,p+1+n,cmpx);
    for(int i=1;i<=n;i++)
    {
        if(i==1||p[i].x!=last) 
        {
            last=p[i].x;
            p[i].x=++cntx;
        }
        else p[i].x=p[i-1].x;
    }
    ans=n;
    l=1;
    r=1;
    for(int i=1;i<=cntx;i++)
    {
        while(r<n&&p[r+1].x==i)
        {
            ans+=addup(p[r+1].y-1)-addup(p[r].y);
            r++;
        }
        for(int j=l;j<=r;j++)
        {
            now[p[j].y]++;
            if(now[p[j].y]==1) change(p[j].y,1);
            if(now[p[j].y]==sum[p[j].y]) change(p[j].y,-1);
        }
        l=r=r+1;
    }
    printf("%d",ans);
}