202403-B-smooth 数(luogu-B3969)
GESP C++ 2024年3月五级真题,数论、埃氏筛思想考点,难度⭐⭐★☆☆,属于五级真题中比较简单的。洛谷难度等级普及−
luogu-B3969 [GESP202403 五级] B-smooth 数
题目要求
题目描述
小杨同学想寻找一种名为 -smooth 数的正整数。
如果一个正整数的最大质因子不超过 ,则该正整数为 -smooth 数。小杨同学想知道,对于给定的 和 ,有多少个不超过 的 -smooth 数。
输入格式
第一行包含两个正整数 和 ,含义如题面所示。
输出格式
输出一个非负整数,表示不超过 的 -smooth 数的数量。
输入输出样例 #1
输入 #1
10 3输出 #1
7说明/提示
数据规模与约定
| 子任务 | 得分 | ||
|---|---|---|---|
对全部的测试数据,保证 。
题目分析
解题思路
核心思想:先筛出所有 的质数,再用“埃氏筛”思想把这些质数的倍数全部标记为 -smooth,最后统计标记个数。
- 先线性筛出所有 的质数,得到“允许质因子”集合。
- 用埃氏筛思想:把上述质数的所有倍数都标记为 -smooth。
- 若某数被 的质数整除,则其最大质因子必 ,故撤销该数的 -smooth 标记。
- 最后统计被标记的个数,并额外加上 (因为 恒为 -smooth)。
示例代码
#include <iostream>
// 标记数组:nums[i]=1 表示 i 是 B-smooth 数,反之不是
int nums[1000005];
// 判断 x 是否为质数
bool is_prime(int x) {
if (x < 2) return false;
for (int i = 2; i * i <= x; ++i)
if (x % i == 0) return false;
return true;
}
int main() {
int n, B;
std::cin >> n >> B;
int count = 0; // 当前 B-smooth 数个数
for (int i = 2; i <= n; ++i) {
if (i <= B && is_prime(i)) {
// i 是 ≤B 的质数,将其倍数全部标记为 B-smooth
for (int j = i; j <= n; j += i) {
if (nums[j] == 0) { // 首次被标记
nums[j] = 1;
count++;
}
}
}
if (i > B && is_prime(i)) {
// i 是 >B 的质数,其倍数都不再是 B-smooth,撤销标记
for (int j = i; j <= n; j += i) {
if (nums[j] == 1) { // 之前被标记过
nums[j] = 0;
count--;
}
}
}
}
// 1 恒为 B-smooth,故结果 +1
std::cout << count + 1 << std::endl;
return 0;
}本文由coderli.com原创,按照CC BY-NC-SA 4.0 进行授权
所有代码已上传至Github:https://github.com/lihongzheshuai/yummy-code
“luogu-”系列题目可在 洛谷题库 在线评测。
“bcqm-”系列题目可在 编程启蒙题库 在线评测。
GESP/CSP认证交流QQ群: 688906745
GESP/CSP 认证学习微信公众号

最后更新于