2020 7.31 暑期最短路练习

相关的模板可以看下我之前博客

dijkstra+heap优化 O(mlogn)

只能处理边权为正数的问题。

spfa算法 O(km)

bellman-ford算法的优化版本,可以处理存在负边权的最短路问题。
最坏情况下的时间复杂度是 O(nm)O(nm),但实践证明spfa算法的运行效率非常高,期望运行时间是 O(km)O(km),其中 kk 是常数。

floyd算法 O(n3)

Floyed算法允许图中有带负权值边,允许有回路,但不允许有带负权值边的回路。

标准弗洛伊德算法,三重循环。循环结束之后 d[i][j]d[i][j] 存储的就是点 ii 到点 jj 的最短距离。需要注意循环顺序不能变:第一层枚举中间点,第二层和第三层枚举起点和终点。

 

 

相关的练习(最后五个题是树的直径)

https://vjudge.net/contest/386791#overview

C - 最短路

思路:可以用dijkstra还可以用floyd 算法

我的这种写法便于改变边 的长度这个也是邻接矩阵的优势

 

code:

#include<iostream>
#include<cstdio>
#include<algorithm>
#include<queue>
#include<vector>
#define INF 0x3f3f3f3f

using namespace std;

const int N=110;
typedef pair<int,int> PII;
int  n,m;
int dis[N];
bool st[N];
int G[N][N]; 

void dijkstra(){
	for(int i=1;i<=n;i++){
		dis[i]=INF;
		st[i]=0;
	}
	priority_queue<PII, vector<PII> , greater<PII> > heap;
	dis[1]=0;
	heap.push({0,1});
	while(heap.size()){
		auto t=heap.top();
		heap.pop();
		int ver=t.second,distance=t.first;
		if(st[ver]) continue;
		st[ver]=1;
		for(int i=1;i<=n;i++){
			if(dis[i]>distance+G[ver][i]){
				dis[i]=distance+G[ver][i];
				heap.push({dis[i],i});
			}
		}
	} 
}
int main(){
	while(cin>>n>>m&&(n+m)){
		for(int i=1;i<=n;i++)
			for(int j=1;j<=n;j++) 
				G[i][j]=INF;
		while(m--){
			int u,v ,w;
			cin>>u>>v>>w;
			G[u][v]=G[v][u]=w;
		}
		dijkstra();
		cout<<dis[n]<<endl;
	}
	 
	 
	return 0;
	
}

 

Frogger

思路:找到所有从1到2 的路中青蛙的最短跳跃距离

 注意这个图是无向图

code

#include<iostream>
#include<cstdio>
#include<cmath>

using namespace std;

const int N=220;
double f[N][N];//记录距离
struct node{
	double x;
	double y;
}point[N]; 
int n;
int main()
{
	int k=0;
	while(cin>>n&&n){
		k++;
		for(int i=1;i<=n;i++) scanf("%lf%lf",&point[i].x,&point[i].y);
		for(int i=1;i<=n-1;i++)
			for(int j=i+1;j<=n;j++) {
				f[i][j]=f[j][i]=sqrt((point[i].x-point[j].x)*(point[i].x-point[j].x)+(point[i].y-point[j].y)*(point[i].y-point[j].y));
			}
		for(int k=1;k<=n;k++)
			for(int i=1;i<=n;i++)
				for(int j=1;j<=n;j++){
					if(f[i][j]>max(f[i][k],f[k][j])){
						f[i][j]=f[j][i]=max(f[i][k],f[k][j]);
					}
				}
		printf("Scenario #%d\n",k);
		printf("Frog Distance = %.3f\n",f[1][2]);
		puts("");
	}
	
	
	return 0;
}

 

 

一个人的旅行

 思路:有多个起点,然后有多个终点,然后将起点的dis都置为0,跑dijkstra,求到哪个终点的距离最近

 

code:

#include<iostream>
#include<queue>
#include<cstdio>
#include<vector>
#include<cstring>
#define INF 0x3f3f3f3f

using namespace std;

typedef pair<int,int> PII;
const int M=2e5+10;
const int N=1100;
int head[N],ne[M],w[M],e[M];
int dis[N];
bool st[N];
vector<int> s;
vector<int> d;
int T,S,D;
int cnt;

