【BZOJ1818】【CQOI2010】【XSY2428】内部白点（树状数组+扫描线）

先把所有点的\(x\)坐标离散化。

然后分别将所有点按\(x\)、\(y\)排序。这里以按\(x\)排序为例，对于\(x\)坐标相同的两个点，我们把它们连成一条线段。那么按\(y\)坐标排序也一样，把\(y\)坐标相同的两个点也连成一条线段。

那么连出来后的图就是这样的：

那么横竖线段的所有交点（图中蓝点）即为可以变dark的点，因为它左右有dark点，上下都有dark点，符合变dark条件。

那么我们怎么维护交点呢？我们考虑用扫描线（如图黑线）。

我们先把横线、竖线的两个端点全部塞进一个数组里，然后把它们按\(y\)坐标排序。

然后从\(y\)坐标的小到大一个一个向上扫描，分3种情况讨论：

1. 如果扫描到一条竖线的下端点\((x,y)\)，那么我们把树状数组里的位置\(x\)加\(1\)，即\(add(x,1)\)。

2. 如果扫描到一条竖线的上端点\((x,y)\)，那么我们把树状数组里的位置\(x\)减\(1\)，即\(add(x,-1)\)。那么对于某一条竖线\(l\)，其上端点为\((x,y_1)\)，下端点为\((x,y_2)\)，我们就能保证当扫描线在区间\([y_1,y_2]\)时树状数组中位置\(x\)会\(+1\)。

3. 如果扫描到一条横线，且两端点横坐标为\(x_1\)、\(x_2\)（\(x_1<x_2\)），我们只要查询\(query(x_2)-query(x_1-1)\)就可以得到这条线段上与其它竖线的交点坐标了，不过由于竖线与横线两端点的交点不能算进去，所以应该查询\(query(x_2-1)-query(x_1)\)。

那么最后答案就是所有交点的个数了。

代码如下：

#include<bits/stdc++.h>
 
#define N 100010
 
using namespace std;
 
struct Point
{
    int x,y;
}p[N];
 
struct Line 
{
    int opt,x,r,y;//1 竖线下端点 0 横线 -1 竖线上端点
}l[N*10];
 
int n,ans,cnt,b[N],c[N];
 
 //树状数组
int lowbit(int x)
{
    return x&-x;
}
 
void add(int x,int y)
{
    for(;x<=n;x+=lowbit(x))c[x]+=y;
}
 
int query(int x)
{
    int ans=0;
    for(;x;x-=lowbit(x))ans+=c[x];
    return ans;
}
 
bool cmpx(Point a,Point b){return a.x==b.x?a.y<b.y:a.x<b.x;}
bool cmpy(Point a,Point b){return a.y==b.y?a.x<b.x:a.y<b.y;}
bool cmpLine(Line a,Line b){return a.y==b.y?a.opt<b.opt:a.y<b.y;}
 
void work()
{
    sort(p+1,p+n+1,cmpx);
    for(int i=1;i<n;i++)
    {
        if(p[i].x==p[i+1].x)
        {
            l[++cnt]=(Line){1,p[i].x,0,p[i].y};
            l[++cnt]=(Line){-1,p[i].x,0,p[i+1].y};
        }
    }
    sort(p+1,p+n+1,cmpy);
    for(int i=1;i<n;i++)
        if(p[i].y==p[i+1].y)
            l[++cnt]=(Line){0,p[i].x,p[i+1].x,p[i].y};
}
 
int main()
{
    scanf("%d",&n);
    for(int i=1;i<=n;i++)
    {
        scanf("%d%d",&p[i].x,&p[i].y);
        b[i]=p[i].x;
    }   
    sort(b+1,b+n+1);//离散化
    int tot=unique(b+1,b+n+1)-b-1;
    for(int i=1;i<=n;i++)
        p[i].x=lower_bound(b+1,b+tot+1,p[i].x)-b;
  	work();//加入线段
    sort(l+1,l+cnt+1,cmpLine);//扫描
    for(int i=1;i<=cnt;i++)
    {
        if(!l[i].opt)ans+=query(l[i].r-1)-query(l[i].x);
        else add(l[i].x,l[i].opt);
    }
    printf("%d\n",ans+n);
    return 0;
}