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

 //欧拉函数,质数
vector<int> euler_range(int n) {
    vector<int> phi(n + 1), prime;
    vector<bool> is_prime(n + 1, true);
    is_prime[1] = 0, phi[1] = 1;
    for (int i = 2; i <= n; i++) {
        if (is_prime[i]) prime.push_back(i), phi[i] = i - 1;
        for (int j = 0; j < (int)prime.size() && i * prime[j] <= n; j++) {
            is_prime[i * prime[j]] = 0;
            if (i % prime[j]) phi[i * prime[j]] = phi[i] * (prime[j] - 1);
            else {
                phi[i * prime[j]] = phi[i] * prime[j];
                break;
            }
        }
    }
    return phi;
}

上一篇[KEYENCE Programming Contest 2023 Autumn（AtCoder Beginner Contest 325） E

下一篇二分图最大匹配模板(匈牙利算法)

本文作者：Ke_scholar

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

posted @ 2023-11-30 17:48 Ke_scholar 阅读(20) 评论(0) 编辑收藏举报

刷新页面返回顶部

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

Ke_scholar

欧拉函数模板

公告

最新随笔

我的标签

合集 (3)

随笔分类 (242)

随笔档案 (228)

阅读排行榜

评论排行榜

推荐排行榜

最新评论

	//欧拉函数,质数
	vector<int> euler_range(int n) {
	vector<int> phi(n + 1), prime;
	vector<bool> is_prime(n + 1, true);
	is_prime[1] = 0, phi[1] = 1;
	for (int i = 2; i <= n; i++) {
	if (is_prime[i]) prime.push_back(i), phi[i] = i - 1;
	for (int j = 0; j < (int)prime.size() && i * prime[j] <= n; j++) {
	is_prime[i * prime[j]] = 0;
	if (i % prime[j]) phi[i * prime[j]] = phi[i] * (prime[j] - 1);
	else {
	phi[i * prime[j]] = phi[i] * prime[j];
	break;
	}
	}
	}
	return phi;
	}