强化阶段 Day 23 算法笔记 10.4 最短路径

目录

1.dijkstra邻接矩阵

2.Emergency

3.Travel Plan

4.bellman ford

5.Emergency(bellman)

6.spfa

7.Floyd


1.dijkstra邻接矩阵

int n,g[maxv][maxv];
int d[maxv];
bool vis[maxv]={false};

void dijkstra(int s){
	fill(d,d+maxv,inf);
	d[s]=0;
	for(int i=0;i<n;i++){
		int u=-1,min=inf;
		for(int j=0;j<n;j++){
			if(vis[j]==false&&d[j]<min){
				u=j;
				min=d[j];
			}
		}
		if(u==-1) return;
		vis[u]=true;
		for(int v=0;v<n;v++){
			if(vis[v]==false&&g[u][v]!=inf&&d[u]+g[u][v]<d[v]){
				d[v]=d[u]+g[u][v];
			}
		}
	}
}

2.Emergency

#include<cstdio>
#include<cstring>
#include<stack>
#include<set>
#include<map>
#include<queue>
#include<algorithm>
#include<iostream>
using namespace std;
typedef long long ll;

const int maxv=1000;
const int nf=1000000000;

int n,m,st,ed,g[maxv][maxv],weight[maxv];
int d[maxv],w[maxv],num[maxv];
bool vis[maxv]={false};

void dijkstra(int s){
	fill(d,d+maxv,nf);
	memset(num,0,sizeof(num));
	memset(w,0,sizeof(w));
	d[s]=0;
	w[s]=weight[s];
	num[s]=1;
	for(int i=0;i<n;i++){
		int u=-1,min_d=nf;
		for(int j=0;j<n;j++){
			if(vis[j]==false&&d[j]<min_d){
				min_d=d[j];
				u=j;
			}
		}
		if(u==-1) return;
		vis[u]=true;
		for(int v=0;v<n;v++){
			if(vis[v]==false&&g[u][v]!=nf){
				if(d[u]+g[u][v]<d[v]){
					d[v]=d[u]+g[u][v];
					w[v]=w[u]+weight[v];
					num[v]=num[u];
				}else if(d[u]+g[u][v]==d[v]){
					if(w[u]+weight[v]>w[v]){
						w[v]=w[u]+weight[v];
					}
					num[v]+=num[u];
				}
			}
		}
	}
}

int main(){
	
	scanf("%d%d%d%d",&n,&m,&st,&ed);
	for(int i=0;i<n;i++){
		scanf("%d",&weight[i]);
	}
	int u,v;
	fill(g[0],g[0]+maxv*maxv,nf);
	for(int i=0;i<m;i++){
		scanf("%d%d",&u,&v);
		scanf("%d",&g[u][v]);
		g[v][u]=g[u][v];
	}
	dijkstra(st);
	printf("%d %d\n",num[ed],w[ed]);
	
	return 0;
}

3.Travel Plan

#include<cstdio>
#include<cstring>
#include<stack>
#include<set>
#include<map>
#include<queue>
#include<algorithm>
#include<iostream>
using namespace std;
typedef long long ll;

const int maxv=505;
const int inf=1000000000;

int n,m,st,ed;
int g[maxv][maxv],cost[maxv][maxv];
int d[maxv],c[maxv];
bool vis[maxv]={false};

vector<int> pre[maxv];
void dijkstra(int st){
	fill(d,d+maxv,inf);
	d[st]=0;
	for(int i=0;i<n;i++){
		int u=-1,min_d=inf;
		for(int j=0;j<n;j++){
			if(vis[j]==false&&d[j]<min_d){
				min_d=d[j];
				u=j;
			}
		}
		if(u==-1) return;
		vis[u]=true;
		for(int v=0;v<n;v++){
			if(vis[v]==false&&g[u][v]!=inf){
				if(d[u]+g[u][v]<d[v]){
					d[v]=d[u]+g[u][v];
					pre[v].clear();
					pre[v].push_back(u);
				}else if(d[u]+g[u][v]==d[v]){
					pre[v].push_back(u);
				}
			}
		}
	}
}

