ZOJ1298 POJ1135 Domino Effect，Dijkstra 算法

最新推荐文章于 2016-08-04 16:37:36 发布

neofung

最新推荐文章于 2016-08-04 16:37:36 发布

阅读量751

点赞数

CC 4.0 BY-SA版权

分类专栏：最短路径 ZOJ POJ 文章标签：算法 path system float email im

本文链接：https://blog.youkuaiyun.com/neofung/article/details/6684841

ZOJ 同时被 3 个专栏收录

125 篇文章

订阅专栏

POJ

105 篇文章

订阅专栏

最短路径

13 篇文章

订阅专栏

本文探讨了Dijkstra算法在解决DominoEffect问题中的应用，详细介绍了如何使用该算法来查找最短路径，并考虑了特殊情况下单一顶点的影响。通过实例演示，展示了算法在实际问题解决中的有效性和效率。

摘要生成于 C知道，由 DeepSeek-R1 满血版支持，前往体验 >

这是很简单Dijkstra算法的应用。把遍历过的顶点记录下来，每一次找到一个顶点，就把最短路径到达这个顶点的边记录下来，在最后找到所有没有记录的边，看看边上是否存在碰撞点。

这题要考虑一种情况，就是一个顶点的情况。这可是贡献了我多个WA啊~~~~~

/*******************************************************************************
 # Author : Neo Fung
 # Email : neosfung@gmail.com
 # Last modified: 2011-08-13 20:08
 # Filename: ZOJ1298 POJ1135 Domino Effect.cpp
 # Description : 
 ******************************************************************************/
// ZOJ1298 POJ1135 Domino Effect.cpp : Defines the entry point for the console application.
//

// #include "stdafx.h"


#include <fstream>
#include <stdio.h>
#include <iostream>
#include <string.h>
#include <string>
#include <vector>
#include <stack>
#include <queue>
#include <deque>
#include <map>
#include <math.h>
#include <algorithm>
#include <numeric>
#include <functional>
#include <memory.h>
#include <limits>


using namespace std;

#define INFINITE numeric_limits<int>::max()

int main(void)
{
// 	ifstream cin("data.txt");
	int n,m;
	int x,y,len;
	int map[502][502];
	int visit[502];
	int lastvisit[502];
	int path[502];
	int ncases(1);

	while(cin>>n>>m &&(n||m))
	{

		if (n==1)
		{
			printf("System #%d\nThe last domino falls after 0.0 seconds, at key domino 1.\n\n",ncases++);
			continue;
		}
		for(int i=1;i<=n;++i)
		{
			for(int j=1;j<=n;++j)
				map[i][j]=INFINITE;
			map[i][i]=0;
		}

		for(int i=0;i<m;++i)
		{
			cin>>x>>y>>len;
			map[x][y]=len;
			map[y][x]=len;
		}

		memset(visit,0,sizeof(visit));
		visit[1]=1;

		memset(lastvisit,0,sizeof(lastvisit));
		memset(path,0,sizeof(path));
		for(int i=1;i<=n;++i)
		{
			path[i]=map[1][i];
			if(path[i]<INFINITE)
				lastvisit[i]=1;
		}


		for(int i=1;i<n;++i)
		{
			int min(INFINITE),minindex(1) ;

			for(int j=1;j<=n;++j)
			{
				if(!visit[j] && path[j]< min)
				{
					min=path[j];
					minindex=j;
				}
			}

			visit[minindex]=1;
			map[lastvisit[minindex]][minindex]=INFINITE;
			map[minindex][lastvisit[minindex]]=INFINITE;
			
			for(int j=1;j<=n;++j)
			{
				if(!visit[j] && map[minindex][j] < INFINITE
						&& path[j]>map[minindex][j]+path[minindex])
				{
					path[j]=map[minindex][j]+path[minindex];
					lastvisit[j]=minindex;
				}
			}
		}

		float total(0.0);
		int total_index;
		for(int i=1;i<=n;++i)
			if(path[i]>total)
			{
				total = path[i];
				total_index = i;
			}

		int flag(0);
		int a,b;
		for(int i=1;i<=n;++i)
			for(int j=1;j<=n;++j)
			{
				if(map[i][j]<INFINITE && (total<(path[i]+path[j]+map[i][j])/2.0))
				{
					total = (path[i]+path[j] + map[i][j])/2.0;
					flag=1;
					a=i;
					b=j;
					map[i][j]=INFINITE;
					map[j][i]=INFINITE;
				}
			}

		cout<<"System #"<<ncases++<<endl;
		if(flag)
		{
			if (a>b)
			{
				swap(a,b);
			}
			printf("The last domino falls after %0.1f seconds, between key dominoes %d and %d.\n",total,a,b);
		}
		else
		{
			printf("The last domino falls after %0.1f seconds, at key domino %d.\n",total,total_index);
		}
		cout<<endl;

	}


	return 0;
}