感觉以前写的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;
}
}