计算几何 - 圆 - 洛谷 P1652

计算几何 - 圆 - 洛谷 P1652

$$\begin{gathered} 给出n个圆,保证任意两个圆都不相交。然后给出两个点(x_1,y_1),(x_2,y_2)， \\ 保证均不在某个圆上，要从(x_1,y_1)->(x_2,y_2)画条曲线，问这条曲线最少穿过多少次圆的边界？ \end{gathered}$$

输入格式

第一行为一个整数 n，表示圆的个数；

第二行是 n 个整数，表示 n 个圆的 x 坐标；

第三行是 n 个整数，表示 n 个圆的 y 坐标；

第四行是 n 个整数，表示 n 个圆的半径 r；

第五行是四个整数 $x_1,y_1,x_2,y_2$

输出格式

仅一个整数，表示最少要穿过多少次圆的边界。

输入输出样例

输入：

7
1 -3 2 5 -4 12 12
1 -1 2 5 5 1 1
8 1 2 1 1 1 2
-5 1 12 1

输出：

3

数据范围

$对于100\%的数据，1\le n\le 50，|x|,|y|\le 1000，1\le r\le 1000。保证圆之间没有公共点。$

---

分析：

$①、首先考虑特殊情况：$

$当两点在同一个小圆内部或者两个点都未被任何圆包含，此时输出0.$

$②、一般情况：$

$当其中一个点被某个圆包围时，我们可将整个圆视作一{个}^{\prime \prime }大{点}^{\prime \prime }，如下图$

$在两{个}^{\prime \prime }大{点}^{\prime \prime }之间连线，无需穿过任何边界，$

$而{对}^{\prime \prime }小{点}^{\prime \prime }A而言，至少穿过2层边界才能与B相连，对小{点}^{\prime \prime }{B}^{\prime \prime }而言，至少穿过1层边界才能与A相连。$

$事实上，对于任意一个圆，若A、B两点在该圆的两侧，则至少要穿过这层边界才能在A、B间连线。$

$于是，问题转化为，统计有多少个圆，将A、B两点划分在圆的内外两侧。$

$我们仅需保证，一个点到圆心的距离小于半径，另一个点到圆心的距离大于半径，即可满足条件。$

代码：

#include<iostream>
#include<cstring>
#include<cstdio>
#include<algorithm>
#include<cmath>

using namespace std;

const int N=55;

struct Point
{
    double x,y,r;
}V[N];

int n;
Point a,b;

double get_dis(Point a,Point b) //两点之间的距离
{
    double dx=a.x-b.x, dy=a.y-b.y;
    return sqrt(dx*dx+dy*dy);
}

int main()
{
    cin>>n;
    for(int i=0;i<n;i++) cin>>V[i].x;
    for(int i=0;i<n;i++) cin>>V[i].y;
    for(int i=0;i<n;i++) cin>>V[i].r;
    cin>>a.x>>a.y>>b.x>>b.y;
    
    int cnt=0;
    for(int i=0;i<n;i++)
        if((get_dis(a,V[i])<V[i].r) ^ (get_dis(b,V[i])<V[i].r)) cnt++;
    
    cout<<cnt<<endl;
    
    return 0;
}