[扫描线]POJ2932 Coneology
题意:有n个圆 依次给了半径和圆心坐标 保证输入的圆不相交(只有 相离 和 内含/外含 的情况)
问 有几个圆 不内含在其他圆中,并分别列出这几个圆的编号(1~n)
(n的范围是[1, 40000])
案例画出来大概是这样的
(那个原点为(50,50)的太远了,就意思一下)
所以答案是3号圆和5号圆 不被包含
好了,若这道题n只有1000,那么只要for两层,每个圆与另外的圆比较, 判断圆心是否在其他圆内即可判断是否包含
这样的复杂度是O($n^2$)
可是现在n有40000,显然不能用O($n^2$)来解决
由这道题,简单粗暴的学习了一下扫描线, 复杂度为O(nlogn)
什么是扫描线呢?
顾名思义,就是一根线,扫描过去。
怎样的一根线,怎么扫过去呢?
①平行于x轴,自上而下/自下而上 扫描
②平行于y轴,自左而右/自右而左 扫描
③绕圆心 逆时针/顺时针 扫描
扫描线要干什么呢?
“扫描线在平面上按给定轨迹移动的同时,不断根据扫描线扫过部分更新信息,从而得到整体所要求的结果”
这道题,可以用上述的①/②
自左而右扫描时,只有当扫描线移动到圆的左右两端时,线与圆的关系才会发生改变,因此先把圆的左右端点取出来排个序
每当遇到某圆的左端点,判断该圆是否包含在其他圆内
因为所有的圆都不相交,因此,每个圆都只可能包含在上下两个与它最相近的圆中
(此处“上下两个”是通过比较 圆心的纵坐标 来确定的)
也就是 每当我们得到一个不包含在其他圆中的圆,我们需要将它存起来,并将这些圆按圆心的纵坐标排序 以便与下一个扫到的圆进行比较
我们可以用set<pair<double, int> >来存 (pair.first是圆心的纵坐标,pair.second是圆的编号) set会自动按pair.first排序
当我们扫到某圆的右端点时,表示该圆的影响范围已经扫完了,后面扫到的圆不可能包含在该圆中,因此可以把该圆从set中去掉
查找“上下两个与它最近的圆”的复杂度为O(logn)
因此遍历n个圆的复杂度为O(nlogn)
const int N=; double x[N], y[N] ,r[N]; bool inside(int i, int j) // i是否在j内
{
double dx=x[i]-x[j], dy=y[i]-y[j];
return dx*dx+dy*dy<=r[j]*r[j];
} int main()
{
int n;
while(~scanf("%d", &n))
{
for(int i=;i<n;i++)
scanf("%lf%lf%lf", &r[i], &x[i], &y[i]);
vector<pair<double, int> > X;
for(int i=;i<n;i++)
{
X.push_back(make_pair(x[i]-r[i], i)); // 左
X.push_back(make_pair(x[i]+r[i], i+n));// 右
}
sort(X.begin(), X.end());
set<pair<double, int> > out;
vector<int> ans;
for(int i=;i<X.size();i++) // 从左到右扫描
{
int id=X[i].second%n;
if(X[i].second<n) // 左
{
set<pair<double, int> >::iterator it=out.lower_bound(make_pair(y[id], id));
if(it!=out.end() && inside(id, it->second)) // id 在 前一个圆内 不加入
continue;
if(it!=out.begin() && inside(id, (--it)->second)) // id 在 后一个圆内 不加入
continue;
ans.push_back(id);
out.insert(make_pair(y[id], id));
}
else // 右
out.erase(make_pair(y[id], id));
}
sort(ans.begin(), ans.end());
printf("%d\n", ans.size());
for(int i=;i<ans.size();i++)
printf("%d%c", ans[i]+, i+==ans.size()? '\n':' ');
}
return ;
}
POJ 2932
[扫描线]POJ2932 Coneology的更多相关文章
- POJ2932 Coneology【圆扫描线】
POJ2932 Coneology 题意: 给出一些不相交的圆,问有多少个圆不被其他圆包围 题解: 扫描线,把所有圆的左边界和右边界放到\(vector\)里排序,遍历到圆左边界的时候判断是否满足条件 ...
- poj2932 Coneology (扫描线)
转载请注明出处: http://www.cnblogs.com/fraud/ ——by fraud Coneology Time Limit: 5000MS Memory Lim ...
- 计算几何值平面扫面poj2932 Coneology
Coneology Time Limit: 5000MS Memory Limit: 65536K Total Submissions: 4097 Accepted: 859 Descript ...
- poj2932 Coneology
地址:http://poj.org/problem?id=2932 题目: Coneology Time Limit: 5000MS Memory Limit: 65536K Total Subm ...
- 刷题总结——coneology(poj2932 扫描线)
题目: Description A student named Round Square loved to play with cones. He would arrange cones with d ...
- poj 2932 Coneology(扫描线+set)
Coneology Time Limit: 5000MS Memory Limit: 65536K Total Submissions: 3574 Accepted: 680 Descript ...
- poj 2932 Coneology (扫描线)
题意 平面上有N个两两不相交的圆,求全部最外层的,即不被其它圆包括的圆的个数并输出 思路 挑战程序竞赛P259页 代码 /* ************************************* ...
- POJ 2932 Coneology(扫描线)
[题目链接] http://poj.org/problem?id=2932 [题目大意] 给出N个两两没有公共点的圆,求所有不包含于其它圆内部的圆 [题解] 我们计算出所有点在圆心所有y位置的x值, ...
- Coneology(POJ 2932)
原题如下: Coneology Time Limit: 5000MS Memory Limit: 65536K Total Submissions: 4937 Accepted: 1086 D ...
随机推荐
- windows下redis服务安装
1.redis简介redis是一个key-value存储系统.和Memcached类似,它支持存储的value类型相对更多,包括string(字符串).list(链表).set(集合).zset(so ...
- .net设计模式之观察者模式
摘要 在今天的设计模式系列中我给大家带来了观察者模式,首先我会以一个生动的故事引入观察者模式的应用的场景,然后描述这个场景中出现的问题,最后我们提出观察者模式的解决方案,并给出C#语言实现的代 ...
- Android知识思维导图
注:图片来源于网络,谢谢分享. 一. 项目目录结构: 布局控件 ListVIew控件 Widget:(窗口小部件) Activity Manager 二. 应用程序的5个模块构成: Activit ...
- 使用DNSSCrypt解决DNS污染问题
本文转自 月光博客,如有需要,请阅读原文. google近期在国内是不能访问了,dropbox这货居然也被DNS污染了,幸好发现DNSCrypt这一神器,防止DNS污染的绝佳工具. 基本原理:DNSC ...
- iOS 非ARC基本内存管理系列 3-循环retain和@class
1.@class 使用场景:对于循环依赖关系来说,比方A类引用B类,同时B类也引用A类: 可以看出Person和Card互相引用,此时如果使用#import编译报错!因此当使用@class在两个类中相 ...
- mysql 运行sql脚本文件
#只运行,不导出 mysql> source /home/user/to_run.sql; #导出 $ mysql -h host -u user -ppassword database < ...
- apache重写
---- 本文旨在提供如何用Apache重写规则来解决一些常见的URL重写方法的问题,通过常见的实例给用户一些使用重写规则的基本方法和线索. 一.为什么需要用重写规则 ---- 网站的生命在于不断地进 ...
- 在Spring中使用cache(EhCache的对象缓存和页面缓存)
Spring框架从version3.1开始支持cache,并在version4.1版本中对cache功能进行了增强. spring cache 的关键原理就是 spring AOP,通过 spring ...
- KnockoutJS(2)-监控属性
本节主要涉及到3个内容: 1. 监控属性 Observables (这个用的比较多,但是其本身使用比较简单,唯一需要注意就是,它监控的对象是一个方法,所以取值和设置值的时候容易混淆) 2. 计算属性 ...
- java 中的异常处理
一. 异常的概念和Java异常体系结构 异常是程序运行过程中出现的错误.本文主要讲授的是Java语言的异常处理.Java语言的异常处理框架, 是Java语言健壮性的一个重要体现. Java把 ...