POJ3268

 注意求单源最短路径,djkstra的复杂度是n^2比弗洛伊德的算法n^3要好很多。这道题目的大意是求单源的来回最短路径之和的最大值，用到矩阵的转置，在输入的时候用数组处理一下就行。注意求最短路径时地图的初始化一定是正无穷。

#include<iostream>
#include<cstdio>
#include<cstring>
#include<algorithm>
using namespace std;
const int inf=0x3f3f3f3f;
const int maxn=1005;
int mapp1[maxn][maxn];
int mapp2[maxn][maxn];
int dis1[maxn];
int dis2[maxn];
int vis1[maxn];
int vis2[maxn];
int temp[maxn];//记录最长时间
int n,m,x,a,b,t,maxt;
void init()
{
    memset(vis1,0,sizeof(vis1));
    memset(vis2,0,sizeof(vis2));
    memset(mapp1,inf,sizeof(mapp1));
    memset(mapp2,inf,sizeof(mapp2));
    memset(temp,0,sizeof(temp));
    memset(dis1,-1,sizeof(dis1));
    memset(dis2,-1,sizeof(dis2));
    maxt=-1;
}
void djkstra1()
{
    for(int i=1; i<=n; i++)
        dis1[i]=mapp1[x][i];
    dis1[x]=0;
    vis1[x]=1;
    for(int i=1; i<=n; i++)
    {
        int minx=inf;
        int p=x;
        for(int j=1; j<=n; j++)
            if(!vis1[j]&&minx>dis1[j])
            {
                minx=dis1[j];
                p=j;
            }
        vis1[p]=1;
        for(int j=1; j<=n; j++)
            if(!vis1[j]&&dis1[j]>dis1[p]+mapp1[p][j])
            {
                dis1[j]=dis1[p]+mapp1[p][j];
            }
    }
    for(int i=1;i<=n;i++)
        temp[i]+=dis1[i];
}
void djkstra2()
{
    for(int i=1; i<=n; i++)
        dis2[i]=mapp2[x][i];
    dis2[x]=0;
    vis2[x]=1;
    for(int i=1; i<=n; i++)
    {
        int minx=inf;
        int p=x;
        for(int j=1; j<=n; j++)
            if(!vis2[j]&&minx>dis2[j])
            {
                minx=dis2[j];
                p=j;
            }
        vis2[p]=1;
        for(int j=1; j<=n; j++)
            if(!vis2[j]&&dis2[j]>dis2[p]+mapp2[p][j])
            {
                dis2[j]=dis2[p]+mapp2[p][j];
            }
    }
    for(int i=1;i<=n;i++)
        temp[i]+=dis2[i];
}
int main()
{
    cin>>n>>m>>x;
    init();
    for(int i=0; i<m; i++)
    {
        cin>>a>>b>>t;
        mapp1[a][b]=mapp2[b][a]=t;
    }
    djkstra1();
    djkstra2();
    for(int i=1;i<=n;i++)
        if(i!=x)maxt=max(maxt,temp[i]);
    cout<<maxt<<endl;
    return 0;
}