ZZULIOJ 1925: 信道安全（Dijkstra，Spfa）

1925: 信道安全

题目描述

Alpha 机构有自己的一套网络系统进行信息传送。情报员 A 位于节点 1，他准备将一份情报
发送给位于节点 n 的情报部门。可是由于最近国际纷争，战事不断，很多信道都有可能被遭到监
视或破坏。
经过测试分析，Alpha 情报系统获得了网络中每段信道安全可靠性的概率，情报员 A 决定选
择一条安全性最高，即概率最大的信道路径进行发送情报。
你能帮情报员 A 找到这条信道路径吗？

输入

第一行： T 表示以下有 T 组测试数据 （ 1≤T ≤8 ）
对每组测试数据：
第一行： n m 分别表示网络中的节点数和信道数 (1<=n<=10000，1<=m<=50000)
接下来有 m 行， 每行包含三个整数 i，j，p，表示节点 i 与节点 j 之间有一条信道，其信
道安全可靠性的概率为 p%。 ( 1<=i, j<=n 1<=p<=100)

输出

每组测试数据，输出占一行，一个实数 即情报传送到达节点 n 的最高概率，精确到小数点后
6 位。

样例输入 Copy

1
5 7
5 2 100
3 5 80
2 3 70
2 1 50
3 4 90
4 1 85
3 1 70

样例输出 Copy

61.200000

先来一个 超时 and 超内存 的 Floyd 算法，深刻理解一下本题的含义

#include<stdio.h>
double e[10010][10010];
int main()
{
    int t;scanf("%d",&t);
    while(t--)
    {
    	int n,m;scanf("%d%d",&n,&m);
        for(int i=1;i<=n;i++) 
            for(int j=1;j<=n;j++)
            {
                if(i==j) e[i][j]=1;
                else e[i][j]=0;
            }
        while(m--)
        {
            int a,b;double c;scanf("%d%d%lf",&a,&b,&c);
            e[a][b]=e[b][a]=c/100;
        }
        for(int k=1;k<=n;k++) 
            for(int i=1;i<=n;i++) 
                for(int j=1;j<=n;j++)
                {
                    if(e[i][j]<e[i][k]*e[k][j])
                    e[i][j]=e[i][k]*e[k][j];
                }
        printf("%.6f\n",100*e[1][n]);
    }
    return 0; 
}

堆优化版的Dijkstra算法

#include<cstdio>
#include<cstring>
#include<queue>
using namespace std;
typedef pair<double,int> PDI; 
const int N = 10010, M = 50010;
int n,m,idx;
int h[N],e[M],w[M],ne[M];
double dist[N];
bool st[N];
void add(int a,int b,int c)  // 添加一条边a->b，边权为c
{
    e[idx]=b,w[idx]=c,ne[idx]=h[a],h[a]=idx++;
}
void dijkstra(int s)  //求 s 号点到  n 号点的最短路距离
{
    memset(dist,0,sizeof dist);
    dist[s]=1; 		//自身到自身的概率为1 
    priority_queue<PDI> heap;	// 默认是大根堆 
    heap.push((PDI){1.0,s});		// first存储距离，second存储节点编号
 
    while(!heap.empty())
    {
        PDI t=heap.top();heap.pop();
        
        int ver=t.second;
        if(st[ver]) continue;
        st[ver]=true;
 
        for(int i=h[ver];~i;i=ne[i])
        {
            int j=e[i];
            if(dist[j]<dist[ver]*w[i]/100)
            {
                dist[j]=dist[ver]*w[i]/100;
                heap.push((PDI){dist[j],j});
            }
        }
    }
}
int main()
{
    int t;scanf("%d",&t);
    while(t--)
    {
        memset(h,-1,sizeof h);
        memset(st,false,sizeof st);
        idx=0;
        scanf("%d%d",&n,&m);
        while(m--)
        {
            int a,b,c;scanf("%d%d%d",&a,&b,&c);
            add(a,b,c);add(b,a,c);
        }
        dijkstra(1);
        printf("%.6f\n",dist[n]*100);
    }
    return 0;
}

Spfa，链式前向星

#include<cstdio>
#include<cstring>
#include<queue>
using namespace std;
 
const int N = 10010, M = 50010;
 
int n,m,idx;
int h[N],e[M],w[M],ne[M];
double dist[N];
bool st[N];
 
void add(int a,int b,int c)  // 添加一条边a->b，边权为c
{
    e[idx]=b,w[idx]=c,ne[idx]=h[a],h[a]=idx++;
}
 
void spfa(int s)
{
    queue<int>Q;
    memset(dist,0,sizeof dist);
    dist[s]=1;st[s]=true;
    Q.push(s);
    
    while(!Q.empty())
	{
        int t=Q.front();Q.pop();
        st[t]=false;//从队列中取出来之后该节点st被标记为false,代表之后该节点如果发生更新可再次入队
        for(int i=h[t];i!=-1;i=ne[i])
		{
            int j=e[i];
            if(dist[j]<dist[t]*w[i]/100)
			{
                dist[j]=dist[t]*w[i]/100;
                if(!st[j])//当前已经加入队列的结点，无需再次加入队列，即便发生了更新也只用更新数值即可，重复添加降低效率
				{
					Q.push(j);
                    st[j]=true;
                }
            }
        }
    }
}
int main()
{
    int t;scanf("%d",&t);
    while(t--)
    {
        memset(h,-1,sizeof h);
        memset(st,false,sizeof st);
        idx=0;
        scanf("%d%d",&n,&m);
        while(m--)
        {
            int a,b,c;scanf("%d%d%d",&a,&b,&c);
            add(a,b,c);add(b,a,c);
        }
        spfa(1);				
        printf("%.6f\n",dist[n]*100);
    }
    return 0;
}