hdoj 1596 find the safest road

题目传送：http://acm.hdu.edu.cn/showproblem.php?pid=1596

分析：Dijkstra变体，最短路径判断计算方式：Safe(P) = s(e1)*s(e2)…*s(ek) e1,e2,ek是P 上的边。 判断方式：lowcost[pre]*cost[pre][j]>lowcost[j]

代码：

//2013-10-30 10:54:06    Accepted    1596    1484MS    9788K    1836 B    C++    空信高手

#include <iostream>

#include <iomanip>

using namespace std;



/*==================================================*\

 | Dijkstra 数组实现O （N^2 ）

| Dijkstra --- 数组实现( 在此基础上可直接改为STL 的Queue实现)

 | lowcost[] --- beg 到其他点的最近距离

| path[] -- beg为根展开的树，记录父亲结点

\*==================================================*/

const int N=1100;

int path[N],vis[N];

double cost[N][N],lowcost[N];

void Dijkstra(int n,int beg)

{

    int i,j;

    double max;

    memset(vis,0,sizeof(vis));

    vis[beg]=1;

    for(i=0; i<n; i++)

    {

        lowcost[i]=cost[beg][i];

        path[i]=beg;

    }

    lowcost[beg]=0;

    path[beg]=-1;

    int pre=beg;

    for(i=1; i<n; i++)

    {

        max=0;

        for(j=0; j<n; j++)

            //下面的加法可能导致溢出，INF不能取太大

            if(vis[j]==0&&lowcost[pre]*cost[pre][j]>lowcost[j])

            {

                lowcost[j]=lowcost[pre]*cost[pre][j];

                path[j]=pre;

            }

        for(j=0; j<n; j++)

            if(vis[j]==0&&lowcost[j]>max)

            {

                max=lowcost[j];

                pre=j;

            }

        vis[pre]=1;

    }

}



int main()

{

 //   freopen("input.txt","r",stdin);

    int n,m,i,j,a,b;

    while(scanf("%d",&n)!=EOF&&n)

    {

        memset(cost,0,sizeof(cost));

        memset(lowcost,0,sizeof(lowcost));

        for(i=0; i<n; i++)

            for(j=0; j<n; j++)

                cin>>cost[i][j];

        cin>>m;

        while(m--)

        {

            cin>>a>>b;

            Dijkstra(n,a-1);

            if(lowcost[b-1]!=0)cout<<setiosflags(ios::fixed)<<setprecision(3)<<lowcost[b-1]<<endl;

            else cout<<"What a pity!"<<endl;

        }

    }

    return 1;

}