【bzoj 2763】2763: [JLOI2011]飞行路线 spfa二维

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

理解一个模板的时候需要更深层次的理解,这题就是一个例子,虽说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;
}


评论
添加红包

请填写红包祝福语或标题

红包个数最小为10个

红包金额最低5元

当前余额3.43前往充值 >
需支付:10.00
成就一亿技术人!
领取后你会自动成为博主和红包主的粉丝 规则
hope_wisdom
发出的红包
实付
使用余额支付
点击重新获取
扫码支付
钱包余额 0

抵扣说明:

1.余额是钱包充值的虚拟货币,按照1:1的比例进行支付金额的抵扣。
2.余额无法直接购买下载,可以购买VIP、付费专栏及课程。

余额充值