理解一个模板的时候需要更深层次的理解,这题就是一个例子,虽说spfa是一个很迷的算法(你能一下子说出为什么spfa的复杂度就是那么低吗?)其实大多数的时候推荐使用dijikstra(复杂度是确定的nlogn,虽然有点常数),以前在cf上的比赛就有过卡spfa的,好像是一个深层菊花什么什么图,在最坏的情况下复杂度可以达到n2(扯得有点远)。回到这道题,事实上spfa也是广义的递推,所以思维不要局限,加了一个条件我们就加一维表示还有多少次免费可以用,每次转移情况只有2,,spfa卡着时限过。
http://blog.csdn.net/wzq_qwq/article/details/48445885 这个博客同样大同小异,但是把k剖成k层图,套dijikstra就快得多
#include<cstdio>
#include<cstring>
#include<iostream>
#include<queue>
#define MK make_pair
#define maxn 50020
using namespace std;
//dis[i][20]
int n,m,k,s,t,head[maxn*2],dis[maxn][12],vis[maxn][12],tot;
struct edge{
int v,w,next;
}e[maxn*2];
void adde(int a,int b,int c){e[tot].v=b,e[tot].next=head[a],e[tot].w=c;head[a]=tot++;}
typedef pair<int,int> pii;
queue<pii>q;
void spfa(){
memset(dis,127,sizeof(dis));
dis[s][k]=0,q.push(MK(s,k));
while(!q.empty()){
int u=q.front() . first,t=q.front() . second;q.pop();
vis[u][t]=0;
for(int i=head[u];i!=-1;i=e[i].next ){
int v=e[i].v ,w=e[i].w ;
if(dis[v][t]>dis[u][t]+w){
dis[v][t]=dis[u][t]+w;
if(!vis[v][t]){
vis[v][t]=1;
q.push(MK(v,t));
}
}
if(t>0&&dis[v][t-1]>dis[u][t]){
dis[v][t-1]=dis[u][t];
if(!vis[v][t-1]){
vis[v][t-1]=1;
q.push(MK(v,t-1));
}
}
}
}
int ans=1e9;
for(int i=0;i<=k;i++){
ans=min(ans,dis[t][i]);
}
printf("%d",ans);
}
int main(){
memset(head,-1,sizeof(head));
scanf("%d%d%d",&n,&m,&k);
scanf("%d%d",&s,&t);
s++,t++;
for(int a,b,c,i=1;i<=m;i++){
scanf("%d%d%d",&a,&b,&c);
a++,b++;
adde(a,b,c),adde(b,a,c);
}
spfa();
return 0;
}

本文深入探讨了SPFA算法的应用及优化技巧,通过一道具体题目解释了如何利用SPFA解决带有限制条件的最短路径问题,并对比了Dijkstra算法的不同之处。

315

被折叠的 条评论
为什么被折叠?



