P1907 设计道路 题解

题目传送门

解法

思路

求“最小不满值”，还给了一些点和边，明显是求最短路。

算法流程

记走 Dirt Road 和 Rome Road 一个单位长度时的不满值为 $D$ 和 $R$。

前置知识：欧几里得距离、Dijkstra。

• 按照输入的建 Rome Road。边权为 $R\times \sqrt{(x_1-x_2{)}^2+(y_1-y_2{)}^2}$。$x_1,y_1,x_2,y_2$ 分别为两点坐标。
• 输入起点（码头）和终点（家）。可以将它们记为 $0$ 号和 $n+1$ 号点。
• 遍历任意两点建 Dirt Road。边权为 $D\times \sqrt{(x_1-x_2{)}^2+(y_1-y_2{)}^2}$。
• 进行 Dijkstra。

代码

#include <iostream>
#include <queue>
#include <vector>
#include <cmath>
using namespace std;

struct node//节点
{
	double x,y;
};
struct ndfp//给堆优化用的
{
	int id;
	double d;
	friend bool operator < (ndfp a,ndfp b)
	{
		return a.d<b.d;
	}
};

//计算两点间欧几里得距离
inline double dist(double x1,double y1,double x2,double y2)
{
	return sqrt((x1-x2)*(x1-x2)+(y1-y2)*(y1-y2));
}

int main()
{
	double D,R;//不满值
	int n;
	scanf("%lf %lf\n%d",&D,&R,&n);
	node p[n+2];
	double e[n+2][n+2];//邻接矩阵
	fill(*e,*e+(n+2)*(n+2),0x3f3f3f3f);
	for(int i=1;i<=n;i++)
	{
		scanf("%lf %lf",&p[i].x,&p[i].y);
	}
	int a,b;
	while(true)
	{
		scanf("%d %d",&a,&b);
		if(!(a||b))
			break;
		e[a][b]=R*dist(p[a].x,p[a].y,p[b].x,p[b].y);//建 Rome Road
		e[b][a]=R*dist(p[a].x,p[a].y,p[b].x,p[b].y);
	}
	double bx,by,ex,ey;
	scanf("%lf %lf %lf %lf",&bx,&by,&ex,&ey);
	p[0].x=bx;
	p[0].y=by;
	p[n+1].x=ex;
	p[n+1].y=ey;
	for(int i=0;i<=n+1;i++)
	{
		for(int j=0;j<=n+1;j++)
		{
			if(i!=j)//建 Dirt Road
				e[i][j]=min(e[i][j],D*dist(p[i].x,p[i].y,p[j].x,p[j].y));
		}
	}
	priority_queue<ndfp> q;//堆
	q.push({0,0});
	ndfp t;
	double dis[n+2];
	fill(dis,dis+(n+2),0x3f3f3f3f);
	dis[0]=0;
	while(!q.empty())//Dijkstra（我写堆优化因为好写）
	{
		t=q.top();
		q.pop();
		if(t.d!=dis[t.id])
		{
			continue;
		}
		for(int i=0;i<=n+1;i++)
		{
			if(dis[i]>dis[t.id]+e[t.id][i])
			{
				dis[i]=dis[t.id]+e[t.id][i];
				q.push({i,dis[i]});
			}
		}
	}
	printf("%.4lf",dis[n+1]);
	return 0;
}