城市间紧急救援 解题报告

https://blog.csdn.net/NCC__dapeng/article/details/84314862

内容都在这篇我之前写的博客里。

这里说一下比赛时WA的地方，就是在计算种类的过程出现了错误，下面给出WA代码与AC代码以便更好的分析。

WA代码：

#include <bits/stdc++.h>
using namespace std;
typedef long long ll;
const int INF=0x3f3f3f3f;
const int maxn=500+10;
int n,m,s,D,res;
int num[maxn];
int cost[maxn][maxn],pre[maxn];
int d[maxn],dn[maxn],zl[maxn];
bool used[maxn];

void init()
{
    scanf("%d %d %d %d",&n,&m,&s,&D);

    res=1;
    memset(num,0,sizeof(num));
    memset(cost,INF,sizeof(cost));
    memset(pre,-1,sizeof(pre));
    memset(used,0,sizeof(used));
    memset(d,INF,sizeof(d));
    memset(dn,0,sizeof(dn));
    memset(zl,0,sizeof(zl));
    for(int i=0;i<n;i++) scanf("%d",&num[i]),dn[i]=num[i];
    for(int i=0;i<m;i++)
    {
        int a,b,c; scanf("%d %d %d",&a,&b,&c);
        cost[a][b]=cost[b][a]=min(cost[a][b],c);
    }
}

void Djk()
{
    d[s]=0;
    zl[s]=1;
    while(true)
    {
        int v=-1;
        for(int u=0;u<n;u++)
        {
            if(!used[u]&&(v==-1||d[u]<d[v])) v=u;
        }

        if(v==-1) break;
        used[v]=true;

        for(int u=0;u<n;u++)
        {
            if(d[u]==d[v]+cost[v][u])
            {
                zl[u]++;
                if(dn[u]<dn[v]+num[u])
                {
                    dn[u]=dn[v]+num[u];
                    pre[u]=v;
                }
            }
            else if(d[u]>d[v]+cost[v][u])
            {
                zl[u]=1;
                dn[u]=dn[v]+num[u];
                d[u]=d[v]+cost[v][u];
                pre[u]=v;
            }
        }
    }



}

void print()
{
    vector<int> ans;
    for(int t=D;t!=-1;t=pre[t]) ans.push_back(t);
    reverse(ans.begin(),ans.end());

    for(int i=0;i<ans.size();i++)
    {
        res*=zl[ans[i]];
    }

    printf("%d %d\n",res,dn[D]);
    for(int i=0;i<ans.size();i++)
    {
        if(i==0) printf("%d",ans[i]);
        else printf(" %d",ans[i]);
    }
    printf("\n");
}


int main()
{
    init();
    Djk();
    print();

    return 0;
}

AC代码：

#include <bits/stdc++.h>
using namespace std;
typedef long long ll;
const int INF=0x3f3f3f3f;
const int maxn=500+10;
int n,m,s,D,res;
int num[maxn];
int cost[maxn][maxn],pre[maxn];
int d[maxn],dn[maxn],zl[maxn];
bool used[maxn];

void init()
{
    scanf("%d %d %d %d",&n,&m,&s,&D);

    res=1;
    memset(num,0,sizeof(num));
    memset(cost,INF,sizeof(cost));
    memset(pre,-1,sizeof(pre));
    memset(used,0,sizeof(used));
    memset(d,INF,sizeof(d));
    memset(dn,0,sizeof(dn));
    memset(zl,0,sizeof(zl));
    for(int i=0;i<n;i++) scanf("%d",&num[i]),dn[i]=num[i];
    for(int i=0;i<m;i++)
    {
        int a,b,c; scanf("%d %d %d",&a,&b,&c);
        cost[a][b]=cost[b][a]=min(cost[a][b],c);
    }
}

void Djk()
{
    d[s]=0;
    zl[s]=1;
    while(true)
    {
        int v=-1;
        for(int u=0;u<n;u++)
        {
            if(!used[u]&&(v==-1||d[u]<d[v])) v=u;
        }

        if(v==-1) break;
        used[v]=true;

        for(int u=0;u<n;u++)
        {
            if(d[u]==d[v]+cost[v][u])
            {
                zl[u]+=zl[v];
                if(dn[u]<dn[v]+num[u])
                {
                    dn[u]=dn[v]+num[u];
                    pre[u]=v;
                }
            }
            else if(d[u]>d[v]+cost[v][u])
            {
                zl[u]=zl[v];
                dn[u]=dn[v]+num[u];
                d[u]=d[v]+cost[v][u];
                pre[u]=v;
            }
        }
    }



}

void print()
{
    vector<int> ans;
    for(int t=D;t!=-1;t=pre[t]) ans.push_back(t);
    reverse(ans.begin(),ans.end());


    printf("%d %d\n",zl[D],dn[D]);
    for(int i=0;i<ans.size();i++)
    {
        if(i==0) printf("%d",ans[i]);
        else printf(" %d",ans[i]);
    }
    printf("\n");
}


int main()
{
    init();
    Djk();
    print();

    return 0;
}

这两种计数方式看似没什么区别，只不过第一种是最后需要处理，而最后一种直接进行求和了而已，但是我们注意第一种的最后的处理方法，他计算的仅仅是最短路劲的点的种类，有一种可能，那就是“殊途同归”想象一下，当所有相同的最短路径，只有终点起点相同，并且其中一条最短路径的一个点出现了2种到达方式，但这并不是最佳的最短路径时，第一种计算方式就会发生严重的错误！