A simple problem

 

题目地址

http://acm.hust.edu.cn/thx/problem.php?id=1425

题目描述：给定m个圆，n个点，求对于每个点，有多少圆包含了它（边界上不算）ans[i] 

对于任何一个点i，ans[i]保证<=100

分析：核心思想是扫描线，对每个圆，拆分成两个点，左边点+和右边点-

当一条扫描线经过时，如果遇到+,那么加入该点，如果遇到-点，删除该点

复杂度分析，如果没有最后一句话，这个题的复杂度是不能保证的，但是由于<=100的限制，也就是说

在扫描线滑动的任何一个时刻，对于一个询问某点的操作，其代价是不会很大的，

因为里面的存在的事件点不会很多的,当然这也只是从一个整体上的观察，实际上任然无法保证复杂度

但是由于其个数的限制，而使得答案存在且很稀疏，如果有更好的算法欢迎与我讨论 

主要过程，得到事件点，然后遍历所有点，吧x坐标<当前点的全都放到或删除在容器中

遍历容器，计数

//题目描述：给定m个圆，n个点，求对于每个点，有多少圆包含了它（边界上不算）ans[i] //对于任何一个点i，ans[i]保证<=100 //分析：核心思想是扫描线，对每个圆，拆分成两个点，左边点+和右边点- //当一条扫描线经过时，如果遇到+,那么加入该点，如果遇到-点，删除该点 //复杂度分析，如果没有最后一句话，这个题的复杂度是不能保证的，但是由于<=100的限制，也就是说 //在扫描线滑动的任何一个时刻，对于一个询问某点的操作，其代价是不会很大的， //因为里面的存在的事件点不会很多的,当然这也只是从一个整体上的观察，实际上任然无法保证复杂度 //但是由于其个数的限制，而使得答案存在且很稀疏，如果有更好的算法欢迎与我讨论 //主要过程，得到事件点，然后遍历所有点，吧x坐标<当前点的全都放到或删除在容器中 //遍历容器，计数 #include<iostream> #include<cstdio> #include<set> using namespace std; const int maxn=50000,maxm=100000; typedef long long LL; int n,m,len,ans[maxn]; struct point{ int x,y,id; point(int xx=0,int yy=0,int ii=0){x=xx;y=yy;id=ii;} bool operator<(const point &p)const{ return x<p.x; } }p[maxn]; struct event{ int x,v,id; event(int xx=0,int vv=0,int ii=0){x=xx;v=vv;id=ii;} bool operator<(const event &p)const{ return x<p.x; } }e[2*maxm]; struct circle{ int x,y,r; circle(int xx=0,int yy=0,int rr=0){x=xx;y=yy;r=rr;} bool in(const point &p){ return (LL)(p.x-x)*(p.x-x)+(LL)(p.y-y)*(p.y-y)<(LL)r*r; } }c[maxm]; set<int>mset; set<int>::iterator ii; int main(){ int test; scanf("%d",&test); for(int caseid=1;caseid<=test;caseid++){ scanf("%d",&m); for(int i=0;i<m;i++)scanf("%d%d%d",&c[i].x,&c[i].y,&c[i].r); scanf("%d",&n); for(int i=0;i<n;i++){scanf("%d%d",&p[i].x,&p[i].y);p[i].id=i;} len=0; for(int i=0;i<m;i++){ e[len++]=event(c[i].x-c[i].r,1,i); e[len++]=event(c[i].x+c[i].r,0,i); } sort(p,p+n); sort(e,e+len); mset.clear(); for(int i=0,j=0;i<n;i++){ for(;e[j].x<p[i].x;j++){ if(e[j].v)mset.insert(e[j].id); else mset.erase(e[j].id); } int k=p[i].id; for(ans[k]=0,ii=mset.begin();ii!=mset.end();ii++){ if(c[*ii].in(p[i]))ans[k]++; } } printf("Case #%d:",caseid); for(int i=0;i<n;i++)printf(" %d",ans[i]); puts(""); } return 0; }