算法数论笔记

---

数论

数论都是针对整数而言

倍数、因数、整除

对于整数$a、b$，存在整数$q$满足a=bq

则称$a$是$b$的倍数、$b$是$a$的因数、$b|a$

最大公因数 记作：$gcd（a,b）$ 或者$（a,b）$

最大的$d$满足 $d|a$和 $d|b$

求最大公因数=>辗转相除法

设$a,b,c$是任意三个不全为0的整数，满足存在整数$q$，$a=bq+c$，那么$(a,b)=(b,c)$。

代码如下 时间复杂度$O(log(min(a，b)))$

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

扩展欧几里得

若$a,b$是任意两个不全为$0$的整数，则存在整数$x,y$，
使得 $ax+by=(a,b)$。

int exgcd(int a,int b,int &x,int &y)
{
    if(b==0) return x=1,y=0,a;
    int r=exgcd(b,a%b,x,y)
    y = y-(a/b)*x;
    return r
}

注意：

1. 满足$ax＋by=(a,b)$的$(x,y)$有无限多组，假设$(x_0,y_0)$是其中一组，则所有解为：$$(x_0+k\frac{b}{(a,b)},y_0-k\frac{a}{(a,b)}),(k\in Z)$$

   口诀 : 前加后减,k倍相反除$gcd$

2. $ax+by=1⟺(a,b)=1$

stein算法(求超过$longlong$范围的最大公因数)

• 如果$a=0$或$b=0$，则$(0,b)=b,(a,0)=a$ （）

• 如果$a,b$ 都是偶数，则$(a,b)=(a/2,b/2)$

• 如果$a$是偶数，b是奇数，则$(a,b)=(a/2,b)$

• 如果$a$是奇数，b是偶数，则$(a,b)=(a,b/2)$

• 如果$a，b$都是奇数，不妨设$a>b$，则$(a，b)=(b，a-b)$

  总结：为0等于其他，偶数除2，奇数取模变相减

**最小公倍数 **

记作：$lcm(a,b)$ 或者 [a,b]

最小的$m$满足 $a|m$ 和 b|m

两个数相乘除最大公因数

int lcm(int a, int b)
{ 
	return a / gcd(a, b) * b; //这里是为了防止溢出 先除后乘
}

互质

两个整数最大公因数是$1$ 即$gcd(a,b)=1$

具有的性质 a、b互质 a| bc 则 a | c

质数（素数）、合数

设p为正整数，p >1，如果p 的正因数只有1及它本身，则称p为质数（素数)，否则称p为合数。

算数基本定理

任意大于1的整数a能够唯一写成 $a={p_1}^{{\alpha }_1}{p_2}^{{\alpha }_2}...{p_k}^{{\alpha }_k}$ ，其中$p_i$是质数，

${\alpha }_1>0,i=1..k$,且$p_i<p_j(i<j)$

埃氏筛法

用质数把质数的倍数筛掉

代码如下 时间复杂度 $O(nlog(log(n)))$非常接近$O(n)$

int pr[N],cnt;
bool st[N]; //st[i] 表示i不是质数
void get_pr(int n)
{
    for(int i=2;i<=n;i++)
    {
        if(st[i]) continue; //如果这个数不是质数跳过
        pr[cnt++]=i; //是的话加入质数数中
        for(int j=i+i;j<=n;j+=i) st[j]=true; //用这个质数筛其他数
    }
}

欧拉筛

每个合数只需被其最小质因子筛掉

代码如下 时间复杂度$O(n)$

int pr[N],cnt;
bool st[N];

for(int i=2;i<=n;i++)
{
    if(!st[i]) pr[cnt++]=i;//如果这个数是质数
    for(int j=0;pr[j]<=n/i;j++)
    {
        st[pr[j]*i]=true;
        if(i%pr[j]==0) break;
    }
}

证明：

假设要筛的数为$x$， $x={p_1}^{{\alpha }_1}{p_2}^{{\alpha }_2}...{p_k}^{{\alpha }_k}$

当枚举到质数$ p_1$ 时候 $i={p_1}^{{\alpha }_1-1}{p_2}^{{\alpha }_2}...{p_k}^{{\alpha }_k}$

当枚举到质数$ p_2$ 时候 $i={p_1}^{{\alpha }_1}{p_2}^{{\alpha }_2-1}...{p_k}^{{\alpha }_k}$ 但是当$i\%pr[j]==0$ 就break掉了不可能存在这种情况 所以只会被最小质因子筛掉

质因数分解

求一个数n的所有因数

vector<int> fac;
for(int i=1;i<=n/i;i++)
{
    if(n%i==0)
    {
        fac.push_back(i);
        if(n/i != i) fac.push_back(n/i);
        
    }
}

求一个数n的所有质因数

vector<int> pfac;
for(int i=1;i<=n/i;i++)
{
    if(n%i==0)
    {
        pfac.push_back(i);
        while(n%i==0) n/=i;
    }
}
if(n>1) pfac.push_back(n)//上面把所有小于等于根号n的质因子筛掉了 
    //最多剩下一个大于根号n的质因子

同余

