欧拉函数模板 筛法 直接求

4 篇文章 0 订阅
2 篇文章 0 订阅

两种方法都是直接根据定义来求,

φ(n)=n(1  - 1/p1) (1 - 1/p2) (1 - 1/p3) (1 - 1/p4)…..(1 - 1/pk)

e[i] == i 与 x % i == 0 均起到了保证 i 为质数的作用,

下面直接看代码


1.

递推:

void euler() {
    e[1] = 0;
    for (int i = 2; i <= maxn; ++i) e[i] = i;
    for (int i = 2; i <= maxn; ++i) {
        if (e[i] == i) {
            for (int j = i; j <= maxn; j += i) {
                e[j] = e[j] / i * (i - 1);
            }
        }
    }
}

2. 

线性筛:

void euler2() {
    int tot = 0;
    e[1] = 1;
    for (int i = 2; i <= maxn; ++i) {
        if (!check[i]) {
            prime[tot++] = i;
            e[i] = i - 1;
        }
        for (int j = 0; j < tot; ++j) {
            if (i * prime[j] > maxn) break;
            check[i * prime[j]] = true;
            if (i % prime[j]) e[i * prime[j]] = e[i] * (prime[j] - 1);
            else {
                e[i * prime[j]] = e[i] * prime[j];
                break;
            }
        }
    }
    e[1] = 0;
}



3. 

直接求:

LL phi(LL x) {
    LL ans = x;
    for (LL i = 2; i * i <= x; ++i) {
        if (x % i == 0) {
            ans -= ans / i;
            while (x % i == 0) x /= i;
        }
    }
    if (x > 1) ans -= ans / x;
    return ans;
}


  • 0
    点赞
  • 0
    收藏
    觉得还不错? 一键收藏
  • 0
    评论

“相关推荐”对你有帮助么?

  • 非常没帮助
  • 没帮助
  • 一般
  • 有帮助
  • 非常有帮助
提交
评论
添加红包

请填写红包祝福语或标题

红包个数最小为10个

红包金额最低5元

当前余额3.43前往充值 >
需支付:10.00
成就一亿技术人!
领取后你会自动成为博主和红包主的粉丝 规则
hope_wisdom
发出的红包
实付
使用余额支付
点击重新获取
扫码支付
钱包余额 0

抵扣说明:

1.余额是钱包充值的虚拟货币,按照1:1的比例进行支付金额的抵扣。
2.余额无法直接购买下载,可以购买VIP、付费专栏及课程。

余额充值