graham模板整理

来源:互联网 发布:大数据应用的典型案例 编辑:程序博客网 时间:2024/06/05 18:03

感觉以前写的graham代码风格太差了,然后又重新整理了一份


struct Node{    double x,y;    int pos;    double val;}a[maxn],s;double across(Node x,Node y,Node z)//z在x,y连线的顺时针为负{    return (x.x - y.x) * (z.y - y.y) - (x.y - y.y) * (z.x - y.x);}bool cmp(Node x,Node y){    if(x.val == y.val)return d2(x,s) < d2(y,s);    return x.val > y.val;}void Init(){for(int i = 1;i <= n;i ++){     scanf("%lf%lf",&a[i].x,&a[i].y);a[i].pos = i;     if(a[i].y < s.y || (a[i].y == s.y && a[i].x < s.x))s = a[i];}for(int i = 1;i <= n;i ++)if(a[i].x == s.x && a[i].y == s.y)a[i].val = 1;else a[i].val = get(a[i],s);sort(a + 1,a + n + 1,cmp);}void graham()  {      Node st[maxn];      int len = 0;      st[++ len] = a[1];      if(n>= 2)      st[++ len] = a[2];      if(n >= 3)      st[++ len] = a[3];      for(int i = 4; i <= n; i ++)      {          while(across(st[len],st[len],a[i]) <= 0)//< 0是可以允许共线,<= 0是不允许          {              len --;          }          st[++ len] = a[i];      }  //周长    double ans = 0;    for(int i = 1; i <= len - 1; i ++)      {          ans += distances(st[i],st[i + 1]);      }      ans += distances(st[1],st[len]) ;//面积:double anss = 0;      for(int i = 1; i <= len - 1; i ++)      {          anss += (st[i].x * st[i + 1].y - st[i + 1].x * st[i].y) / 2;      }    }  


原创粉丝点击