void add(int u,int v,int val ){
	w[cnt]=val,e[cnt]=v,ne[cnt]=head[u],head[u]=cnt++;
}

void dijkstra()
{
	for(int i=1;i<=N;i++) dis[i]=INF,st[i]=0;
	priority_queue<PII, vector<PII> ,greater<PII> > heap;
	for(int i=0;i<s.size();i++){
		dis[s[i]]=0;
		heap.push({0,s[i]});
	}
//	heap.push({0,1});
	while(heap.size()){
		auto t=heap.top();
		heap.pop();
		int ver=t.second,distance=t.first;
		if(st[ver]) continue;
		st[ver]=1;
		for(int i=head[ver];i!=-1;i=ne[i]){
			int v=e[i];
			if(dis[v]>distance+w[i]){
				dis[v]=distance +w[i];
				heap.push({dis[v],v});
			}
		}
	}
}
int main()
{
	while(~scanf("%d%d%d",&T,&S,&D)){
		s.clear();
		d.clear();
		int a,b,c;
		memset(head,-1,sizeof(head));
		
		cnt=0;
		while(T--){
			scanf("%d%d%d",&a,&b,&c);
			add(a,b,c);
			add(b,a,c);
		}
		int x;
		for(int i=1;i<=S;i++){
			scanf("%d",&x);
			s.push_back(x);	
		} 
		for(int i=1;i<=D;i++){
			scanf("%d",&x);
			d.push_back(x);
		}
		dijkstra();
		int ans=INF;
		for(int i=0;i<d.size();i++){
			ans=min(ans,dis[d[i]]);
		}
//		cout<<dis[3]<<endl;
		printf("%d\n",ans);
	}
	return 0;
}

D - HDU Today

 

思路:这个题需要注意,起点和终点可能是一样的,而且这是个无向图,真是恶心,将string变为int可以采用unordered_map

需要用到<unordered_map>这个头文件

code:

#include<iostream>
#include<cstdio>
#include<unordered_map>
#include<queue>
#include<cstring>
#define INF 0x3f3f3f3f

using namespace std;

typedef pair<int,int> PII;
const int N=2e4+10;
int head[200],ne[N],w[N],e[N];
bool st[200];
int dis[N];
int sx,ex;
int n;
int cnt;
int index;

void add(int u,int v,int val){
	e[cnt]=v,w[cnt]=val,ne[cnt]=head[u],head[u]=cnt++;
}

void dijkstra()
{
	for(int i=1;i<=index;i++) dis[i]=INF,st[i]=0;
	priority_queue<PII, vector<PII>, greater<PII> > heap;
	dis[sx]=0;
	heap.push({0,sx});
	while(heap.size()){
		auto t=heap.top();
		heap.pop();
		int ver=t.second,distance=t.first;
		if(st[ver]) continue;
		st[ver]=1;
		for(int i=head[ver];i!=-1;i=ne[i]){
			int v=e[i];
			if(dis[v]>distance+w[i]){
				dis[v]=distance+w[i];
				heap.push({dis[v],v});
			}
		}
	}
	
}
int main()
{
	while(cin>>n&&n!=-1){
		memset(head,-1,sizeof head);
		cnt=0;
		index=0;
		unordered_map<string,int> ma;
		string ss,sss;
		cin>>ss>>sss;
		if(!ma.count(ss)){
			ma[ss]=++index;
			sx=index;
		}
		if(!ma.count(sss)){
			ma[sss]=++index;
			ex=index; 
		} 
		else ex=ma[sss];
		for(int i=1;i<=n;i++){ 
			int a,b,c;
			cin>>ss>>sss>>c;
			if(!ma.count(ss)){
				ma[ss]=++index;
				a=index;
			}
			else{
				a=ma[ss];
			}
			
			if(!ma.count(sss)){
				ma[sss]=++index;
				b=index;
			}
			else{
				b=ma[sss];
			}
			add(a,b,c);
			add(b,a,c);
		}
		dijkstra();
//		cout<<dis[2]<<endl;
		if(dis[ex]==INF) cout<<"-1"<<endl;
		else cout<<dis[ex]<<endl;
	} 
		
	
	return 0;
}

 

