数论学习笔记 欧拉函数（一些性质和运用）+杜教筛

本文应sharpland男神要求转载一下这位大佬的欧拉函数博客！！
原文传送门

定义

在数论中，对正整数n，欧拉函数是小于等于n的数中与n互质的数的数目。并且用符号$\varphi(n)$表示一个整数的欧拉函数。例如$\varphi(8)=4$。$\varphi(1)=1$。

一些欧拉函数的性质

性质一

对于一个质数n，$\varphi(n)=n-1$。
证明：因为n是质数。

性质二

若$n=p^k$，则$\varphi (n)=p^k-p^{k-1}=(p-1)\times p^{k-1}$
证明：
因为除了p的倍数以外，其他数都跟n互质。

性质三

当$gcd(n,m)=1$时，$\varphi (nm)=\varphi(n)*\varphi(m)$
证明：$\varphi(n)$是积性函数。

性质四

设$n=p_1^{k_1}\times p_2^{k_2}\times …… \times p_m^{k_m}$，则$\varphi (n)=n\times (1-\frac{1}{p_1})\times (1-\frac{1}{p_2})\times ……\times (1-\frac{1}{p_m})$
证明：
根据性质二

$$\begin{align} \ \varphi(n)& =\prod (p_i-1)p_i^{k_i-1}(p_i|n)\\ & = n\prod (1-\frac 1{p_1})\\ \end{align}$$
当然也可以用容斥的想法去理解。

性质五

欧拉定理：对于互质的整数$a,m$有$a^{\varphi (m)}\equiv 1(mod m)$
证明：
这小于n且与n互质的集合时$Z$显然$|Z|=\varphi(n)$，$Z={p_1,p_2…p_{\varphi(n)}}$。
令集合$S=\{a\times p_1 \mod n,a\times p_2 \mod n… a\times p_{\varphi (n)}\mod n\}$。
因为：

1. 因为$a$与$n$互质，$p_i$与$n$互质，所以$a\times p_i$与$n$互质，所以$a\times p_1\mod n ∈Z$
2. 若$i≠j$，那么$a\times p_i \mod n≠a\times p_j \mod n$
   反证：
   假如$a\times p_i \mod n=a\times p_j \mod n$，设$a\times p_i=k_i\times n +b$
   那么
   $$\begin{align} a\times p_i & = k_i\times n+b=k_j\times n+b\\ a\times (p_i-p_j) & =n\times (k_i-k_j)\\ \end{align}$$
   因为$a$与$n$互质，即$n|(p_i-p_j)$，不成立。
   所以$S=Z$
   由此我们可以列出等式：
   $$\begin{align} a\times p_1 \times a \times p_2…\times a \times p_{\varphi(n)}& \equiv p_1\times p_2\times …p_{\varphi(n)}(\mod n)\\ a^{\varphi (n)}& \equiv 1(\mod n) \end{align}$$
   延伸：
   费马定理：如果正整数$a$与$p$互质，则$a^{p-1}\equiv 1(\mod p)$。
   证明：由性质一可得$\varphi (p)=p-1$，代入欧拉定理即可。

性质六

设小于n的所有与n互质的数的和为$Sum，Sum=n\times \frac{\varphi (n)}{2}$
证明：
首先证明一个结论：如果$gcd(n,i)=1$则$gcd(n,n-i)=1$。
反证法:
如果存在$k≠1$使$gcd(n,n-i)=k$那么
$(n-i)\mod k=0$，$n \mod k=0$
可以得到$i \mod k=0$，即$gcd(n,i)=k$，也就是说如果$gcd(n,i)=1$，那么$gcd(n,n-i)$就不能大于1。
那么就可以得知与n互质的数都是成对存在的，并且和为n，那么就可以得出
$Sum=n\times \frac{\varphi(n)}{2}$的公式。

性质七

首先$p$是个质数。如果$i\mod p=0$，那么$\varphi (i\times p)=p\times \varphi(i)$（结论一），否则$\varphi (i\times p)=\varphi(i)\times (p-1)$（结论二）。
证明：
对于第一个结论我们只需证明$gcd(n,m)=1$可以得出$gcd(n,m+n)=1$即可。
反证法：
假设$gcd(n,m+n)=b(b≠1)$。设$n+m=k_1b,m=k_2b$

