数论-欧拉函数

欧拉函数

定义：

从$1$~$n$中与$n$互质的数的个数，$a,b$互质即$gcd(a,b)=1$

公式：

$$N=p_1^{\phi 1}\ast p_2^{\phi 2}\ast ......\ast p_k^{\phi k}$$
$$\phi (n)=N\ast (1-\frac{1}{p_1})\ast (1-\frac{1}{p_2})\ast ......\ast (1-\frac{1}{p_k})$$
证明：基于容斥原理，详细太复杂可看别的博客

复杂度：

单次$\sqrt{n}$,可线性求1-n的欧拉函数

应用：

欧拉定理：
若$a$与$n$互质，则$a^{\phi (n)}\equiv 1(\mathrm{m}\mathrm{o}\mathrm{d}n)$
推论(费马定理)：
若p为质数时
$$a^{\phi (p)}\equiv 1(\mathrm{m}\mathrm{o}\mathrm{d}p)$$
$$a^{p-1}\equiv 1(\mathrm{m}\mathrm{o}\mathrm{d}p)$$

板子：

//单次
#include<bits/stdc++.h>
using namespace std;
int main(){
    int t;cin>>t;while(t--){
        int n;
        cin>>n;
        int res=n;
        for(int i=2;i<=n/i;i++){
            if(n%i==0){
                res=res/i*(i-1);
                while(n%i==0) n/=i;
            }
        }
        if(n>1) res=res/n*(n-1);
        
        cout<<res<<endl;
    }
}

//线性求1-n的欧拉函数
#include<bits/stdc++.h>
using namespace std;
const int N=1e6+10;
int primes[N], cnt;     // primes[]存储所有素数
int euler[N];           // 存储每个数的欧拉函数

bool st[N];         // st[x]存储x是否被筛掉（合数标记）
void get_eulers(int n)
{
    euler[1] = 1;
    for (int i = 2; i <= n; i ++ )
    {
        if (!st[i])
        {
            // 质数
            primes[cnt ++ ] = i;
            euler[i] = i - 1;
        }
        for (int j = 0; primes[j] <= n / i; j ++ )
        {
            int t = primes[j] * i;
            st[t] = true;
            if (i % primes[j] == 0)
            {
                euler[t] = euler[i] * primes[j];
                break;
            }
            euler[t] = euler[i] * (primes[j] - 1);
        }
    }
}

int main(int argc, char const *argv[]) {
    get_eulers(1e6);
    int n;
    cin>>n;
    long long res=0;
    for(int i=1;i<=n;i++){
        res+=euler[i];
    }
    cout<<res<<endl;
    return 0;
}