BNU 29379斩

来源：互联网发布：geo数据库使用方法编辑：程序博客网时间：2024/04/29 17:16

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

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

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

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

31 1 3 3 -1 1 01 1 3 4 -1 1 01 2 3 4 -1 1 0

Sample Output

2.0002.0000.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;}

0 0