POJ 2932 扫描线

第一次做扫描线的题，这应该还算是道比较简单的。。。

看题的时候没有注意到input里说没有两个圆是相交的，纠结半天。

因为只要判断一个圆是否被其他圆所包含。那么我们只需要找当前这个圆，左端点<=当前圆的左端点 ，那些圆中y在它附近的圆是否完全包含了它，就能判断它是否被别的圆所包含。

#include<cstdio>
#include<cstring>
#include<algorithm>
#include<set>
#define maxl 40010
#define eps 1e-8
using namespace std;

typedef pair<double,int> p;

int n,cnt,ans;
double x[maxl],y[maxl],r[maxl];
bool in[maxl];
p a[maxl*2];
set<p> s;
set<p> :: iterator it;

inline void prework()
{
	scanf("%d",&n);
	for(int i=1;i<=n;i++)
		scanf("%lf%lf%lf",&r[i],&x[i],&y[i]);
	for(int i=1;i<=n;i++)
	{
		a[++cnt]=make_pair(x[i]-r[i],i);
		a[++cnt]=make_pair(x[i]+r[i],i+n);
	}
}

inline bool ainb(int a,int b)
{
	return (x[a]-x[b])*(x[a]-x[b])+(y[a]-y[b])*(y[a]-y[b])<=r[b]*r[b];
}

inline void mainwork()
{
	sort(a+1,a+1+cnt);
	for(int i=1;i<=cnt;i++)
	{
		int ind=a[i].second;
		if(ind<=n)
		{
			it=s.lower_bound(make_pair(y[ind],ind));
			if(it!=s.end() && ainb(ind,it->second))
				continue;
			if(it!=s.begin())
			{
				--it;
				if(ainb(ind,it->second))
					continue;
			}
			in[ind]=true;ans++;
			s.insert(make_pair(y[ind],ind));		
		}
		else
			s.erase(make_pair(y[ind-n],ind-n));
	}
}

inline void print()
{
	printf("%d\n",ans);
	for(int i=1;i<=n;i++)
	if(in[i])
		printf("%d ",i);
}

int main()
{
	prework();
	mainwork();
	print();
	return 0;
}