程序设计思维与实践 Week7 作业 (1/2/智能班)

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

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

解题思路

  • 本题可转换为求多源最短路问题
  • A胜过B,B胜过C,则A也能胜过C 可转换为有向图中有A到B的路,B到C的路,则A可直接到C。
  • 用二维数组表示边,对输入的所有胜负关系,A胜过B则dis[A][B]=1,dis[B][A]=0;
  • 遍历所有点,如果dis[l][i]=1并且dis[i][j]=1,则dis[l][j]也置为1;
  • 最后对每个点,遍历它“后面”的所有点,如果dis[i][j]和dis[j][i]一方等于1则count++,否则无法判断;最后用所有情况的数目减去count即可

注意事项
如不进行剪枝,则复杂度为 O(n³),会超时。所以在遍历所有点确定最短路时,只有当dis[l][i]=1时再往下判断i到其他点的输赢。

完整代码

#include<iostream>
#include<cstring>
using namespace std;
int dis[501][501];
int n,m;

int Floyd()
{
	for(int i=1;i<=n;i++)
		for(int l=1;l<=n;l++)
			if(dis[l][i]==1)
			for(int j=1;j<=n;j++)
			{
			if(dis[i][j]==1)
				dis[l][j]=1;}
	int count=0;
	for(int i=1;i<=n;i++)
		for(int j=i+1;j<=n;j++)
			{
			if(dis[i][j]==1||dis[j][i]==1)
				count++;
				 
			}
		return count;
}
int main()
{	
	int number;
	cin>>number;
	while(number--)
	{cin>>n>>m;
	int a,b;
	memset(dis,0,sizeof(dis));
	for(int i=1;i<=m;i++)
	{	cin>>a>>b;
		dis[a][b]=1;
		dis[b][a]=0;	
	}
	int answer=0;
	for(int i=n-1;i>0;i--)
		answer=answer+i;
	answer=answer-Floyd();
	cout<<answer<<endl;
	}	

}

B - TT 的旅行日记

问题描述

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

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

输入

输入包含多组数据。每组数据第一行为 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 行是商业线路段的描述,格式同经济线。

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

输出

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

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

输入样例

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

输出样例

1 2 4
2
5

解题思路

  • 本题为单源最短路问题,可使用Dijkstra算法
  • 如果没有商业线,则设从起点到所有点的距离为无穷大,再用最小堆遍历所有能到达的点,每次取堆顶元素u,对每个邻接点v,根据dis[v]和dis[u]+w(u,v)比较进行松弛。
  • 考虑到商业线,则可以以起点为源点求单源最短路(dis1),再以终点为源点求单源最短路(dis2),枚举商业线(u, v, w),取min{dis1[u]+dis2[v]+w, dis1[v]+dis2[u]+w},最终再与不走商业线的答案取min即可

完整代码

#include<cstdio>
#include<vector>
#include<cstring>
#include<queue>
using namespace std;

struct Edge
{
	int to; 
	int w;
}; 

vector <vector< Edge > > m;//二维数组,e[i][j]表示与i节点 相邻的第j个节点  
priority_queue<pair<int,int> > q;
int reach[510],dis1[510],dis2[510],pre1[510],pre2[510];
int X,Y,Z,M,K,n,s,e,u,v,ans;
int x1,y1,w1;


void dijkstra1()
{	memset(reach,0,sizeof(reach));
	pair<int,int>pa;
	while(!q.empty()) 
		q.pop();
	for(int i=1;i<=n;i++)
		dis1[i]=5000100;
	dis1[s]=0;
	pa.first=0;
	pa.second=s;
	q.push(pa);
	while(!q.empty())
	{int x2=q.top().second;
	q.pop();
	if(reach[x2]==1)
	 	continue;
	reach[x2]=1;
	for(int i=0;i<m[x2].size();i++)//访问x2所有临边 
		{int y2=m[x2][i].to;
		if(dis1[y2]>dis1[x2]+m[x2][i].w)
		{	dis1[y2]=dis1[x2]+m[x2][i].w;
				pa.first=-dis1[y2];
				pa.second=y2;
				pre1[y2]=x2;
				q.push(pa);}
				}
	}


}

