POJ2932 Coneology

题目：
http://poj.org/problem?id=2932
《挑战程序设计竞赛》P258

题意：
xy平面内给一堆圆，不相交，可以内含；
问所有不被包含的圆是哪些；

分析：
计算几何扫描线例题
将圆左右端点的x坐标及圆序号封装到结构体中，以x坐标排序；
按x从小到大扫描；
开一个set：outers记录目前所有最外层的圆；
扫描到圆的左端点，如果不被任何圆包含，则加入outers中；
扫描到圆的右端点，如果在outers中，则在outers中删除此圆；

那么，怎么判断一个圆是否被包含呢？
题目规定：所有圆不相交；
对于一个圆O，假设在其上方有两个圆在outers中，圆心y坐标分别是y1,y2(y1< y2)，那么O一定不被圆y2包含
反证法证明如下：

扫描线扫描到圆O的左端点，y=y1只能在y0~y2之间动；
那么圆y1要么和圆y2相交，要么被圆y2包含；
然而，题目保证所有圆不相交；
而且若圆y1、y2均存在于outers中，则不存在包含关系；
因此，对于一个圆，只可能被y坐标第一个比它大的(outers中的)圆包含；
同理，只能被y坐标第一个比他小的(outers中的)圆包含；
那么每次在outer中二分查找就行啦；

代码：

#include <cstdio>
#include <algorithm>
#include <vector>
#include <cstring>
#include <set>
using namespace std;
const int Tmax=40005;
int n;
double r[Tmax],x[Tmax],y[Tmax];
vector<pair<double,int> > events;  //左右端点的x坐标 
vector<int> ans;                 //最外圆列表，即答案列表 
set<pair<double,int> > outers;   //与扫描线相交的最外层圆集合,存其圆心y坐标 
set<pair<double,int> >::iterator it;
bool inside(int i,int j)
{
    double dx=x[i]-x[j],dy=y[i]-y[j];
    return dx*dx+dy*dy<=r[j]*r[j];
}
void work()
{
    int i,id;
    for(i=0;i<events.size();i++)
    {
        id=events[i].second %n;
        if(events[i].second<n)   //扫描左端点 
        {
            it=outers.lower_bound(make_pair(y[id],id));   //找到第一个大于等于y[id]的圆
            if(it!=outers.end()&&inside(id,it->second)==true) continue;
            if(it!=outers.begin()&&inside(id,(--it)->second)==true) continue;
            ans.push_back(id);
            outers.insert(make_pair(y[id],id)); 
        }
        else                     //扫到右端点 
          outers.erase(make_pair(y[id],id));
    }
    return;
}
int main()
{
    int i;
    scanf("%d",&n);
    for(i=0;i<n;i++)
    {
        scanf("%lf%lf%lf",&r[i],&x[i],&y[i]);
        events.push_back(make_pair(x[i]-r[i],i));
        events.push_back(make_pair(x[i]+r[i],i+n));
    }
    sort(events.begin(),events.end());
    work();
    sort(ans.begin(),ans.end());
    printf("%d\n",ans.size());
    for(i=0;i<ans.size();i++)
      printf("%d ",ans[i]+1);
    return 0;
}