BNUOJ--29379 斩

最新推荐文章于 2019-08-22 17:23:52 发布

for_further

最新推荐文章于 2019-08-22 17:23:52 发布

阅读量721

点赞数

分类专栏： bnuoj 我要学数学！QAQ 文章标签：几何

本文链接：https://blog.csdn.net/fad_further/article/details/23048461

版权

bnuoj 同时被 2 个专栏收录

10 篇文章 0 订阅

订阅专栏

我要学数学！QAQ

9 篇文章 0 订阅

订阅专栏

一条直线分割一个矩形，求分割后面积小的那部分的面积。

http://www.bnuoj.com/bnuoj/problem_show.php?pid=29379

直线方程：ax + by + c = 0

首先特判直线平行于坐标轴的情况。。。也就是最简单的情况（虽然比赛的时候我就是折在这上面。。。QuQ

我们将剩下的情况分为红色线的情况和蓝色线的情况

红色线的情况就是直线交在矩形相邻的两条边上，即求红色区域的面积

蓝色线的情况是直线交在对边上，先求三角形的面积，在加上更小的矩形的面积

#include <cmath>
#include <cstdio>
#include <cstring>
#include <algorithm>
using namespace std;
const double  eps = 1e-5;
int t;
double xl ,yl, xr, yr, a, b, c;
bool A, B, C, D;
int main(){
    scanf("%d", &t);
    while(t--){
        scanf("%lf%lf%lf%lf%lf%lf%lf", &xl, &yl, &xr, &yr, &a, &b, &c);
        A = B = C = D = false;
        if(fabs(a) < eps){
            double y = (yl + yr) / 2.0;
            if(-c / b < y)
                printf("%.3f\n", (xr - xl) * (-c / b - yl));
            else
                printf("%.3f\n", (xr - xl) * (yr + c / b));
            continue;
        }
        if(fabs(b) < eps){
            double x = (xl +xr) / 2;
            if(-c / a < x)
                printf("%.3f\n", (yr - yl) * (-c / a - xl));
            else
                printf("%.3f\n", (yr - yl) * (xr + c / a));
            continue;
        }
        double xb = (-c / a) - (b / a) * yr;
        if((xb - xl) >= 0.0 && (xr - xb) >= 0.0) B = true;
        double xd = (-c / a) - (b / a) * yl;
        if((xd - xl) >= 0.0 && (xr - xd) >= 0.0) D = true;
        double ya = (-a / b) * xl - c / b;
        if((ya - yl) >= 0.0 && (yr - ya) >= 0.0) A = true;
        double yc = (-a / b) * xr - c / b;
        if((yc - yl) >= 0.0 && (yr - yc) >= 0.0) C = true;
        if(A && B && !C && !D){
            printf("%.3f\n", 0.5 * (yr - ya) * (xb - xl));
        }
        else if(A && D && !B && !C){
            printf("%.3f\n", 0.5 * (ya - yl) * (xd - xl));
        }
        else if(C && B && !A && !D){
            printf("%.3f\n", 0.5 * (yr - yc) * (xr - xb));
        }
        else if(C && D && !A && !B){
            printf("%.3f\n", 0.5 * (yc - yl) * (xr - xd));
        }
        else if(A && C){
            if((ya - (yl + yr) / 2.0) > eps){
                double ans = 0.5 * (xr - xl) * (ya - yc);
                if(((yr - ya) - (yc - yl)) > eps) ans += (xr - xl) * (yc - yl);
                else ans += (xr - xl) * (yr - ya);
                printf("%.3f\n", ans);
            }
            else{
                double ans = 0.5 * (xr - xl) * (yc - ya);
                if(((yr - yc) - (ya - yl)) > eps) ans += (xr - xl) * (ya - yl);
                else ans += (xr - xl) * (yr - yc);
                printf("%.3f\n", ans);
            }
        }
        else if(B && D){
            if((xb - (xr + xl) / 2.0) >eps){
                double ans = 0.5 * (yr - yl) * (xb - xd);
                if(((xr - xb) - (xd - xl)) >eps) ans += (yr - yl) * (xd - xl);
                else ans += (yr - yl) * (xr - xb);
                printf("%.3f\n", ans);
            }
            else {
                double ans = 0.5 * (yr - yl) * (xd - xb);
                if(((xr - xd) - (xb - xl)) >eps) ans += (yr - yl) * (xb - xl);
                else ans += (yr - yl) * (xr - xd);
                printf("%.3f\n", ans);
            }
        }
    }
    return 0;
}