hdu 1385 Minimum Transport Cost

/*用floyd算出每两点间的最短路，用一个theta矩阵记录由 i 到 j 的中间节点中编号最高的一个，据此可以不断缩小寻找的范围，找到theta[i][j]=0时，就直接输出i-->j。print函数和path函数是类似的，就是找由I到j的路径。至于字典序的前后，只要把两条路径经过的节点用path（）找出来，再比较一下即可。*/

#include<iostream>
using namespace std;
    int i,j,k,map[1001][1001],tax[1001],theta[1001][1001],t;
    int ch1[1001],ch2[1001];
int max(int a,int b,int c)
{
    if(a>b && a>c)
     return c;
     else if(b>c)
     return b;
     else
     return c;
}
void print(int a,int b)
{
     if(theta[a][b]==0){
         printf("%d-->",a);
         return ;
     }
     print(a,theta[a][b]);
     print(theta[a][b],b);
}

void path(int i,int j)
{
     if(theta[i][j]==0)
     {
         ch1[t++]=i;
         return;
     }
     path(i,theta[i][j]);
     path(theta[i][j],j);
}

bool find(int i,int j,int k)
{
     int a;
     t=0;
     path(i,j);
     ch1[t++]=j;
     for(a=0;a<t;a++)
     ch2[a]=ch1[a];
     int t2=t;
    
     t=0;
     path(i,k);
     ch1[t++]=k;
     int t1=t;
    
     for(a=0;;a++){
                   if(a>t1)
                   return true;
                   if(a>t2)
                   return false;
                   if(ch1[a]>ch2[a])
                   return false;
                   if(ch1[a]<ch2[a])
                   return true;
     }
}
int main()
{
    int n,a,b;
    while(scanf("%d",&n)&& n!=0)
    {
        for(i=1;i<=n;i++)
            for(j=1;j<=n;j++)
                scanf("%d",&map[i][j]);
        memset(theta,0,sizeof(theta));
        for(i=1;i<=n;i++)
        scanf("%d",&tax[i]);
        for(k=1;k<=n;k++)
            for(i=1;i<=n;i++)
                for(j=1;j<=n;j++)
                    if(map[i][k]!=-1&&map[k][j]!=-1 && i!=k && j!=k)
                    {
                        int s=map[i][k]+map[k][j]+tax[k];
                        if(map[i][j]==-1 || map[i][j]>s
                        || map[i][j]==s && find(i,j,k)){//题目没说是来回路径不一样，但实际上有这种情况，调了很久。
                             map[i][j]=s;
                            theta[i][j]=max(theta[i][k],theta[k][j],k);
                            }
                    }
        while(scanf("%d%d",&a,&b)&&a>=0)
        {
            printf("From %d to %d :/nPath: ",a,b);
            if(a!=b)
            print(a,b);
            printf("%d/nTotal cost : %d/n/n",b,map[a][b]);
        }
    }
    return 0;
}