【解题报告】PKU 2826 An Easy Problem?!

原题链接：http://poj.org/problem?id=2826

一题很蛋疼的一题。目前为止还有一个问题我没搞清楚，问题注在代码中。

题目大意：

外面下雨了，农民Johnoson的bull(??是什么)Ben要接点雨水去浇花。Ben拿两个木板在墙上，如图所示，墙上的两块木板可以看作是空中的两条线段，它们有一定的宽度。

你的任务是计算两块木板能接多少水。

最后要计算两块木板能装水的部分的面积。但是要注意一种情况，雨是垂直往下滴的，如果装水的部分没有雨能够进去的话，也是接不到水的。如图：

#include<stdio.h>
#include<math.h>
#define MAX(x,y) ((x)>(y)?(x):(y))
#define MIN(x,y) ((x)<(y)?(x):(y))
#define ABS(x) (((x)>0)?(x):(-(x)))
#define EPS 0.0000001/*精度控制*/
/*坐标的定义*/
typedef double coo;/*int*/
/*判相等*/
int is_equel(double a,double b)
{
double c=ABS(a-b);
if(c<=EPS) return 1;/*相等*/
else return 0;/*不相等*/
}
/*点、向量*/
typedef struct POINT
{
coo x,y;
}point,vector;
/*线段*/
typedef struct SEGMENT
{
point p1,p2;/*p[2];*/
}segment,line;
typedef struct LINE2/*一般式*/
{
coo A,B,C;
}line2;
/*向量的减法p1-p2*/
vector vector_minus(vector p1,vector p2)
{
vector p;
p.x=p1.x-p2.x;
p.y=p1.y-p2.y;
return p;
}
/*向量叉乘*/
double cross_product(vector p1,vector p2)
{/*x1y2-x2y1*/
return p1.x*p2.y-p1.y*p2.x;
}
int line_location(line l1,line l2)
{
if((l2.p2.y-l2.p1.y)*(l1.p2.x-l1.p1.x)==(l1.p2.y-l1.p1.y)*(l2.p2.x-l2.p1.x))
{
if((l2.p2.y-l1.p1.y)*(l1.p2.x-l1.p1.x)==(l1.p2.y-l1.p1.y)*(l2.p2.x-l1.p1.x))
return -1;/*重合*/
else return 1;/*平行*/
}
else return 0;/*相交*/
}

/*判断线段是否相交*/
int segment_intersect(segment l1,segment l2)
{
/*快速排斥试验*/
coo min,max;
min=MIN(l1.p1.x,l1.p2.x);
max=MAX(l2.p1.x,l2.p2.x);
if(min>max) return 0;
min=MIN(l2.p1.x,l2.p2.x);
max=MAX(l1.p1.x,l1.p2.x);
if(min>max) return 0;
min=MIN(l1.p1.y,l1.p2.y);
max=MAX(l2.p1.y,l2.p2.y);
if(min>max) return 0;
min=MIN(l2.p1.y,l2.p2.y);
max=MAX(l1.p1.y,l1.p2.y);
if(min>max) return 0;
/*跨立试验*/
if(cross_product(vector_minus(l1.p1,l2.p1),vector_minus(l2.p2,l2.p1))*cross_product(vector_minus(l1.p2,l2.p1),vector_minus(l2.p2,l2.p1))>0) return 0;
if(cross_product(vector_minus(l2.p1,l1.p1),vector_minus(l1.p2,l1.p1))*cross_product(vector_minus(l2.p2,l1.p1),vector_minus(l1.p2,l1.p1))>0) return 0;
return 1;
}
line2 line_line2(line l)
{
line2 l2;
l2.A=l.p2.y-l.p1.y;
l2.B=l.p1.x-l.p2.x;
l2.C=l.p2.x*l.p1.y-l.p1.x*l.p2.y;
return l2;
}
/*直线相交的交点（两点式）*//*未处理平行和重合的情况*/
point line_intersection(line a,line b)
{
line2 l1,l2;
point i={0,0};
l1=line_line2(a);
l2=line_line2(b);
i.x= (l2.B*l1.C-l1.B*l2.C)/(l2.A*l1.B-l1.A*l2.B);
i.y=-(l2.A*l1.C-l1.A*l2.C)/(l2.A*l1.B-l1.A*l2.B);
return i;
}
void scan(line *l)
{
scanf("%lf%lf%lf%lf",&l->p1.x,&l->p1.y,&l->p2.x,&l->p2.y);
}
int main()
{
int t,i,j;
segment l1,l2;
point P,p2;
/*freopen("data.txt","r",stdin);
freopen("ouput2.txt","w",stdout);*/
scanf("%d",&t);
while(t--)
{
scan(&l1);
scan(&l2);
if(segment_intersect(l1,l2)==0) printf("0.00\n");
else if(l1.p1.y==l1.p2.y||l2.p1.y==l2.p2.y) printf("0.00\n");
else if(line_location(l1,l2)!=0) printf("0.00\n");
else
{
point p[4];
P=line_intersection(l1,l2);
p[0]=l1.p1;
p[1]=l1.p2;
p[2]=l2.p1;
p[3]=l2.p2;
for(i=1;i<4;i++)
for(j=0;j<i;j++)
{
if(p[i].y>p[j].y) {p2=p[i];p[i]=p[j];p[j]=p2;}
}
{
segment s1,s2,s0;/*| -*/
double S;
s0.p1=p[0];s0.p2=P;
s1.p1=p[1];s1.p2.x=p[1].x;s1.p2.y=p[1].y-1;
s2.p1=p[1];s2.p2.x=p[1].x-1;s2.p2.y=p[1].y;
p[2]=line_intersection(s1,s0);
if(p[2].y<=p[0].y&&p[2].y>p[1].y) printf("0.00\n");
else
{
p[3]=line_intersection(s2,s0);
S=cross_product(vector_minus(p[1],P),vector_minus(p[3],P))/2;
S=ABS(S);
printf("%.2lf\n",S+EPS);/*问题：为什么最后这里不加EPS就过不了呢？？？求大神解答。*/
}
}
}
}
return 0;
}

View Code
不知道为什么调试过程中最后一步输出S之前已经取了绝对值，可是有次还是输出了-0.00。是什么情况？？？