Weekly7作业——Floyd/Dijkstra/带负环的SPFA

最新推荐文章于 2022-12-09 11:47:22 发布

明天可以有猫吗

最新推荐文章于 2022-12-09 11:47:22 发布

阅读量178

点赞数

分类专栏：程序设计思维与实践

本文链接：https://blog.csdn.net/weixin_43793842/article/details/105571648

版权

程序设计思维与实践专栏收录该内容

20 篇文章 0 订阅

订阅专栏

A - TT 的魔法猫

题目

众所周知，TT 有一只魔法猫。

这一天，TT 正在专心致志地玩《猫和老鼠》游戏，然而比赛还没开始，聪明的魔法猫便告诉了 TT 比赛的最终结果。TT 非常诧异，不仅诧异于他的小猫咪居然会说话，更诧异于这可爱的小不点为何有如此魔力？

魔法猫告诉 TT，它其实拥有一张游戏胜负表，上面有 N 个人以及 M 个胜负关系，每个胜负关系为 A B，表示 A 能胜过 B，且胜负关系具有传递性。即 A 胜过 B，B 胜过 C，则 A 也能胜过 C。

TT 不相信他的小猫咪什么比赛都能预测，因此他想知道有多少对选手的胜负无法预先得知，你能帮帮他吗？

Input

第一行给出数据组数。

每组数据第一行给出 N 和 M（N , M <= 500）。

接下来 M 行，每行给出 A B，表示 A 可以胜过 B。

Output

对于每一组数据，判断有多少场比赛的胜负不能预先得知。注意 (a, b) 与 (b, a) 等价，即每一个二元组只被计算一次。

Example

Sample Input
3
3 3
1 2
1 3
2 3
3 2
1 2
2 3
4 2
1 2
3 4
Sample Output
0
0
4

思路

因为胜负关系具有传递性，因此可以用Floyd算法求出任意两点的胜负关系（传递闭包）
- dis[a][b]=1 表示a比b强
- dis[a][b]=0 表示a与b的胜负关系不明
- dis[a][b]=0 且 dis[b][a]=0 表示a与b的胜负关系无法预先判断
由于Floyd算法的时间复杂度为O（n^3），所以进行剪枝，即若dis[i][k]=0,直接continue
注意结果的去重

总结

Floyd算法
- 用于求取图中任意两点之间的关系
- 多源最短路，求任意两点的距离关系
- 图上的传递闭包，任意两点的连通关系
- 时间复杂度较高，可能需要剪枝

实现代码

#include<iostream>
#include<cmath>
#include<cstdio>
#include<string.h>
#include<algorithm>
using namespace std;

int N,m,num,dis[505][505]={0};
int a,b;
int res=0;

void floyd(int n)
{
	for(int k=1;k<=n;k++)
	{
		for(int i=1;i<=n;i++)
		{
			if(dis[i][k]==0) continue;
			for(int j=1;j<=n;j++)
			{
				dis[i][j]=max(dis[i][j],dis[i][k]&dis[k][j]);
			}	
		}
	}
	return;
} 

int main()
{
	cin>>num;
	for(int x=0;x<num;x++)//数据组数 
	{
		res=0;
		cin>>N>>m;
		for(int i=1;i<=N;i++)
		{
			for(int j=1;j<=N;j++)
				dis[i][j]=0;
		}
		for(int y=0;y<m;y++)
		{
			cin>>a>>b;
			dis[a][b]=1;
		}
		floyd(N);
		for(int i=1;i<=N;i++)
		{
			for(int j=1;j<=N;j++)
			{
				if(dis[i][j]==0&&dis[j][i]==0) res++;
			}
		}
		cout<<(res-N)/2<<endl;
	}
	return 0;
}

B - TT 的旅行日记

题目

众所周知，TT 有一只魔法猫。

今天他在 B 站上开启了一次旅行直播，记录他与魔法猫在喵星旅游时的奇遇。 TT 从家里出发，准备乘坐猫猫快线前往喵星机场。猫猫快线分为经济线和商业线两种，它们的速度与价钱都不同。当然啦，商业线要比经济线贵，TT 平常只能坐经济线，但是今天 TT 的魔法猫变出了一张商业线车票，可以坐一站商业线。假设 TT 换乘的时间忽略不计，请你帮 TT 找到一条去喵星机场最快的线路，不然就要误机了！

Input

输入包含多组数据。每组数据第一行为 3 个整数 N, S 和 E (2 ≤ N ≤ 500, 1 ≤ S, E ≤ 100)，即猫猫快线中的车站总数，起点和终点（即喵星机场所在站）编号。

下一行包含一个整数 M (1 ≤ M ≤ 1000)，即经济线的路段条数。

接下来有 M 行，每行 3 个整数 X, Y, Z (1 ≤ X, Y ≤ N, 1 ≤ Z ≤ 100)，表示 TT 可以乘坐经济线在车站 X 和车站 Y 之间往返，其中单程需要 Z 分钟。

下一行为商业线的路段条数 K (1 ≤ K ≤ 1000)。

接下来 K 行是商业线路段的描述，格式同经济线。

