Пересечение двух окружностей

#include <iostream>
#include <fstream>
#include <iomanip>
 
using namespace std;
ifstream in("circles.in");
ofstream out("circles.out");
 
struct Point {
    double x, y;
};
 
struct Circle {
    double x, y, r;
};
 
struct Segment {
    double a, b, c;
};
 
void IntersectCircle(double cx0, double cy0, double radius0, double cx1, double cy1, double radius1) {
    double dx = cx0 - cx1;
    double dy = cy0 - cy1;
    double dist = sqrt(dx * dx + dy * dy);
 
    if (dist > radius0 + radius1) {
        out << "0" << endl;
    }
    else if (dist < abs(radius0 - radius1)) {
        out << "0" << endl;
    }
    else if ((dist == 0) && (radius0 == radius1)) {
        out << "-1" << endl;
    }
    else {
        double a = (radius0 * radius0 - radius1 * radius1 + dist * dist) / (2 * dist);
        double h = sqrt(radius0 * radius0 - a * a);
 
        double cx2 = cx0 + a * (cx1 - cx0) / dist;
        double cy2 = cy0 + a * (cy1 - cy0) / dist;
 
        if (dist == radius0 + radius1) {
            out << "1" << endl;
            out << setprecision(9) << (cx2 + h * (cy1 - cy0) / dist) << " " << setprecision(9) << (cy2 - h * (cx1 - cx0) / dist) << endl;
        }
        else {
            out << "2" << endl;
            out << setprecision(9) << (cx2 + h * (cy1 - cy0) / dist) << " " << setprecision(9) << (cy2 - h * (cx1 - cx0) / dist) << " ";
            out << setprecision(9) << (cx2 - h * (cy1 - cy0) / dist) << " " << setprecision(9) << (cy2 + h * (cx1 - cx0) / dist) << endl;
        }
    }
}
 
int main() {
    int n;
 
    double r1, r2, x1, x2, y1, y2;
 
    in >> n;
 
    for (int i = 0; i < n; i++) {
        in >> x1 >> y1 >> r1;
        in >> x2 >> y2 >> r2;
        
        IntersectCircle(x1, y1, r1, x2, y2, r2);
        out << endl;
        
    }
 
    return 0;
}