HDU 6354 Everything Has Changed（计算几何）

最新推荐文章于 2019-07-25 20:18:00 发布

ExRoc

最新推荐文章于 2019-07-25 20:18:00 发布

阅读量293

点赞数

本文链接：https://blog.csdn.net/CSDNjiangshan/article/details/81749218

版权

题目链接：Everything Has Changed

题意

在平面直角坐标系内有一个圆心在原点，半径为 $R$ 的初始圆，以及 $m$ 个圆心在 $(x_i,y_i)$ ，半径为 $r_i$ 的圆，这 $m$ 个圆两两互不相交且不会包含整个初始圆，这 $m$ 个小圆将对初始圆进行切割，问最终初始圆的周长。

输入

第一行为一个整数 $T~(1\leq T\leq1000)$ ，接下去有 $T$ 组数据，每组数据第一行为两个整数 $m,R~(1\leq m\leq100)$ ，接下去 $m$ 行每行三个整数 $x_i,y_i,r_i~(-1000\leq x_i,y_i\leq1000,1\leq R,r_i\leq1000)$ 。

输出

输出最终圆的周长，误差在 $10^{-6}$ 以内都认为答案正确。

样例

输入
1 4 10 6 3 5 10 -4 3 -2 -4 4 0 9 1
输出
81.62198908430238475376
提示
切割圆的方式如下：其中红线为圆的周长。

题解

对于两圆相离或者外切的情况，对周长的贡献都为 $0$ ，对于两圆内含的情况，贡献也为零，因此只要考虑两圆相交与内切的情况，先考虑第一种相交（两圆圆心在交点连线的两侧）：

设 $AG=h,O_1O_2=d,O_1G=a,O_2G=b$ ，则可以联立方程：

$⎧ ⎩ ⎨ ⎪ ⎪ a 2 + h 2 = r 21 b 2 + h 2 = r 22 a + b = d$ $\left\{ \begin{array}{c} a^2+h^2=r_1^2\\ b^2+h^2=r_2^2\\ a+b=d \end{array}\right.$ 解方程得到 $b=\frac{r_2^2-r_1^2+d^2}{2d}$ ，而对于两圆圆心在交点连线同侧的情况：

只需要将上面方程组的第 $3$ 个方程改为 $a-b=d$ 即可，得到 $b=-\frac{r_2^2-r_1^2+d^2}{2d}$ ，发现和上式差一个负号，要求圆 $O_2$ 的贡献，就是求圆 $O_2$ 在圆 $O_1$ 内部的圆弧长度，因此如果直接对 $\frac{b}{r_2}$ 取 $\arccos$ ，第二种情况的负号正好使得圆 $O_2$ 的圆心角计算得到大于 $90°$ 的部分。
最后还需要减掉圆 $O_1$ 在圆 $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;
}