zoj Area

WA了n久最终曲网上找了个程序，一块一块的替换。

 

1.count++ 放在了continue 的后面，以至于Impossible 的时候count没有+1

 

2.向量aXb=0 以为就在同一直线上，只判断了x是否在线段中间。其实其中一个向量可能为0

 

3.Figuer 1: 0.75 中间的空格没有加

 

代码如下:

#include <stdio.h> #define MAXLEN 1010 #define PRECISION 0.000001 typedef struct{ double x,y; } Point; Point p[MAXLEN]; int n; double area; int dblcmp(double x) { if(x<-PRECISION) return -1; if(x>PRECISION) return 1; return 0; } /* * 向量叉乘: p1p0 X p2p0 */ double det(Point p0,Point p1,Point p2) { return (p1.x-p0.x)*(p2.y-p0.y) - (p2.x-p0.x)*(p1.y-p0.y); } double Max(double a,double b) { return (a>b) ? a:b; } double Min(double a,double b) { return (a<b) ? a:b; } /* * 判断在同一直线上的三点a,b,c. * a是不是在线段<b,c>上 */ int on_segment(Point p0,Point p1,Point p2) { double xMax = Max(p1.x,p2.x); double xMin = Min(p1.x,p2.x); double yMax = Max(p1.y,p2.y); double yMin = Min(p1.y,p2.y); if(p0.x >= xMin && p0.x <= xMax && p0.y >= yMin && p0.y <= yMax) return 1; else return 0; } int judge_cross() { int i,j,start; double d[4]; if(n<3) return 1; for(i=3;i<=n;i++) { start= (i==n)? 2:1; /**判断 边<i-1,i> 和 边<j-1,j> 是否相交**/ for(j=start;j<i-1;j++) { /*判断边<i-1,i> 和 边<j-1,j> 是否正规相交*/ d[0] = det(p[j-1],p[i],p[i-1]); d[1] = det(p[j],p[i],p[i-1]); d[2] = det(p[i-1],p[j],p[j-1]); d[3] = det(p[i],p[j],p[j-1]); if( d[0]*d[1]<-PRECISION && d[2]*d[3]<-PRECISION) return 1; /*判断点j-1或j是否在线段<i-1,i>上*/ if( dblcmp(d[0])==0 && on_segment(p[j-1],p[i],p[i-1]) ) return 1; if( dblcmp(d[1])==0 && on_segment(p[j],p[i],p[i-1]) ) return 1; /*判断点i-1或i是否在线段<j-1,j>上*/ if( dblcmp(d[2])==0 && on_segment(p[i-1],p[j-1],p[j]) ) return 1; if( dblcmp(d[3])==0 && on_segment(p[i],p[j-1],p[j]) ) return 1; } /*判断边<i-1,i> 和 边<i-2,i-1> 是否有重合 即点i是否在边<i-2,i-1>上*/ d[0] = det(p[i],p[i-1],p[i-2]); if( dblcmp(d[0])==0 && on_segment(p[i],p[i-1],p[i-2]) ) return 1; } return 0; } void compute() { int i; area = 0.0; for(i=0;i<n;i++) { area += det(p[0],p[i-1],p[i]); } if(dblcmp(area)==-1) area = -area; area /= 2.0; } int main() { int i,count=1; while(scanf("%d",&n)&&n) { for(i=0;i<n;i++) scanf("%lf%lf",&p[i].x,&p[i].y); p[n] = p[0]; printf("Figure %d: ",count++); if(judge_cross() == 1) printf("Impossible/n"); else { compute(); printf("%.2lf/n", area); } } return 0; }