设$a,b,m$为三个整数，其中$m$ >0，如果满足$a\%m=b\%m$，则称a与b对模m同余，记作$a\equiv b(modm)$。

等价类

令集合 i = { x ∣ x ≡ i ( m o d    m ) } i = \{x | x \equiv i (\mod m) \} i={x∣x≡i(modm)} ,$(i\in \mathbb{Z})$

假设$j=i+m$,则$j=i$

等价类上的运算

在${\mathbb{Z}}_m=0,1,...,m-1$上定义两种运算$+,\times$,其中

$A+B=C$当且仅当$\forall a\in A,b\in B,a+b\in C$

$A\times B=C$ 当且仅当$\forall a\in A,b\in B,a\times b\in C$

逆元

如果$a,b\in {\mathbb{Z}}_m$,满足$ab=1$,则称$b$是$a$的逆元,记作$a^{-1}$

判定在${\mathbb{Z}}_m$里，a是否有逆元 如果有，如何求$a^{-1}$

• $ab=km+1$ ,扩展欧几里得算法
• $a$在${\mathbb{Z}}_m$内有逆元，当且仅当$(a,m)=1$。
• 如果$m$为质数，则任意$0<a<m$，$a\in {\mathbb{Z}}_m$有逆.

欧拉函数

定义欧拉函数$$\phi (n)$$为正整数$n$与序列$1,2,...,n-1$,$n$中互质的数的个数

求欧拉函数

设$n={p_1}^{{\alpha }_1}{p_2}^{{\alpha }_2}...{p_k}^{{\alpha }_k}$ ,则$\phi (n)=(p_1-1){p_1}^{{\alpha }_1-1}(p_2-1){p_2}^{{\alpha }_2-1}...(p_k-1){p_k}^{{\alpha }_k-1}$

ll get_phi(int n)
{
	ll phi =1;
	for(int i=2;i<=n/i;i++)
	{
		if(n%i==0)
		{
			phi *= (i-1);
			n/=i;
			while(n%i==0) phi *=i,n/=i;
		}
	}
	if(n>1) phi*=n;
	return phi;
}

求$$\phi (1)..\phi (n)$$ 欧拉筛的基础上做修改

int pr[N],cnt;
int phi[N];
bool st[N];

for(int i=2;i<=n;i++)
{
    if(!st[i]) pr[cnt++]=i;phi[i]=i-1//如果这个数是质数
    for(int j=0;pr[j]<=n/i;j++)
    {
        st[pr[j]*i]=true;
        if(i%pr[j]==0) 
        {
        	phi[i*pr[j]]=phi[i]*pr[j];
        	break;
        }
        phi[i*pr[j]] = phi[i]*(pr[j]-1);
    }
}

欧拉定理

如果$(a,m)=1$，那么$a^{\phi (m)}=1(\mathrm{m}\mathrm{o}\mathrm{d}m)$。

扩展欧拉定理

$a^c=a^{c\%\phi (m)+\phi (m)}(\mathrm{m}\mathrm{o}\mathrm{d}m)$ 当$c>=\phi (m)$

费马小定理

设$p$为质数，如果$p\times a$，那么$a^{P-1}=1(modp)$。

设$m$为正整数，$a\in {\mathbb{Z}}_m,(a,m)=1$，逆元$a^{-1}=a^{\phi (m)}-1$

设p为质数，$a\in {\mathbb{Z}}_p,0<a<p，$逆元$a^{-1}=a^{p-2}$

同余方程组

设$m_1,m_2,..,m_k$是两两互质的$k$个正整数,求方程组$\{\begin{array}{rl}& x\equiv a_1(\mathrm{m}\mathrm{o}\mathrm{d}{m}_1)\\ & x\equiv a_2(\mathrm{m}\mathrm{o}\mathrm{d}{m}_2)\\ .& \\ .& \\ & x\equiv a_k(\mathrm{m}\mathrm{o}\mathrm{d}{m}_k)\end{array}$ 的解$x$

先假设$k=2$的情况下 $x=y_1{m}_1+a_1=y_2{m}_2+a_2$

根据扩展欧几里得可得$x=x_0+k\frac{m_1{m}_2}{[m_1,m_2]}$

有解当且仅当$(m_1,m_2)|a_1-a_2$

解的形式为$x\equiv x_0\mathrm{m}\mathrm{o}\mathrm{d}[m_1,m_2]$

然后每两项合并,$x\equiv x_0\mathrm{m}\mathrm{o}\mathrm{d}[m_1,m_2..m_k]$

中国剩余定理(孙子定理)

$x\equiv {\sum }_{i=1}^k{M_i}^{-1}{M}_i{a}_i(\mathrm{m}\mathrm{o}\mathrm{d}M)$

其中$M=m_1{m}_2..m_k,M_i=M/m_i$

计数问题解决方法

• 直接计算
• dp
• 容斥
• 生成函数
• 其他的定理

加法原理和乘法原理

加法原理:做某件事情有几种选择，每种选择的方案数之和就是做这件事情的方案数

乘法原理:做某件事情分为几步，每步的方案数是独立的，则它们的积就是做这件事情的方案数