乐师理工积分赛7.18

CodeForces - 483A Counterexample 【暴力】

传送门:CodeForces483A Counterexample

题目大意

  给定区间 [l,r],问是否能在区间中找到3个不同的数a,b,c满足a与b互质,b与c互质,且a与c不互质。找到则输出3个数,否则输出-1。

解题思路

   根据题意 r - l ≤ 50,暴力枚举即可。

AC代码

#include<stdio.h>
#include<string.h>
#include<algorithm>
using namespace std;
typedef long long ll;
const int INF = 0x3f3f3f3f;
int main()
{
	ll l,r;
	scanf("%lld%lld",&l,&r);
	for(ll i = l;i < r; i++)
	{
		for(ll j = i + 1; j < r; j++)
		{
			for(ll k = j + 1; k <= r; k++)
			{
				if(__gcd(i,j) == 1 && __gcd(j,k) == 1 && __gcd(i,k) != 1)
				{
					printf("%lld %lld %lld\n",i,j,k);
					return 0;
				}
			}
		}
	}
	printf("-1\n");
	return 0;
}

CodeForces - 485B Valuable Resources 【贪心】

传送门:CodeForces485B Valuable Resources

题目大意

  给定n个坐标,问:能将n个坐标包含在内的最小的正方形面积(坐标可以在正方形边上)。

解题思路

  确定最长的一条边即可。

AC代码

#include<stdio.h>
#include<string.h>
#include<algorithm>
using namespace std;
typedef long long ll;
const int INF=0x3f3f3f3f;
ll x,y;
ll minx = INF,miny = INF,maxx = -INF,maxy = -INF;
int main()
{
	int n;
	scanf("%d",&n);
	while (n--)
	{
		scanf("%lld%lld",&x,&y);
		minx = min(x,minx);
		miny = min(y,miny);
		maxx = max(x,maxx);
		maxy = max(y,maxy);
	}
	ll ans = max(maxx-minx,maxy-miny);
	printf("%lld\n",ans * ans);
	return 0;
}

CodeForces - 483B Friends and Presents 【二分 + 容斥】

传送门:CodeForces - 483B Friends and Presents

题目大意

  你有两个朋友,现在要向每个朋友分配一些正整数。你希望将cnt1个数字分给第一个朋友,将cnt2个数字分给第二个朋友。所有分配的数字都是不同的,这也意味着不应该同时向两个朋友分配同一个数字。另外,第一个朋友不喜欢能被质数x除尽的数,第二个朋友不喜欢能被质数y除尽的数,当然,你不会给你朋友他们不喜欢的数。
  找到最小的数字v,以便可以使用集合1,2,…,v中的数字来分配。

解题思路

在这里插入图片描述

AC代码

#include<bits/stdc++.h>
using namespace std;
typedef long long LL;
const int MAXN = 1e5;
const int INF = 0x3f3f3f3f;
LL cnt1, cnt2, x, y;
bool check(LL v)
{
	if (v - v / x < cnt1 || v - v / y < cnt2)
	{
		return false;
	}
	LL lcm = x / __gcd(x, y) * y;
	// if (v - v / x / y < cnt1 + cnt2)
	if (v - v / lcm < cnt1 + cnt2)
	{
		return false;
	}
	return true;
}
int main()
{
	scanf("%lld%lld%lld%lld",&cnt1, &cnt2, &x, &y);
	LL l = 0, r = 1e18, mid;
	// 二分枚举 v
	while (l <= r)
	{
		mid = (l + r) >> 1;
		if (check(mid))
		{
			r = mid - 1;
		}
		else
		{
			l = mid + 1;
		}
	}
	printf("%lld\n",r + 1);
	return 0;
}

CodeForces - 485C Bits 【位或】

传送门:CodeForces - 485C Bits

题目大意

  n个询问,对于每个询问的[l,r]区间,找出区间内的一个x,要求x的二进制位上的1要尽量多,当有多个x二进制位上1的个数相同时,要最小的x。

解题思路

  在 l <= r 的前提下将 l 从右到左的每一位尽量变为1。因为是从右到左变的,所以最终找到的数必定是最小的。

AC代码

#include<bits/stdc++.h>
using namespace std;
typedef long long LL;
const int MAXN = 1e5;
const int INF = 0x3f3f3f3f;
int main()
{
	int n;
	scanf("%d", &n);
	while (n--)
	{
		LL l, r, t = 1;
		scanf("%lld%lld", &l, &r);
		while ((l | t) <= r)
		{
			l = l | t;
			t <<= 1;
		}
		printf("%lld\n",l);
	}
	return 0;
}

CodeForces485D Maximum Value【思维 + 二分 / 预处理 】

传送门:CodeForces485D Maximum Value

题目大意

  给定一个n个数的整数序列,问: 在1 ≤ i, j ≤ n 且 ai ≥ aj 的情况下,ai % aj 的最大值为多少。

解题思路

  对于一个数x,它的最大余数是x - 1。
  对于ai ≥ aj,找出小于 k * aj 最大的 ai(k >= 2),此时ai % aj 为一个可能的答案,对所有答案取最大即可。

AC代码1【二分】

#include<stdio.h>
#include<algorithm>
using namespace std;
typedef long long ll;
const int MAXN = 2e5 + 5;
const int INF = 0x3f3f3f3f;
int a[MAXN];
int main()
{
	int n, num;
	scanf("%d",&n);
	for (int i = 0; i < n; i++)
	{
		scanf("%d",&a[i]);
	}
	sort(a, a + n);
	// 去重
	n = unique(a, a + n) - a;
	int ans = 0;
	for (int i = 0; i < n; i++)
	{
		int k = a[n - 1] / a[i] +1;
		for (int j = 2; j <= k; j++)
		{
			int pos = lower_bound(a, a + n, a[i] * j) - a - 1;
			ans = max(ans, a[pos] % a[i]);
			// 优化,当 ans大于取余 a[i]的最大余数,就没必要对 a[i]处理了
			if (ans >= a[i] - 1)
			{
				break;
			}
		}
	}
	printf("%d\n",ans);
	return 0;
}

AC代码2【预处理】

#include<bits/stdc++.h>
using namespace std;
typedef long long ll;
const int MAXN = 2e6 + 5;
const int INF = 0x3f3f3f3f;
// 数组要开 2e6,因为下标对应数据,后续要处理 a[i] * k(k >= 2)
int a[MAXN]; 
int main()
{
	int n, num;
	scanf("%d",&n);
	for (int i = 0; i < n; i++)
	{
		scanf("%d",&num);
		a[num] = num;
	}
	// 预处理
	for (int i = 0; i < MAXN; i++)
	{
		if (a[i] != i)
		{
			a[i] = a[i - 1];
		}
	}
	int ans = 0;
	for (int i = 2; i < MAXN; i++)
	{
		if (a[i] == i)
		{
			for (int j = 2; i * j < MAXN; j++)
			{
				// a[i * j - 1] 是比 i 的 j 倍小的最大的数 
				ans = max(ans,a[i * j - 1] % i);
				// 优化,当 ans大于取余 i 的最大余数,就没必要对 i 处理了
				if (ans >= i - 1)
				{
					break;
				}
			}
		}
	}
	printf("%d\n",ans);
	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、付费专栏及课程。

余额充值