Биссектриса угла по 3 точкам

#include <iostream>
#include <vector>
#include <cmath>
#include <iomanip>
 
#define PI acos(-1.0)
 
#define fast_cin() std::ios_base::sync_with_stdio(false); std::cin.tie(NULL); std::cout.tie(NULL)
#define precise_cout std::cout << std::setprecision(12)
 
class Point {
public:
    explicit Point(double _x = 0, double _y = 0) : x(_x), y(_y) {};
 
    void read_point() { std::cin >> x >> y; }
 
    void print_point() const {
        std::cout << "(" << x << "," << y << ")" << " ";
    }
 
    Point operator+(const Point &other) const {
        return Point(x + other.x, y + other.y);
    }
 
    Point operator-(const Point &other) const {
        return Point(x - other.x, y - other.y);
    }
 
    Point operator/(double num) const { // На 0 поделим, получим NaN - ибо нефиг так делить
        return Point(x / num, y / num);
    }
 
    Point operator*(double num) const {
        return Point(x * num, y * num);
    }
 
 
public:
    double x, y;
};
 
class Line {
public:
    explicit Line(double _a = 0, double _b = 0, double _c = 0) : A(_a), B(_b), C(_c) {}; // По коэффициентам
    Line(Point a, Point b) : A(b.y - a.y), B(a.x - b.x), C(-(A * a.x + B * a.y)) {}; // Две точки
    Line(Point a, Vector v) : A(v.b.x), B(v.b.y), C(-(A * a.x + B * a.y)) {}; // Точка и вектор нормали к этой прямой
 
    [[nodiscard]] bool is_point_on_the_line(Point p) const {
        return A * p.x + B * p.y + C == 0;
    }
 
    [[nodiscard]] bool are_two_points_on_the_same_side_of_line(Point p1, Point p2) const {
        return (A * p1.x + B * p1.y + C) * (A * p2.x + B * p2.y + C) > 0;
    }
 
    [[nodiscard]] Line get_perpendicular_line(Point p) const {
        double a = B, b = -A, c = -a * p.x - b * p.y;
        return Line(a, b, c);
    }
 
    [[nodiscard]] Line get_line_on_dist(double dist) const {
        return Line(A, B, dist * sqrt(pow(A, 2) + pow(B, 2)) + C);
    }
 
    [[nodiscard]] double get_dist_to_point(Point p) const {
        return std::abs(A * p.x + B * p.y + C) / sqrt(pow(A, 2) + pow(B, 2));
    }
 
    void print_line_equation_coefficients() const {
        precise_cout << A << " " << B << " " << C << "\n";
    }
 
    void read_line_equation_coefficients() {
        std::cin >> A >> B >> C;
    }
 
public:
    double A, B, C;
};
 
double sq(Line l) {
    return sqrt(pow(l.A, 2) + pow(l.B, 2));
}
 
int main() {
    fast_cin();
 
    Point x, y, z;
    x.read_point();
    y.read_point();
    z.read_point();
 
    Line XY(x, y), XZ(x, z);
 
    double t1 = sq(XY), t2 = sq(XZ);
 
    auto b1 = Line(XY.A * t2 + XZ.A * t1, XY.B * t2 + XZ.B * t1, XY.C * t2 + XZ.C * t1),
            b2 = Line(XY.A * t2 - XZ.A * t1, XY.B * t2 - XZ.B * t1, XY.C * t2 - XZ.C * t1);
 
    if (!b1.are_two_points_on_the_same_side_of_line(y, z)) {
        b1.print_line_equation_coefficients();
    } else { b2.print_line_equation_coefficients(); }
    return 0;
}