数据结构-pta-最短路径（Dijkstra）+输出路径

 

①这道题与边无关，不适用邻接表，甚至只需一维数组不用邻接矩阵

②多源最短路的Dijkstra写法，多了一些辅助变量判断，为保证第一跳最近，排序→贪心

③Dijkstra结构→while→初始化→for找最小值→for更新邻接、vis

在老电影“007之生死关头”（Live and Let Die）中有一个情节，007被毒贩抓到一个鳄鱼池中心的小岛上，他用了一种极为大胆的方法逃脱 —— 直接踩着池子里一系列鳄鱼的大脑袋跳上岸去！（据说当年替身演员被最后一条鳄鱼咬住了脚，幸好穿的是特别加厚的靴子才逃过一劫。）

设鳄鱼池是长宽为100米的方形，中心坐标为 (0, 0)，且东北角坐标为 (50, 50)。池心岛是以 (0, 0) 为圆心、直径15米的圆。给定池中分布的鳄鱼的坐标、以及007一次能跳跃的最大距离，你需要给他指一条最短的逃生路径 —— 所谓“最短”是指007要跳跃的步数最少。

输入格式：

首先第一行给出两个正整数：鳄鱼数量 N（≤）和007一次能跳跃的最大距离 D。随后 N 行，每行给出一条鳄鱼的 ( 坐标。注意：不会有两条鳄鱼待在同一个点上。

输出格式：

如果007有可能逃脱，首先在第一行输出007需要跳跃的最少步数，然后从第二行起，每行给出从池心岛到岸边每一步要跳到的鳄鱼的坐标 (。如果没可能逃脱，就在第一行输出 0 作为跳跃步数。如果最短路径不唯一，则输出第一跳最近的那个解，题目保证这样的解是唯一的。

输入样例 1：

17 15
10 -21
10 21
-40 10
30 -50
20 40
35 10
0 -10
-25 22
40 -40
-30 30
-10 22
0 11
25 21
25 10
10 10
10 35
-30 10

输出样例 1：

4
0 11
10 21
10 35

输入样例 2：

4 13
-12 12
12 12
-12 -12
12 -12

输出样例 2：

0

#include <bits/stdc++.h>
using namespace std;

bool vis[105]; 
int shortd[105];
int Path[105];//Dijkstra算法的辅助数组 
int res=INT_MAX;//最短步数 
int cnt[105];//最短路径 

struct point1  //点集 
{
	double x;
	double y;
};

struct fj1    //每个点到原点的距离和相应下标，用于判断第一次跳 
{
	double d;
	int index;
};

struct Graph
{
	point1 point[105];
	fj1 fj[105];
	int c[105][105];//邻接矩阵，判断两点之间是否可跳 
	int n,edge;
	double D;     //最大可跳距离 
};

bool cmp(fj1 a,fj1 b)//最小在前面 
{
	return a.d<b.d;
}

void Dijkstra(Graph G,int begin)
{
    int midest,midindex,cur;//最小值，最小值所在下标 
    
    shortd[begin]=1;//初始化 
    while(1)
    {
    	midindex=-1;
    	midest=INT_MAX;
	    	for(int i=0;i<G.n;i++)//找最小值，下一个要加入的顶点 
	    {
	    	if(vis[i]==0&&shortd[i]<midest)
	    	{
			     	midindex=i;   
			     	midest=shortd[i];
			} 
	    		
	    		         
		}

		if(midindex==-1)break;//已经全部遍历过了 
	
    	cur=midindex;
		vis[cur]=1;//此处已被访问 
    
		for(int i=0;i<G.n;i++)//更新周围数据 
		{
			
			if(G.c[cur][i]==1&&shortd[cur]+1<shortd[i])//可达并且最短，更新 
			{
			    	shortd[i]=shortd[cur]+1; 
			    	Path[i]=cur;//记录路径 
			}
			    	
		}
	}
	
    
}

void firstjump(Graph G)//判断并求出对每个第一次可以跳到的点，作为最短路算法起点， 比较每个起点得出的最短路径结果 
{
	
	for(int i=0;i<G.n;i++)
{
	int begin=G.fj[i].index;
	
	for(int j=0;j<G.n;j++)
	{
			shortd[j]=INT_MAX;
			Path[j]=-1;
	}
	
		memset(vis,0,sizeof(vis));
	//以上为初始化操作 
    if(G.fj[i].d<=G.D+7.5)//满足条件 
    {
    
    		Dijkstra(G,begin);
	
	for(int j=0;j<G.n;j++)//判断终点 
	{

		if(abs(G.point[j].x)+G.D>=50||(abs(G.point[j].y)+G.D>=50))
		{

		if(res>shortd[j])//比较不同起点的结果 
		{
			//输出最短步数 
			res=shortd[j];
			//输出路径 
			int size=0;
			int temp=j;

			for(int z=0;z<res;z++)
			{

				cnt[size++]=temp;
				if(temp!=-1)
				temp=Path[temp];
			}
		}
		
		}
	}
	}

}
}

int main()
{
	Graph G;
cin>>G.n>>G.D;

for(int i=0;i<G.n;i++)
{

	cin>>G.point[i].x>>G.point[i].y;
	G.fj[i].d= sqrt((G.point[i].x)*(G.point[i].x)+(G.point[i].y)*(G.point[i].y));
	G.fj[i].index=i; //因为要排序，下标会打乱，先记录 
	
}
sort(G.fj,G.fj+G.n,cmp);//按题意离出发点最近的优先 
for(int i=0;i<G.n;i++)
{
	for(int j=0;j<G.n;j++)
	{
	    if(sqrt(pow(G.point[i].x-G.point[j].x,2)+pow(G.point[i].y-G.point[j].y,2))<=G.D)
		G.c[j][i]=G.c[i][j]=1;//两点可达，邻接矩阵为1
		else
		G.c[j][i]=G.c[i][j]=0;	
	}
}
if(G.D<42.5)//一次跳不出去 
firstjump(G);
else
res=0;//一次跳出去 

if(res!=INT_MAX)//找得到 
{
	cout<<res+1<<endl;//算上第一次跳的 
	for(int i=res-1;i>=0;i--)//Path是倒过来存的，所以逆向输出 
	{
		cout<<G.point[cnt[i]].x<<" "<<G.point[cnt[i]].y<<endl;
	}
}

else
cout<<"0"<<endl;//找不到 
return 0;
}