$$k_2b+n=k_1b$$
$$n=(k_1-k_2)b$$
所以 $gcd(n,m)$至少等于 $b$。
得证。
对于第二个结论，我们可知由于$gcd(i,p)=1，\varphi (i\times p)=\varphi(p)\times \varphi(i)$，并且 $\varphi(p)=p-1$（通过性质一和性质三），得证。

性质八

直接给式子吧…
$n=\sum_{d|n}\varphi(d)$
根据上面那条式子可以继续推点显而易见的东西 。

$$\sum_{i=1}^{n}i=\sum_{i=1}^n\sum_{d|i}\varphi(d)=\sum_{d=1}^n\varphi(d)\times \lfloor \frac{n}{d}\rfloor$$
反演一下， $\varphi(n)=\sum_{d|n}\mu \times \frac{n}{d}$

信息学中的应用

应用一：线筛$\varphi$函数

这个算法可以让我们从$O(n)$的时间得出$\varphi(1)$到$\varphi(n)$的值。
首先我们要用到几条上面提到的：性质一，性质七。
那么根线筛素数的原理一样（积性函数）我们只需根据欧拉函数的性质来进行线筛。

//YxuanwKeith
void Getphi(int Max) {
    phi[1] = 1;
    for (int i = 2; i <= Max; i ++) {
        if (!Flag[i]) {
            phi[i] = i - 1; // i是质数，第一种情况。
            pri[++ pri[0]] = i;
        }
        for (int j = 1; j <= pri[0]; j ++) {
            if (1ll * i * pri[j] > Max) break;
            Flag[i * pri[j]] = 1; //筛质数。
            if (i % pri[j] == 0) {
                phi[i * pri[j]] = phi[i] * pri[j];//（i % pri[j] = 0），第二种情况
                break;
            }
            phi[i * pri[j]] = phi[i] * (pri[j] - 1); //（i % pri[j] ！= 0），第三种情况
        }
    }
}

应用二：$O(\sqrt n)$得到$\varphi(n)$

直接根据性质四，$O(\sqrt n)$的枚举$n$的所有质因子，然后直接套用公式算，如果$n$特别大的话还可以选择用$millerrabin+pollarrho$来找质数。

//YxuanwKeith
long long Getphi(long long n) {
    long long phi = n;
    for (long long i = 2; i * i <= n; i ++) {
        if (n % i == 0) {
            phi /= i;
            phi *= i - 1;
            while (n % i == 0) n /= i;
        }
    }
    if (n != 1) phi /= n, phi *= n - 1;
    return phi;
}

应用三：杜教筛求$\varphi$的前缀和

即求$\sum_{i=1}^n\varphi(i)$，$n\leq10^{10}$
51nod题目链接
由于本文着重讨论欧拉函数的性质，所以就不细讲杜教筛相关内容，有兴趣的可以通过下面的blog了解：
杜教筛
基本原理：假设我们要计算$S(n)=\sum^n_{i=1}f(i)$

