【数论】欧拉函数（基本性质、递推法、公式法、线性筛法）

欧拉函数

• 欧拉函数 $\phi (n)$ ：不超过 $n$ 且与 $n$ 互素的正整数的个数，$n\in N^{\ast }$
• $p$ 是素数$\iff \phi (p)=p-1$
• $p$ 是素数，$a\in N^{\ast }\Rightarrow \phi \left(p^a\right)=p^a-p^{a-1}$
  • 证明：与$p^a$ 不互素的只有$p,2p,3p,\cdots ,p^{k-1}p$ ，等比求和减去即可
• 任意正整数 $n$ 的素数幂分解 $n=p_1^{a_1}\cdot p_2^{a_2}\cdots p_s{}^{a_s}\Rightarrow \phi (n)=n\cdot \left(1-\frac{1}{p_1}\right)\cdot \left(1-\frac{1}{p_2}\right)\cdots \left(1-\frac{1}{p_s}\right)$
  • 证明：上一定理还可写成$\phi \left(p^a\right)=p^a(1-\frac{1}{p})$ ，由欧拉函数是积性函数可得上式。
  • 推论：当 $n$ 为奇数时, 有 $\phi (2n)=\phi (n)$
• $n>2,n\in N^{\ast }\Rightarrow \phi (n)$ 是偶数
• $n\in N^{\ast },{\sum }_{d\mid n}\phi (d)=n$
• 欧拉定理： $m$ 是一正整数, $a$ 是一个整数且 $(a,m)=1$, 那么 $a^{\phi (m)}\equiv$ $1(\mathrm{m}\mathrm{o}\mathrm{d}m)$
• $x\mathrm{m}\mathrm{o}\mathrm{d}p=0\Rightarrow \phi (xp)=\phi (x)\cdot p$ ，证明如下：
  • $[1,x]$ 中与$x$ 不互质的数$y$ ，$y+x$ 与$x$ 仍不互质，$y+cx,c\le p$ 仍与$x$ 不互质，即$[1,x]$ 区间中与$x$ 不互质的数在$+cx$ 后仍不互质。前面素数性质可证$[1,x]$ 中与$x$ 互质的数在$+cx$ 后仍互质。所以可得上述结论。

注：以下均来自kuangbin的模板

分解质因数法

int Euler(int n)
{
    getFactors(n); // 前面提过的分解质因数
    int ret = n;
    for (int i = 0; i < fatCnt; i++)
        ret = ret / factor[i][0] * (factor[i][0] - 1);
    return ret;
}

递推法

将$\phi (i)$ 先置为$i$ ，把$ki$ 的原始$\phi$ 乘上这个$i$ 做出的贡献。

过程中，$\phi$ 没有赋值说明是素数

int euler[N];
void getEuler()
{
    memset(euler, 0, sizeof(euler));
    euler[1] = 1;
    for (int i = 2; i < N; i++)
        if (!euler[i])
            for (int j = i; j < N; j += i)
            {
                if (!euler[j])
                    euler[j] = j;
                euler[j] = euler[j] / i * (i - 1);
            }
}

求单个欧拉函数

找素因子，用公式

int euler(int n)
{
    int ans = n;
    for (int i = 2; i * i <= n; i++)
        if (n % i == 0)
        {
            ans -= ans / i;
            while (n % i == 0)
                n /= i;
        }
    if (n > 1)
        ans -= ans / n;
    return ans;
}

线性筛

bool check[N + 10];
int phi[N + 10];
int prime[N + 10];
int tot; //素数的个数
void phi_and_prime_table(int N)
{
    memset(check, false, sizeof(check));
    phi[1] = 1;
    tot = 0;
    for (int i = 2; i <= N; i++)
    {
        if (!check[i])
        {
            prime[tot++] = i;
            phi[i] = i - 1;
        }
        for (int j = 0; j < tot; j++)
        {
            if (i * prime[j] > N)
                break;
            check[i * prime[j]] = true;
            if (i % prime[j] == 0)
            {
                phi[i * prime[j]] = phi[i] * prime[j]; // 利用性质
                break;
            }
            else
                phi[i * prime[j]] = phi[i] * (prime[j] - 1); // 积性函数
        }
    }
}