CF2C 计算几何

题目大意：

三个不相交的圆，且三个圆心不在同一直线上，圆外找一点，到三个圆的切线角相等（多个，则要求切线角最大），如果没有则不输出

解题思路：

• 寻思着二分切线角，判断是否三个圆有交集从而往大分往小分，但失败了。。。有个测试点，发现怎样都不对，于是就改了别的方法
• 设$R_i=dis(p,o_i)$为圆外点到第$i$个圆圆心距离，题目要求的切线角相等，其实就是要求$arcsin(\frac{r_1}{R_1})=arcsin(\frac{r_2}{R_2})=arcsin(\frac{r_3}{R_3})$
• 也就是要求$\frac{r_1}{R_1}=\frac{r_2}{R_2}=\frac{r_3}{R_3}$
• 那么就是$(\frac{r_1}{R_1}-\frac{r_2}{R_2}{)}^2+(\frac{r_3}{R_3}-\frac{r_2}{R_2}{)}^2+(\frac{r_1}{R_1}-\frac{r_3}{R_3}{)}^2=0$
• 对于这类求平面代价最小题目，很多就可以用模拟退火的思想来做
• 太久没做计算几何的题目了，思维都锈了。。。早该想到的

AC代码

#include <bits/stdc++.h>
#define ft first
#define sd second
#define pb push_back
#define IOS ios::sync_with_stdio(false), cin.tie(0), cout.tie(0) //不能跟puts混用
#define seteps(N) fixed << setprecision(N)
#define endl "\n"
const int maxn = 4;
using namespace std;
typedef long long ll;
typedef double db;
typedef pair<int, int> pii;
const ll mod = 1e9 + 7;
const db eps = 1e-5;
const db pi = acos((db)-1.0);
int sgn(db x) {
    if (fabs(x) < eps) return 0;
    else return x < 0 ? -1 : 1;
}
struct Point {
    db x, y, w;
    Point operator - (Point p) {return {x - p.x, y - p.y};}
    Point operator + (Point p) {return {x + p.x, y + p.y};}
    db operator ^ (Point p) {return x * p.y - y * p.x;}
    db dis(Point p) {return hypot(x - p.x, y - p.y);}
} ans;
struct Circle {
    Point p;
    db r;
} c[4];
db energy(Point tmp) {
    db d[4];
    for (int i = 1; i <= 3; i++) d[i] = c[i].p.dis(tmp) / c[i].r;
    db res = 0;
    for (int i = 1; i <= 3; i++) res += (d[i] - d[i % 3 + 1]) * (d[i] - d[i % 3 + 1]);
    return res;
}
int main() {
    for (int i = 1; i <= 3; i++) cin >> c[i].p.x >> c[i].p.y >> c[i].r, ans.x += c[i].p.x, ans.y += c[i].p.y;
    ans.x /= 3.0, ans.y = 3.0;
    ans.w = energy(ans);
    db T = 1;
    while (T > eps) {
        int s = 0;
        ans.w = energy(ans);
        if (energy({ans.x + T, ans.y}) < ans.w) {
            ans.x = ans.x + T;
            s = 1;
        }
        else if (energy({ans.x - T, ans.y}) < ans.w) {
            ans.x = ans.x - T;
            s = 1;
        }
        else if (energy({ans.x, ans.y + T}) < ans.w) {
            ans.y = ans.y + T;
            s = 1;
        }
        else if (energy({ans.x, ans.y - T}) < ans.w) {
            ans.y = ans.y - T;
            s = 1;
        }
        if (!s) T /= 2.0;    //说明可以在当前步长继续走
    }
    ans.w = energy(ans);
    if (ans.w < eps) cout << seteps(5) << ans.x << " " << ans.y << endl;
    return 0;
}
/*
0 0 1
10 10 1
20 0 2
*/