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;
}