Invitation Cards

思路:这个题开始是从单点到多点求最短路,然后是多点到单点求最短路,单点到多点的好求,直接dijkstra,多点到单点的话需要反向建图

注意答案ans需要long long 类型啊printf("%lld) 

我总是写%d写顺手了

code:

 

#include<iostream>
#include<cstdio>
#include<cstring>
#include<queue>
#define INF 0x3f3f3f3f3f3f3f3f

using namespace std;

typedef long long ll;
typedef pair<ll,int> PII;
const int N=1e6+10; 
int cnt;
int head[N],ne[N],w[N],e[N];
ll dis[N];
bool st[N];
int n;
int p,q;
int a[N],b[N],c[N];

void add(int u,int v,int val){
	e[cnt]=v,w[cnt]=val,ne[cnt]=head[u],head[u]=cnt++;
}

ll dijkstra()
{
	for(int i=1;i<=p;i++) dis[i]=INF,st[i]=0;
	priority_queue<PII,vector<PII> ,greater<PII> >heap;
	dis[1]=0;
	heap.push({0,1});
	while(heap.size()){
		PII t= heap.top();
		heap.pop();
		ll ver= t.second, distance =t.first;
		if(st[ver]) continue;
		st[ver]=1;
		for(int i=head[ver];i!=-1;i=ne[i]){
			int v=e[i];
			if(dis[v]>distance+w[i]){
				dis[v]=distance+w[i];
				heap.push({dis[v],v});
			}
		} 
	}
	ll ans=0;
	for(int i=1;i<=p;i++) ans+=dis[i];	
	return ans; 
}
int main()
{
	scanf("%d",&n);
	while(n--){
		scanf("%d%d",&p,&q);
		for(int i=1;i<=p;i++) head[i]=-1;
		cnt=0;
		for(int i=1;i<=q;i++){
			scanf("%d%d%d",&a[i],&b[i],&c[i]);
			add(a[i],b[i],c[i]);
		}
		ll ans=0;
		ans=dijkstra();
		cnt=0;
		for(int i=1;i<=p;i++) head[i]=-1;
		for(int i=1;i<=q;i++) add(b[i],a[i],c[i]);
		ans+=dijkstra();
		printf("%lld\n",ans); //注意这个输出的是longlong类型 
	}
	
	
	return 0;
}

F - 畅通工程续

板子

code:

#include<iostream>
#include<queue>
#include<cstdio>
#include<cstring>
#define INF 0x3f3f3f3f

using namespace std;

typedef pair<int,int> PII;
int n,m;
int cnt;
int head[200],ne[2200],w[2200],e[2200];
int sx,ex; 
int dis[200];
bool st[200];

void add(int u,int v,int val){
	e[cnt]=v,w[cnt]=val,ne[cnt]=head[u],head[u]=cnt++;
}

void dijkstra()
{
	for(int i=0;i<=n;i++) dis[i]=INF,st[i]=0;
	priority_queue<PII,vector<PII>,greater< PII> > heap;
	dis[sx]=0;
	heap.push({0,sx});
	while(heap.size()){
		auto t=heap.top();
		heap.pop();
		int ver= t.second, distance = t.first;
		if(st[ver]) continue;
		st[ver]=1;
		for(int i=head[ver];i!=-1;i=ne[i]){
			int v=e[i];
			if( dis[v]>distance+w[i]){
				dis[v]=distance+w[i]; 
				heap.push({dis[v],v});
			}
		}
	}
	
}
int main()
{
	while(~scanf("%d%d",&n,&m)){
		for(int i=0;i<=n;i++) head[i]=-1;
		cnt=0;
		int a,b,c;
		while(m--){
			scanf("%d%d%d",&a,&b,&c);
			add(a,b,c);
			add(b,a,c);
		}
		scanf("%d%d",&sx,&ex);
		dijkstra();
		if(dis[ex]==INF) printf("-1\n");
		else  printf("%d\n",dis[ex]);
	}
	
	
	return 0;
 } 

 

评论
添加红包

请填写红包祝福语或标题

红包个数最小为10个

红包金额最低5元

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

抵扣说明:

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

余额充值