POJ 1106 极角排序

题意：

给定半圆的圆心和半径和n个点，求能在半圆内的最多的点数。

 

题解：

先把所有的在圆外的点删除，然后求出所有的点关于圆心的极角，然后排序，因为是环形的，极角复制一份加在原来的后面，维护双指针即可~

 

View Code

 1 #include <iostream>

 2 #include <cstdlib>

 3 #include <cstdio>

 4 #include <cstring>

 5 #include <algorithm>

 6 #include <cmath>

 7 

 8 #define PI 3.141592653589

 9 #define EPS 1e-7

10 #define N 600

11 

12 using namespace std;

13 

14 struct PO

15 {

16     double x,y,angle;

17 }p[N];

18 

19 int n,ans;

20 double sr;

21 

22 inline bool cmp(const PO &a,const PO &b)

23 {

24     return a.angle<b.angle;

25 }

26 

27 inline int doublecmp(double x)

28 {

29     if(x>EPS) return 1;

30     else if(x<-EPS) return -1;

31     return 0;

32 }

33 

34 inline double get_dis2(PO &a,PO &b)

35 {

36     return (a.x-b.x)*(a.x-b.x)+(a.y-b.y)*(a.y-b.y);

37 }

38 

39 inline void read()

40 {

41     n=0;

42     int tmp;scanf("%d",&tmp);

43     for(int i=1;i<=tmp;i++)

44     {

45         ++n;

46         scanf("%lf%lf",&p[n].x,&p[n].y);

47         if(doublecmp(get_dis2(p[n],p[0])-sr*sr)==1) n--;

48         else p[n].angle=atan2(p[n].y-p[0].y,p[n].x-p[0].x);

49     }

50 }

51 

52 inline void go()

53 {

54     ans=0;

55     sort(p+1,p+1+n,cmp);

56     for(int i=1;i<=n;i++) p[i+n].angle=p[i].angle+2*PI;

57     for(int g1=1,g2=1;g1<=n&&g2<=n+n;g2++)

58         if(doublecmp(p[g2].angle-p[g1].angle-PI)==1)

59         {

60             ans=max(g2-g1,ans);

61             g1++;

62         }

63     printf("%d\n",ans);

64 }

65 

66 int main()

67 {

68     while(scanf("%lf%lf%lf",&p[0].x,&p[0].y,&sr),doublecmp(sr)>=0) read(),go();

69     return 0;

70 }