欢迎您访问程序员文章站本站旨在为大家提供分享程序员计算机编程知识!
您现在的位置是: 首页

Cows POJ - 3348 (平面凸包graham算法+多边形面积)

程序员文章站 2024-03-16 19:37:46
...

传送门

题解:使用平面凸包graham算法+多边形面积

平面凸包——>传送门

多边形面积——>传送门

附上代码:


#include<iostream>
#include<cstdio>
#include<algorithm>

using namespace std;

const int maxn=10010;

struct point{
    int x,y;
};
point A[maxn],result[maxn];

int dist(point a,point b)
{
    return (a.x-b.x)*(a.x-b.x)+(a.x-b.x)*(a.x-b.x);
}

int cp(point p1,point p2,point p3)
{
    return (p3.x-p1.x)*(p2.y-p1.y)-(p3.y-p1.y)*(p2.x-p1.x);
}

bool cmp(point a,point b)
{
    int ans=cp(A[0],a,b);
    if(ans==0){
        return dist(A[0],a)-dist(A[0],b)<=0;
    }else{
        return ans>0;
    }
}

int main()
{
    int n;
    while(scanf("%d",&n)!=EOF){
        int pos=0;
        for(int i=0;i<n;i++){
            scanf("%d%d",&A[i].x,&A[i].y);
            if(A[pos].y>=A[i].y){
                if(A[pos].y==A[i].y){
                    if(A[pos].x>A[i].x){
                        pos=i;
                    }
                }
            }else{
                pos=i;
            }
        }
        if(n<3){
            printf("0\n");
            continue;
        }
        swap(A[pos],A[0]);
        sort(A+1,A+n,cmp);
        int top=1;
        result[0]=A[0],result[1]=A[1];
        for(int i=2;i<n;i++){
            while(cp(result[top-1],result[top],A[i])<0){
                top--;
            }
            result[++top]=A[i];
        }
        double s=0;
        for(int i=0;i<=top;i++){
            double area=(result[(i+1)%(top+1)].x*result[i].y-result[(i+1)%(top+1)].y*result[i].x)/2.0;
			s+=area;
        }
        printf("%d\n",(int)(s/50.0));
    }
    return 0;
}

 

相关标签: computer geometry