(Dijkstra+bellon_ford+floyd) Silver Cow Party (P3268)

首先向这篇文章学习并比较了四种求最短路的基本做法：http://www.cnblogs.com/XBWer/archive/2012/08/15/2640990.html

题意：给出N个点，M条单向路。求每个点到指定的一个点去开party，

            求所有来回最短路中最长的一条

解法：求两次party点对所有点的最短路，一次是开party点每一个点的最短距离。

           再把路径的方向倒一下再求一次即可。

网址：http://poj.org/problem?id=3268

（1）DIJKSTRA： 104MS

#include<iostream>
using namespace std;

#define N 1001
#define inf 0x3f3f3f3f
int map[N][N];
int n,m,x;
int hi[N],a[N],b[N];
bool vis[N];

int main()
{
	freopen("in","r",stdin);
	int i,j,k;
	cin>>n>>m>>x;
	memset(map,0x3f,sizeof(map));
	while (m--) cin>>i>>j>>k,map[i][j]=k;
 	for (i=1;i<=n;i++) a[i]=map[x][i],b[i]=map[i][x];
	vis[x]=true,k=n-1,j=x,a[x]=0;
	while (k--)
	{
		int big=inf;
		for (i=1;i<=n;i++) if (!vis[i] && big>a[i])
			big=a[i],j=i;
		if (vis[j]) break;
		vis[j]=true;
		for (i=1;i<=n;i++) if (!vis[i])
			a[i]=min(a[i],a[j]+map[j][i]);
	}
	memset(vis,false,sizeof(vis));
	k=n-1,j=x,vis[x]=true,b[x]=0;
	while (k--)
	{
		int big=inf;
		for (i=1;i<=n;i++) if (!vis[i] && big>b[i])
			big=b[i],j=i;
		if (vis[j]) break;
		vis[j]=true;
		for (i=1;i<=n;i++)  if (!vis[i])
			b[i]=min(b[i],b[j]+map[i][j]);
	}
	for (k=0,i=1;i<=n;i++) k=max(a[i]+b[i],k);
	cout<<k<<endl;
	return 0;
}

同理有：

（2）BELLON_FORD    94MS

#include<iostream>
using namespace std;

#define N 1010
struct point
{
	int from,to,v;
}go1[100010],go2[100010];
int n,m,x;
int dis1[N],dis2[N];

void bellon_ford(point go[],int dis[])
{
	int i,j,k;
	for (i=1;i<=n;i++) dis[i]=1000000000;
	for (i=0,dis[x]=0;i<n;i++)
	{
		bool isok=false;
		for (j=0;j<m;j++) if(dis[go[j].to]>dis[go[j].from]+go[j].v)
			dis[go[j].to]=dis[go[j].from]+go[j].v,isok=true;			
		if (!isok) break;
	}
}

int main()
{
	freopen("in","r",stdin);
	int i,j,k;
	cin>>n>>m>>x;
	for (i=0;i<m;i++)
		cin>>go1[i].from>>go1[i].to>>go1[i].v,go2[i]=go1[i],swap(go2[i].from,go2[i].to);
	bellon_ford(go1,dis1);
	bellon_ford(go2,dis2);
	for (i=1,k=0;i<=n;i++) k=max(k,dis1[i]+dis2[i]);
	cout<<k<<endl;
}

SPFA: 优化前：150MS左右 ；  优化后：110MS

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

#define N 1001
int n,m,x;
struct my
{
	int to,v;
}cur;
vector<my> v1[N],v2[N];
int dis1[N],dis2[N];

void spfa(vector<my> v[],int dis[])
{
	int i,j,k;
	for (i=1;i<=n;i++) dis[i]=1000000000;
	bool vis[N];
	memset(vis,false,sizeof(vis));
	deque<int> q;
	q.push_back(x),vis[x]=true,dis[x]=0;
	while (!q.empty())
	{
		int d=q.front();
		q.pop_front();
		vis[d]=false;
		for (i=0;i<v[d].size();i++)
		{
			cur=v[d][i];
			if (dis[cur.to]<=cur.v+dis[d]) continue;
			dis[cur.to]=cur.v+dis[d];
			if (!vis[cur.to])
			{
				vis[cur.to]=true;
				if (q.empty() || q.front()>=dis[cur.to])
					q.push_front(cur.to);
				else
					q.push_back(cur.to);
			}
		}
	}
}

int main()
{
	freopen("in","r",stdin);
	int i,j,k;
	cin>>n>>m>>x;
	while (m--) cin>>i>>j>>k,cur.to=j,cur.v=k,v1[i].push_back(cur),cur.to=i,v2[j].push_back(cur);
	spfa(v1,dis1);
	spfa(v2,dis2);
	for (i=1,k=0;i<=n;i++) k=max(k,dis1[i]+dis2[i]);
	cout<<k<<endl;
}

还有FLOYD算法，不过在这个题不能实现。因为N最大为1000，时间复杂度太高。