AcWing数学知识

数学知识

质数

试除法$O(\sqrt{n})$

public static boolean isPrime(int x) {
   
    if (x < 2) return false;
    for (int i = 2; i <= x / i; i ++) {
   
        if (x % i == 0) return false;
    }
    return true;
}

分解质因数$O(\sqrt{n})$

public static void divide(int x) {
   
    for (int i = 2; i <= x / i; i ++) {
   
        if (x % i == 0) {
   
            int s = 0;
            while (x % i == 0) {
   
                x /= i;
                s ++;
            }
            System.out.println(i + " " + s);
        }
    }
    if (x > 1) System.out.println(x + " 1");
    System.out.println();
}

筛质数

// 埃氏筛法 $O(nloglogn)$
public static int getPrimes(int x) {
   
    for (int i = 2; i <= x; i ++) {
   
        if (!st[i]) {
   
            prime[++ cnt] = i;
            for (int j = i + i; j <= x; j += i) {
   
                st[j] = true;
            }
        }
    }
    return cnt;
}

// 线性筛法 $O(n)$， n = 1e7时比埃氏筛法快一倍
// 算法核心：x仅会被其最小质因子筛去
import java.util.*;

class Main {
   
    static final int N = 1000010;
    static int[] primes = new int[N];
    static boolean[] st = new boolean[N];
    static int cnt = 0;
    
    public static int getPrimes(int x) {
   
        for (int i = 2; i <= x; i ++) {
    // 循环遍历所有数
            if (!st[i]) {
   
                primes[cnt ++] = i;
            }
            for (int j = 0; primes[j] <= x / i; j ++) {
   
                // 对于任意一个合数x，假设pj为最小质因子，当i<x/pj时，pj*i一定会被筛掉
                st[primes[j] * i] = true;
                if (i % primes[j] == 0) break;
            	/*
            	1. i%pj == 0，pj一定为i的最小质因子，pj也一定为i*pj的最小质因子
            	2. i%pj != 0, pj一定小于i的所有质因子，所以pj也为pj*i的最小质因子
            	*/
            }
        }
        return cnt;
    }
    
    public static void main(String[] args) {
   
        Scanner sc = new Scanner(System.in);
        int n = sc.nextInt();
        int res = getPrimes(n);
        System.out.println(res);
    }
}

约数

试除法求约数

class Main {
   
    
    public static void main(String[] args) {
   
        Scanner sc = new Scanner(System.in);
        int n = sc.nextInt();
        List<Integer> temp;
        while (n -- > 0) {
   
            int x = sc.nextInt();
            temp = getDivisors(x);
            for(int i : temp) {
   
                System.out.print(i + " ");
            }
            System.out.println();
        }
    }
    
    public static List<Integer> getDivisors(int n) {
   
        ArrayList<Integer> res = new ArrayList<Integer>();
        
        for (int i = 1; i <= n / i; i ++) {
   
            if (n % i == 0) {
   
                res.add(i);
                if (i != n / i) {
   
                    res.add(n / i);
                }
            }
        }
        
        Collections.sort(res);
        return res;
    }
}

约数个数