$$\begin{align} \sum_{i=1}^n(f\times g)(i)& =\sum_{i=1}^n\sum_{d|i}g(d)f(\frac{i}{d}) \\ &= \sum_{d=1}^n g(d)\sum_{1\leq i\leq n,d|i}f(\frac id) \\ &=\sum_{d=1}^ng(d)\sum_{1\leq i\leq \lfloor \frac nd\rfloor}f(i)\\ &=\sum_{d=1}^ng(d)S(\lfloor \frac nd\rfloor) \end{align}$$
$$∴g(1)S(n)=\sum_{i=1}^n(f\times g)(i)-\sum_{i=2}^ng(i)S(\lfloor \frac ni\rfloor )$$
根据杜教筛的应用，设 $f(i)=\varphi(i)$，我们要找到一个函数 $g(i)$使得 $g(i)$的前缀和和 $f(i)\times g(i)$的前缀和都很好求。根据性质八，不难发现 $f(i)\times l=id(l函数满足各项都为1，id函数满足id(i)=i)$
设 $\phi(n)=\sum_{i=1}^n\varphi(i)$，根据杜教筛:
$$\begin{align} \phi(n)&=\sum_{i=1}^n\varphi(i)\\ &=\sum_{i=1}^n\biggl( i-\sum_{d|i,d<i}\varphi(d)\biggl)\\ &=\frac{n\times(n+1)}{2}-\sum_{i=2}^n\sum_{d|i,d<i}\varphi(d)\\ &=\frac{n\times(n+1)}{2}-\sum_{\frac id=2}^n\sum_{d=1}^{\lfloor \frac{n}{\frac id}\rfloor}\varphi(d)\\ &=\frac{n\times(n+1)}{2}-\sum_{i=2}^n\sum_{d=1}^{\lfloor \frac ni\rfloor}\varphi(d)\\ &=\frac{n\times(n+1)}{2}-\sum_{i=2}^n\phi(\lfloor \frac ni\rfloor) \end{align}$$
然后递推进去算就可以了，直接这样算是 $O(n^{\frac34})$但是由于 $\varphi$是积性函数，可以预处理出前 $n^{\frac23}$的 $\phi$然后再做杜教筛，复杂度就优化成了 $O(n^{\frac23})$

//YxuanwKeith
//51nod1239
#include <cstring>
#include <cstdio>
#include <algorithm>
#include <cmath>

using namespace std;

typedef long long LL;

const int MAXN = 5e6 + 5, MAXM = 1e6 + 5;
const int Mo = 1e9 + 7;

LL n, que[MAXM];
int m, num, pri[MAXN], sum[MAXN], phi[MAXN], ids[MAXM], idl[MAXM], s[MAXM], inv[MAXN];
bool flag[MAXN];

int power(int x, int y) {
    int ans = 1;
    for (; y; y >>= 1, x = 1ll * x * x % Mo)
        if (y & 1) ans = 1ll * ans * x % Mo;
    return ans;
}

void prepare() {
    phi[1] = 1;
    for (int i = 2; i < MAXN; i ++) {
        if (!flag[i]) pri[++ pri[0]] = i, phi[i] = i - 1;
        for (int j = 1; j <= pri[0] && 1ll * i * pri[j] < MAXN; j ++) {
            int x = i * pri[j];
            flag[x] = 1;
            if (i % pri[j] == 0) {
                phi[x] = phi[i] * pri[j];
                break;
            }
            phi[x] = phi[i] * (pri[j] - 1);
        }
    }
    for (int i = 1; i < MAXN; i ++) sum[i] = (sum[i - 1] + phi[i]) % Mo;
}

int main() {
    scanf("%lld", &n);
    prepare();
    m = sqrt(n);
    for (LL l = 1; l <= n; l ++) {
        LL d = n / l, r = n / d;
        que[++ num] = d;
        if (d <= m) ids[d] = num; else idl[l] = num;
        l = r;
    }
    int inv = power(2, Mo - 2);
    for (int i = num; i; i --) {
        LL x = que[i];
        if (x < MAXN) s[i] = sum[x]; else {
            int y = x % Mo;
            s[i] = 1ll * y * (y + 1) % Mo * inv % Mo;
            for (LL l = 2; l <= x; l ++) {
                LL p = x / l, r = x / p;
                p = (p <= m) ? ids[p] : idl[n / p];
                (s[i] -= 1ll * (r - l + 1) * s[p] % Mo) %= Mo;
                l = r;
            }

        }
    }
    printf("%d\n", (s[1] + Mo) % Mo);
}

另外还有一些欧拉函数的变形也可以通过杜教筛来算：
1.$f(i)=\varphi(i)\times i$，$f(i)\times id=id^2$

$$\sum_{d|n}d\times \varphi(d) \times \frac nd=n\times \sum_{d|n}\varphi(d)=n^2$$
这个可以扩展到 $f(i)=\varphi(i)\times i^k$， $f(i)\times id^k=id^{k+1}$