筛选法求欧拉函数

可以借鉴之前的求质因数分解的线性求法。对于要求1~N上的欧拉函数的和
欧拉函数定义如下:
在这里插入图片描述
那么φ(pi * N)等于多少呢?
由定义可以知道,欧拉函数和质因子的次方无关。因此这时候要分类讨论
1.pi 如果是N的质因子,那么φ(pi * N) = pi * φ(N)
2.如果不是,φ(pi * N) = pi * φ(N)* (pi - 1) / pi

代码实现如下:

#include<iostream>

using namespace std;

const int N = 1000010;
typedef long long LL;
int euler[N], prime[N], cnt;
bool st[N];

//借助筛选欧拉函数


LL get_euler(int n)
{
    
    
    euler[1] = 1;   //根据定义
    for(int i = 2; i <= n; i++){
    
        //挑选质数
        if(!st[i]){
    
    
            prime[cnt ++] = i;
            euler[i] = i - 1;
        }
        
        for(int j = 0; prime[j] <= n / i; j ++){
    
    
            st[prime[j] * i] = true;
            if(i % prime[j] == 0){
    
    
                euler[prime[j] * i] = prime[j] * euler[i];
                break;
            }
            euler[prime[j] * i] = (prime[j]  - 1 )* euler[i];
        }
    }
    
    LL res = 0;
    for(int i = 1; i <= n ; i++) res += euler[i];
    return res;
}
int main(){
    
    
    int n;
    cin >> n;
    LL ans = get_euler(n);
    cout << ans << endl;
    return 0;
}

猜你喜欢

转载自blog.csdn.net/qq_44879626/article/details/108083724
今日推荐