BNU 29379斩

在电影里面，我们经常可以看到武士们拔出刀，然后一刀斩下去，结果………………………………一张纸片被砍成了两半，囧…………

而在本题中，我们需要计算一下被斩下去较小的那一部分的面积。

我们假设纸片是矩形的，平行于坐标轴的，武士砍纸片的轨迹是一条直线。

Input

第一行一个整数N(2<=N<=350)，表示数据组数。

接下来一行，每行7个整数，x_l,y_l,x_r,y_r,a,b,c ，表分别表示矩形左下角坐标(x_l,y_l)，右上角坐标(x_r,y_r)，以及轨迹方程ax+by+c=0。整数的绝对值均小于200。

Output

对于每一组数据，输出一个三位小数，表示砍后较小那部分的面积。数据保证这个面积大于0.001。

Sample Input

3
1 1 3 3 -1 1 0
1 1 3 4 -1 1 0
1 2 3 4 -1 1 0

Sample Output

2.000
2.000
0.500

这是计算几何的题目，题目保证有切割的，故分情况讨论：

1、直线与矩形边平行时。 此时a=0 或者b=0；

2、不平行时刻找出四个点x_1,y_1,x_2,y_2，分别是直线与矩形下，上，左，右四条边的交点的横和纵坐标。

      然后讨论直线斜率看点的位置。我是以y_1，y_2为准。

#include <iostream>
#include <cstdio>
#include <cstring>
#include <algorithm>
#include <cmath>

using namespace std;

int main(){
    int cas;
    scanf("%d",&cas);
    int xl,yl,xr,yr,a,b,c;

    while(cas--){
        scanf("%d%d%d%d%d%d%d",&xl,&yl,&xr,&yr,&a,&b,&c);
        if(a==0){
            double y=-1*(double)c/b;
            if(y-yl > yr-y)
                printf("%.3lf\n",(double)(yr-y)*(xr-xl));
            else printf("%.3lf\n",(double)(y-yl)*(xr-xl));
            continue;
        }
        if(b==0){
            double x=-1*(double)c/a;
            if(x-xl > xr-x) printf("%.3lf\n",(double)(yr-yl)*(xr-x));
            else printf("%.3lf\n",(double)(yr-yl)*(x-xl));
            continue;
        }
        double x_1=(double)(-1*c-b*yl)/a;
        double x_2=(double)(-1*c-b*yr)/a;
        double y_1=(double)(-1*c-a*xl)/b;
        double y_2=(double)(-1*c-a*xr)/b;
        double k=(double)(-1*a)/b;
        double S=(double)(yr-yl)*(xr-xl);
        if(k>0){
            if(y_2>=yr && y_1<=yl){
                double s=(double)(xr-x_2+xr-x_1)*(yr-yl)/2;
                printf("%.3lf\n",min(s,S-s));
            }
            else if((y_1>=yl && y_1<=yr) && y_2>=yr)
                printf("%.3lf\n",(double)(x_2-xl)*(yr-y_1)/2);
            else if((y_2<=yr && y_2>=yl) && y_1<=yl)
            printf("%.3lf\n",(double)(xr-x_1)*(y_2-yl)/2);
            else if((y_1>=yl && y_1<=yr) && (y_2<=yr && y_2>=yl)){
                double s=(double)(yr-y_2+yr-y_1)*(xr-xl)/2;
                printf("%.3lf\n",min(s,S-s));
            }
            continue;
        }
        else{
            if(y_1>=yr && y_2<=yl){
                double s=(double)(x_2-xl+x_1-xl)*(yr-yl)/2;
                printf("%.3lf\n",min(s,S-s));
            }
            else if((y_1<=yr && y_1>=yl) && y_2<=yl)
                printf("%.3lf\n",(double)(y_1-yl)*(x_1-xl)/2);
            else if((y_2>=yl && y_2<=yr) && y_1>=yr)
                printf("%.3lf\n",(double)(xr-x_2)*(yr-y_2)/2);
            else{
                double s=(double)(yr-y_1+yr-y_2)*(xr-xl)/2;
                printf("%.3lf\n",min(s,S-s));
            }
            continue;
        }
    }
    return 0;
}