BZOJ 1818 Cqoi2010 内部白点 树状数组

题目大意：给定平面上的一些黑点，其它位置都是白点，一个白点如果上下左右都有黑点就会变成黑点，求最终会有多少个黑点

语文题？

总之我们现在得到了一些横向和纵向的线段(注意线段不能包含两端点！)，求出交点数后+n就是答案

我们可以将横向线段看做修改，纵向线段看做询问，那么一个询问可以拆成两个

然后我们离线做，将询问按y坐标排序，对于每个询问将y坐标小于等于这个询问的修改都加入树状数组，然后查询这个询问的x坐标即可

懒得离散化，第一次尝试用Hash表写树状数组，效果还不错= =

时间复杂度O(nlogn)

#include <cstdio>#include <cstring>#include <iostream>#include <algorithm>#define M 100100using namespace std;struct Point{int x,y;friend istream& operator >> (istream &_,Point &p){scanf("%d%d",&p.x,&p.y);p.x+=1000000007;p.y+=1000000007;return _;}}points[M];struct Interval{int l,r,y;Interval() {}Interval(int _,int __,int ___):l(_),r(__),y(___) {}bool operator < (const Interval &i) const{return y < i.y;}}intervals[M],*I=intervals;struct Query:Point{int flag;Query() {}Query(int _,int __,int ___):flag(___){x=_;y=__;}}queries[M<<1],*Q=queries;int n;long long ans;bool Compare1(const Point &p1,const Point &p2){if(p1.x!=p2.x)return p1.x<p2.x;return p1.y<p2.y;}bool Compare2(const Point &p1,const Point &p2){if(p1.y!=p2.y)return p1.y<p2.y;return p1.x<p2.x;}namespace BIT{struct Hash_Set{struct List{int hash,val;List *next;List(int _,List *__):hash(_),val(0),next(__){}}*head[1001001];int& operator [] (int hash){int pos=hash%1001001;List *temp;for(temp=head[pos];temp;temp=temp->next)if(temp->hash==hash)return temp->val;return (head[pos]=new List(hash,head[pos]))->val;}}c;void Update(int x,int y){for(;x;x-=x&-x)c[x]+=y;}int Get_Ans(long long x){int re=0;for(;x<=2000000007;x+=x&-x)re+=c[x];return re;}}int main(){using namespace BIT;int i;cin>>n;for(i=1;i<=n;i++)cin>>points[i];sort(points+1,points+n+1,Compare2);for(i=2;i<=n;i++)if(points[i].y==points[i-1].y&&points[i].x>=points[i-1].x+2)new (I++)Interval(points[i-1].x+1,points[i].x-1,points[i].y);sort(points+1,points+n+1,Compare1);for(i=2;i<=n;i++)if(points[i].x==points[i-1].x&&points[i].y>=points[i-1].y+2){new (Q++)Query(points[i].x,points[i-1].y,-1);new (Q++)Query(points[i].x,points[i].y-1,1);}sort(intervals,I);sort(queries,Q,Compare2);Interval *_i=intervals;Query *_q=queries;for(;_q<Q;_q++){for(;_i<I&&_i->y<=_q->y;_i++)Update(_i->l-1,-1),Update(_i->r,1);ans+=Get_Ans(_q->x)*_q->flag;}cout<<ans+n<<endl;return 0;}