求欧拉函数

欧拉函数

一、欧拉函数简介

    欧拉函数：对于正整数n，欧拉函数是小于或者等于(等于好像没有用哎)n的正整数中与n互质的数的数目，即为$\varphi (n)$，其中，$\varphi (1)$ = 1 (人为规定)。

二、某数的欧拉函数

公式一 ：$\varphi (m\ast n)=\varphi (m)\ast \varphi (n);条件：m与n互质,即gcd(m,n)=1$

根据唯一分解定理：

$$n=p_1^{\alpha 1}\ast p_2^{\alpha 2}....\ast p_k^{\alpha k};$$

因此 $\varphi (n)=\varphi (p_1^{\alpha 1})\ast \varphi (p_2^{\alpha 2})\ast ...\ast \varphi (p_k^{\alpha k})$

   从定义出发，$\varphi (p_k^{\alpha k})$ = 小于等于$p_k^{\alpha k}$的整数与$p_k^{\alpha k}$互质的数的数目，从1 ~$p_k^{\alpha k}$中有$p_k^{\alpha k}$个数字，其中与$p_k^{\alpha k}$不互质的有$p_k,2p_k,...p_k^{\alpha k-1}{p}_k$，一共$p_k^{\alpha k-1}$项，所以有：

$$\varphi (p_k^{\alpha k})=p_k^{\alpha k}-p_k^{\alpha k-1}=p_k^{\alpha k}\ast (1-\frac{1}{p_k});$$

进一步，有：
$$\varphi (n)=n\ast \sum _{i=1}^k(1-\frac{1}{p_i});$$
代码如下：

#include <bits/stdc++.h>

using namespace std;

int main()
{
    int a; cin >> a;
    int res = a;
    for (int i = 2; i <= a / i; ++ i) 
    {
        if (a % i == 0) 
        {
            while (a % i == 0) a /= i;
            res = res / i * (i - 1);
        }
    }
    if (a > 1) res = res / a * (a - 1);
    cout << res <<  '\n';
    
    system("pause");
    
    return 0;
}

---

三、1~n的欧拉函数和

结合线性筛质数，不懂线性筛质数的请移步：

1. 如果p为质数，则$\varphi (p)=n-1$，因为质数p与$\le$ p的数均是互质的

2. 若m1为质数：

   2.1 m2 % m1 = 0，则$\varphi (m1\ast m2)=m1\ast \varphi (m2)$;

   (其实2.1中m1不为质数，等式也成立，但2.2不行)

   2.2 m2 % m1 $\ne 0,则\varphi (m1\ast m2)=\varphi (m1)\ast \varphi (m2)$;

3. 步骤：
   第一步，取质数，得到质数的欧拉函数$\varphi (p)=n-1$
   第二步将得到的质数当成筛选因子，得到$\varphi (i\ast p)$的值(根据2中的法则，将m2为i，m1为p)

总结公式：
$$\begin{gathered} \varphi (p)=p-1,p为质数 \\ \varphi (i\ast p)=\{\begin{array}{rl}& p\ast \varphi (i),i\%p=0\\ & \varphi (i)\ast \varphi (p),i\%p\ne 0\end{array} \end{gathered}$$
代码如下：

#include <bits/stdc++.h>

using namespace std;

typedef long long LL;

const int N = 1000010;
int primes[N], cnt; 
int phi[N];
bool st[N];
LL get_eulers(int n)
{
    phi[1] = 1;
    for (int i = 2; i <= n; i ++) 
    {
        if (!st[i]) 
        {
            primes[cnt ++] = i;
            phi[i] = 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);
        }
    }
    LL res = 0;
    for (int i = 1; i <= n; i ++) res += phi[i];
    return res;
}

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

    cout << get_eulers(n) << endl;
    return 0;
}

---

THE END…