BZOJ 2763 飞行路线

2763: [JLOI2011]飞行路线

Time Limit: 10 Sec  Memory Limit: 128 MB
Submit: 3113  Solved: 1168
[ 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

8

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.

也是一个类似于DP的图论题，设计状态为f[i][k]表示从s点到i点，使用k次免费机会的最小花费。考虑转移，无法从前面推过来，故从当前状态往i+1方向推，不难看出f[i][k]可以向f[j][k]和f[j][k+1]转移（当i与j连通时）。

#include<cstdio>
#include<cstring>
#include<string>
#include<algorithm>
#include<iostream>
#include<cmath>
#include<queue>
#define pa pair<int,int>
#define pai pair<pa,int>
#define INF 50000010

using namespace std;

const int MAXN=10010;
const int MAXM=100010;
const int MAXK=15;
struct edge
{
	int to,next,val;
};
edge e[MAXM];
int head[MAXN],cnt;
int f[MAXN][MAXK];
void addedge(int u,int v,int w)
{
	e[++cnt].next=head[u];
	e[cnt].to=v;
	e[cnt].val=w;
	head[u]=cnt;
}

priority_queue < pai , vector<pai> , greater<pai> > q;
bool visit[MAXN][MAXK];

int main()
{
	int n,m,k;
	scanf("%d%d%d",&n,&m,&k);
	int s,t;
	scanf("%d%d",&s,&t);
	s++,t++;
	for (int i=1;i<=n;i++)
	{
		for (int j=0;j<=k;j++)
		{
			f[i][j]=INF;	
		}
	}
	for (int i=1;i<=m;i++)
	{
		int u,v,w;
		scanf("%d%d%d",&u,&v,&w);
		u++,v++;
		addedge(u,v,w);
		addedge(v,u,w);
	}
	memset(visit,false,sizeof(visit));
	q.push(make_pair(make_pair(0,s),0));
	f[s][0]=0;
	while (!q.empty())
	{
		int x=q.top().first.second;
		int y=q.top().second;
		q.pop();
		if (visit[x][y])
		{
			continue;
		}
		visit[x][y]=true;
		for (int i=head[x];i;i=e[i].next)
		{
			int v=e[i].to;
			if (f[v][y]>f[x][y]+e[i].val)
			{
				f[v][y]=f[x][y]+e[i].val;
				q.push(make_pair(make_pair(f[v][y],v),y));
			}
			if (y+1>k)
			{
				continue;
			}
			if (f[v][y+1]>f[x][y])
			{
				f[v][y+1]=f[x][y];
				q.push(make_pair(make_pair(f[v][y+1],v),y+1));
			}
		}
	}	
	cout<<f[t][k];
    return 0;
}