题目链接:传送门
题目大意:
给出两条平行线和 个圆,计算从一条线走到另一条线所消耗的体力。同时,在线上、圆上、园内走不消耗体力,求体力的最小花费。
解题思路:
一开始没想到是最短路,后来仔细一想,除了最短路也没啥其他的了。将线和圆都想象成一个个点,每个点之间的距离满足以下关系:
到
之间连边权值
线
与圆
之间连边权值
圆
与圆
之间连边权值
最后求最短路即可
代码思路:
注意要用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]); //输出最短路结果
return 0;
}