void dijkstra2()
{	
	memset(reach,0,sizeof(reach)); 
	pair<int,int>pa;
	while(!q.empty()) 
		q.pop();
	for(int i=1;i<=n;i++)
		dis2[i]=5000100;
	dis2[e]=0;
	pa.first=0;
	pa.second=e;
	q.push(pa);
	while(!q.empty())
	{int x2=q.top().second;
	q.pop();
	if(reach[x2]==1)
	 	continue;
	reach[x2]=1;
	for(int i=0;i<m[x2].size();i++)
		{int y2=m[x2][i].to;
		
		if(dis2[y2]>dis2[x2]+m[x2][i].w)
		{	dis2[y2]=dis2[x2]+m[x2][i].w;
				pa.first=-dis2[y2];
				pa.second=y2;
		   		q.push(pa);
				pre2[y2]=x2;}
				}
	}


}




void out1(int begin,int end)
{
	
	if(end==begin)
	{printf("%d",end);	return;}
	out1(begin,pre1[end]);
	printf(" %d",end);
		
}

void out2(int begin,int end)
{  
	int ss=begin;
	while(ss!=0&&ss!=end){
		printf("%d ",ss);
		ss=pre2[ss];
	}
	printf("%d",ss);
}

int main()
{	int t=0;
	while(scanf("%d %d %d",&n,&s,&e)!=EOF)
		{
		if(t!=0)
			printf("\n");
		Edge ed;
		pre1[s]=0;
		pre2[e]=0;
		scanf("%d",&M);
		m.clear();
		m.resize(M+10); 
		
		for(int i=0;i<M;i++)
		{	scanf("%d %d %d",&x1,&y1,&w1);
			ed.to=y1;
			ed.w=w1;
			m[x1].push_back(ed);
			ed.to=x1;
			m[y1].push_back(ed);
		
			}
			
		dijkstra1();
		dijkstra2();
		
		scanf("%d",&K);
		ans=5000100; 
		for(int i=0;i<K;i++)
		{	scanf("%d %d %d",&x1,&y1,&w1);
			if(dis1[x1]+dis2[y1]+w1<ans)
			{
				u=x1,v=y1,ans=dis1[x1]+dis2[y1]+w1;
			}  
			
			if(dis1[y1]+dis2[x1]+w1<ans)
			{
				u=y1,v=x1,ans=dis1[y1]+dis2[x1]+w1;
			} 
			 
			}	
			if(ans>dis1[e])//不走商业线
			{
			out1(s,e);
			printf("\nTicket Not Used\n");
			printf("%d\n",dis1[e]);  	 
			}
			else//先u,再u到v,再v到e 
			{
				out1(s,u);
				printf(" ");
				out2(v,e);
			printf("\n%d\n",u);
			printf("%d\n",ans);
			}
			t++;
		}
	


}

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

解题思路

  • 本题为最短路问题,但是税可能是负数,即边权可能为负,所以不能用Dijkstra来求解。
  • 本题采用SPFA算法。使用队列对所有边进行松弛,每次从队首取点并松弛其邻接点,松弛成功则将其加入队列,直到队空。
  • 如果存在负环,则税收可以无限小,最短路不存在。
  • 可以通过判断最短路的边数来发现负环,如果到某一点的最短路的边数超过了n-1 则说明有负环。通过一遍bfs,将所有与负环上的点的可达的点全部标记。
  • 输出‘?’:
    • 输出时标记了的点
    • 距离小于3的点
    • 未到达的点
  • 其余正常输出最短路。

完整代码

#include<iostream>
#include<vector>
#include<queue>
#include<cmath>
#define inf 10000000
using namespace std;

struct Edge
{
	int u;
	int v;
	int w;
};

vector <vector< Edge > > e; 
int T,N,M,Q,P,A,B;
int reach1[210],reach2[210],dis[210],cnt[210],a[210];
queue <int> q;

void initial()   //初始化
{
	for(int i=0;i<=N;i++)
	{
		dis[i] = inf;
		reach1[i] = 0;
		reach2[i] = 0;
		cnt[i] = 0;
	}
	while(!q.empty()) q.pop();
	e.clear();
	e.resize(N+10); 
}

