BZOJ2763[JLOI2011]飞行路线 [分层图最短路]

2763: [JLOI2011]飞行路线

Time Limit: 10 Sec Memory Limit: 128 MB
Submit: 2523 Solved: 946
[Submit][Status][Discuss]

Description

Alice和Bob现在要乘飞机旅行，他们选择了一家相对便宜的航空公司。该航空公司一共在n个城市设有业务，设这些城市分别标记为0到n-1，一共有m种航线，每种航线连接两个城市，并且航线有一定的价格。Alice和Bob现在要从一个城市沿着航线到达另一个城市，途中可以进行转机。航空公司对他们这次旅行也推出优惠，他们可以免费在最多k种航线上搭乘飞机。那么Alice和Bob这次出行最少花费多少？

Input

数据的第一行有三个整数，n,m,k，分别表示城市数，航线数和免费乘坐次数。

第二行有两个整数，s,t，分别表示他们出行的起点城市编号和终点城市编号。(0<=s,t<n)

接下来有m行，每行三个整数，a,b,c，表示存在一种航线，能从城市a到达城市b，或从城市b到达城市a，价格为c。(0<=a,b<n,a与b不相等，0<=c<=1000)

Output

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

Sample Input

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

Sample Output

HINT

对于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.

d[i][j]表示到节点i用了j次免费的最小花费

//

//  main.cpp

//  bzoj2763

//

//  Created by Candy on 9/26/16.

//  Copyright © 2016 Candy. All rights reserved.

//

#include <iostream>

#include <cstdio>

#include <algorithm>

#include <queue>

#include <cstring>

using namespace std;

const int N=1e4+,M=5e4+,K=;

inline int read(){

    char c=getchar();int x=,f=;

    while(c<''||c>''){if(c=='-')f=-;c=getchar();}

    while(c>=''&&c<=''){x=x*+c-'';c=getchar();}

    return x*f;

}

int n,m,k,s,t,a,b,c;

struct edge{

    int v,w,ne;

}e[M<<];

int h[N],cnt=;

inline void ins(int u,int v,int w){

    cnt++;

    e[cnt].v=v;e[cnt].w=w;e[cnt].ne=h[u];h[u]=cnt;

    cnt++;

    e[cnt].v=u;e[cnt].w=w;e[cnt].ne=h[v];h[v]=cnt;

}

struct hn{

    int u,d,f;

    hn(int a=,int b=,int c=):u(a),d(b),f(c){}

    bool operator<(const hn &rhs)const{return d>rhs.d;}

};

int d[N][K],done[N][K];

void bfs(){

    memset(d,,sizeof(d));

    priority_queue<hn> q;

    q.push(hn(s,,));

    d[s][]=;

    while(!q.empty()){

        hn x=q.top();q.pop();

        int u=x.u,dis=x.d,f=x.f;

        if(done[u][f]) continue;

        done[u][f]=;

        if(u==t){printf("%d",dis);return;}

        for(int i=h[u];i;i=e[i].ne){

            int v=e[i].v,w=e[i].w;

            if(f<k&&!done[v][f+]&&d[v][f+]>dis){

                d[v][f+]=dis;

                q.push(hn(v,dis,f+));

            }

            if(!done[v][f]&&d[v][f]>dis+w){

                d[v][f]=dis+w;

                q.push(hn(v,d[v][f],f));

            }

        }

    }

}

int main(int argc, const char * argv[]) {

    n=read();m=read();k=read();s=read();t=read();

    for(int i=;i<=m;i++){a=read();b=read();c=read();ins(a,b,c);}

    bfs();

    return ;

}

巴特西