HDU ——1385 FLOYD

题意就是在给定你n个城市之间的路费和税费后，问你任意两点的费用，然后输出路径，就酱。

从我的描述大概也能猜出来是使用Floyd算法吧，因为dij算法是求单源最短路，而floyd是求多源最短路，在这里并不是说dij算法不能用，可以，只是floyd更方便而已，在后面我会放出自己dij算法超时的代码，大佬有兴趣帮忙优化一下？QAQ

如果对floyd算法不太了解的话，在这里做简单的描绘，其实floyd算法是对dij算法的一种扩展，因为dij算法是以其中一个点为固定点，然后找他到剩余所有点的最近距离，而floyd算法则是通过第三个点来作为中介去求的i到j的最近距离，也是通过k的不断变化来更新i到j的最近距离

下面放上floyd的代码：

#include<cstdio>
#include<cstring>

using namespace std;

const int maxn=110;
const int inf=0x3f3f3f3f;
int map[maxn][maxn];
int path[maxn][maxn];
int dis[maxn][maxn];
int tax[maxn];
int n;

void floyd()
{
    for(int k=1;k<=n;++k)
    {
        for(int i=1;i<=n;++i)
        {
            for(int j=1;j<=n;++j)
            {
                if(w[i][k]!=inf&&w[k][j]!=inf)
                {
                    int temp=dis[i][k]+dis[k][j]+tax[k];
                    if(w[i][j]>temp)
                    {
                        dis[i][j]=temp;
                        path[i][j]=path[i][k];
                    }
                    else if(w[i][j]==temp&&path[i][j]>path[i][k])
                        path[i][j]=path[i][k];
                }
            }
        }
    }
}

void ini()
{
    for(int i=0;i<=n;i++)
    {
        for(int j=0;j<=n;j++)
        {
            path[i][j]=j;
        }
    }
    memset(dis,0x3f,sizeof(dis));
}

int main()
{
    while(scanf("%d",&n)&&n)
    {
        ini();
        for(int i=1;i<=n;i++)
        {
            for(int j=1;j<=n;j++)
            {
                scanf("%d",&dis[i][j]);
                if(dis[i][j]==-1) dis[i][j]=inf;
            }
        }
        for(int i=1;i<=n;i++)
            scanf("%d",&tax[i]);
        floyd();
        int a,b;
        while(scanf("%d %d",&a,&b)&&a!=-1&&b!=-1)
        {
            printf("From %d to %d :\n",a,b);
            printf("Path: ");
            printf("%d",a);
            int point=a;
            while(point!=b)
            {
                printf("->%d",path[point][b]);
                point=path[point][b];
            }
            puts("");
            printf("Total cost : %d\n",dis[a][b]);
            puts("");
        }
    }
    return 0;
}

然后默默的挂上自己的超时的dij算法，

#include<cstdio>
#include<cstring>
#include<queue>

using namespace std;

const int maxn=1e6+10;
const int inf=0x3f3f3f3f;
typedef pair <int , int> pii;
int path[maxn];
int head[maxn];
int dis[maxn];
int tax[maxn];
int re[maxn];
int n,m;
int cnt=1;

struct Edge
{
    int to,next,val;
    Edge(){}
    Edge(int to,int next,int val)
    {
        to=to;
        next=next;
        val=val;
    }
}edge[maxn];

void add(int u,int v,int val)
{
    edge[cnt].to=v;
    edge[cnt].next=head[u];
    edge[cnt].val=val;
    head[u]=cnt++;
}

void dij(int st,int fa)
{
    memset(dis,0x3f,sizeof(dis));
    dis[st]=0;
    priority_queue<pii, vector<pii>, greater<pii> > que;
    //priority_queue< pii,vector<pii>, greater<pii> >que;
    que.push(make_pair(0,st));
    while(!que.empty())
    {
        pii p=que.top();
        que.pop();
        int v=p.second;
        if(dis[v]<p.first) continue;
        for(int i=head[v];~i;i=edge[i].next)
        {
            Edge e=edge[i];
            if(e.to==fa)
            {
                if(dis[e.to]>dis[v]+e.val)
                {
                    path[e.to]=v;
                    dis[e.to]=dis[v]+e.val;
                    que.push(make_pair(dis[e.to],e.to));
                }
            }
            else
            {
                if(dis[e.to]>dis[v]+e.val+tax[e.to])
                {
                    path[e.to]=v;
                    dis[e.to]=dis[v]+e.val+tax[e.to];
                    que.push(make_pair(dis[e.to],e.to));
                }
            }
        }
    }
    re[0]=fa;
    int point=1;
    printf("Path: ");
    for(int i=fa;~i;i=path[i])
    {
        if(path[i]==-1) break;
        re[point++]=path[i];
    }
    point--;
    int flag=0;
    for(int i=point;i>=0;i--)
    {
        if(flag==0)
        {
            flag=1;
            printf("%d",re[i]);
        }
        else printf("->%d",re[i]);
    }
    printf("\n");
    printf("Total cost : %d\n",dis[fa]);
    puts("");
}

void ini()
{
    memset(head,-1,sizeof(head));
    memset(tax,0,sizeof(tax));
    memset(path,-1,sizeof(path));
    memset(re,0,sizeof(re));
}

int main()
{
    while(scanf("%d",&n)&&n)
    {
        ini();
        for(int i=1;i<=n;i++)
        {
            for(int j=1;j<=n;j++)
            {
                int x;
                scanf("%d",&x);
                if(x==-1) x=inf;
                add(i,j,x);
            }
        }
        for(int i=1;i<=n;i++)
            scanf("%d",&tax[i]);
        int st,fa;
        while(scanf("%d %d",&st,&fa)&&st!=-1&&fa!=-1)
        {
            printf("From %d to %d :\n",st,fa);
            dij(st,fa);
            memset(path,-1,sizeof(path));
        }
    }
    return 0;
}