void dfs(int x)
{
	reach2[x] = 1;
	for(int i=0;i<e[x].size();i++)
	{
		x = e[x][i].v;
		if(reach2[x] == 1)  continue;
		dfs(x);
	}
}

void SPFA()
{	
	dis[1]=0;
	q.push(1);
	while(!q.empty())
	{
		int node=q.front();
		q.pop(); 
		reach1[node]=0;
		for(int i=0;i<e[node].size();i++)
			{int u=e[node][i].u;
			int v=e[node][i].v;
			int w=e[node][i].w;
			if(dis[v]>dis[u]+w)
				{
				dis[v]=dis[u]+w;
				cnt[v]=cnt[u]+1;
				if(cnt[v]>=N)//负环
				 	dfs(v);
				
			if(reach1[v]==0&&reach2[v]==0)	
			{	q.push(v);
				reach1[v]=1;}
			}}
		
	
	
	}	
	


}
int main()
{
	cin>>T;
	for(int times=1;times<=T;times++)
	{
	cin>>N;
	initial();
	for(int i=1;i<=N;i++)
		cin>>a[i];
	cin>>M;
	while(M--)
	{cin>>A>>B;
	Edge edg;
	edg.u=A;
	edg.v=B;
	edg.w=pow(a[B]-a[A],3);
	e[A].push_back(edg);}
	SPFA();
	cin>>Q;
	cout<<"Case "<<times<<":"<<endl;
		while(Q--)
	{cin>>P;

	if(reach2[P]==1 || dis[P] == inf || dis[P]<3) 
		{cout<<"?"<<endl;
		continue;}
	cout<<dis[P]<<endl;	
	}
	}


}
  • 0
    点赞
  • 0
    收藏
    觉得还不错? 一键收藏
  • 0
    评论
在信号处理领域,DOA(Direction of Arrival)估计是一项关键技术,主要用于确定多个信号源到达接收阵列的方向。本文将详细探讨三种ESPRIT(Estimation of Signal Parameters via Rotational Invariance Techniques)算法在DOA估计中的实现,以及它们在MATLAB环境中的具体应用。 ESPRIT算法是由Paul Kailath等人于1986年提出的,其核心思想是利用阵列数据的旋转不变性来估计信号源的角度。这种算法相比传统的 MUSIC(Multiple Signal Classification)算法具有较低的计算复杂度,且无需进行特征值分解,因此在实际应用中颇具优势。 1. 普通ESPRIT算法 普通ESPRIT算法分为两个主要步骤:构造等效旋转不变系统和估计角度。通过空间平移(如延时)构建两个子阵列,使得它们之间的关系具有旋转不变性。然后,通过对子阵列数据进行最小二乘拟合,可以得到信号源的角频率估计,进一步转换为DOA估计。 2. 常规ESPRIT算法实现 在描述中提到的`common_esprit_method1.m`和`common_esprit_method2.m`是两种不同的普通ESPRIT算法实现。它们可能在实现细节上略有差异,比如选择子阵列的方式、参数估计的策略等。MATLAB代码通常会包含预处理步骤(如数据归一化)、子阵列构造、旋转不变性矩阵的建立、最小二乘估计等部分。通过运行这两个文件,可以比较它们在估计精度和计算效率上的异同。 3. TLS_ESPRIT算法 TLS(Total Least Squares)ESPRIT是对普通ESPRIT的优化,它考虑了数据噪声的影响,提高了估计的稳健性。在TLS_ESPRIT算法中,不假设数据噪声是高斯白噪声,而是采用总最小二乘准则来拟合数据。这使得算法在噪声环境下表现更优。`TLS_esprit.m`文件应该包含了TLS_ESPRIT算法的完整实现,包括TLS估计的步骤和旋转不变性矩阵的改进处理。 在实际应用中,选择合适的ESPRIT变体取决于系统条件,例如噪声水平、信号质量以及计算资源。通过MATLAB实现,研究者和工程师可以方便地比较不同算法的效果,并根据需要进行调整和优化。同时,这些代码也为教学和学习DOA估计提供了一个直观的平台,有助于深入理解ESPRIT算法的工作原理。
评论
添加红包

请填写红包祝福语或标题

红包个数最小为10个

红包金额最低5元

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

抵扣说明:

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

余额充值