POJ 2826 An Easy Problem?! <计算几何>

题目

题目大意：墙上有两个等宽的木板（宽度为一个单位），问它们能够接到多少雨水。

分析：就是求能够接到雨水的区域的截面面积，本来这道题没什么思考上的难度的，但是情况有很多种，一开始我就忽略掉了几种能够接到雨水的情况。。代码写了详细的注释，我感觉自己的代码写得还是挺清晰简洁的-.-

代码：

#include <iostream>#include <cmath>#include <algorithm>#include <iomanip>using namespace std;const double EPS=1e-8;struct Point;typedef Point Vec;int dblcmp(double x){    return fabs(x)<EPS?0:(x>0?1:-1);}struct Point{    double x,y;    Point(){}    Point(double xx,double yy):x(xx),y(yy){}    Vec operator -(Point p){    //两点相减得向量        return Vec(x-p.x,y-p.y);    }    double operator ^(Vec v){   //叉积        return x*v.y-y*v.x;    }};struct Segment{    Point p1,p2;    Segment(){}    Segment(Point pp1,Point pp2):p1(pp1),p2(pp2){}    bool isParallel(Segment s){ //是否共线        return dblcmp((p2-p1)^(s.p1-p1))==0&&                dblcmp((p2-p1)^(s.p2-p1))==0;    }    bool isCross(Segment s,double& s1,double &s2){  //线段相交判断，允许端点相交，但不允许延长线相交        return max(p1.x,p2.x)>=min(s.p1.x,s.p2.x)&&                max(p1.y,p2.y)>=min(s.p1.y,s.p2.y)&&                max(s.p1.x,s.p2.x)>=min(p1.x,p2.x)&&                max(s.p1.y,s.p2.y)>=min(p1.y,p2.y)&&                dblcmp(s1=((p2-p1)^(s.p1-p1)))*dblcmp(s2=((p2-p1)^(s.p2-p1)))<=0&&                dblcmp((s.p2-s.p1)^(p1-s.p1))*dblcmp((s.p2-s.p1)^(p2-s.p1))<=0;    }    Point intersectionPoint(double s1,double s2){//求交点，定点分比法        double x=(p1.x*s2-s1*p2.x)/(s2-s1);        double y=(p1.y*s2-s1*p2.y)/(s2-s1);        return Point(x,y);    }};Segment s[2];void solve(){    double s1,s2;    Point p[2];    int c=0;    if(s[0].isParallel(s[1])) cout<<"0.00"<<endl;    else if(s[0].isCross(s[1],s1,s2)){        Point o=s[1].intersectionPoint(s1,s2);        for(int i=0;i<2;++i){   //找出在交点上方的点            if(s[i].p1.y>o.y) p[c++]=s[i].p1;            if(s[i].p2.y>o.y) p[c++]=s[i].p2;        }        if(c!=2) cout<<"0.00"<<endl; //交点上方不是两个点，一定不能接到雨水        else{            Vec v1=p[0]-o,v2=p[1]-o;            int top=0;            if(p[1].y>p[0].y) top=1;            Point pl=p[top^1],pt=p[top]; //pt是较高的点            int flag1=dblcmp(fabs(pt.x-o.x)-fabs(pl.x-o.x));//水平方向上，较高点是否离交点更远            int flag2=dblcmp((pt-o)^(pl-o)); //较高点与交点形成的木板，与较低点与交点形成的木板之间的相对位置（逆时针还是顺时针）            int flag3=dblcmp(pt.x-o.x); //当两点位于交点同侧是，是在右侧还是左侧            int flag=dblcmp(flag1*flag2*flag3); //综合影响。画一画图就会发现，当两点位于交点同侧时，以上三个判断的结果相乘大于0才能接到雨水            if(dblcmp(v1.x*v2.x)==1&&flag<=0) //只有当两点位于交点同侧时，才需判断flag                cout<<"0.00"<<endl;            else{                double yg=pt.y,yl=pl.y;                double r=(yl-o.y)/(yg-o.y); //能接到雨水的区域所占的比例                cout<<fixed<<setprecision(2)<<fabs((v1^v2)*r/2.0)<<endl;            }        }    }    else cout<<"0.00"<<endl; //没有交点}int main(){    ios::sync_with_stdio(false);    int n;    cin>>n;    while(n--){        double x1,y1,x2,y2;        for(int i=0;i<2;++i){            cin>>x1>>y1>>x2>>y2;            s[i]=Segment(Point(x1,y1),Point(x2,y2));        }        solve();    }    return 0;}