vector<int> temppath,path;
int mincost=inf;
void dfs(int ed){
	if(ed==st){
		temppath.push_back(ed);
		int tempcost=0;
		for(int i=temppath.size()-1;i>0;i--){
			int id=temppath[i],idnext=temppath[i-1];
			tempcost+=cost[id][idnext];
		}
		if(tempcost<mincost){
			mincost=tempcost;
			path=temppath;
		}
		temppath.pop_back();
		return;
	}
	temppath.push_back(ed);
	for(int i=0;i<pre[ed].size();i++){
		dfs(pre[ed][i]);
	}
	temppath.pop_back();
}

int main(){
	scanf("%d%d%d%d",&n,&m,&st,&ed);
	fill(g[0],g[0]+maxv*maxv,inf);
	
	int u,v;
	for(int i=0;i<m;i++){
		scanf("%d %d",&u,&v);
		scanf("%d%d",&g[u][v],&cost[u][v]);
		g[v][u]=g[u][v];
		cost[v][u]=cost[u][v];
	}
	dijkstra(st);
	dfs(ed);
	for(int i=path.size()-1;i>=0;i--){
		printf("%d ",path[i]);
	}
	printf("%d %d",d[ed],mincost);
	
	
	
	return 0;
	
}
#include<cstdio>
#include<cstring>
#include<stack>
#include<set>
#include<map>
#include<queue>
#include<algorithm>
#include<iostream>
using namespace std;
typedef long long ll;

const int maxv=505;
const int inf=1000000000;

int n,m,st,ed;
int g[maxv][maxv],cost[maxv][maxv];
int d[maxv],c[maxv];
bool vis[maxv]={false};
int pre[maxv];

void dijkstra(int s){
	fill(d,d+maxv,inf);
	fill(c,c+maxv,inf);
	for(int i=0;i<n;i++) pre[i]=i;
	d[s]=0;
	c[s]=0;
	for(int i=0;i<n;i++){
		int u=-1,min_d=inf;
		for(int j=0;j<n;j++){
			if(vis[j]==false&&d[j]<min_d){
				min_d=d[j];
				u=j;
			}
		}
		if(u==-1) return;
		vis[u]=true;
		for(int v=0;v<n;v++){
			if(vis[v]==false&&g[u][v]!=inf){
				if(d[u]+g[u][v]<d[v]){
					d[v]=d[u]+g[u][v];
					c[v]=c[u]+cost[u][v];
					pre[v]=u;
				}else if(d[u]+g[u][v]==d[v]){
					if(c[v]>c[u]+cost[u][v]){
						c[v]=c[u]+cost[u][v];
						pre[v]=u;
					}
				}
			}
		}
	}
}

void dfs(int v){
	if(v==st){
		printf("%d ",v);
		return;
	}
	dfs(pre[v]);
	printf("%d ",v);
}

int main(){
	
	scanf("%d%d%d%d",&n,&m,&st,&ed);
	int u,v;
	fill(g[0],g[0]+maxv*maxv,inf);
	fill(cost[0],cost[0]+maxv*maxv,inf);
	for(int i=0;i<m;i++){
		scanf("%d%d",&u,&v);
		scanf("%d%d",&g[u][v],&cost[u][v]);
		g[v][u]=g[u][v];
		cost[v][u]=cost[u][v];
	}
	dijkstra(st);
	dfs(ed);
	printf("%d %d",d[ed],c[ed]);
	return 0;
	
}

4.bellman ford

struct node{
	int v,dis;
};
const int maxv=5005;
const int inf=100000000;
vector<node> adj[maxv];

int n;
int d[maxv];

