P3897 [湖南集训]Crazy Rabbit
2024-09-28 17:23:23
\(\color{#0066ff}{ 题目描述 }\)
兔子们决定在自己的城堡里安排一些士兵进行防守。
给出 n 个点的坐标,和城堡里一个圆心在原点的圆形的障碍,兔子们希望从中选出 k 个兔子,使得它们两两所在的直线都不与圆相交。
兔子们希望知道最多能选出多少兔子。
\(\color{#0066ff}{输入格式}\)
第一行两个整数 N 和 R, 表示兔子的个数和圆的半径接下来 N 行,每行两个整数 xi 和 yi ,表示第 i 只兔子的坐标
保证每只兔子都严格在障碍外部,且两两的所在的直线不与圆相切。
\(\color{#0066ff}{输出格式}\)
输出一行一个整数, 表示最多能选出多少兔子
\(\color{#0066ff}{输入样例}\)
6 3
0 6
-7 -4
-3 -2
7 -5
-2 3
8 -3
\(\color{#0066ff}{输出样例}\)
4
\(\color{#0066ff}{数据范围与提示}\)
选择第 1, 2, 6, 4 只兔子即可。
- 对于 10% 的测试数据, 1 ≤ n ≤ 20.
- 对于 30% 的测试数据, 1 ≤ n ≤ 100.
- 对于 100% 的测试数据, 1 ≤ n ≤ 2000; 1 ≤ R; xi; yi ≤ 5000.
\(\color{#0066ff}{ 题解 }\)
过每一个点做圆的两条切线,每个点占据了圆上的一段劣弧
不难发现,两个点的连线与圆有交,当且仅当两段劣弧不相离且互补包含,即普通的相交
于是我们开始想办法获得这个劣弧
我们该怎么存储这种东西呢
考虑把环断开,以极角作为链上的数值,强制让线段范围为\([-\pi,\pi]\),对于不在这个范围上的端点,可以直接\(+-2\pi\),显然没影响
那么怎么获得这两个端点值(角度)呢
∠AOB=∠AOC, 因此只要让极角∠AOD加减∠AOB即可获得端点值
然后我们把所有区间按左端点排序
然后枚举每个区间作为选择的起始区间,枚举之后选择的区间,那之后选择的区间要满足左端点小于起始区间右端点,右端点大于起始区间右端点,这样才是相交
我们把这些区间放在数组里,左端点是递增的且一定小于起始区间右端点
那么可以发现,最多选择的数量就是这些区间右端点的最长上升子序列!
于是\(O(n^2logn)\)即可解决本题
#include<bits/stdc++.h>
#define LL long long
LL in() {
char ch; LL x = 0, f = 1;
while(!isdigit(ch = getchar()))(ch == '-') && (f = -f);
for(x = ch ^ 48; isdigit(ch = getchar()); x = (x << 1) + (x << 3) + (ch ^ 48));
return x * f;
}
const int maxn = 2050;
const double pi = acos(-1);
using std::pair;
using std::make_pair;
int n;
double R;
pair<double, double> e[maxn];
int num, ans;
double f[maxn], t[maxn];
int query() {
int len;
if(!num) return 0;
f[len = 1] = t[1];
for(int i = 2; i <= num; i++) {
if(t[i] > f[len]) f[++len] = t[i];
else *std::upper_bound(f + 1, f + len + 1, t[i]) = t[i];
}
return len;
}
int main() {
n = in(), R = in();
double x, y;
for(int i = 1; i <= n; i++) {
x = in(), y = in();
double bt = atan2(y, x);
double len = sqrt(x * x + y * y);
double af = acos(R / len);
x = bt - af, y = bt + af;
while(x <= -pi) x += 2.0 * pi;
while(y > pi) y -= 2.0 * pi;
if(x > y) std::swap(x, y);
e[i] = make_pair(x, y);
}
std::sort(e + 1, e + n + 1);
for(int i = 1; i <= n; i++) {
num = 0;
for(int j = i + 1; j <= n; j++)
if(e[j].first <= e[i].second && e[j].second >= e[i].second)
t[++num] = e[j].second;
ans = std::max(ans, query());
}
printf("%d", ans + 1);
return 0;
}
最新文章
- [Java] 特殊正则-替换字符串
- JAVA 设计模式 享元模式
- eclipse/intellij Idea集成jetty
- Spring学习 Ioc篇(二 )
- 标准盒模型与ie盒模型
- Linux之kernal分析与启动20160610
- Red5 1.0.5安装过程记录
- 为什么子线程不能做UI操作
- FileUtils类应该有哪些内容
- 如果有需要确解MD5的,可以尝试这个网站。
- openNebula rgister img instance vms error collections
- Kmeans算法与KNN算法的区别
- Dynamics CRM 2015-Auto Save
- C# Excel写入数据及图表
- 7-EL表达式和JSTL表达式
- 7.ASP.NET MVC 5.0中的Routing【路由】
- python安装pandas和lxml
- http链接的性能测试工具httping
- P/Invoke 光标的操作
- QtAV的编译方法