UVALive 6195 —— The Dueling Philosophers Problem

原题:https://icpcarchive.ecs.baylor.edu/index.php?option=com_onlinejudge&Itemid=8&page=show_problem&problem=4206


题意:有n个点,给出m条限制关系,即a必须在b前面;将这n个点进行先后排序,如果只有一种排序方法,则输出 1 ;如果有多种,则输出 2 ;如果无法给出排序(即存在矛盾关系),则输出 0 ;

思路:如果存在强连通分量,那么就输出0,因为强连通分量中的任意两点可以相互到达,既要满足a 在 b 前面,又要满足 b 在 a 前面,显然矛盾,无法给出排序;

进行拓扑排序,当遇到入度 = 0 的点超过一个时,就输出 2,因为此时这些点之间不存在偏序关系,可以把 a 放 b 前面,也可以把b 放 a前面,有多种方式;

注意:题目里的 m 的边界值少了个 0;


#include<bits/stdc++.h>
using namespace std;
const int maxn = 1100;
const int maxm = 500005;
int n, m;
int head[maxn], edgenum;
int Time, taj, top;
int DFN[maxn], Low[maxn], Stack[maxn];
bool Instack[maxn];
int in[maxn];

struct Edge
{
	int from, to, next;
}edge[maxm];

void add(int u, int v)
{
	edge[edgenum].from = u;
	edge[edgenum].to = v;
	edge[edgenum].next = head[u];
	head[u] = edgenum++;
}

vector<int>bcc[maxn];
void Tarjan(int u)
{
	DFN[u] = Low[u] = ++Time;
	Stack[top++] = u;
	Instack[u] = true;
	for(int i = head[u];i != -1;i = edge[i].next)
	{
		int v = edge[i].to;
		if(DFN[v] == -1)
		{
			Tarjan(v);
			Low[u] = min(Low[u], Low[v]);
		}
		else if(Instack[v] && Low[u] > DFN[v])
			Low[u] = DFN[v];
	}
	if(Low[u] == DFN[u])
	{
		taj++;
		bcc[taj].clear();
		while(1)
		{
			int now = Stack[--top];
			Instack[now] = false;
			bcc[taj].push_back(now);
			if(now == u)	break;
		}
	}
}

void init()
{
	memset(head, -1, sizeof head);
	edgenum = 0;
	memset(DFN, -1, sizeof DFN);
	memset(Instack, false, sizeof Instack);
	Time = taj = top = 0;
}

queue<int>Q;
bool toposort()
{
	while(!Q.empty())	Q.pop();
	int cnt = 0;
	for(int i = 1;i<=n;i++)
	{
		if(in[i] == 0)
		{
			cnt++;
			in[i]--;
			Q.push(i);
		}
	}
	if(cnt > 1)	return false;
	while(!Q.empty())
	{
		int u = Q.front();
		Q.pop();
		cnt = 0;
		for(int i = head[u];i != -1;i = edge[i].next)
		{
			int v = edge[i].to;
			in[v]--;
			if(in[v] == 0)
			{
				in[v]--;
				Q.push(v);
				cnt++;
			}
		}
		if(cnt > 1)	return false;
	}
	return true;
}

int main()
{
	while(~scanf("%d%d", &n, &m))
	{
		if(n == 0 && m == 0)	break;
		init();
		memset(in, 0, sizeof in);
		while(m--)
		{
			int u, v;
			scanf("%d%d", &u, &v);
			add(u, v);
			in[v]++;
		}
		for(int i = 1;i<=n;i++)
		{
			if(DFN[i] == -1)
			Tarjan(i);
		}
		if(taj != n)
		{
			printf("0\n");
			continue;
		}
		if(toposort())	printf("1\n");
		else	printf("2\n");
		
	}
	return 0;
}


  • 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、付费专栏及课程。

余额充值