bool bellman(int s){
	fill(d,d+maxv,inf);
	d[s]=0;
	for(int i=0;i<n-1;i++){
		for(int u=0;u<n;u++){
			for(int j=0;j<adj[u].size();j++){
				int v=adj[u][j].v;
				int dis=adj[u][j].dis;
				if(d[u]+dis<d[v]){
					d[v]=d[u]+dis;
				}
			}
		}
	}
	for(int u=0;u<n;u++){
		for(int j=0;j<adj[u].size();u++){
			int v=adj[u][j].v;
			int dis=adj[u][j].dis;
			if(d[u]+dis<d[v]){
				return false;
			}
		}
	}
	return true;
}

5.Emergency(bellman)

#include<cstdio>
#include<cstring>
#include<stack>
#include<set>
#include<map>
#include<queue>
#include<algorithm>
#include<iostream>
using namespace std;
typedef long long ll;

const int maxv=505;
const int inf=0x3fffffff;

struct node{
	int v,dis;
	node(int _v,int _dis) : v(_v),dis(_dis) {}
};
vector<node> adj[maxv];

int n,m,st,ed,weight[maxv];
int d[maxv],w[maxv],num[maxv];
set<int> pre[maxv];

void bellman(int s){
	fill(d,d+maxv,inf);
	memset(num,0,sizeof(num));
	memset(w,0,sizeof(w));
	d[s]=0;
	w[s]=weight[s];
	num[s]=1;
	for(int i=0;i<n-1;i++){
		for(int u=0;u<n;u++){
			for(int j=0;j<adj[u].size();j++){
				int v=adj[u][j].v;
				int dis=adj[u][j].dis;
				if(d[u]+dis<d[v]){
					d[v]=d[u]+dis;
					w[v]=w[u]+weight[v];
					num[v]=num[u];
					pre[v].clear();
					pre[v].insert(u);
				}else if(d[u]+dis==d[v]){
					if(w[u]+weight[v]>w[v]){
						w[v]=w[u]+weight[v];
					}
					pre[v].insert(u);
					num[v]=0;
					set<int>::iterator it;
					for(it=pre[v].begin();it!=pre[v].end();it++){
						num[v]+=num[*it];
					}
				}
			}
		}
	}
}

int main(){
	scanf("%d%d%d%d",&n,&m,&st,&ed);
	for(int i=0;i<n;i++){
		scanf("%d",&weight[i]);
	}
	int u,v,wt;
	for(int i=0;i<m;i++){
		scanf("%d%d%d",&u,&v,&wt);
		adj[u].push_back(node(v,wt));
		adj[v].push_back(node(u,wt));
	}
	bellman(st);
	printf("%d %d\n",num[ed],w[ed]);
	return 0;
}

6.spfa

bool spfa(int s){
	fill(d,d+maxv,inf);
	memset(inq,false,sizeof(inq));
	memset(num,0,sizeof(num));
	
	queue<int> q;
	q.push(s);
	inq[s]=true;
	num[s]++;
	d[s]=0;
	while(!q.empty()){
		int u=q.front();
		q.pop();
		inq[u]=false;
		for(int j=0;j<adj[u].size();j++){
			int v=adj[u][j].v;
			int dis=adj[u][j].dis;
			if(d[u]+dis<d[v]){
				d[v]=d[u]+dis;
				if(inq[v]==false){
					q.push(v);
					inq[v]=true;
					num[v]++;
					if(num[v]>n) return false;
				}
			}
		}
	}
	return true;
}

7.Floyd

void floyd(){
	for(int k=0;k<n;k++){
		for(int i=0;i<n;i++){
			for(int j=0;j<n;j++){
				if(dis[i][k]!=inf&&dis[k][j]!=inf){
					if(dis[i][k]+dis[k][j]<dis[i][j]){
						dis[i][j]=dis[i][k]+dis[k][j];
					}
				}
			}
		}
	}
}

  • 0
    点赞
  • 0
    收藏
    觉得还不错? 一键收藏
  • 0
    评论
评论
添加红包

请填写红包祝福语或标题

红包个数最小为10个

红包金额最低5元

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

抵扣说明:

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

余额充值