MPI Maelstrom——最短路练习（Dijkstra、Bellman-Ford、SPFA、Floyd）

MPI Maelstrom

题意：程序的输入将是A的下部三角形部分。也就是说，输入的第二行将包含一个条目A（2,1）。下一行将包含两个条目A（3,1）和A（3,2），依此类推。
您的程序应输出从第一个处理器向所有其他处理器广播消息所需的最短通信时间。
分析：找到最短路径没有负边，不存在环不环的问题。可采用Dijkstra。

Dijkstra代码

#include <iostream>
#include <math.h>
#include <algorithm>
#include <cmath>
#include <string.h>
#include <stdio.h>
#include<queue>
using namespace std;
#define INF 0x3f3f3f3f
typedef pair<int,int>P;
int n,ans=-1;
int f[1005][1005];
int used[1005],d[10005];
void dijkstra()
{
    for(int i=0; i<=n; i++)
    {
        used[i]=0;
        d[i]=INF;
    }
    d[1]=0;
    for(int i=1; i<=n; i++)
    {
        int tmp=INF,pos;
        for(int j=1; j<=n; j++)
        {
            if(d[j]<tmp&&!used[j])
            {
                tmp=d[j];
                pos=j;
            }
        }
        used[pos]=1;
        for(int j=1; j<=n; j++)
        {
            if(!used[j]&&d[pos]+f[pos][j]<d[j])
            {
                d[j]=d[pos]+f[pos][j];
            }
        }
    }
    for(int j=1; j<=n; j++)
    {
        ans=max(ans,d[j]);
    }
    printf("%d\n",ans);
}
int main()
{
    scanf("%d",&n);
    for(int i=0; i<=n; i++)
    {
        for(int j=0; j<=n; j++)
        {
            if(i==j)
                f[i][j]=0;
            else
                f[i][j]=INF;
        }
    }
    for(int i=2; i<=n; i++)
    {
        for(int j=1; j<i; j++)
        {
            char s[100];
            scanf("%s",s);
            if(s[0]!='x')
            {
                f[i][j]=f[j][i]=atoi(s);
            }
        }
    }
    dijkstra();
    return 0;
}

used，d，f。

Bellman-Ford代码

#include <iostream>
#include <math.h>
#include <algorithm>
#include <cmath>
#include <string.h>
#include <stdio.h>
#include<queue>
using namespace std;
#define INF 0x3f3f3f3f
typedef pair<int,int>P;
struct B
{
    int s,e,c;
} b[10005];
int d[10005];
int ans=-1,n,k;
void Bellman_Ford()
{
    for(int i=0; i<=n; i++)
        d[i]=INF;
    d[1]=0;
    for(int i=1; i<n; i++)
    {
        for(int j=0; j<k; j++)
        {
            if(d[b[j].e]>d[b[j].s]+b[j].c)
            {
                d[b[j].e]=d[b[j].s]+b[j].c;
            }
        }
    }
    for(int i=1; i<=n; i++)
        ans=max(ans,d[i]);
    printf("%d\n",ans);
}
int main()
{
    scanf("%d",&n);
    k=0;
    for(int i=2; i<=n; i++)
    {
        for(int j=1; j<i; j++)
        {
            char s[100];
            scanf("%s",s);
            if(s[0]!='x')
            {
                b[k].s=i;
                b[k].e=j;
                b[k++].c=atoi(s);
                b[k].s=j;
                b[k].e=i;
                b[k++].c=atoi(s);
            }
        }
    }
    Bellman_Ford();
    return 0;
}

SPFA代码（Bellman-Ford的优化版）邻接矩阵

#include <iostream>
#include <math.h>
#include <algorithm>
#include <cmath>
#include <string.h>
#include <stdio.h>
#include<queue>
using namespace std;
#define INF 0x3f3f3f3f
typedef pair<int,int>P;
int n,ans=-1;
int f[1005][1005];
int used[1005],d[10005];
void spfa()
{
    for(int i=0; i<=n; i++)
    {
        used[i]=0;
        d[i]=INF;
    }
    queue<int>q;
    d[1]=0;
    used[1]=1;
    q.push(1);
    while(!q.empty())
    {
        int pos=q.front();
        q.pop();
        used[pos]=0;
        for(int j=1; j<=n; j++)
        {
            if(d[pos]+f[pos][j]<d[j])
            {
                d[j]=d[pos]+f[pos][j];
                if(!used[j])
                {
                    used[j]=1;
                    q.push(j);
                }
            }
        }
    }
    for(int j=1; j<=n; j++)
    {
        ans=max(ans,d[j]);
    }
    printf("%d\n",ans);
}
int main()
{
    scanf("%d",&n);
    for(int i=0; i<=n; i++)
    {
        for(int j=0; j<=n; j++)
        {
            if(i==j)
                f[i][j]=0;
            else
                f[i][j]=INF;
        }
    }
    for(int i=2; i<=n; i++)
    {
        for(int j=1; j<i; j++)
        {
            char s[100];
            scanf("%s",s);
            if(s[0]!='x')
            {
                f[i][j]=f[j][i]=atoi(s);
            }
        }
    }
    spfa();
    return 0;
}

Floyd代码

#include <iostream>
#include <math.h>
#include <algorithm>
#include <cmath>
#include <string.h>
#include <stdio.h>
#include<queue>
using namespace std;
#define INF 0x3f3f3f3f
typedef pair<int,int>P;
int n,ans=-1;
int f[105][105];
void floyd()
{
    for(int t=1; t<=n; t++)
    {
        for(int i=1; i<=n; i++)
        {
            for(int j=1; j<=n; j++)
            {
                f[i][j]=min(f[i][t]+f[t][j],f[i][j]);
            }
        }
    }
    for(int j=1; j<=n; j++)
    {
        ans=max(ans,f[1][j]);
    }
    printf("%d\n",ans);
}
int main()
{
    scanf("%d",&n);
    for(int i=0; i<=n; i++)
    {
        for(int j=0; j<=n; j++)
        {
            if(i==j)
                f[i][j]=0;
            else
                f[i][j]=INF;
        }
    }
    for(int i=2; i<=n; i++)
    {
        for(int j=1; j<i; j++)
        {
            char s[100];
            scanf("%s",s);
            if(s[0]!='x')
            {
                f[i][j]=f[j][i]=atoi(s);
            }
        }
    }
    floyd();
    return 0;
}