题目链接:Everything Has Changed
题意
在平面直角坐标系内有一个圆心在原点,半径为 R R 的初始圆,以及 个圆心在 (xi,yi) ( x i , y i ) ,半径为 ri r i 的圆,这 m m 个圆两两互不相交且不会包含整个初始圆,这 个小圆将对初始圆进行切割,问最终初始圆的周长。
输入
第一行为一个整数 T (1≤T≤1000) T ( 1 ≤ T ≤ 1000 ) ,接下去有 T T 组数据,每组数据第一行为两个整数 ,接下去 m m 行每行三个整数 。
输出
输出最终圆的周长,误差在 10−6 10 − 6 以内都认为答案正确。
样例
输入 |
---|
1 4 10 6 3 5 10 -4 3 -2 -4 4 0 9 1 |
输出 |
81.62198908430238475376 |
提示 |
切割圆的方式如下: 其中红线为圆的周长。 |
题解
对于两圆相离或者外切的情况,对周长的贡献都为 0 0 ,对于两圆内含的情况,贡献也为零,因此只要考虑两圆相交与内切的情况,先考虑第一种相交(两圆圆心在交点连线的两侧):
设 ,则可以联立方程:
⎧⎩⎨⎪⎪a2+h2=r21b2+h2=r22a+b=d { a 2 + h 2 = r 1 2 b 2 + h 2 = r 2 2 a + b = d解方程得到 b=r22−r21+d22d b = r 2 2 − r 1 2 + d 2 2 d ,而对于两圆圆心在交点连线同侧的情况:
只需要将上面方程组的第 3 3 个方程改为 即可,得到 b=−r22−r21+d22d b = − r 2 2 − r 1 2 + d 2 2 d ,发现和上式差一个负号,要求圆 O2 O 2 的贡献,就是求圆 O2 O 2 在圆 O1 O 1 内部的圆弧长度,因此如果直接对 br2 b r 2 取 arccos arccos ,第二种情况的负号正好使得圆 O2 O 2 的圆心角计算得到大于 90° 90 ° 的部分。
最后还需要减掉圆 O1 O 1 在圆 O2 O 2 内部的圆弧长度,计算方式与第一种情况相同。
过题代码
#include <iostream>
#include <cstdio>
#include <cstdlib>
#include <cmath>
#include <climits>
#include <cstring>
#include <string>
#include <vector>
#include <list>
#include <queue>
#include <stack>
#include <map>
#include <set>
#include <bitset>
#include <algorithm>
#include <functional>
#include <iomanip>
using namespace std;
#define LL long long
const double eps = 1e-8;
const double PI = acos(-1.0);
const int maxn = 100 + 100;
struct Circle {
double x, y;
double r, val;
};
int T, n;
double R;
Circle c[maxn];
int sign(const double &x) {
if(fabs(x) < eps) {
return 0;
}
if(x < 0) {
return -1;
}
return 1;
}
void get_val(Circle &c, int Index) {
double dis = c.x * c.x + c.y * c.y;
int sgn = sign(dis - (R + c.r) * (R + c.r));
if(sgn > 0) {
c.val = 0;
return ;
}
sgn = (R - c.r) * (R - c.r) - dis;
if(sgn > 0) {
c.val = 0;
return ;
}
double a = (c.r * c.r - R * R + dis) / (2 * sqrt(dis));
double theta = acos(a / c.r) * 2;
c.val = c.r * theta;
double b = (R * R - c.r * c.r + dis) / (2 * sqrt(dis));
theta = acos(b / R) * 2;
c.val -= R * theta;
}
bool Set(Circle &a, Circle &b) {
double dis = (a.x - b.x) * (a.x - b.x) + (a.y - b.y) * (a.y - b.y);
int sgn = sign(dis - (a.r + b.r) * (a.r + b.r));
if(sgn > 0) {
return false;
}
sgn = (a.r - b.r) * (a.r - b.r) - dis;
if(sgn > 0) {
return false;
}
return true;
}
int main() {
#ifdef LOCAL
freopen("test.txt", "r", stdin);
// freopen("out.txt", "w", stdout);
#endif // LOCAL
ios::sync_with_stdio(false);
scanf("%d", &T);
while(T--) {
scanf("%d%lf", &n, &R);
for(int i = 1; i <= n; ++i) {
scanf("%lf%lf%lf", &c[i].x, &c[i].y, &c[i].r);
get_val(c[i], i);
}
double ans = 2 * PI * R;
for(int i = 1; i <= n; ++i) {
ans += c[i].val;
}
printf("%.10f\n", ans);
}
return 0;
}