zoj 1655

整理一下昨晚做的一道不难但是让我纠结好久的题目。

/*
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;
}