所有路段都是双向的，但有可能必须使用商业车票才能到达机场。保证最优解唯一。

Output

对于每组数据，输出3行。第一行按访问顺序给出 TT 经过的各个车站（包括起点和终点），第二行是 TT 换乘商业线的车站编号（如果没有使用商业线车票，输出"Ticket Not Used"，不含引号），第三行是 TT 前往喵星机场花费的总时间。

本题不忽略多余的空格和制表符，且每一组答案间要输出一个换行

Example

思路

以起点为源点求单源最短路，得到dis1数组
再以终点为源点求单源最短路，得到dis2数组
枚举商业线（u,v,w）取min{dis1[u]+dis2[v]+w,dis1[v]+dis2[u]+w},最终再与不走商业线的答案取最小值
输出格式比较麻烦：若使用了商业线，则从u到起点进行递归，从v到终点进行递归，若没有使用商业线，则从起点到终点递归即可（~~空格的问题，主要卡在这儿了_(:з)∠)_~~ ）

总结

dijkstra算法主要解决图中没有负边的单元最短路问题，时间复杂度为O((n+m)logn)
这道题也可以利用类似于动态规划的思想解决，记录dis[u][0/1]，数组的第二维表示当前是否经过商业线，多条商业线也可以采用这种方法。

实现代码

#include<iostream>
#include<cmath>
#include<cstdio>
#include<string.h>
#include<algorithm>
const int N = 5*1e6+100;
const int M = 5*1e6+100;
const int inf = 5*1e8;
using namespace std;
 
struct edge{
	int to,next,w;
}e[M];
int head[N],tot,n,m,vis[N],dis1[N],dis2[N],pre1[N],pre2[N],S,E;
priority_queue< pair<int,int> >q;
int a,b,c;
int num;
bool shangye=false;
int s,l;//若使用商业线，记录两端
bool first=true;

void add(int x,int y,int w)
{
	e[++tot].to=y,e[tot].next=head[x];
	e[tot].w=w,head[x]=tot;
}

void dijkstra(int s,int *dis,int *pre)
{
	while(q.size()) q.pop();
	memset(vis,0,sizeof(vis));
	memset(pre,0,sizeof(pre));
	for(int i=1;i<=n;i++) dis[i]=inf;
	dis[s]=0;
	q.push(make_pair(0,s));
	while(q.size())
	{
		int x=q.top().second;
		q.pop();
		if(vis[x]) continue;
		vis[x]=1;
		for(int i=head[x];i;i=e[i].next)
		{
			int y=e[i].to,w=e[i].w;
			if(dis[y]>dis[x]+w)
			{
				pre[y]=x;
				dis[y]=dis[x]+w;
				q.push(make_pair(-dis[y],y));
			}
		}
	}
}

void shows(int x)
{
	if(x==S)
	{
		cout<<x;
		return;
	}
	shows(pre1[x]);
	cout<<" "<<x;
}

void showe(int x)
{
	if(x==E)
	{
		cout<<" "<<x;
		return;
	}
	cout<<" "<<x;
	showe(pre2[x]);
}

int main()
{
	while(scanf("%d%d%d",&n,&S,&E)!=EOF)
	{
		shangye=false;
		cin>>m;
		for(int i=0;i<N;i++)
		{
			head[i]=-1;
		}
		for(int i=0;i<m;i++)
		{
			scanf("%d%d%d",&a,&b,&c);
			add(a,b,c);
			add(b,a,c);
		}
		dijkstra(S,dis1,pre1);
		dijkstra(E,dis2,pre2);
		int Min=dis1[E];
		cin>>num;
		for(int i=0;i<num;i++)
		{
			scanf("%d%d%d",&a,&b,&c);
			if(dis1[a]+dis2[b]+c<Min)
			{
				Min=dis1[a]+dis2[b]+c;
				shangye=true;s=a,l=b;
			}
			if(dis1[b]+dis2[a]+c<Min)
			{
				Min=dis1[b]+dis2[a]+c;
				shangye=true;s=b,l=a;
			}
		}
		if(first==false) cout<<endl;
		if(shangye==true)
		{
			shows(s),showe(l);
			cout<<endl;
			cout<<s<<endl;
			cout<<Min<<endl;
		}
		else
		{
			shows(E);
			cout<<endl;
			cout<<"Ticket Not Used"<<endl;
			cout<<Min<<endl; 
		}
		first=false;
	}
	return 0;
}

C - TT 的美梦

题目

这一晚，TT 做了个美梦！

在梦中，TT 的愿望成真了，他成为了喵星的统领！喵星上有 N 个商业城市，编号 1 ～ N，其中 1 号城市是 TT 所在的城市，即首都。

喵星上共有 M 条有向道路供商业城市相互往来。但是随着喵星商业的日渐繁荣，有些道路变得非常拥挤。正在 TT 为之苦恼之时，他的魔法小猫咪提出了一个解决方案！TT 欣然接受并针对该方案颁布了一项新的政策。

具体政策如下：对每一个商业城市标记一个正整数，表示其繁荣程度，当每一只喵沿道路从一个商业城市走到另一个商业城市时，TT 都会收取它们（目的地繁荣程度 - 出发地繁荣程度）^ 3 的税。

