WEEK7作业题

TTT的魔法猫

众所周知,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) 等价,即每一个二元组只被计算一次。

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算法,floyd算法的时间复杂度为O(n³)
算法本身的开销很大,在使用时要注意数据范围我,如果数据太大将不能使用该方法,本题目数据量为500,可以使用
代码如下

#include<cstdio>
#include<iostream>
#include<cstring>
using namespace std;
int n;
int a;
int b;
int g[510][510];

void floyd()
{
	for(int k = 1;k<=n;k++)
		for(int i = 1; i<=n;i++)
		{
			if(i==k) continue;
			if(!g[i][k])continue;
			for(int j = 1;j<=n;j++)
			{
				if(i==j||j==k) continue;
				if(g[i][k]&&g[k][j]) g[i][j]=1;
			}	
		}


}

int main()
{
	int cnt;
	cin>>cnt;
	int ans;
	while(cnt--)
	{
		memset(g,0,sizeof(g));
		ans = 0;
		int m;
		cin>>n>>m;
		while(m--)
		{
			cin>>a>>b;
			g[a][b] = 1;		
		}
		floyd();
		for(int i = 1;i <= n;i++)
		{
			for(int j = 1;j<i;j++)
			{
				if(g[i][j]==0&&g[j][i]==0) ans++;
			}
		 } 
		 cout<<ans<<endl;
	}
}

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 前往喵星机场花费的总时间。
本题不忽略多余的空格和制表符,且每一组答案间要输出一个换行

Sample input:

4 1 4
4
1 2 2
1 3 3
2 4 4
3 4 5
1
2 4 3

Sample output:

1 2 4
2
5

题目分析:
这道题主要掌握djiskra方法,djiskra方法中每个点访问一次,可以通过松弛迭代求出最短路径或者得带最小的解,是用这种方法最重要的一点是题目是否能被松弛且贪心,如果能被松弛且贪心那么说明可以采用djiskra的方法,如果不能被贪心,例如求得最长路径的问题,首先无法使用贪心,因为每个点只访问一遍,局部最优解不是全局的最优解,这种题目通常采用SPFA方法。
在这里我们从两边分头进行最短路径的求解,再用不走商业线的花费和走商业线的花费进行对比
代码如下;

#include<iostream>
#include<cstdio>
#include<cstring>
#include<queue>
using namespace std;
const int M = 10200;
struct edge{
	int to,next,w;
};
edge e[M*2];
int Head[1001],tot,n,S,E;
int vis[1001],dis1[1001],dis2[1001],pre1[1001],pre2[1001],a[1001];

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

priority_queue<pair<int,int> > q;

void dijkstra(int s,int* dis,int* pre)
{
	while(q.size()) q.pop();
	for(int i = 1;i<=n;i++)
	{
		vis[i] = 0;dis[i] = 1e9;
	}
	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)
			{
				dis[y] = dis[x] +w;
				pre[y] = x;
				q.push(make_pair(-dis[y],y));
			}
		}
	}
 } 
 int main()
 {
 	int x,y,z, t = 0;
 	while(cin>>n>>S>>E)
 	{
 		if(t!=0)
		{
			printf("\n");
		}
		t++;
		tot=0;
		for(int i=1;i<=n;i++) pre1[i]=pre2[i]=0;
		int m;cin>>m;
		for(int i=1;i<=2*m;i++) Head[i]=0;
		for(int i=1;i<=m;i++)
		{
			
			scanf("%d%d%d",&x,&y,&z);
			add(x,y,z);add(y,x,z);
		}
		dijkstra(S,dis1,pre1);
	//	for(int i = 1;i<=n;i++)cout<<dis1[i]<<" ";
		dijkstra(E,dis2,pre2);
	//	for(int i = 1;i<=n;i++)cout<<dis2[i]<<" ";
	 	int k;
	 	cin>>k;
	 	int mn = dis1[E];
	 	int x1=0,y1=0,cut = 0;
	 	while(k--)
	 	{
	 		cin>>x>>y>>z;
	 		if(dis1[y]+dis2[x]+z<mn)
			{
				mn = dis1[y]+dis2[x]+z;
				x1=y;y1 = x;
			}
			if(dis1[x]+dis2[y]+z<mn)
			{
				mn = dis1[x]+dis2[y]+z;
				x1 = x;y1 = y;
			}
		}
		if(!x1)
			{
				int s=E;
				while(s){
					a[++cut]=s;
					if(s==S) break;
					s=pre1[s];
				}
				for(int i=cut;i>=2;i--) printf("%d ",a[i]);
				printf("%d\n",a[1]);
				printf("Ticket Not Used\n");
			}
			else{
				int s=y1;
				while(s){
					a[++cut]=s;
					if(s==E) break;
					s=pre2[s];
				}
				for(int i=1;i<=cut/2;i++) swap(a[i],a[cut-i+1]);
				s=x1;
				while(s){
					a[++cut]=s;
					if(s==E) break;
					s=pre1[s];
				}
				for(int i=cut;i>=2;i--) printf("%d ",a[i]);
				printf("%d\n",a[1]);
				printf("%d\n",x1);
			}
			cout<<mn<<endl;
	}
 }
 
 

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:
?
?

