zoj 1010 Area

给你一些点（按顺序），先判断这些点是否组成一个多边形，题目已经很清楚了，如果不是一个多边形，那么它的一条边和其他非直接相邻的边相交了。如果是多边形，输出多边形面积。

 

n = 1 和2的时候必然不是多边形，需要特判。

 

这个求多边形是把一个多边形分解为多个三角形，还是利用叉积来计算的，手工模拟太痛苦了><。。。不弄了。

 

#include <queue> #include <stack> #include <math.h> #include <stdio.h> #include <stdlib.h> #include <iostream> #include <limits.h> #include <string.h> #include <algorithm> using namespace std; const int MAX = 1010; const double eps = 1e-6; struct point{ double x,y; }; point p[MAX]; bool dy(double x,double y) { return x > y + eps; } // x > y bool xy(double x,double y) { return x < y - eps; } // x < y bool dyd(double x,double y) { return x > y - eps; } // x >= y bool xyd(double x,double y) { return x < y + eps; } // x <= y bool dd(double x,double y) { return fabs( x - y ) < eps; } // x == y double crossProduct(point a,point b,point c)//向量 ac 在 ab 的方向 { return (c.x - a.x)*(b.y - a.y) - (b.x - a.x)*(c.y - a.y); } double area_polygon(point p[],int n) { double s1 = 0.0,s2 = 0.0; for(int i=0; i<n; i++) { s1 += p[(i+1)%n].y * p[i].x; s2 += p[(i+1)%n].y * p[(i+2)%n].x; } return fabs(s1 - s2)/2.0; } bool onSegment(point a, point b, point c) { double maxx = max(a.x,b.x); double maxy = max(a.y,b.y); double minx = min(a.x,b.x); double miny = min(a.y,b.y); if( dd(crossProduct(a,b,c),0.0) && dyd(c.x,minx) && xyd(c.x,maxx) && dyd(c.y,miny) && xyd(c.y,maxy) ) return true; return false; } bool segIntersect(point p1,point p2, point p3, point p4) { double d1 = crossProduct(p3,p4,p1); double d2 = crossProduct(p3,p4,p2); double d3 = crossProduct(p1,p2,p3); double d4 = crossProduct(p1,p2,p4); if( xy(d1 * d2,0.0) && xy( d3*d4,0.0 ) ) return true; if( dd(d1,0.0) && onSegment(p3,p4,p1) ) return true; if( dd(d2,0.0) && onSegment(p3,p4,p2) ) return true; if( dd(d3,0.0) && onSegment(p1,p2,p3) ) return true; if( dd(d4,0.0) && onSegment(p1,p2,p4) ) return true; return false; } bool is_polygon(point p[],int n) { p[n] = p[0]; for(int i=0; i<n; i++) for(int k=i+1; k<n; k++) if( i+1 != k && (k+1)%n != i && segIntersect(p[i],p[i+1],p[k],p[k+1]) ) return false; return true; } int main() { int n; int ind = 1; while( ~scanf("%d",&n) && n ) { if( ind != 1 ) printf("/n"); for(int i=0; i<n; i++) scanf("%lf %lf",&p[i].x,&p[i].y); printf("Figure %d: ",ind++); if( n <= 2 ) { printf("Impossible/n"); continue; } if( is_polygon(p,n) ) printf("%.2lf/n",area_polygon(p,n)); else printf("Impossible/n"); } return 0; }