首页 > 代码库 > POJ - 2252 Frogger(Dijkstra变形题)

POJ - 2252 Frogger(Dijkstra变形题)

题意:

      题目撰写者的英语真是艰难晦涩,看了别人题解,才知道这题题意。

      两个forger 一个froger 要蹦到另外一个froger处,他们的最短距离是这样定义的 :

     The frog distance (humans also call it minimax distance) between two stones therefore is defined as the minimum necessary jump range

     over all possible paths between the two stones. 

     即为 :两个石头之间的frog distance就是在这两个石头之间的所有路径中最大跳(necessary jump range)最小的。 (以上引自aowarmen‘s blog)

题目要我们求解的frog distance和dijksta算法中的最短路径距离不是同一个概念。

分析:

理解了题意之后,我们把Dijsktra中的松弛条件改成:

dist[i] = max( dist[u],edge[u][i])  //u为某一中间节点,dist[i]表示源点到结点i的距离  (以上摘自Eucalyptus)

同时,我们可以优化距离的开根号运算,两点之间的距离定义改成:欧几里德距离的平方和。

这样一来,我们只需要对结果开根号即可,节省中间的运算时间。

提交时请不要用g++,要用c++编译器,同样的代码,我用G++提交时总是WA

 

代码:

优先队列实现的Dijkstra模板,只要修改对应的松弛条件即可。

#include<iostream>#include<cstdio>#include<string.h>#include<queue>#include<cmath>#include<vector>using namespace std;#define maxn 201#define inf 0x3f3f3f3f#define eps 0.0001typedef pair<double,int> P;double max(double a,double b){	return (a>b?a:b);}struct points {	int x, y;	points(int xi,int yi){ x=xi,y=yi;}	points(){x=0,y=0;}};points nodes[maxn];double map[maxn][maxn];double dist[maxn];double cal_dist(points a,points b){	return (a.x-b.x)*(a.x-b.x)+(a.y-b.y)*(a.y-b.y);}void init(int n){	for(int i=1;i<=n;i++){		for(int j=1;j<=n;j++){			map[i][j]=cal_dist(nodes[i],nodes[j]);            map[j][i]=map[i][j];		}	}}void dijkstra(int s,int n){	priority_queue<P,vector<P>,greater<P> > Q;	bool visited[maxn];	memset(visited,0,sizeof(visited));	for(int i=1;i<=n;i++)		dist[i]=inf;	dist[s]=0.0;	Q.push(P(0.0,s));	while(!Q.empty()){		int v=Q.top().second;		Q.pop();		if(visited[v]) continue;		visited[v]=true;		for(int i=1;i<=n;i++){			double local_max=max(dist[v],map[v][i]);			if(dist[i]-local_max>eps){				dist[i]=local_max;				Q.push(P(local_max,i));			}		}	}}int main(){	//freopen("in.txt","r",stdin);	int cases=0,n;	while(scanf("%d",&n)!=EOF && n){		cases++;		for(int i=1;i<=n;i++)			scanf("%d %d",&nodes[i].x,&nodes[i].y);		init(n);		dijkstra(1,n);		printf("Scenario #%d\n",cases);		printf("Frog Distance = %.3lf\n\n",sqrt(dist[2]));	}}

 

普通方法实现的Dijkstra算法

#include<iostream>#include<cstdio>#include<string.h>#include<queue>#include<cmath>#include<vector>using namespace std;#define maxn 210#define inf 0x3f3f3f3ftypedef pair<double,int> P;double max(double a,double b){	return (a>b?a:b);}struct point {	double x, y;	point(double xi,double yi){ x=xi,y=yi;}	point(){x=0,y=0;}};point nodes[maxn];double map[maxn][maxn];double dist[maxn];double cal_dist(point a,point b){	return sqrt((a.x-b.x)*(a.x-b.x)+(a.y-b.y)*(a.y-b.y));}void init(int n){	for(int i=1;i<=n;i++){		for(int j=1;j<=n;j++){			if(i==j) map[i][j]=0.0;			else				map[i][j]=cal_dist(nodes[i],nodes[j]);			map[j][i]=map[i][j];		}	}}void dijkstra(int s,int n){	bool visited[maxn];	memset(visited,0,sizeof(visited));	visited[s]=true;	for(int i=1;i<=n;i++)		dist[i]=map[s][i];	dist[s]=0.0;	for(int i=1;i<=n;i++){		double min=inf;		int u=0;		for(int j=1;j<=n;j++)			if(!visited[j] && dist[j]<min){				min=dist[j];				u=j;			}		visited[u]=true;		if(min==inf)  break;		for(int j=1;j<=n;j++){			double tmp=max(dist[u],map[u][j]);			if(!visited[j]&&dist[j]>tmp)				dist[j]=tmp;		}	}}int main(){	//freopen("in.txt","r",stdin);	int cases=0,n;	while(scanf("%d",&n)!=EOF && n){		cases++;		for(int i=1;i<=n;i++)			scanf("%lf %lf",&nodes[i].x,&nodes[i].y);		init(n);		dijkstra(1,n);		printf("Scenario #%d\n",cases);		printf("Frog Distance = %.3lf\n\n",dist[2]);	}}

POJ - 2252 Frogger(Dijkstra变形题)