POJ 2187 凸包旋转卡壳

题意：

求平面最远点对。输出最远距离的平方。

参考：

http://www.cnblogs.com/Booble/archive/2011/04/03/2004865.html

http://www.cppblog.com/staryjy/archive/2009/11/19/101412.html

View Code

 1 #include <iostream>
 2 #include <cstdlib>
 3 #include <cstdio>
 4 #include <cstring>
 5 #include <algorithm>
 6 #include <cmath>
 7 
 8 #define N 50050
 9 
10 using namespace std;
11 
12 struct PO
13 {
14     int x,y;
15 }p[N];
16 
17 int stk[N],top,n;
18 
19 inline bool cmp(const PO &a,const PO &b)
20 {
21     if(a.x==b.x) return a.y<b.y;
22     else return a.x<b.x;
23 }
24 
25 inline int cross(const PO &o,const PO &a,const PO &b)
26 {
27     return (a.x-o.x)*(b.y-o.y)-(b.x-o.x)*(a.y-o.y);
28 }
29 
30 inline int get_dis2(const PO &a,const PO &b)
31 {
32     return (a.x-b.x)*(a.x-b.x)+(a.y-b.y)*(a.y-b.y);
33 }
34 
35 inline void read()
36 {
37     for(int i=1;i<=n;i++) scanf("%d%d",&p[i].x,&p[i].y);
38     sort(p+1,p+1+n,cmp);
39 }
40 
41 inline void graham()
42 {
43     sort(p+1,p+1+n,cmp);
44     top=0;
45     stk[++top]=1; stk[++top]=2;
46     for(int i=3;i<=n;i++)
47     {
48         while(top>=2&&cross(p[stk[top-1]],p[stk[top]],p[i])<=0) top--; 
49         stk[++top]=i;
50     }
51     int num=top;
52     for(int i=n-1;i>=1;i--)
53     {
54         while(top>num&&cross(p[stk[top-1]],p[stk[top]],p[i])<=0) top--;
55         stk[++top]=i;
56     }
57 }
58 
59 inline int rotating_calipers()
60 {
61     int ans=0,q=2;
62     for(int i=1;i<top;i++)
63     {
64         while(cross(p[stk[i+1]],p[stk[q+1]],p[stk[i]])>cross(p[stk[i+1]],p[stk[q]],p[stk[i]]))
65         {
66             q=(q+1)%top;
67             if(q==0) q++;
68         }
69         ans=max(ans,max(get_dis2(p[stk[i]],p[stk[q]]),get_dis2(p[stk[i+1]],p[stk[q+1]])));
70     }
71     return ans;
72 }
73 
74 inline void go()
75 {
76     graham();
77     printf("%d\n",rotating_calipers());
78 } 
79 
80 int main()
81 {
82     while(scanf("%d",&n)!=EOF) read(),go();
83     return 0;
84 }

表示和dyf神牛讲的一点都不一样。。

目测用旋转法解会恶心死。。