bzoj 4206 最大团 几何+lis

最大团

Time Limit: 10 Sec  Memory Limit: 256 MB
Submit: 142  Solved: 65
[Submit][Status][Discuss]

Description

给出平面上N个点的坐标，和一个半径为R的圆心在原点的圆。对于两个点，它们之间有连边，当且仅当它们的连线与圆不相交。

求此图的最大团。

 

Input

第一行两个整数N和R, 表示点数和圆的半径。

接下来N 行，每行两个整数xi 和yi,表示第i个点的坐标

保证每个点都严格在园外，且两两直线不与圆相切。

 

Output

输出一个整数：最大团的大小。

 

Sample Input

6 3
0 6
-7 -4
-3 -2
7 -5
-2 3
8 -3

Sample Output

4

HINT

对于100%的数据，1≤N≤2000,|xi|,|yi|,R≤5000

 #include<cstring>
 #include<cstdio>
 #include<algorithm>
 #include<iostream>
 #include<cmath>
 #include<queue>
 #include<map>

 #define N 2007
 #define pi acos(-1)
 #define inf 1000000007
 using namespace std;
 inline int read()
 {
     int x=,f=;char ch=getchar();
     while(!isdigit(ch)){if(ch=='-')f=-;ch=getchar();}
     while(isdigit(ch)){x=(x<<)+(x<<)+ch-'';ch=getchar();}
     return x*f;
 }

 int n,tot,st[N];
 double r;
 struct Node
 {
     double l,r;
     Node()
     {
         l=r=;
     }
     friend bool operator<(Node x,Node y)
     {
         return x.l<y.l;
     }
 }a[N];

 int getlis(int x)
 {
     st[tot=]=x,a[].r=-inf;
     for (int i=x+;i<=n&&a[i].l<a[x].r;i++)
         if (a[i].r>a[st[tot]].r) st[++tot]=i;
         else
         {
             int l=,r=tot,ans=tot;
             while(l<=r)
             {
                 int mid=(l+r)>>;
                 if (a[st[mid]].r>=a[i].r) ans=min(mid,ans),r=mid-;
                 else l=mid+;
             }
             if (ans==) continue;
             if (a[i].r<a[st[ans]].r) st[ans]=i;
         }
     return tot;
 }
 double get_dis(double x,double y)
 {
     return sqrt(x*x+y*y);
 }
 int main()
 {
     n=read(),r=read();
     for (int i=;i<=n;i++)
     {
         double x=read(),y=read();
         if (get_dis(x,y)<r){i--,n--;continue;}
         a[i].l=atan2(y,x)-acos(r/get_dis(x,y));
         a[i].r=atan2(y,x)+acos(r/get_dis(x,y));
         if (a[i].r>pi) a[i].r-=*pi,swap(a[i].l,a[i].r);
         if (a[i].l<-pi) a[i].l+=*pi,swap(a[i].l,a[i].r);
     }
     sort(a+,a+n+);

     int ans=;
     for (int i=;i<=n;i++)
         ans=max(ans,getlis(i));
     printf("%d\n",ans);
 }