Uva11722 joining with friend 概率 几何计算

/*
    题目描述：已知甲在s1到s2时间内各时刻到达车站的概率相等，乙在t1到t2时间内各时刻到达车站的概率相等，
                    火车将在车站停留w时间，问在火车停靠时间内甲乙见面的概率是多少
                    
    方法：典型的几何概型，高中时必修三最常考的，画出x= s1 , x = s2 , y = t1 ,y = t2所围成的区域，再做出直线
              y = x + w ,  直线 y = x - w，矩形区域夹在两线间的面积与矩形面积比即为答案
              直线所夹面积求法：cal（w）函数的含义是求出矩形在y = x + w 下方部分的面积，那么所夹面积可表示为
              cal（w）- cal(-w)
*/
#include<cstdio>
#include<cstring>
#include<algorithm>
#define mem(a,x) memset(a,x,sizeof(a))
using namespace std;
int s1,s2,t1,t2,w;
double sum;
double cal(int w)               //作用：计算y = x + w 直线下方图形的面积
{
    //以下四个变量为直线是否与矩形的四条边有交点
    bool left , right ,up , down;
    //以下四个变量分别为过矩形四个顶点的斜率为1直线的截距
    int lu = t2 - s1 , ld = t1 - s1 , ru = t2 - s2 , rd = t1 - s2 ;
    double ret;
    left = (w >= ld && w <= lu);
    right = (w >= rd && w <= ru);
    up = (w <= lu && w >= ru);
    down = (w <= ld && w >= rd);

    if(left && up)       ret = sum - 0.5*(t2 - w - s1)*(t2 - w - s1);
    else if(left && right)      ret = 0.5 * ( 2* w + s1 + s2 - 2 * t1)*(s2 - s1);
    else if(up && down)      ret = 0.5 * (2 * w + 2*s2 - t1 - t2) * (t2 - t1);
    else if(down && right)      ret = 0.5 * (s2 - t1 + w)*(s2 - t1 + w);
    if(w>=lu)      ret = sum;
    if(w<=rd)      ret = 0;
    return ret;
}
int main()
{
//    freopen("D:\\professional debugging\\in9.txt","r",stdin);
//    freopen("D:\\professional debugging\\out92.txt","w",stdout);
    int T , kase = 0;
    scanf("%d",&T);
    while(T--){
        scanf("%d %d %d %d %d",&s1,&s2,&t1,&t2,&w);
        sum = 1.0*(s2 - s1) * (t2 - t1);
        double area = cal(w) - cal(-w);
        double ans = area / sum;
        printf("Case #%d: %.8lf\n",++kase,ans);
    }
    return 0;
}