分析;
题目存在负环问题,所以不能采用的diiskra方法,SPFA 方法原理是通过多次的遍历最后逼近结果的数值,党江所有的可能性遍历之后,剩下的结果就是最终的方案,在本题中因为存在负环的问题,所亦需要采用SPFA方法,因为最长的路径为n-1,如果一个点的访问的位数是n则证明有负环,在这里走了一个圈,即采用DFS镜可以走道的地方全部找出来,这些地方都会一次一次的被更新,将他们标记以后不在访问
代码如下;

#include<iostream>
#include<cstdio>
#include<cstring>
#include<algorithm>
#include<queue>
#include<cmath> 
using namespace std;
const int inf = 1e9;
int T,N,M,K,a[1001],dis[1001],tot,head[1001],cnt[1001];
bool vis[1001],vis2[1001];

struct edge{
	int to,next,w;
}e[400001];

void add(int x,int y,int w){
	e[++tot].to = y;e[tot].next = head[x];head[x] = tot;e[tot].w= w;
}
void dfs(int u)
{
	vis2[u]=true;
	for(int i=head[u];i;i=e[i].next)
	{
		int v=e[i].to;
		if(vis2[v]) continue;
		dfs(v);
	}
	return ;
} 

queue<int> q; 
void spfa(int s)
{
	while(q.size()) q.pop();
	for(int i = 0; i <= N; i++)
		vis2[i] = vis[i] = cnt[i] = 0,dis[i] = inf;
	dis[s] = 0 ; vis[s] = 1 ; q.push(s);
	while(q.size())
	{
		//将现在的点弹出,用新的距离更新现在的最短 
		int x = q.front(); q.pop();vis[x] = 0;
		//if(vis2[x] == 1) continue;
		for(int i = head[x]; i; i = e[i].next){
			int y = e[i].to;
			if(dis[y]>dis[x]+e[i].w){
				dis[y] = dis[x]+e[i].w;
				cnt[y]++;
				if(cnt[y]>N){
					dfs(y);
					continue;
				}
				if(!vis[y]) q.push(y),vis[y] = 1;
			}
		}
	}
}


int main()
{
	cin>>T;
	int x,y;
	for(int t = 1;t<=T;t++)
	{
		printf("Case %d:\n",t);
		
	scanf("%d",&N);
		tot = 0;
		for(int i = 1;i <= N; i++) head[i] = 0;
		for(int i = 1;i<=N;i++) scanf("%d",&a[i]);
		cin>>M;
		for(int i = 1;i<=M;i++){
			scanf("%d%d",&x,&y);
			add(x,y,pow((a[y]-a[x]),3));
			//cout<<(a[B]-a[A])*(a[B]-a[A])*(a[B]-a[A])<<endl;
		}
		spfa(1);
		int Q,p;
		cin>>Q;
		while(Q--)
		{
			scanf("%d",&p);
			if(vis2[p]==1||dis[p]<3||dis[p]==inf) cout<<"?"<<endl;
			else{
				cout<<dis[p]<<endl;
			}
		}
	}
 } 
  • 0
    点赞
  • 0
    收藏
    觉得还不错? 一键收藏
  • 0
    评论
评论
添加红包

请填写红包祝福语或标题

红包个数最小为10个

红包金额最低5元

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

抵扣说明:

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

余额充值