bzoj2763: [JLOI2011]飞行路线

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

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

分层图最短路。

一直感觉这东西云里雾里的,做完这道题才发现其实也就那样。

有一点DP的思想,其实还是将状态全部枚举全就可以了。

这道题除了到达哪一个点这个状态之外,还会有用了多少次免费机会的这个状态。\(k\leq10\),看到这个我就放心了(想不出了k大一点这道题要怎么做。。。)

直接设状态为:\(dis[i][j]\)表示到达i这个点y用了j次机会的最短路。

那么我们用类似背包的思想去处理第二维状态。

也就是说我们到达每一个点,都会有两种状态:

一种是当前u到v的这条边我们不用免费条件,可以得到式子:

\[dis(v)(ks)=min(dis(v)(ks),dis(u)(ks)+edge(i).dis)
\]

那么相应地,就会有用掉免费条件的情况:

\[dis(v)(ks+1)=min(dis(v)(ks+1),dis(u)(ks))
\]

再跑最短路的时候,到达每一个点都用着两种情况更新即可。

另外,bzoj数据不卡SPFA,洛谷的数据卡SPFA,实测。

#include<iostream>
#include<cstdio>
#include<queue>
#include<cstring>
using namespace std;
const int wx=50017;
inline char get_char(){
static char buf[1000001],*p1=buf,*p2=buf;
return p1==p2&&(p2=(p1=buf)+fread(buf,1,1000000,stdin),p1==p2)?EOF:*p1++;
}
#define short long long
inline short read(){
short num=0;
char c;
while(!isdigit(c=get_char()));
for(num=c-48;isdigit(c=get_char());num=((num+(num<<2))<<1)+c-48);
return num;
}
int n,m,k,s,t,num;
int head[wx],vis[wx][11],dis[wx][11];
struct e{
int nxt,to,dis;
}edge[wx*2];
void add(int from,int to,int dis){
edge[++num].nxt=head[from];
edge[num].to=to;
edge[num].dis=dis;
head[from]=num;
}
struct node{
int u,pos,d;
friend bool operator < (const node & a,const node & b){
return a.d>b.d;
}
};
priority_queue<node > q;
void Dij(){
memset(dis,0x3f,sizeof dis);
dis[s][0]=0;q.push((node){s,0,0});
while(q.size()){
int u=q.top().u;int ks=q.top().pos;
q.pop();
if(vis[u][ks])continue;vis[u][ks]=1;
for(int i=head[u];i;i=edge[i].nxt){
int v=edge[i].to;
if(dis[v][ks]>dis[u][ks]+edge[i].dis){
dis[v][ks]=dis[u][ks]+edge[i].dis;
q.push((node){v,ks,dis[v][ks]});
}
if(ks+1<=k){
if(dis[v][ks+1]>dis[u][ks]){
dis[v][ks+1]=dis[u][ks];
q.push((node){v,ks+1,dis[v][ks+1]});
}
}
}
}
}
int main(){
n=read();m=read();k=read();s=read();t=read();s++;t++;
for(int i=1;i<=m;i++){
int x,y,z;
x=read();y=read();z=read();
x++;y++;
add(x,y,z);add(y,x,z);
}
Dij();
printf("%d\n",dis[t][k]);
return 0;
}

最新文章

  1. MySQL存储过程动态SQL语句的生成
  2. 在MVC中处理异常的总结
  3. mongdb查询与排序
  4. C#中 导入和导出Excel的方法
  5. array,vertor,arraylist,hashable,hashmap等几个易混淆概念的区别
  6. C# random(number)
  7. linux动态库默认搜索路径设置的三种方法
  8. Mybatis上路_06-使用Java自动生成
  9. C++实现一个限制对象实例个数的类
  10. apache2部署django以及静态文件
  11. SQL复习一(基础知识)
  12. windows下pip安装python模块时报错
  13. [Codeforces]663E Binary Table
  14. requests库使用:通过cookie跳过验证码登录,并用Session跨请求保持cookie
  15. Django的auth【认证】模块简介
  16. 智能卡操作系统COS概述
  17. Drawing Graphs using Dot and Graphviz
  18. 使用awk分割字符串并且获取分割后的最后一个字符串
  19. memcached整理の编译
  20. 【java】java获取对象属性类型、属性名称、属性值

热门文章

  1. MySQL--Basic(二)
  2. python实现Deque
  3. 浏览器全屏的JS代码实现
  4. 苹果公司CEO乔布斯在斯坦福大学毕业典礼上的演讲
  5. cloudera上面安装Spark2.0
  6. bzoj 3998 弦论 —— 后缀自动机
  7. 51单片机的TXD、 RXD 既接了 232 又接了 485芯片 ,会导致通信失败!
  8. linux日常管理-vmstat命令
  9. Jmeter测试接口简单使用教程
  10. require用法及源码解析