[number theory]数论

一些数论函数

约数和函数（数论积性）
$f(n)=\sum_{d|n}d$
若$n=p_1^{a_1}p_2^{a_2}p_3^{a_3}...p_k^{a_k}$
则$f(n)=(1+p_1^2+p_1^3...+p_1^{a_1})(1+p_2^2+p_2^3...+p_2^{a_2})...(1+p_k^2+p_k^3...+p_k^{a_k})$
约数个数函数（数论积性）
$f(n)=\sum_{d|n}1$
若$n=p_1^{a_1}p_2^{a_2}p_3^{a_3}...p_k^{a_k}$
则$f(n)=(1+{a_1})(1+{a_2})...(1+a_k)$

分解质因数

朴素算法$O(\sqrt{N})$

for(LL i = 2; i * i <= n; i++)
    while (n % i == 0) n /= i;

Pollard Rho算法$O(n^{1/4})$

LL pollard_rho(LL n, LL c)  
{  
    LL i = 1, k = 2;  
    LL x = rand() % (n - 1) + 1;  
    LL y = x;  
    while(true)  
    {  
        i++;  
        x = (multi(x, x, n) + c) % n;  
        LL d = gcd((y - x + n) % n, n);  
        if(1 < d && d < n) return d;  
        if(y == x) return n;  
        if(i == k)  
        {  
            y = x;  
            k <<= 1;  
        }  
    }  
}  

组合数取模$\binom{n}{m}\mod p$

$Lucas$定理（$p$为质数）

$n=sp+q,m=tp+r$

$$\binom{n}{m}=\binom{sp+q}{tp+r}=\binom{s}{t}\binom{q}{r}\mod p$$
$O(log_{p}n*p)$计算组合数取模

$p$不为质数

首先我们有$p=p_1^{a_1}p_2^{a_2}p_3^{a_3}...p_k^{a_k}$
那我们将模数变为$p_{i}^{a_{i}}$再用中国剩余定理合并就可以了
问题转化为

$$\binom{n}{m}=\frac{n!}{m!(n-m)!} \mod p^M$$
令$n! = T * p^k$，其中$p不能整除T$，如果能够成功地将$n!$分解成这个样子，那么
$$\binom{n}{m}=\frac{n!}{m!(n-m)!} = T_{n}*T_{m}^{-1}*T_{n-m}^{-1}*p^{k_{n}-k_{m}-k_{n-m}} \mod p^M$$
而$T与p$互质，能求逆元，问题就可做了。
如何将$n!分解为T*p^k$呢？
首先我们可以将阶乘写成这个样子
$n!=1*2*...p*...*2p*...*kp*....n$ 然后我们就有
$n!=(1*2*...*(p-1)*(p+1)*...*(2p-1)*(2p+1)*...*(kp-1)*(kp+1)*...*n)*k!*p^k$
前面一坨东西算出来，$k!$是个子问题，我们在这一层成功提取了$p^k$，解决子问题之后累积起来就行了

中国剩余定理（$CRT$）

$x=a_{1} \mod m_{1}$
$x=a_{2} \mod m_{2}$
…
$x=a_{n} \mod m_{n}$

假设整数$m_{1},m_{2},...,m_{n}$互质，则对于任意$a_{1},a_{2},...,a_{n}$，方程组有解
$设M=m_{1}*m_{2}*...*m_{n}, M_{i}=M/m_{i}$
$设T_{i}=M_{i}^{-1}\mod m_{i}$
则$x=a_{1}t_{1}M_{1}+a_{2}t_{2}M_{2}+...a_{n}t_{n}M_{n}+kM$