KaTeX parse error: No such environment: align at position 8: \begin{̲a̲l̲i̲g̲n̲}̲ N = &p_1^{\al…

import java.util.*;

class Main {
   
    public static void main(String[] args) {
   
        Scanner sc = new Scanner(System.in);    
        int n = sc.nextInt();
        HashMap<Integer, Integer> map = new HashMap<Integer, Integer>();
        long ans = 1;
        long mod = 1000000007;
        while (n -- > 0) {
   
            int x = sc.nextInt();
            for (int i = 2; i <= x / i; i ++) {
   
                while(x % i == 0) {
   
                    x /= i;
                    map.put(i, map.getOrDefault(i, 0) + 1);
                }
                
            }
            if (x > 1) {
   
                map.put(x, map.getOrDefault(x, 0) + 1);
            }
        }
        
        for (int i : map.values()) {
   
            ans = ans * (i + 1) % mod;    
        }
        System.out.println(ans);
    }
}

约数之和

$$\begin{gathered} N=p_1^{c_1}\times p_2^{c_2}\times ...\times p_k^{c_k} \\ 约数个数：(c_1+1)\times (c_2+1)\times ...\times (c_k+1) \\ 约数之和：({p_1}^0+{p_1}^1+...{p_1}^{c_1})\times ...\times ({p_k}^0+{p_k}^1+...{p_k}^{c_k}) \end{gathered}$$

while (pow -- > 0) {
   
    temp = (temp * key + 1) % mod
}

$t=t\ast p+1$
$t=1$
$t=p+1$
$t=p^2+p+1$
$t=p^3+p^2+p+1$
$...$
$t=p^{pow}+p^{pow-1}+...+1$

import java.util.*;

class Main {
   
    public static void main(String[] args) {
   
        Scanner sc = new Scanner(System.in);
        int n = sc.nextInt();
        HashMap<Integer, Integer> map = new HashMap<Integer, Integer>();
        
        while (n -- > 0) {
   
            int x = sc.nextInt();
            for (int i = 2; i <= x / i; i ++) {
   
                while (x % i == 0) {
   
                    x /= i;
                    map.put(i, map.getOrDefault(i,0) + 1);
                }
            }
            if (x > 1) {
   
                map.put(x, map.getOrDefault(x, 0) + 1);
            }
        }
        
        long res = 1L;
        long mod = 1000000007L;
        for (int key : map.keySet()) {
   
            int pow = map.get(key);
            long temp = 1;
            for (int i = 0; i < pow; i ++) {
   
                temp = (temp * key + 1) % mod;
            }
            res = res * temp % mod;
        }
        
        System.out.println(res);
    }
}

最大公约数

$a$和$b$的最大公约数，记为$\mathrm{g}\mathrm{c}\mathrm{d}(a,b)$，$a$和$b$的最小公倍数，记为$\mathrm{l}\mathrm{c}\mathrm{m}(a,b)$
$$\forall a,b\in \mathbb{N},\mathrm{g}\mathrm{c}\mathrm{d}(a,b)\ast \mathrm{l}\mathrm{c}\mathrm{m}(a,b)=a\ast b$$

设$d=\mathrm{g}\mathrm{c}\mathrm{d}(a,b),a_0=a/d,b_0=b/d$。根据最大公约数的定义，有$\mathrm{g}\mathrm{c}\mathrm{d}(a_0,b_0)=1$（反证法证明，假设$\mathrm{g}\mathrm{c}\mathrm{d}(a_0,b_0)=x>1,y=a_0/x,z=b_0/x$，则存在$d_0=d\ast x$使得$y=a/d_0,z=b/d_0,d_0=d\ast x>d$，与$d$是$a,b$的最大公因数相矛盾）。根据最下公倍数的定义，有$\mathrm{l}\mathrm{c}\mathrm{m}(a_0,b_0)=a_0\ast b_0$
于是$\mathrm{l}\mathrm{c}\mathrm{m}(a,b)=\mathrm{l}\mathrm{c}\mathrm{m}(a_0\ast d,b_0\ast d)=\mathrm{l}\mathrm{c}\mathrm{m}(a_0,b_0)\ast d=a_0\ast b_0\ast d=a\ast b/d$。

更相减损术

$\forall a,b\in \mathbb{N},a\ge b$，有$\mathrm{g}\mathrm{c}\mathrm{d}(a,b)=\mathrm{g}\mathrm{c}\mathrm{d}(b,a-b)=\mathrm{g}\mathrm{c}\mathrm{d}(a,a-b)$

对于$a,b$的任意公约数$d$，因为$d|a,d|b$，所以$d|(a-b)$。因此$d$也是$b,a-b$的公约数。

欧几里得算法（辗转相除法）

$$\forall a,b\in \mathbb{N},b\ne 0,\mathrm{g}\mathrm{c}\mathrm{d}(a,b)=\mathrm{g}\mathrm{c}\mathrm{d}(b,a\mathrm{m}\mathrm{o}\mathrm{d}b)$$

若$a<b$，则${\rm gcd}(b, a\ {\rm mod}\ b) = {\rm gcd} (b,a) = {\rm gcd}(a,b) $，命题成立若$a\ge b$，不妨设$a=q\ast b + r$，其中$0\le r<b$。显然$r=a\mathrm{m}\mathrm{o}\mathrm{d}b$。对于$a,b$的任意公约数$d$，因为$d|a,d|q\ast b$，故$d|(a-qb)$，即$d$也是$b,r$的公约数。反之亦成立。因此$a,b$的公约数集合与$b,a\mathrm{m}\mathrm{o}\mathrm{d}b$的公约数集合相同。
于是它们的最大公约数自然也相同。

任何数都可以整除$0$，因此$b=0$时，$\mathrm{g}\mathrm{c}\mathrm{d}(a,b)=a$。

import java.util.*;

class Main {
   
    static int n;
    
    public static int gcd(int a, int b) {
   
        return b != 0 ? gcd(b, a % b) : a;
    }
    
    public static void main(String[] args) {
   
        Scanner sc = new Scanner(System.in);
        n = sc.nextInt();
        while (n -- > 0) {
   
            int a = sc.nextInt();
            int b = sc.nextInt();
            System.out.println(gcd(a, b));
        }
    }
}

欧拉函数

欧拉函数

欧拉函数
对于正整数$n$,欧拉函数是小于或等于$n$的正整数中与$n$互质的数的数目，记作$\phi (n)$。$\phi (1)=1$

求n的欧拉值
首先， 欧拉函数是一个积性函数，当$m,n$互质时，$\phi (mn)=\phi (m)\ast \phi (n)$根

据唯一分解定理知$n=p_1^{a_1}\ast p_2^{a_2}\ast ...\ast p_x^{a_x}$

因此$\phi (n)=\phi (p_1^{a_1})\ast ...\ast \phi (p_x^{a_x})$

对于任意一项$\phi (p_s^{a_s})=p_s^{a_s}-p_s^{a_s-1}$

从定义出发$\phi (p_s^{a_s})$等于小于或等于$p_s^{a_s}$的正整数中与$p_s^{a_s}$互质的数的数目

从$1$到$p_s^{a_s}$中共有$p_s^{a_s}$个数字

其中与$p_s^{a_s}$不互质的有$p_s,2p_s,...,p_s^{a_s-1}\ast p_s$共$p_s^{a_s-1}$项

所以 $\phi (p_s^{a_s})=p_s^{a_s}-p_s^{a_s-1}=$