POJ2932

又是烦死人的几何题……

题意是给出平面上若干个没有公共点的圆，问哪些圆不被其他圆包含……

做法是平面扫描法，将所有圆左右端点拿出来排序，然后从左到右用map维护现在包括的圆的所有圆心纵坐标和对应的圆的编号。每次更新只要找与现在添加的圆纵坐标最接近的圆判断是否包含即可。注意最开始排序的时候是先排左端点再排右端点，因为有可能出现半径为0的情况……

给出代码吧这样说真乱……

#include<stdio.h>
#include<vector>
#include<set>
#include<map>
#include<math.h>
#include<string.h>
#include<algorithm>
using namespace std;
#define sqr(x) (x)*(x)

int n;
double x[100005],y[100005],r[100005];

typedef pair<double,int> point;

map<double,int> ma;

int bh(int i,int j)
{
	if(!i||!j) return 0;
	if(i==n+1||j==n+1) return 0;
	double dis=sqrt(sqr(x[i]-x[j])+sqr(y[i]-y[j]));
	if(dis>r[i]+r[j]+1e-8) return 0;
	if(dis<r[i]-r[j]+1e-8) return 1;
	return -1;	
}
bool cmp(point t1,point t2)
{
	if(fabs(t1.first-t2.first)>1e-8) return t1.first<t2.first;
	return t1.second>t2.second;
}

int vis[500005];

point c[1000005];

int main()
{
	scanf("%d",&n);
	x[0]=y[0]=-1e12;x[n+1]=y[n+1]=1e12;r[0]=r[n+1]=0;
	for(int i=1;i<=n;i++)
	{
		scanf("%lf %lf %lf",r+i,x+i,y+i);
		//r[i]=i*i%103+134;x[i]=y[i]=i*i%100;
		c[i+i-1]=point(x[i]-r[i],i);
		c[i+i]=point(x[i]+r[i],-i);
	}
	memset(vis,-1,sizeof(vis));
	sort(c+1,c+n+n,cmp);
	ma.insert(make_pair(-1e12,0));ma.insert(make_pair(1e12,n+1));
	for(int i=1;i<=n+n;i++)
	{
		int nut=c[i].second;
		if(nut<0) {nut=-nut;if(vis[nut]) ma.erase(y[nut]);continue;}
		if(ma.find(y[nut])!=ma.end())
		{
			if(r[ma[y[nut]]]>r[nut]) vis[nut]=0;
			else {vis[ma[y[nut]]]=0;ma.erase(y[nut]);ma[y[nut]]=nut;}
			continue;
		}
		ma[y[nut]]=nut;
		map<double,int>::iterator ite=ma.find(y[nut]),it2=ite;
		--it2;
		int ttt=bh(nut,it2->second);
		if(ttt==-1)
		{
			point tst=*it2;
			vis[nut]=0;
			ma.erase(y[nut]);
			continue;
		}
		if(ttt==1)
		{
			while(ttt==1)
			{
				vis[it2->second]=0;
				if(it2->second>0&&it2->second<n+1)
				{
					ma.erase(it2->first);
				}
				else break;
				ite=it2=ma.find(y[nut]);--it2;
				ttt=bh(nut,it2->second);
			}
		}
		ite=it2=ma.find(y[nut]);++it2;
		ttt=bh(nut,it2->second);
		if(ttt==-1)
		{
			point tst=*it2;
			vis[nut]=0;
			ma.erase(y[nut]);
			continue;
		}
		if(ttt==1)
		{
			while(ttt==1)
			{
				vis[it2->second]=0;
				if(it2->second>0&&it2->second<n+1)
				{
					ma.erase(it2->first);
				}
				else break;
				ite=it2=ma.find(y[nut]);++it2;
				ttt=bh(nut,it2->second);
			}
		}
	}
	int sz=0;
	for(int i=1;i<=n;i++) if(vis[i]) sz++;
	printf("%d\n",sz);
	for(int i=1,j=0;i<=n;i++) if(vis[i]) {if(j)printf(" ");++j;printf("%d",i);}
	printf("\n");
//	scanf("%d",&n);
	return 0;
}