Graham算法模板
1 #include <iostream>
2 #include <cstdio>
3 #include <cstring>
4 #include <cmath>
5 #include <algorithm>
6 using namespace std;
7 typedef long long ll;
8 const int maxn=1100;
9 const double pi=acos(-1.0);
10 struct node{
11 double x,y;
12 }p[maxn],a[maxn];
13 int n,tot;
14
15 //算距离
16 double dis(node a, node b){
17 return hypot(a.x-b.x,a.y-b.y);
18 }
19
20 //return 正:p2在向量p0p1的左侧;return 负:p2在向量p0p1的右侧;return 0:三点共线
21 double multi(node p0,node p1,node p2){
22 return (p1.x-p0.x)*(p2.y-p0.y)-(p2.x-p0.x)*(p1.y-p0.y);
23 }
24 //极角排序：根据坐标系内每一个点与x轴所成的角，逆时针比较，按照角度从小到大排序
25 int cmp(node a, node b){
26 int x=multi(a,b,p[0]);
27 if( x>0 || (x==0&& dis(a,p[0])<dis(b,p[0]))) return 1;
28 return 0;
29 }
30
31 void Graham(){
32 int k=0;
33 for(int i=1;i<n;i++){
34 if( p[i].y<p[k].y || (p[i].y==p[k].y&&p[i].x<p[k].x)) k=i;
35 }
36 swap(p[0],p[k]);
37 sort(p+1,p+n,cmp);
38 tot=2;
39 a[0]=p[0];
40 a[1]=p[1];
41 for(int i=2;i<n;i++){
42 while( tot>1 && multi(p[i],a[tot-1],a[tot-2])>=0 ){//a[tot-2]在p[i]a[tot-1]的左侧或三点共线,则p[i]在a[tot-2]a[tot-1]的右侧,则a[tot-1]不是极点
43 tot--;
44 }
45 a[tot++]=p[i];
46 }
47 }
48
49 int main()
50 {
51 int t;
52 scanf("%d",&t);
53 while( t-- ){
54 double r;
55 scanf("%d%lf",&n,&r);
56 for(int i=0;i<n;i++) scanf("%lf%lf",&p[i].x,&p[i].y);
57 Graham();
58 }
59 return 0;
60 }