数论篇-------欧拉函数（简记）

欧拉函数

特别的，1的欧拉函数是1

思路:你会发现，利用到算术基本定理，然后再乘以每个质因数的式子即可

例题：给定n个正整数ai，请你求出每个数的欧拉函数。

# include <iostream>
# include <algorithm>
using namespace std;
int main(void)
{
    int t;
    cin>>t;
    while(t--){
        int x;
        cin>>x;
        long long res = x;
        for(int i=2;i*i<=x;++i){
            if(x%i==0){
                //res = res*(1-1/i); 这里拆开，主要是避免小数的存在
                res = (res/i)*(i-1);
                while(x%i==0) x/=i;
            }
        }
        if(x>1) res = (res/x)*(x-1);
        cout<<res<<endl;
    }
    return 0;
}

那求连续的数的欧拉函数之和怎么求？

最简单的做法就是直接用for循环加上面的代码，但是这样的时间复杂度是O(n * 根号n)，这里用到欧拉线性筛素数法来优化，以达到O(n)。

（下面的图主要是给自己看的）

#include<bits/stdc++.h>
using namespace std;
typedef long long ll;
const int N=1e6+10;

int primes[N],cnt;
int phi[N];   //存欧拉函数值
bool st[N];


int get_eulers(int n){
    phi[1]=1;
    for(int i=2;i<=n;i++){
        if(!st[i]){
            primes[cnt++]=i;
            phi[i]=i-1; //从定义出发，这个数i是质数，那么1到i-1都与它互质
        }
        for(int j=0;primes[j]<=n/i;j++){
            st[primes[j]*i]=true;
            if(i%primes[j]==0){  //对应上面第一张图
                phi[primes[j]*i]=phi[i]*primes[j];
                break;
            }
            //对应上面第二张图
            else phi[primes[j]*i]=phi[i]*(primes[j]-1);
        }
    }
}

int main(){
    int n;
    cin>>n;

    get_eulers(n);

    ll res=0;
    for(int i=1;i<=n;i++){
        res+=phi[i];
    }
    cout<<res<<endl;
}