NOIP模拟赛 bzoj2934业务（好题）

题解：

这道题是一道最短路，但光凭最短路是不行的，还需要有一定的dp思想

相信Floyd这个算法对于各位来说都并不陌生，而这道题所用到的dp思想

也正与Floyd中的类似，事实上，至今为止，我见过的图论中能用的dp似

乎也就这个

我们定义状态ans[i][j]表示从i到j的最小费用

dis[i][0]表示从k出发到j，费用最大的城市是k的最小费用

dis[i][1]表示从k出发到j，费用最次大的城市是k的最小费用

这里其实是省略了一维k的，详情见代码

然后是转移：

dis[e][0]=min(dis[s][0]+arr[s][e]); (cost[s]<=cost[k])

dis[e][1]=min(dis[s][0]+arr[s][e]); (cost[s]>cost[k])

dis[e][1]=min(dis[s][1]+arr[s][e]); (cost[s]<=cost[k])

处理出这个后，ans就可以转移了

ans[i][j]=min(dis[i][0]+dis[j][1],dis[i][1]+dis[j][0])+cost[k];

code：

#include<iostream>
#include<cstdio>
#include<cstring>
#include<cmath>
#include<cstdlib>
#include<algorithm>
#include<queue>
 
const int MAXN=255;
using namespace std;
int arr[MAXN][MAXN],cost[MAXN];
int n,m,a,b,K,x,dis[MAXN][2];
int ans[MAXN][MAXN];
queue<int>Q; bool vis[MAXN];
 
inline void Read(int &Ret){
    char ch;bool flag=0;
    for(;ch=getchar(),ch<'0'||ch>'9';)if(ch=='-')flag=1;
    for(Ret=ch-'0';ch=getchar(),'0'<=ch&&ch<='9';Ret=Ret*10+ch-'0');
    flag&&(Ret=-Ret);
}
inline void spfa(int k){
    memset(dis,0x3f,sizeof dis);
    dis[k][1]=dis[k][0]=0; Q.push(k);
    while(!Q.empty())
    {
        int i=Q.front(); vis[i]=0; Q.pop();
        for(int j=1;j<=n;j++)
        {
            if(arr[i][j]==arr[0][0]) continue;
            if(cost[j]>cost[k])
            {
                if(dis[j][1]>dis[i][0]+arr[i][j])
                {
                    dis[j][1]=dis[i][0]+arr[i][j];
                    if(!vis[j]) Q.push(j),vis[j]=1;
                }
            }
            else
            {
                if(dis[j][1]>dis[i][1]+arr[i][j])
                {
                    dis[j][1]=dis[i][1]+arr[i][j];
                    if(!vis[j]) Q.push(j),vis[j]=1;
                }
                if(dis[j][0]>dis[i][0]+arr[i][j])
                {
                    dis[j][0]=dis[i][0]+arr[i][j];
                    if(!vis[j]) Q.push(j),vis[j]=1;
                }
            }
        }
    }
    return;
}
int main()
{
    Read(n); Read(m); Read(K);
    for(int i=1;i<=n;i++) Read(cost[i]);
    memset(arr,0x3f,sizeof arr);
    memset(ans,0x3f,sizeof ans);
    for(int i=1;i<=m;i++)
    {
        Read(a); Read(b); Read(x);
        arr[a][b]=min(arr[a][b],x);
        arr[b][a]=arr[a][b];
    }
    for(int k=1;k<=n;k++)
    {
        spfa(k);
        for(int i=1;i<=n;i++)
        for(int j=1;j<=n;j++)
        {
            int tmp=min(dis[i][0]+dis[j][1],dis[j][0]+dis[i][1]);
            tmp+=cost[k]; ans[i][j]=min(ans[i][j],tmp);
        }
    }
    for(int i=1;i<=K;i++)
    {
        Read(a); Read(b);
        printf("%d",ans[a][b]);
        putchar(10);
    }
    return 0;
}