【GESP】C++五级真题(数论、埃氏筛思想考点) luogu-B3969 [GESP202403 五级] B-smooth 数

### 解析 在 **GESP2403 五级编程考试** 中,关于 **B-smooth ** 的题目要求是:给定两个正整 $ n $ 和 $ B $,统计不超过 $ n $ 的所有 **B-smooth ** 的量。所谓 **B-smooth ** 是指一个正整的最大质因子不超过 $ B $。 #### 算法思路 为了解决这个问题,可以采用 **筛法** 来计算每个的最大质因子,并据此判断是否为 B-smooth 。具体步骤如下: 1. **线性筛法(欧拉筛)**: - 使用线性筛法来遍历从 $ 2 $ 到 $ n $ 的所有整- 对于每个 $ i $,如果它是质,则将其最大质因子标记为其本身;否则,在筛的过程中记录其最大的质因子。 2. **筛选最大质因子**: - 在筛法过程中,通过维护一个质列表 `prime` 和组 `mx_prime_factor`,其中 `mx_prime_factor[i]` 表示 $ i $ 的最大质因子。 - 当处理合时,利用当前的质 $ p $ 和当前 $ i $ 构造新的合 $ i \times p $,并更新该合的最大质因子为 $ \max(mx\_prime\_factor[i], p) $。 3. **统计结果**: - 遍历从 $ 1 $ 到 $ n $ 的所有,检查其最大质因子是否小于等于 $ B $,如果是,则计加一。 #### 时间复杂度 该算法的时间复杂度为 $ O(n \log \log n) $,因为使用了线性筛法来高效地处理所有[^3]。 #### 示例代码 以下是解决问题的 C++ 实现代码: ```cpp #include <bits/stdc++.h> using namespace std; int main() { int n, B; cin >> n >> B; vector<bool> vis(n + 5, false); // 标记是否已处理 vector<int> mx_prime_factor(n + 5, 0); // 记录每个的最大质因子 vector<int> prime; // 存储质列表 mx_prime_factor[1] = 1; // 特殊情况处理:1 的最大质因子设为 1 for (int i = 2; i <= n; i++) { if (!vis[i]) { mx_prime_factor[i] = i; // 如果是质,最大质因子是自身 prime.push_back(i); } for (int p : prime) { if ((long long)p * i > n) break; vis[i * p] = true; mx_prime_factor[i * p] = max(mx_prime_factor[i], p); if (i % p == 0) break; } } int ans = 0; for (int i = 1; i <= n; i++) { if (mx_prime_factor[i] <= B) ans++; } cout << ans << endl; return 0; } ``` #### 样例解析 以样例输入 `n = 10` 和 `B = 3` 为例: - 不超过 $ 10 $ 的所有整中,最大质因子不超过 $ 3 $ 的包括:`{1, 2, 3, 4, 6, 8, 9}`,共 $ 7 $ 个[^1]。 --- ###
评论
添加红包

请填写红包祝福语或标题

红包个数最小为10个

红包金额最低5元

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

抵扣说明:

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

余额充值