飞行路线
这道题是在海亮时的一道练习题,但是当时没做完,打包带了回来(一定不要在bzoj上打注释,会挂)
#include<iostream>
#include<cstring>
#include<string.h>
#include<cstdio>
#include<cctype>
#include<algorithm>
#include<queue>
using namespace std;
const int maxn = 110000;
const int maxxn = 12000;
const int inf = 0x3f3f3f;
queue<int > q;
int ans = 214748364;
int k;
int st,e;
int n,m,tot = 0;
int d[maxxn][20];//用第二位表示用了几次免费航行
bool v[maxxn][20];//由于用不同次数的免费造成的这里也要开二位
int head[maxn];
struct node {
int fr,ver,value,Next;
};
node edge[maxn * 2];
void spfa(int s) {
memset(d,inf,sizeof(d));
memset(v,0,sizeof(v));
q.push(s);
q.push(0);
v[s][0] = 1;
while (!q.empty()) {
int li = q.front();
q.pop();
int vv = q.front();
q.pop();//这里用了两次pop,实现了用一个以int为内容的queue来实现一个以pair为内容的queue所做的事;
v[li][vv] = 0;
d[s][vv] = 0;
for (int i = head[li];i;i = edge[i].Next) {
int y = edge[i].ver;
if (d[y][vv] > d[li][vv] + edge[i].value) {
d[y][vv] = d[li][vv] + edge[i].value;//在同一层中跑最短路
if (!v[y][vv]) {
q.push(y);
q.push(vv);
v[y][vv] = 1;
}
}
if (vv < k) {//这里就是用来进行分层图最短路的不同,一个状态的继承 ,由于下面是对 vv + 1 进行更新且vv为整数,所以 vv < k就是
//vv + 1 <= k;
if (d[y][vv + 1] > d[li][vv]) {//可以知道,这个时候d[y][ vv + 1]还未遍历,因为上面都不涉及vv + 1 的遍历 ,但是之所以还要写
//这一行,是因为这是spfa,所以一个点可以被便利多次;
d[y][vv + 1] = d[li][vv];//一个继承操作,用一次免费航行后从li到y,所以两个相等
if (!v[y][vv + 1]) {
q.push(y);
q.push(vv + 1);
v[y][vv + 1] = 1;
}
}
}
}
}
}
inline int Read() {
int xx = 0;
int ff = 1;
char ch = getchar();
while (ch > '9'||ch < '0') {
if (ch == '-') ff = -1;
ch = getchar();
}
while (ch <= '9'&&ch >= '0') {
xx = (xx << 1) + (xx << 3) + ch - '0';
ch = getchar();
}
return xx * ff;
}
inline void add (int x,int y,int z) {
edge[++tot].fr = x; edge[tot].ver = y; edge[tot].value = z; edge[tot].Next = head[x]; head[x] = tot;
}
int main () {
n = Read(); m = Read(); k = Read();
st = Read(); e = Read();
memset(head,0,sizeof(head));
memset(edge,0,sizeof(edge));
for (int i = 1;i <= m;i++) {
int a = Read(); int b = Read(); int c = Read();
add (a,b,c);
add (b,a,c);
}
spfa(st);
for (int i = 0;i <= k;i++) {
ans = min(ans,d[e][i]);
}
printf("%d\n",ans);
return 0;
}