Нужно найти площадь пересечения двух окружностей (формулы взял тут):
S = S1 + S2
где,
- R1 — радиус первой окружности;
- R2 — радиус второй окружности;
- D — расстояние между центрами окружностей.
Код:
float find_area(float x1, float y1, float r1, float x2, float y2, float r2) {
// Расстояние между центрами окружностей
float distance = sqrt((x1 - x2) * (x1 - x2) + (y1 - y2) * (y1 - y2));
float f1, f2;
float s1, s2;
if(distance >= r1 + r2) {
// Если не пересекаются (можно дополнить вычислением площади внутренней окружности)
return 0;
}else if(distance <= fabs(r1 - r2)) {
// Если окружность внутри другой
return 0;
}
else {
// Если пересекаются
f1 = 2 * acos((r1*r1 - r2*r2 - distance*distance) / (2 * r1 * distance));
f2 = 2 * acos((r2*r2 - r1*r1 - distance*distance) / (2 * r2 * distance));
s1 = (r1*r1 * (f1 - sin(f1))) / 2;
s2 = (r2*r2 * (f2 - sin(f2))) / 2;
return s1 + s2;
}
return 0;
}
Площадь пересечения окружностей
(x1 = 0; y1 = 0; r1 = 2) ∩ (x2 = 3; y2 = 0; r2 = 4)
должна быть равна примерно 9.57019, но у меня выводит что-то непонятное:
-1.#IND00
+D^2
, а у вас минус.double
вместоfloat
в подобных расчётах, если нет специальных ограничений на этот счёт.