TT 打算测试一下这项政策是否合理，因此他想知道从首都出发，走到其他城市至少要交多少的税，如果总金额小于 3 或者无法到达请悄咪咪地打出 ‘?’。

Input

第一行输入 T，表明共有 T 组数据。（1 <= T <= 50）

对于每一组数据，第一行输入 N，表示点的个数。（1 <= N <= 200）

第二行输入 N 个整数，表示 1 ～ N 点的权值 a[i]。（0 <= a[i] <= 20）

第三行输入 M，表示有向道路的条数。（0 <= M <= 100000）

接下来 M 行，每行有两个整数 A B，表示存在一条 A 到 B 的有向道路。

接下来给出一个整数 Q，表示询问个数。（0 <= Q <= 100000）

每一次询问给出一个 P，表示求 1 号点到 P 号点的最少税费。

Output

每个询问输出一行，如果不可达或税费小于 3 则输出 ‘?’。

Sample Input
2
5
6 7 8 9 10
6
1 2
2 3
3 4
1 5
5 4
4 5
2
4
5
10
1 2 4 4 5 6 7 8 9 10
10
1 2
2 3
3 1
1 4
4 5
5 6
6 7
7 8
8 9
9 10
2
3 10

Sample Output
Case 1:
3
4
Case 2:
?
?

思路

经过分析可知，这是一个含负权边的单源最短路问题，采用SPFA算法解决
判断是否存在负环：使用数组cnt[x]表示x当前最短路上的边数，若cnt[x]≥n，则说明存在负环
找到负环后，将负环所在的连通块里的点打上标记即可
当不可达/最短路小于3时，输出“？”

总结

~~数组越界这种错误……~~
注意题目的要求，两地的繁华程度不同→可能出现负权边
SPFA通常用来处理带负权边的单源最短路问题

实现代码

#include<iostream>
#include<cmath>
#include<cstdio>
#include<string.h>
#include<algorithm>
const int N = 1e5+100;
const int M = 1e5+100;
const int inf = 1e5;
using namespace std;
 
struct edge{
	int to,next,w;
}e[M];
int head[N],tot,dis[N],cnt[N];
bool vis[N],fu[N];
queue<int>q;
int T,n,m,a[N],A,B,Q,P,C;

void add(int x,int y,int w)
{
	e[++tot].to=y,e[tot].next=head[x];
	e[tot].w=w,head[x]=tot;
}

void dfs(int u)
{
	for(int i=head[u];i!=-1;i=e[i].next)
	{
		if(!fu[e[i].to])
		{
			fu[e[i].to]=true;
			dfs(e[i].to);
		}
	}
}

void spfa(int s)
{
	for(int i=1;i<=n;i++)
	{
		dis[i]=inf;
		vis[i]=false;
		cnt[i]=0;
	}
	dis[s]=0;
	vis[s]=true;
	q.push(s);
	while(!q.empty())
	{
		int u=q.front();q.pop();
		vis[u]=false;
		for(int i=head[u];i!=-1;i=e[i].next)
		{
			int v=e[i].to;
			if(dis[v]>dis[u]+e[i].w)
			{
				cnt[v]=cnt[u]+1;
				if(cnt[v]>=n)
				{
					dfs(v);
					continue;
				}
				dis[v]=dis[u]+e[i].w;
				if(!vis[v])
				{
					q.push(v);
					vis[v]=true;
				}
			}
		}
	}
}

int main()
{
	cin>>T;
	for(int x=1;x<=T;x++)
	{
		tot=0;
		for(int i=0;i<N;i++)
		{
			head[i]=-1;
			fu[i]=false;
			a[i]=0;
		}
		cin>>n;
		for(int i=1;i<=n;i++)
			cin>>a[i];
		cin>>m;
		for(int i=1;i<=m;i++)
		{
			scanf("%d%d",&A,&B);
			C=(a[B]-a[A])*(a[B]-a[A])*(a[B]-a[A]);
			add(A,B,C);
		}
		spfa(1);
		cin>>Q;
		cout<<"Case "<<x<<":"<<endl;
		for(int i=1;i<=Q;i++)
		{
			cin>>P;
			if(dis[P]<3||fu[P]||dis[P]==inf)
				cout<<"?"<<endl;
			else
				cout<<dis[P]<<endl;
		}
	}
	return 0;
}

明天可以有猫吗

关注

0
点赞
踩
1

收藏

觉得还不错? 一键收藏
0
评论
Weekly7作业——Floyd/Dijkstra/带负环的SPFA

A - TT 的魔法猫题目众所周知，TT 有一只魔法猫。这一天，TT 正在专心致志地玩《猫和老鼠》游戏，然而比赛还没开始，聪明的魔法猫便告诉了 TT 比赛的最终结果。TT 非常诧异，不仅诧异于他的小猫咪居然会说话，更诧异于这可爱的小不点为何有如此魔力？魔法猫告诉 TT，它其实拥有一张游戏胜负表，上面有 N 个人以及 M 个胜负关系，每个胜负关系为 A B，表示 A 能胜过 B，且胜负关系具...
复制链接

扫一扫