区间素数筛模板

1.并查集模板2023-10-31 2.求质因数模板2023-10-31 3.二分图最大匹配模板(匈牙利算法)2023-12-02 4.欧拉函数模板2023-11-30 5.ST表模板2023-11-26 6.快速幂模板2023-10-31 7.字典树模板2023-12-04 8.矩阵模板2023-12-05 9.Dijkstra单源最短路模板2023-12-09 10.最近公共祖先模板(LCA)2023-12-12 11.拓扑排序模板2023-12-12

12.区间素数筛模板2023-12-16

13.Kruskal和Prim模板2023-12-24 14.树状数组模板2023-12-26 15.二维坐标离散化模板2024-03-07 16.单点修改区间查最值-树状数组模板2024-03-11 17.KMP模板2024-03-11 18.二叉搜索树模板2024-03-30 19.DIjkstra进阶模板路径记录按权重(结点数最小等)记录2024-03-31 20.判断负环模板2024-04-28 21.Exgcd 模板2024-05-03 22.压位高精度模板2024-05-06 23.线段树模板2024-05-12 24.扫描线模板2024-05-18 25.莫队模板2024-07-17 26.带修莫队模板2024-07-17 27.SCC缩点模板2024-07-18 28.取模+组合数2024-07-24 29.FFT 高精度乘法模板2024-07-24 30.字符串自然溢出哈希/单哈希/双哈希模板2024-07-25 31.树模板2024-07-27 32.dsu on tree 模板2024-07-28 33.线段树模板重制2024-08-03 34.主席树模板2024-08-07 35.大数质因数分解模板2024-08-08 36.线段树合并模板2024-08-09 37.int128输入输出流2024-08-09 38.Meissel_Lehmer模板2024-08-09 39.浮点高精度2024-08-13 40.自适应辛普森法2024-10-02 41.unordered_map随机底数种子2024-10-07

例题素数密度

 template<typename T>
struct segment_sieve {
	vector<bool> is_prime, is_prime_small;
	vector<T> prime;
	segment_sieve() {
		is_prime.resize(1000010);
		is_prime_small.resize(1000100);
	}
 
	//对区间[a,b]内的整数执行筛法，is_prime[i-a]=true  ---  表示i是素数 注意这里下标偏移了a，所以从0开始。
	void get_num(T a, T b) {
		if (a == 1) a = 2;
		for (T i = 0; i * i <= b; ++i)
			is_prime_small[i] = true; //对[2,sqrt(b))的初始化全为质数
		for (T i = 0; i <= b - a; ++i)
			is_prime[i] = true; //对下标偏移后的[a,b)进行初始化
 
		for (T i = 2; i * i <= b; ++i) {
			if (is_prime_small[i]) {
				for (T j = 2 * i; j * j <= b; j += i)
					is_prime_small[j] = false; //筛选[2,sqrt(b));
				//(a+i-1)/i得到最接近a的i的倍数，最低是i的2倍，然后筛选
				for (T j = max((T)2, (a + i - 1) / i) * i; j <= b; j += i)
					is_prime[j - a] = false;
			}
		}
		for (T i = 0; i <= b - a; ++i) //统计个数
			if (is_prime[i])
				prime.push_back(i + a);
	}
};

上一篇素数判断优化

下一篇Kruskal和Prim模板

本文作者：Ke_scholar

本文链接：https://www.cnblogs.com/Kescholar/p/17904984.html

posted @ 2023-12-16 16:27 Ke_scholar 阅读(9) 评论(0) 编辑收藏举报

刷新页面返回顶部

登录后才能查看或发表评论，立即登录或者逛逛博客园首页

Ke_scholar

区间素数筛模板

公告

最新随笔

我的标签

合集 (3)

随笔分类 (242)

随笔档案 (228)

阅读排行榜

评论排行榜

推荐排行榜

最新评论

	template<typename T>
	struct segment_sieve {
	vector<bool> is_prime, is_prime_small;
	vector<T> prime;
	segment_sieve() {
	is_prime.resize(1000010);
	is_prime_small.resize(1000100);
	}

	//对区间[a,b]内的整数执行筛法，is_prime[i-a]=true --- 表示i是素数注意这里下标偏移了a，所以从0开始。
	void get_num(T a, T b) {
	if (a == 1) a = 2;
	for (T i = 0; i * i <= b; ++i)
	is_prime_small[i] = true; //对[2,sqrt(b))的初始化全为质数
	for (T i = 0; i <= b - a; ++i)
	is_prime[i] = true; //对下标偏移后的[a,b)进行初始化

	for (T i = 2; i * i <= b; ++i) {
	if (is_prime_small[i]) {
	for (T j = 2 * i; j * j <= b; j += i)
	is_prime_small[j] = false; //筛选[2,sqrt(b));
	//(a+i-1)/i得到最接近a的i的倍数，最低是i的2倍，然后筛选
	for (T j = max((T)2, (a + i - 1) / i) * i; j <= b; j += i)
	is_prime[j - a] = false;
	}
	}
	for (T i = 0; i <= b - a; ++i) //统计个数
	if (is_prime[i])
	prime.push_back(i + a);
	}
	};