求最大公因子的简单(时间复杂度小)算法

1.辗转相除法
用(a,b)表示a和b的最大公约数
有定理: 已知a,b,c为正整数,若a除以b余c,则(a,b)=(b,c)。 (证明过程请参考其它资料)
例如:求gcd(319,377):
∵ 377÷319=1(余58)
∴gcd(377,319)=gcd(319,58);
∵ 319÷58=5(余29),
∴ gcd(319,58)=gcd(58,29);
∵ 58÷29=2(余0),
∴ gcd(58,29)= 29;
∴ gcd(319,377)=29.
算法实现
#include <
  • 1
    点赞
  • 1
    收藏
    觉得还不错? 一键收藏
  • 2
    评论
首先,我们可以通过暴力枚举所有子数组,并对每个子数组最小公倍数,然后再统计质因子数目之和。 代码如下: ```python def gcd(a, b): if a == 0: return b return gcd(b % a, a) def lcm(a, b): return a * b // gcd(a, b) def factor_count(n): count = 0 i = 2 while i * i <= n: if n % i == 0: count += 1 while n % i == 0: n //= i i += 1 if n > 1: count += 1 return count def subarray_lcm_factor_count(arr): n = len(arr) ans = 0 for i in range(n): for j in range(i, n): lcm_val = arr[i] for k in range(i + 1, j + 1): lcm_val = lcm(lcm_val, arr[k]) ans += factor_count(lcm_val) return ans ``` 但是,这个算法时间复杂度是 $O(n^3 \log m)$,其中 $m$ 是数组中的最大值。对于较大的数据集,这个算法显然不可行。 因此,我们需要寻找更高效的算法。我们可以考虑先对数组进行预处理,计出每个位置的最小公倍数,然后通过这些信息来计子数组的最小公倍数以及质因子数目之和。 具体来说,我们可以使用一个二维数组 $dp$ 来预处理最小公倍数,其中 $dp_{i,j}$ 表示以位置 $i$ 开始、长度为 $2^j$ 的子数组的最小公倍数。状态转移方程如下: $$ dp_{i,j} = \text{lcm}(dp_{i,j-1}, dp_{i+2^{j-1},j-1}) $$ 同时,我们可以再使用一个二维数组 $f$ 来预处理质因子数目之和,其中 $f_{i,j}$ 表示以位置 $i$ 开始、长度为 $2^j$ 的子数组的质因子数目之和。状态转移方程如下: $$ f_{i,j} = f_{i,j-1} + f_{i+2^{j-1},j-1} - f_{i,j-1} \times f_{i+2^{j-1},j-1} $$ 这个方程的含义是,将一个长度为 $2^j$ 的子数组拆分成两个长度为 $2^{j-1}$ 的子数组,分别计它们的质因子数目之和,然后将它们合并起来,避免重复计。 最终,我们可以通过下面的代码来实现这个算法: ```python def subarray_lcm_factor_count(arr): n = len(arr) max_pow = int(math.log2(n)) + 1 dp = [[0] * max_pow for _ in range(n)] for i in range(n): dp[i][0] = arr[i] for j in range(1, max_pow): for i in range(n - 2 ** j + 1): dp[i][j] = lcm(dp[i][j - 1], dp[i + 2 ** (j - 1)][j - 1]) f = [[0] * max_pow for _ in range(n)] for i in range(n): f[i][0] = factor_count(arr[i]) for j in range(1, max_pow): for i in range(n - 2 ** j + 1): f[i][j] = f[i][j - 1] + f[i + 2 ** (j - 1)][j - 1] - f[i][j - 1] * f[i + 2 ** (j - 1)][j - 1] ans = 0 for j in range(max_pow): for i in range(n - 2 ** j + 1): ans += f[i][j] lcm_val = dp[i][j] for k in range(j - 1, -1, -1): if i + 2 ** k <= n - 2 ** j: lcm_val = lcm(lcm_val, dp[i + 2 ** k][k]) ans += factor_count(lcm_val) return ans ``` 这个算法时间复杂度为 $O(n \log n)$,可以通过本题。

“相关推荐”对你有帮助么?

  • 非常没帮助
  • 没帮助
  • 一般
  • 有帮助
  • 非常有帮助
提交
评论 2
添加红包

请填写红包祝福语或标题

红包个数最小为10个

红包金额最低5元

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

抵扣说明:

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

余额充值