题目链接:传送门
题目大意:
给出两条平行线和 n n n个圆,计算从一条线走到另一条线所消耗的体力。同时,在线上、圆上、园内走不消耗体力,求体力的最小花费。
解题思路:
一开始没想到是最短路,后来仔细一想,除了最短路也没啥其他的了。将线和圆都想象成一个个点,每个点之间的距离满足以下关系:
L
L
L
1
1
1 到
L
L
L
2
2
2 之间连边权值
∣
C
1
−
C
2
∣
A
2
+
B
2
\frac{| C_1- C_2 |} {\sqrt{A^2 + B^2}}
A2+B2∣C1−C2∣
线
L
L
L与圆
i
i
i之间连边权值
m
a
x
(
0
,
d
(
O
i
,
L
1
)
−
r
i
)
max(0, d(O_i,L_1)-r_i)
max(0,d(Oi,L1)−ri)
圆
i
i
i与圆
j
j
j之间连边权值
m
a
x
(
0
,
d
(
O
i
,
O
j
)
−
r
i
−
r
j
)
max(0, d(O_i,O_j)-r_i-r_j)
max(0,d(Oi,Oj)−ri−rj)
最后求最短路即可
代码思路:
注意要用double,double的绝对值是fabs函数
核心:将线和圆看成一个个点来建图,灵活地运用最短路知识
#include<bits/stdc++.h>
using namespace std;
const int N = 1005;
const double INF= 1e10+10;
int m, n, vis[N];
double dis[N], mp[N][N];
struct node {
double x, y, r;
} cc[1005];
void init() {
for (int i=0; i<=n+1; i++)
for (int j=0; j<=n+1; j++) {
if (i == j) mp[i][j] = 0;
else mp[i][j] = INF;
}
}
void creatgraph() {
double a, b, c1, c2;
scanf("%lf%lf%lf%lf", &a, &b, &c1, &c2);
for(int i=1; i<=n; i++)
scanf("%lf%lf%lf", &cc[i].x, &cc[i].y, &cc[i].r);
mp[0][n+1] = mp[n+1][0] = fabs(c1-c2) / sqrt(a*a+b*b);
for(int i=1; i<=n; i++) {
double tp = fabs(cc[i].x*a+cc[i].y*b+c1) / sqrt(a*a+b*b) - cc[i].r;
if(tp<1e-4) tp=0;
mp[0][i] = mp[i][0] = tp;
tp = fabs(cc[i].x*a+cc[i].y*b+c2) / sqrt(a*a+b*b) - cc[i].r;
if(tp<1e-4) tp=0;
mp[n+1][i] = mp[i][n+1] = tp;
}
for(int i=1; i<=n; i++)
for(int j=1; j<=n; j++) {
if(i==j) continue;
double tp = fabs(cc[i].x-cc[j].x)*fabs(cc[i].x-cc[j].x) + \
fabs(cc[i].y-cc[j].y)*fabs(cc[i].y-cc[j].y);
tp = sqrt(tp) - cc[i].r - cc[j].r;
if(tp <= 1e-4) tp = 0;
mp[i][j] = mp[j][i] = tp;
}
}
void dijkstra(int start) {
for (int i=0; i<=n+1; i++)
dis[i] = mp[start][i];
memset(vis, 0, sizeof(vis));
vis[start] = 1;
for (int i=1; i<=n; i++) {
/*找出离起点最近的点*/
double minn = INF;
int k = -1;
for (int j=0; j<=n+1; j++) {
if (!vis[j] && dis[j]<minn) {
minn = dis[j];
k = j;
}
}
vis[k] = 1;
for (int j=0; j<=n+1; j++) { //松弛操作,找到媒介使得出现新的最短路
if (!vis[j] && mp[k][j] < INF)
dis[j] = min(dis[j], dis[k] + mp[k][j]);
}
}
}
int main() {
scanf("%d", &n); //n个顶点
init(); //初始化地图
creatgraph(); //建图
dijkstra(0); //核心迪杰斯特拉算法
printf("%.6f", dis[n+1]); //输出最短路结果
}