算法分析 千里走单骑

民间流行着有关三国时期的许多故事，其中最有名的就数千里走单骑。关羽保护着两位嫂嫂，从许都出发，过五关，斩六将，终于与刘备、张飞重聚于古城。

现在，你的任务是替关羽寻找一条新的路线，以找到他那比兔子跑得还快的哥哥。假定关羽知道从许都到古城各关卡间的所有道路、以及每个关卡平安通过的概率。请帮助关羽找到一条具有最大概率的有效路线。注意每个关卡他最多只能经过一次。

输入：

第一行有一个整数 T (1 <= T <= 10) ，表示测试情况的个数。

对于每种情况，第一行只有一个整数， N (3 <= N <= 1000) 。接下来的 N 行是一个矩阵，表示从许都到古城的 N – 2 个关卡之间的道路。如果第 i 行第 j 列处是 1 ，则从 i 到 j 有一条道路；而 0 表示从 i 到 j 没有道路。每种情况的最后一行包含 N – 2 个实数，其间以一个空格分开，分别表示关羽平安通过各关卡的概率 P (0 <= P <= 1) 。

输出：由许都出发平安到达古城的最大概率，要求精确到四位小数。如果结果小于 0.0001 ，则输出“ Cannot reach!”。

比如输出Case 1: Cannot reach!

#include<iostream>  
using namespace std;  
int m;  
double road[1001][1001];  
double pro[1001];  
bool mark[1001];  
double Pro[1001];  
void dijkstra()  
{  
for(int i=0;i<m;i++)  
{  
    mark[i] = 0;  
    Pro[i] = 0;  
}  
Pro[0] = 1;  
int vir = 0;  
int max;  
for (int i = 0;i < m;i++)  
{  
    for (int j = 0;j < m;j++)  
        if (!mark[j] && Pro[j] > max)  
        {  
            max = Pro[j];  
            vir = j;  
        }  
    mark[vir] = 1;  
    for (int j = 0;j < m;j++)  
    {  
        if (!mark[j] && Pro[j] < Pro[i] * road[i][j])  
            Pro[j] = Pro[i] * road[i][j];  
    }  
}  
}  
int main()  
{  
    int n;  
    cin>>n;  
    for(int k=1;k<=n;k++)  
    {  
        cin >> m;  
        for (int i = 0;i < m;i++)  
            for (int j = 0;j < m;j++)  
                cin >> road[i][j];  
        for (int i = 1;i < m - 1;i++)  
            cin >> pro[i];  
        pro[m - 1] = 1;  
        for (int i = 0;i < m;i++)  
            for (int j = 0;j < m;j++)  
                if (road[i][j] == 1)  
                    road[i][j] = pro[j];      
        dijkstra();  
        if(Pro[m-1]>=0.0001)  
        {  
        cout<<"Case "<<k<<": ";  
        printf("%.4f\n",Pro[m-1]);  
        }  
        else  
        cout<<"Case "<<k<<": Cannot reach!"<<endl;  
    }  
}