P1621 集合
1.P1077 [NOIP2012 普及组] 摆花2.P3842 [TJOI2007] 线段3.P1164 小A点菜4.P2392 kkksc03考前临时抱佛脚5.P1102 A-B 数对6.P1002 [NOIP2002 普及组] 过河卒7.P163 银行贷款8.P1182 数列分段 Section II9.[NOIP2015 提高组] 跳石头10.P1734 最大约数和11.P8772 [蓝桥杯 2022 省 A] 求和12.P1216 [USACO1.5] [IOI1994]数字三角形 Number Triangles13.P5602 小 E 与美食14.P3817 小A的糖果15.P1130 红牌16.P1196 [NOI2002] 银河英雄传说17.P1955 [NOI2015] 程序自动分析
18.P1621 集合
19.P5250 【深基17.例5】木材仓库20.P2184 家谱21.[USACO16DEC] Cities and States S22.P4653 [CEOI2017] Sure Bet23.P1578 奶牛浴场24.P1387 最大正方形25.小鸟的设备26.P1439 【模板】最长公共子序列27.P1091 [NOIP2004 提高组] 合唱队形28.P1020 [NOIP1999 提高组] 导弹拦截29.P2340 [USACO03FALL] Cow Exhibition G30.P2285 [HNOI2004] 打鼹鼠31.P4995 跳跳!32.P3817 小A的糖果33.P5019 [NOIP2018 提高组] 铺设道路34.P1969 [NOIP2013 提高组] 积木大赛35.P1094 [NOIP2007 普及组] 纪念品分组36.P4447 [AHOI2018初中组] 分组集合
题目描述
Caima 给你了所有 $[a,b]$ 范围内的整数。一开始每个整数都属于各自的集合。每次你需要选择两个属于不同集合的整数,如果这两个整数拥有大于等于 $p$ 的公共质因数,那么把它们所在的集合合并。
重复如上操作,直到没有可以合并的集合为止。
现在 Caima 想知道,最后有多少个集合。
输入格式
一行,共三个整数 $a,b,p$,用空格隔开。
输出格式
一个数,表示最终集合的个数。
样例 #1
样例输入 #1
10 20 3
样例输出 #1
7
提示
样例 1 解释
对于样例给定的数据,最后有 ${10,20,12,15,18},{13},{14},{16},{17},{19},{11}$ 共 $7$ 个集合,所以输出应该为 $7$。
数据规模与约定
- 对于 $80%$ 的数据,$1 \leq a \leq b \leq 10^3$。
- 对于 $100%$ 的数据,$1 \leq a \leq b \leq 10^5,2 \leq p \leq b$。
算法1
(筛法+并查集)
Q.什么是公共质因数?
该数为质数,同时是公共因子。
比如,3是6,9的公共质因子。
整体思路
运用并查集思想,一开始每一个数是一个独立的集合。
1.枚举所有大于等于p,且小于等于b的质数
判断质数的板子
bool prime(int x)
{
bool flag = true;
for(int i = 2; i <= x/i ; i++){
if(x % i == 0) {
flag = false;
break;
}
}
if(flag) return true;
else return false;
}
- 当x是一个很大的数的时候,每次调用sqrt(x),速度慢
- 当i是一个接近于int的最大范围的时候 i*i 会溢出
2.找出该质素的所有倍数 - 合并为一个集合
3.当选取的两个数集合不同时,合并
4.优化: 标记所有已经属于某一个集合的数
5.最后,枚举所有根节点是自己的集合
C++ 代码
#include <bits/stdc++.h>
using namespace std;
const int N = 1e5+10;
int a,b,p;
int f[N];
bool mark[N];
int ans;
bool prime(int x)
{
bool flag = true;
if(x == 2) return true;
for(int i = 2; i <= x/i; i++){
if(x % i == 0) {
flag = false;
break;
}
}
if(flag) return 1;
else return 0;
}
int find(int x)
{
if(x!=f[x]) f[x] = find(f[x]);
return f[x];
}
int main(){
scanf("%d%d%d",&a,&b,&p);
for(int i = 1; i <= b; i ++) f[i] = i;
//从p 开始枚举所有大于等于p且小于等于b的质数
for(int i = p; i <= b; i++){
//已经属于某一个集合,则跳过
if(mark[i]) continue;
if(prime(i)){
mark[i] = true;
//找出所有p的倍数
for(int j = 2; i * j <=b; j++){
if(i * j < a) continue;
mark[i * j] = true;
int f1 = find(i), f2 = find(i*j);
if(f1 != f2) f[f1] = f2; //合并
}
}
}
for(int i = a; i <= b;i ++){
if(i == f[i]) ans++;
}
cout << ans << endl;
return 0;
}
【推荐】国内首个AI IDE,深度理解中文开发场景,立即下载体验Trae
【推荐】编程新体验,更懂你的AI,立即体验豆包MarsCode编程助手
【推荐】抖音旗下AI助手豆包,你的智能百科全书,全免费不限次数
【推荐】轻量又高性能的 SSH 工具 IShell:AI 加持,快人一步
· 分享一个免费、快速、无限量使用的满血 DeepSeek R1 模型,支持深度思考和联网搜索!
· 基于 Docker 搭建 FRP 内网穿透开源项目(很简单哒)
· ollama系列01:轻松3步本地部署deepseek,普通电脑可用
· 25岁的心里话
· 按钮权限的设计及实现