几何 Crazy Rabbit

问题 H: Crazy Rabbit
时间限制: 2 Sec 内存限制: 512 MB
提交: 37 解决: 12
[提交][状态][讨论版]
题目描述
兔子们决定在自己的城堡里安排一些士兵进行防守。给出 n 个点的坐标，和城堡里一个圆心在原点的圆形的障碍
，兔子们希望从中选出 k 个兔子，使得它们两两所在的直线都不与圆相交。兔子们希望知道最多能选出多少兔子
输入
第一行两个整数 N 和 R, 表示兔子的个数和圆的半径接下来 N 行，每行两个整数 xi 和 yi ，表示第 i 只兔子
的坐标保证每只兔子都严格在障碍外部，且两两的所在的直线不与圆相切。
对于 100% 的测试数据, 1 <= n <= 2000; 1 <= R, xi, yi <= 5000
输出
输出一行一个整数, 表示最多能选出多少兔子
样例输入
6 3
0 6
-7 -4
-3 -2
7 -5
-2 3
8 -3
样例输出
4
【样例1解释】
选择第 1, 2, 6, 4 只兔子即可。

我开始想到了二分个数。。。貌似check只能2^n…
实际上从一个兔子作圆的两条切线，作为兔子的视线范围，那么就有几个性质。
1.如果两个兔子的视线范围相交，那么这两个兔子的连线就不会接触到圆。
2.如果两个兔子视线范围互相覆盖，那么连线所在的直线过圆。
3.如果不相交，则连线过圆。
这样问题就变成了在圆上选一些区间，让他们互相相交而不包含。
但这样还是很麻烦QAQ，其实还有一个性质，和线段覆盖的是优弧还是劣弧无关。现在把连线过圆的两种情况视为一种，如果当前被全覆盖，那么换成劣弧就成了不相交。怎么换都是连线过圆的情况。同理，如果当前两线段相交而不互相覆盖，说明线段2覆盖了线段1无法覆盖的一段地方，而线段1改成他的补集后，就覆盖了刚才他无法覆盖的一段。
所以现在问题转化成线段上找区间使其相交而不包含。
可以发现所有区间的左端点都在最靠左区间内部，而区间右端点都在之外。所以我们只要给区间按左端点排序，逐个枚举区间，选出满足这个的区间，然后让区间右端点跑最长上升子序列，更新答案即可。
其实把兔子转化成区间的过程用到了计几。
%%%%%达哥。

#include <cstdio>
#include <cstdlib>
#include <cstring>
#include <iostream>
#include <algorithm>
#include <cmath>
#define N 2005
#define fir first
#define sec second
#define d double
using namespace std;
int read()
{
    int sum=0,f=1;char x=getchar();
    while(x<'0'||x>'9'){if(x=='-')f=-1;x=getchar();}
    while(x>='0'&&x<='9'){sum=(sum<<1)+(sum<<3)+x-'0';x=getchar();}
    return sum*f;
}
int n,ans,r;
pair<d,d> a[N]; 
d b[N],c[N],P=acos(-1.0);
int check(int x)
{
    if(x==0)return 0;
    int len=0;c[++len]=b[1];
    for(int i=2;i<=x;i++)
    {
        if(b[i]>c[len]){c[++len]=b[i];continue;}
        int l=upper_bound(c+1,c+len+1,b[i])-c;
        c[l]=b[i];
    }
    return len;
}
int main()
{
    n=read();r=read();
    for(int i=1,x,y;i<=n;i++)
    {
        x=read();y=read();
        d k=sqrt(x*x+y*y),jj=atan2(y,x),Cos=acos(r/k);
        d l=jj-Cos,r=jj+Cos;
        while(l<=-P)l+=2.0*P;
        while(r>P)r-=2.0*P;
        if(l>r)swap(l,r);
        a[i]=make_pair(l,r);
    }
    sort(a+1,a+n+1);
    for(int i=1,len;i<=n;i++)
    {
        len=0;
        for(int j=i+1;j<=n;j++)
            if(a[j].fir<=a[i].sec&&a[j].sec>a[i].sec)
                b[++len]=a[j].sec;
        ans=max(ans,check(len)+1);
    }
    printf("%d\n",ans);
}