第九届河南省程序设计大赛
C 信道安全
内存限制:64MB 时间限制:1s Special Judge: No
题目描述:
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
此道题我第一次用的是dfs,当时忘了考虑时间的问题,代码就是下边的代码,我居然考虑到了空间的问题,没有使用邻接矩阵,而是使用的邻接表(正确AC的代码在最下方)
#include<iostream>
#include<cstdio>
#include<algorithm>
#include<cstring>
#include<string>
using namespace std;
int t,n,m;
int book[10006]={0};
int x[100010],y[100010];
double w[100010];
int first1[10006];
int next1[100010];
double sum=0;
void dfs(int i,double sum1)
{
if(sum1<=sum)
return ;
if(i==n){
sum=sum1;
return ;
}
int k1=first1[i];
while(k1!=-1)
{
int y1=y[k1];
if(book[y1]==0){
book[y1]=1;
dfs(y1,sum1*w[k1]*0.01);
book[y1]=0;
}
k1=next1[k1];
}
return ;
}
int main()
{
cin>>t;
while(t--)
{
cin>>n>>m;
memset(book,0,sizeof(book));
memset(first1,-1,sizeof(first1));
memset(next1,-1,sizeof(next1));
for(int i=1;i<=m;i++){
cin>>x[i]>>y[i]>>w[i];
next1[i]=first1[x[i]];
first1[x[i]]=i;
}
for(int i=m+1;i<=2*m;i++){
w[i]=w[i-m];
y[i]=x[i-m];
x[i]=y[i-m];
next1[i]=first1[x[i]];
first1[x[i]]=i;
}
int k=first1[1];
book[1]=1;
sum=0;
while(k!=-1)
{
book[y[k]]=1;
dfs(y[k],0.01*w[k]);
book[y[k]]=0;
k=next1[k];
}
printf("%.6f\n",sum*100);
}
return 0;
}
在nyoj以及郑轻oj上都超时了,我才考虑到应该换一种思路,然后我就使用了Bellman-Ford算法,最后提交也就正确了
代码如下:
#include<iostream>
#include<cstdio>
using namespace std;
int t,n,m,k;
int x[100010],y[100010];
double w[100010];
double dis[10009];
int main()
{
cin>>t;
while(t--)
{
cin>>n>>m;
for(int i=0;i<n+5;i++)
dis[i]=-1;
for(int i=1;i<=m;i++){
cin>>x[i]>>y[i]>>w[i];
if(x[i]==1)
dis[y[i]]=w[i]/100;
}
for(int i=m+1;i<=2*m;i++){
w[i]=w[i-m];
y[i]=x[i-m];
x[i]=y[i-m];
if(x[i]==1)
dis[y[i]]=w[i]/100;
}
dis[1]=0;
for(int i=0;i<n-1;i++){
k=0;
for(int j=1;j<=2*m;j++){
if(dis[y[j]]<dis[x[j]]*w[j]/100){
k=1;
dis[y[j]]=dis[x[j]]*w[j]/100;
}
}
if(k==0)
break;
}
printf("%.6f\n",dis[n]*100);
}
return 0;
}