当前位置: 代码迷 >> 综合 >> poj 2932 Coneology 扫面线法
  详细解决方案

poj 2932 Coneology 扫面线法

热度:91   发布时间:2024-01-19 05:34:32.0

题意:

给n个无交点的圆,求这n个圆中不被其它圆包含的圆。

分析:

扫面线法,用二叉树(set+lowerbound方法)维护最外圆的集合。

代码:


#include <iostream>
#include <algorithm>
#include <vector>
#include <set>
using namespace std;
const int maxN=40012;
double r[maxN],x[maxN],y[maxN];
int n;
vector<pair<double,int> > events;
set<pair<double,int> > outers;
vector<int> res;bool inside(int i,int j)
{double dx=x[i]-x[j],dy=y[i]-y[j],dr=r[i]-r[j];return r[i]<=r[j]&&dx*dx+dy*dy<=dr*dr;	
}int main()
{scanf("%d",&n);for(int i=0;i<n;++i)scanf("%lf%lf%lf",&r[i],&x[i],&y[i]);for(int i=0;i<n;++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());for(int i=0;i<events.size();++i){int ids=events[i].second%n;if(events[i].second<n){set<pair<double,int> >::iterator it=outers.lower_bound(make_pair(y[ids],ids));if(it!=outers.end()&&inside(ids,it->second)) continue;if(it!=outers.begin()&&inside(ids,(--it)->second)) continue;res.push_back(ids);outers.insert(make_pair(y[ids],ids));}elseouters.erase(make_pair(y[ids],ids));			}sort(res.begin(),res.end());printf("%d\n",res.size());for(int i=0;i<res.size();++i)printf("%d ",res[i]+1);return 0;	
}