线性筛（时间复杂度o(n)）+欧拉函数

dnsn7

已于 2023-05-09 23:12:03 修改

阅读量77

点赞数 1

分类专栏：算法文章标签：算法 c++ 数据结构

于 2023-05-04 23:08:33 首次发布

本文链接：https://blog.csdn.net/qq_61594973/article/details/130497212

版权

算法专栏收录该内容

1 篇文章 0 订阅

订阅专栏

文章介绍了使用线性筛法来找出一定范围内的素数，并计算欧拉函数值。在筛法中，对于每个合数i，它能被其最小质因子p1整除，因此只需更新不大于p1的质数倍数的欧拉函数值。这种方法避免了时间复杂度过高的问题。同时，文章提到了欧拉函数的一些性质。

摘要由CSDN通过智能技术生成

const int N = 1e8 + 10;
bool is_prime[N];
ll prime[N],cnt=0;
void init(ll n) {
    for (int i = 0; i <= n; ++i)is_prime[i] = true;
    for (int i = 2; i <= n; ++i) {
        if (is_prime[i])prime[++cnt]=i;
        for (ll j = 1; j <= cnt;++j) {
            if (prime[j] * i > n)break;
            is_prime[prime[j] * i] = 0;
            if (i % prime[j] == 0)break;
        }
    }
}

用vector容易TLE

思路：

这个代码有一个关键点就是上面的herehere 这个意义就是对于一个合数mm可以分解为m=pr11∗...∗prnnm=p1r1∗...∗pnrn其中
pipi为质数，那么我们筛mm的时候之前把p1p1筛掉了，所以在枚举ii的时候。

如果ii为素数没问题，直接向后继续推（因为筛出的质数都类似m=p1∗p2m=p1∗p2的形式，所以不可能重复）。
如果为合数，那么ii可以分解成i=p1r1∗...∗pnrni=p1r1∗...∗pnrn形式其中p1−pnp1−pn是递增的，
那么p1p1是最小的那个质数。imodp1=0imodp1=0的时候，就不用继续枚举了，所以我们就只能筛出不大于p1p1的质数∗i∗i。

O(n)求莫比乌斯函数，直接在线性筛里加

欧拉函数的性质(x为质数）
① ϕ(x) = x - 1;
②如果i mod x == 0, 那么 ϕ(i * x) = x * ϕ(i);
③如果i mod x != 0，那么 ϕ(i * x) = (p - 1) * ϕ(i);

const int N = 1e7 + 10;
ll prime[N],phi[N];
bool vis[N];int cnt = 0;
#vis用于标记该点是否为非素数，prime存储素数，phi代表欧拉函数值
void init(ll n) {
	
	//phi[1] = 1;
	for (ll i = 2; i <= n; ++i) {
		if (!vis[i]) {
			prime[++cnt] = i;
			phi[i] = i - 1;
		}
		for (ll j = 1; j <= cnt && i * prime[j] <= n; ++j) {
			vis[prime[j] * i] = 1;
			if (i % prime[j]==0) {
				phi[i * prime[j]] = phi[i] * prime[j];
				break;
			}
		    phi[i * prime[j]] = phi[i] * (prime[j] - 1);
		}
	}
}