整理一下昨晚做的一道不难但是让我纠结好久的题目。
/*
zoj_1655 最短路
简单dijkstra
最近两天感觉状态实在是烂。。老是理解错题目。。。
方法:想到每个城市的货物都单独考虑,然后只需计算每个城市到captial的最大比率乘积。用源点为captial
的dijk很容易就出来。
注意:两点之间路径可能多条(这个在以后做题时不管题目有没有说最好都判断一下)
让我纠结的是100个点的比率乘积精度不会丢失吗?所以我觉得上面那种方法总感觉不怎么安全。。但是换别的
思路写起来又实在不简单。百度了一下,发现几乎所有人用的都是1-rate 的乘积。个人觉得这是因为rate数据
基本都小于0.5才没有wa的。
举个例子:
5 6
10
10
10
10
1 3 0.01
1 4 0.002
2 3 0.0003
2 4 0.00004
3 5 0.000005
4 5 0.0000006
用rate的话数据出来40.00(见错误代码),用1-rate的话出来39.98
但是我想,万一存在
5 6
10
10
10
10
1 3 0.99
1 4 0.998
2 3 0.9997
2 4 0.99996
3 5 0.999995
4 5 0.9999994
这样的数据1-rate 精度不是一样丢失了吗,事实证明,这样的数据是不存在的。。
虽然过了,还是感觉没有完美解决问题。
*/
#include <iostream>
#include <cstdio>
#include <limits.h>
#include <string.h>
using namespace std;
double map[110][110];
double weight[110];
double dist[110];
void dijkstra( int sta,int n )
{
int i,j,mark;
double maxi;
bool flag[110];
memset( flag,0,sizeof(flag) );
memset( dist,0,sizeof(dist) ); //注意初始不能为1因为rate存在0的数据
flag[sta]=1;
dist[sta]=1;
for( i=1;i<=n;i++ )
{
maxi=-1;
for( j=1;j<=n;j++ )
{
if( !flag[j] )
{
if( map[sta][j]!=-1 && dist[j]<dist[sta]*map[sta][j] )
dist[j]=dist[sta]*map[sta][j];
else if( map[sta][j]!=-1 && dist[j]==0 )
dist[j]=dist[sta]*map[sta][j];
if( maxi<dist[j] && dist[j]!=0 )
{
mark=j;
maxi=dist[j];
}
}
}
if( maxi==-1 ) break;
sta=mark;
flag[sta]=1;
}
}
void output( int n )
{
int i;
double sum;
sum=0.0;
for( i=1;i<n;i++ )
{
//printf( "%lf %lf\n",dist[i],weight[i] );
sum+=dist[i]*weight[i];
}
printf( "%.2f\n",sum );
}
int main()
{
//freopen( "a.txt","r",stdin );
int n,m,a,b,i,j;
double rate;
while( scanf( "%d%d",&n,&m )!=EOF )
{
for( i=0;i<=n;i++ )
for( j=0;j<=n;j++ )
map[i][j]=-1;
for( i=1;i<n;i++ )
scanf( "%lf",&weight[i] );
for( i=1;i<=m;i++ )
{
scanf( "%d%d%lf",&a,&b,&rate );
if( 1-rate>map[a][b] )
map[b][a]=map[a][b]=1-rate;
}
dijkstra( n,n );
output( n );
}
return 0;
}
//错误代码
#include <iostream>
#include <cstdio>
#include <limits.h>
#include <string.h>
#define inf 999999.0
using namespace std;
double weight[110];
double map[110][110];
double dist[110];
double sum;
void dijkstra( int sta,int n )
{
int i,j,mark;
double mini;
bool flag[110];
memset( flag,0,sizeof(flag) );
for( i=0;i<=n;i++ )
dist[i]=1;
flag[sta]=1;
for( i=1;i<=n;i++ )
{
mini=inf;
for( j=1;j<=n;j++ )
{
if( !flag[j] )
{
if( map[sta][j]!=-1 && dist[j]>dist[sta]*map[sta][j] )
dist[j]=dist[sta]*map[sta][j];
else if( map[sta][j]!=-1 && dist[j]==1 )
dist[j]=dist[sta]*map[sta][j];
if( mini>dist[j] && dist[j]!=1 )
{
mark=j;
mini=dist[j];
}
}
}
if( mini==inf ) break;
sta=mark;
flag[sta]=1;
}
}
void output( int n )
{
int i;
for( i=1;i<n;i++ )
sum-=dist[i]*weight[i];
printf( "%.2f\n",sum );
}
int main()
{
freopen( "a.txt","r",stdin );
int n,m,a,b,i,j;
double rate;
while( scanf( "%d%d",&n,&m )!=EOF )
{
for( i=0;i<=n;i++ )
for( j=0;j<=n;j++ )
map[i][j]=-1;
sum=0;
for( i=1;i<n;i++ )
{
scanf( "%lf",&weight[i] );
sum+=weight[i];
}
for( i=1;i<=m;i++ )
{
scanf( "%d%d%lf",&a,&b,&rate );
if( map[a][b]==-1 || map[a][b]>1-rate )
{
map[a][b]=rate;
map[b][a]=rate;
}
}
dijkstra( n,n );
output( n );
}
return 0;
}