信道安全
时间限制:1000 ms | 内存限制:65535 KB
难度:2
描述
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 位。
样例输入
1
5 7
5 2 100
3 5 80
2 3 70
2 1 50
3 4 90
4 1 85
3 1 70
样例输出
61.200000
解题思路:一个简单最短路算法,只要把我们以前那种找最短的换成找最长的就可以了,最开始是想直接用dijkstra写的结果内存超了,所以选择用了邻接链表的SPFA最短路径算法。下面我将给出两种SPFA算法,一个是自己写的链表,一个是借用STL中vector。
自己写的邻接链表:
#include <iostream>
#include <bits/stdc++.h>
#define INF 0x3f3f3f3f
using namespace std;
struct Node
{
int v,next;//v用来存起始的顶点,next用来存下一条边,w存权重;
double w;
};
Node edge[50005];
int n,m,num;
double dist[10005];
int vis[10005],head[10005];//head是用来看是否当前链表是否遍历结束
void add_edge(int u,int v,double w)//添加边构造邻接链表
{
edge[num].v=v;
edge[num].w=w;
edge[num].next=head[u];
head[u]=num++;
}
void SPFA()
{
int u,v;//v代表当前节点。u代表当前节点的前一个节点
queue<int> Q;//初始化队列,借助它来用来保存待优化的结点
for(int i=1;i<=n;i++)
dist[i]=-INF;
memset(vis,0,sizeof(vis));
dist[1]=1.0;//保存最安全的路径
vis[1]=1;
Q.push(1);//起始节点入队
while(!Q.empty())
{
u=Q.front();Q.pop();//取出队首节点
vis[u]=0;//标明此节点没有在队里出现过
for(int i=head[u];i!=-1;i=edge[i].next)//对于队首节点相连的节点进行松弛操作
{
v=edge[i].v;
double w = edge[i].w;//当前节点的权重
if(dist[v]<dist[u]*w)//如果当前节点的权重与前一个节点的最安全度之积大于当前节点的最安全度就更新当前节点的最安全度
{
dist[v]=dist[u]*w;
if(vis[v]==0)//如果这个节点没有在队列中就加入队列
{
vis[v]=1;
Q.push(v);
}
}
}
}
}
int main()
{
int T;scanf("%d",&T);
while(T--)
{
num=0;
int x,y;
double p;
memset(head,-1,sizeof(head));
scanf("%d %d",&n,&m);
for(int i=0;i<m;i++)
{
scanf("%d %d %lf",&x,&y,&p);
add_edge(x,y,p/100);//无向图,所以两边都要加
add_edge(y,x,p/100);
}
SPFA();
printf("%lf\n",dist[n]*100);
}
return 0;
}
借用STL中vector的SPFA
#include <iostream>
#include <bits/stdc++.h>
#define INF 0x3f3f3f3f
using namespace std;
struct Node
{
int v;
double w;
};
vector<Node> edge[50005];
int n,m,num;
double dist[10005];
int vis[10005],head[10005];
void SPFA()
{
int u,v;
queue<int> Q;
for(int i=1;i<=n;i++)
dist[i]=-INF;
memset(vis,0,sizeof(vis));
dist[1]=1.0;
vis[1]=1;
Q.push(1);
while(!Q.empty())
{
u=Q.front();Q.pop();
for(int i=0;i<edge[u].size();i++)
{
v=edge[u][i].v;
double w = edge[u][i].w;
if(dist[v]<dist[u]*w)
{
dist[v]=dist[u]*w;
if(vis[v]==0)
{
vis[v]=1;
Q.push(v);
}
}
}
vis[u]=0;
}
}
int main()
{
int T;scanf("%d",&T);
while(T--)
{
num=0;
int x,y;
Node t;
double p;
scanf("%d %d",&n,&m);
for(int i = 1; i <=n;i++ )
{
edge[i].clear();
}
for(int i=0;i<m;i++)
{
scanf("%d %d %lf",&x,&y,&p);
t.v=x;t.w=p/100;
edge[y].push_back(t);
t.v=y;t.w=p/100;
edge[x].push_back(t);
}
SPFA();
printf("%lf\n",dist[n]*100);
}
return 0;
}