数论基础知识

质数

1.质数的判定——试除法

$O(n)$

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

性质1：如果 $d|n$ ，那么$\frac{n}{d}|n$

因为一个数的约数是成对出现的，因此在判断约数的时候只需要枚举较小的约数。

$O(\sqrt{n})$

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

2.分解质因数——试除法

$O(n)$

for (int i = 2; i <= n; i ++)
    if (n % i == 0)
    {
        int s = 0;
        while (n % i == 0)
        {
            n /= i;
            s ++;
		}
        cout << i << ' ' << s << endl;
    }

优化：$n$中最多只包含一个大于$\sqrt{n}$的质因子
$O(\sqrt{n})$

for (int i = 2; i <= n / i; i ++)
    if (n % i == 0)
    {
        int s = 0;
        while (n % i == 0)
        {
            n /= i;
            s ++;
		}
        cout << i << ' ' << s << endl;
    }
if (n > 1) cout << n << ' ' << 1 << endl;

3.筛法

埃氏筛 $O(n\log \log n)$

质数定理：$[1,n]$ 中有$\frac{n}{\ln n}$ 个质数

int primes[N], cnt;     // primes[]存储所有素数
bool vis[N];         // st[x]存储x是否被筛掉

void get_primes(int n)
{
    for (int i = 2; i <= n; i ++ )
    {
        if (vis[i]) continue;
        primes[cnt ++ ] = i;
        for (int j = i + i; j <= n; j += i)
            vis[j] = true;
    }
}

欧拉筛(线性筛) $O(n)$

int primes[N], cnt;     // primes[]存储所有素数
bool vis[N];         // st[x]存储x是否被筛掉

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

约数

1.试除法求一个数的所有约数

vector<int> get_divisors(int x)
{
    vector<int> res;
    for (int i = 1; i <= x / i; i ++ )
        if (x % i == 0)
        {
            res.push_back(i);
            if (i != x / i) res.push_back(x / i);
        }
    sort(res.begin(), res.end());
    return res;
}

2.约数个数

若$x=p_1^{a_1}\ast p_2^{a_2}\ast p_3^{a_3}...\ast p_k^{a_k}$，则其约数个数为：$n=(a_1+1)\ast (a_2+1)\ast ...\ast (a_k+1)$

3.约数之和

若$x=p_1^{a_1}\ast p_2^{a_2}\ast p_3^{a_3}...\ast p_k^{a_k}$，则其约束之和为：$sum=(p_1^0+p_1^2+...+p_1^{a_1})\ast (p_2^0+...+p_2^{a_2})\ast ...\ast (p_k^0+...+p_k^{a_k})$

4.欧几里得算法(辗转相除法)

原理：$gcd(a,b)=gcd(b,a\%b)$

证明：$a\%b=a-\lfloor \frac{a}{b}\rfloor \ast b$， 令$c=\lfloor \frac{a}{b}\rfloor$，因此$a\%b=a-c\ast b$，即只需要证明$gcd(a,b)=gcd(b,a-c\ast b)$

①：左式，对于$\forall d，有d|a且d|b，则d|(a-c\ast b)$

②：右式，对于$\forall d，有d|b且d|(a-c\ast b)，则d|b$

因此$gcd(a,b)=gcd(b,a\%b)$

int gcd(int a, int b)
{
    return b ? gcd(b, a%b) : a;
}

欧拉函数

欧拉定理

求欧拉函数

定义

$\phi (n)$表示在$[1,n-1]$中与n互质的数的个数(互质：$gcd(a,b)=1$)

求法

设$x=p_1^{a_1}\ast p_2^{a_2}\ast ...\ast p_k^{a_k}$，则$\phi (x)=x\ast (1-\frac{1}{p_1})\ast (1-\frac{1}{p_2})...\ast (1-\frac{1}{p_k})$

证明

利用容斥原理证明：

模板

int phi(int x)
{
    int res = x;
    for (int i = 2; i <= x / i; i ++ )
        if (x % i == 0)
        {
            res = res / i * (i - 1);//先除再乘防止溢出
            while (x % i == 0) x /= i;
        }
    if (x > 1) res = res / x * (x - 1);
    return res;
}

筛法求欧拉函数(基于线性筛)

证明：

模板

int primes[N], cnt;     // primes[]存储所有素数
int euler[N];           // 存储每个数的欧拉函数
bool st[N];         // st[x]存储x是否被筛掉
void get_eulers(int n)
{
    euler[1] = 1;
    for (int i = 2; i <= n; i ++ )
    {
        if (!st[i])
        {
            primes[cnt ++ ] = i;
            euler[i] = i - 1;
        }
        for (int j = 0; primes[j] <= n / i; j ++ )
        {
            int t = primes[j] * i;
            st[t] = true;
            if (i % primes[j] == 0)
            {
                euler[t] = euler[i] * primes[j];
                break;
            }
            euler[t] = euler[i] * (primes[j] - 1);
        }
    }
}

快速幂

模板

求 m^k mod p，时间复杂度 O(logk)。

int qmi(int m, int k, int p)
{
    int res = 1 % p, t = m;
    while (k)
    {
        if (k&1) res = res * t % p;
        t = t * t % p;
        k >>= 1;
    }
    return res;
}

快速幂求逆元

乘法逆元定义

费马定理

参考：同余式，整除

扩展欧几里得算法

裴蜀定理

定义

对于任意正整数$a,b$，一定存在非零整数$x,y$使得 $ax+by=gcd(a,b)$

证明

扩展欧几里得算法

在扩展欧几里得算法中，我们得到的$x,y$不唯一，但是得到一组$x,y$的值之后，就可以表示出所有的$(x,y)$

// 求x, y，使得ax + by = gcd(a, b)
int exgcd(int a, int b, int &x, int &y)
{
    if (!b)
    {
        x = 1; y = 0;
        return a;
    }
    int d = exgcd(b, a % b, y, x);
    y -= (a/b) * x;
    return d;
}

高斯消元