质数

质数的定义

质数又称作素数.
质数是除了1和本身没有其他因数的数.
与之相反的是合数.
1不属于质数或合数.

质数的判断

假设我们要判断$n$是不是质数.
从2枚举$n-1$,若能整除$n$,$n$就不是质数.
时间复杂度$O(n)$.
但因数是成对的(除了$\sqrt{n}$)
令$n=c_1{c}_2$,假设$c_1\le c_2$
则$c_1\le \sqrt{n},c_2\ge \sqrt{n}$
所以只需枚举到$\sqrt{n}$.

bool isprime(long long x) {
    for(int i=2; i<=sqrt(x); i++) {
        if(x%i==0) return false;
    }
    return true;
}

时间复杂度$O(\sqrt{n})$.

质数筛法

现在我们要求$1$到$n$中全部质数.

暴力

最简单的办法就是判断每个数是否为质数.
时间复杂度$O(n\sqrt{n})$.

埃氏筛法

从$1$到$\sqrt{n}$,若这个数没有被标记,即为质数.
只要发现了一个质数$t$,就把从$1$到$n$中全部$t$的倍数标记.
时间复杂度$O(n\log n)$

欧拉线性筛法

对于每个$i(1<i<n)$,筛去所有小于$i$的质数与$i$的乘积.
不难发现,有些合数被重复筛去.
设所有小于$i$的质数分别为$p_1,p_2,...,p_j$.
当$i\%p_k=0$时,停止筛$i$的倍数,即可避免重复.
因为当$l>k$时,所有$i\ast p_l$都会被比$i$大的数与$p_k$的乘积筛去.

for(int i=2;i<=n;i++){
    if(!vis[i]) {
        cnt++;
        pri[cnt]=i;
    }
    for(int j=1;i*pri[j]<=n;j++){
        vis[i*pri[j]]=true;
        if(i%pri[j]==0) break;
    }
}

时间复杂度$O(n)$

质因数分解

一个数的质因数必定小于等于$\sqrt{n}$.
只需把$i$从$2$枚举到$\sqrt{n}$,从$n$中去掉所有$i$并做记录,就能求出$n$的质因数.

for(int i=2;i<=sqrt(n);i++) {
    while(n%i==0) {
        n/=i;
        cnt[i]++;
    }
}

时间复杂度$O(\sqrt{n})$

欧拉函数

欧拉函数是从$1$到$n-1$的数中与$n$互质的数的个数.
$\phi (1)=\phi (2)=1$.

性质

1. 对于质数$n$,$\phi (n)=n-1$.
2. 对于质数$p$,$\phi (p^k)=(p-1)\ast p^{k-1}$.
3. 对于互质的$n,m$,$\phi (n)\ast \phi (m)=\phi (n\ast m)$.
4. 对于奇数$n$,$\phi (2\ast n)=\phi (n)$.
5. 对于$n(n>1)$,${\sum }_{k=1}^{n}k(gcd(n,k)=1)=\phi (n)\ast n/2$.
6. 若$k\%p=0$,$\phi (k\ast p)=\phi (k)\ast p$.
7. 若$k\%p!=0$,$\phi (k\ast p)=\phi (k)\ast (p-1)$.
8. 对于互质的$a,m$,$a^{\phi (p)}\equiv 1(\mathrm{m}\mathrm{o}\mathrm{d}p)$.
9. ${\sum }_{k\mid n}^n\phi (k)=n$
10. $\phi (n)={\sum }_{k\mid n}^n\phi (k)\ast \frac{n}{k}$.

计算公式及证明

设$n$的质因数分解为$p_1^{c_1}\ast p_2^{c_2}\ast$…$\ast p_i^{c_i}$.
则$\phi (n)=n\ast \frac{p_1-1}{p_1}\ast \frac{p_2-1}{p_2}\ast$…$\ast \frac{p_i-1}{p_i}$
证明:设$p,q$是$n$的质因子.
则$1$-$n$中,$p$的倍数有$\frac{n}{p}$个,$q$的倍数有$\frac{n}{q}$个,$p\ast q$的倍数有$\frac{n}{p+q}$个.
所以$1$到$n-1$中与$p,q$互质的数的个数为
KaTeX parse error: No such environment: align* at position 7: \begin{̲a̲l̲i̲g̲n̲*̲}̲ n-\frac{n}{p}-…
其他质因子同理.

算法实现

对于求单个欧拉函数,分解质因数即可.

for(long long i=2;i<=sqrt(num);i++) {
    bool div=false;
    while(num%i==0) {
        div=true;
        num/=i;
    }
    if(div) ans=ans*(i-1)/i;
}
if(num>1) ans=ans*(num-1)/num;

时间复杂度与质因数分解一样.
但如果要求从$1$到$n$的数的欧拉函数,就可以通过打表法来求解.

for(int i=2;i<=MAXN;i++) {
    if(!phi[i]) {
        for(int j=i;j<=MAXN;j+=i) {
            if(!phi[j]) phi[j]=j;
            phi[j]=phi[j]/i*(i-1);
        }
    }
}

或是在欧拉筛时同时求出欧拉函数.

for(int i=2;i<=n;i++){
    if(!vis[i]) {
        cnt++;
        pri[cnt]=i;
        phi[cnt]=i-1;
    }
    for(int j=1;i*pri[j]<=n;j++){
        vis[i*pri[j]]=true;
        if(i%pri[j]==0) {
            phi[i*pri[j]]=phi[i]*pri[j];
            break;
        }
        phi[i*pri[j]]=phi[i]*(pri[j]-1);
    }
}