洛谷P4568 [JLOI2011]飞行路线（分层图最短路）

题目描述

Alice和Bob现在要乘飞机旅行，他们选择了一家相对便宜的航空公司。该航空公司一共在nn个城市设有业务，设这些城市分别标记为00到n-1n−1，一共有mm种航线，每种航线连接两个城市，并且航线有一定的价格。

Alice和Bob现在要从一个城市沿着航线到达另一个城市，途中可以进行转机。航空公司对他们这次旅行也推出优惠，他们可以免费在最多kk种航线上搭乘飞机。那么Alice和Bob这次出行最少花费多少？

输入格式

数据的第一行有三个整数，n,m,kn,m,k，分别表示城市数，航线数和免费乘坐次数。
第二行有两个整数，s,ts,t，分别表示他们出行的起点城市编号和终点城市编号。
接下来有m行，每行三个整数，a,b,ca,b,c，表示存在一种航线，能从城市aa到达城市bb，或从城市bb到达城市aa，价格为cc。

输出格式

只有一行，包含一个整数，为最少花费。

输入输出样例

输入 #1复制

5 6 1
0 4
0 1 5
1 2 5
2 3 5
3 4 5
2 3 3
0 2 100

输出 #1复制

8

说明/提示

对于30%的数据,2≤n≤50,1≤m≤300,k=0;
对于50%的数据,2≤n≤600,1≤m≤6000,0≤k≤1;
对于100%的数据,2≤n≤10000,1≤m≤50000,0≤k≤10,0≤s,t<n,0≤a,b<n,a≠b,0≤c≤1000

2018.12.10 增加一组 hack 数据

题解

今天学了一个神奇的最短路DP，于是洛谷一搜，名曰分层图最短路

即：在一张图上，有k次机会可以通过一条边而不需要计算权值（免费过路），求从起点到终点的最短路径

大概意思就是 dis[u][use]表示用了use次免费机会，到达节点u的最短路

对于路径u->v

1.如果我们用了免费机会
dis[v][use] = min{dis[u][use - 1]}

即u->v的路径上使用一次免费机会，这条路的边权就不要加了

(程序里就是dis[v][use+1]>dis[u][use]就更新答案)
2.如果我们没用免费机会
dis[v][use] = min{dis[u][use] + w[u->v]}

然后考场上杯具了：dijkstra写挂了

#include<queue>
#include<cstdio>
#include<cstring>
#include<iostream>
#include<algorithm>
using namespace std;
const int N=100005;
const int oo=0x3f3f3f3f;
int n,m,k,S,T;
inline int read()
{
    int s=0,f=1;char ch=getchar();
    while(ch<'0'||ch>'9'){if(ch=='-')f=-1;ch=getchar();}
    while(ch>='0'&&ch<='9'){s=s*10+ch-'0';ch=getchar();}
    return s*f;
}
int head[N],nxt[N],adj[N],tot,w[N];
inline void add(int u,int v,int ww)
{
    tot++;
    nxt[tot]=head[u];
    head[u]=tot;
    adj[tot]=v;
    w[tot]=ww;
}
int dis[N][15],vis[N][15];
struct node
{
    int u;
    int d;
    int use;//用了几次免费的线路 
};
bool operator<(node p,node q)
{
    return p.d>q.d;
}
void dijkstra(int s)
{
    int i;
    memset(dis,0x3f,sizeof(dis));
    priority_queue<node>q;
    q.push((node){s,0,0});
    dis[s][0]=0;
    while(!q.empty())
    {
        int u=q.top().u,use=q.top().use,d=q.top().d;
        q.pop();
        if(vis[u][use])continue;
        vis[u][use]=1;
        for(i=head[u];i;i=nxt[i])
        {
            int v=adj[i];
            if(!vis[v][use]&&dis[u][use]+w[i]<dis[v][use])//不使用免费次数 
            {
                dis[v][use]=dis[u][use]+w[i];
                q.push((node){v,dis[v][use],use});
            }
            if(!vis[v][use+1]&&use+1<=k&&dis[v][use+1]>dis[u][use])//使用免费次数，但要记得判断use+1<=k 
            {
                dis[v][use+1]=dis[u][use];
                q.push((node){v,dis[v][use+1],use+1});
            }
        }
    } 
}
inline int mn(int x,int y)
{
    return x<y?x:y;
}
int main()
{
    freopen("flight.in","r",stdin);
    freopen("flight.out","w",stdout);
    int i,x,y,z;
    n=read();m=read();k=read();
    S=read()+1;T=read()+1;
    for(i=1;i<=m;i++)
    {
        x=read();y=read();z=read();
        x++;y++;
        add(x,y,z);add(y,x,z);
    }
    dijkstra(S);
    int ans=oo;
    for(i=0;i<=k;i++)
        ans=mn(ans,dis[T][i]);
    printf("%d\n",ans);
    return 0;
}

最后强调一句：SPFA真的会死人的--不管花了大力气加了循环队列优化啦，SLF优化啦，都麻烦而没有什么用