【算法竞赛学习笔记】数论（一）-数学提升计划

---

title : LINNO’s数论笔记（一）
date : 2021-8-4
tags : ACM,数论
author : Linno

---

整除

性质

$$\begin{gathered} 1.如果a|b且b|c，那么a|c。 \\ 2.a|b且a|c等价于任意的整数x和y,有a|(b\ast x+c\ast y) \\ 3.设m\ne 0,那么a|b等价于(m\ast a)|(m\ast b) \\ 4.设整数x和y满足下式：a\ast b+b\ast y=1,且a|n、b|n，那么(a\ast b)|n \end{gathered}$$

奇妙的性质

1.如果2能整除a的末一位，则2|a

2.如果4能整除a的末两位，则4|a

3.如果8能整除a的末三位，则8|a

4.如果3能整除a的各位数字之和，则3|a

5.如果11能整除a的偶数位数字之和与奇数位数字之和的差，则11|a

6.如果一个数的末三位与末三位前面的数字组成的数之差能被7、11、13整除，那么这个数就能够被7、11、13整除。

二元一次不定方程

一般形式$ax+by=c$，此方程有整数解的充要条件是GCD(a,b)|c
$$\begin{gathered} 设x_0,y_0是该方程的一组整数解，那么该方程的所有整数解可表示为: \\ x=x_0+\frac{b}{GCD(a,b)}t,y=y_0-\frac{a}{GCD(a,b)}t \end{gathered}$$

同余

性质

自反性：$a\equiv a(\mathrm{m}\mathrm{o}\mathrm{d}m)$

对称性：$若a\equiv b(\mathrm{m}\mathrm{o}\mathrm{d}m)，则b\equiv a(\mathrm{m}\mathrm{o}\mathrm{d}m)$

传递性：$若a\equiv b(\mathrm{m}\mathrm{o}\mathrm{d}m),b\equiv c(\mathrm{m}\mathrm{o}\mathrm{d}m)，则a\equiv c(\mathrm{m}\mathrm{o}\mathrm{d}m)$

同加性：$若a\equiv b(\mathrm{m}\mathrm{o}\mathrm{d}m),则a+c\equiv b+c(\mathrm{m}\mathrm{o}\mathrm{d}m)$

同乘性：$若a\equiv b(\mathrm{m}\mathrm{o}\mathrm{d}m)，则a\ast c\equiv b\ast c(\mathrm{m}\mathrm{o}\mathrm{d}m)，若a\equiv b(\mathrm{m}\mathrm{o}\mathrm{d}m)，c\equiv d(\mathrm{m}\mathrm{o}\mathrm{d}m),则a\ast c\equiv b\ast d(\mathrm{m}\mathrm{o}\mathrm{d}m)$

同幂性：$若a\equiv b(\mathrm{m}\mathrm{o}\mathrm{d}m)，则a^n\equiv b^n(\mathrm{m}\mathrm{o}\mathrm{d}m)$

推论1：$a\ast b\mathrm{m}\mathrm{o}\mathrm{d}k=(a\mathrm{m}\mathrm{o}\mathrm{d}k)\ast (b\mathrm{m}\mathrm{o}\mathrm{d}k)\mathrm{m}\mathrm{o}\mathrm{d}k$

推论2：$$\begin{gathered} 若a\mathrm{m}\mathrm{o}\mathrm{d}p=x,a\mathrm{m}\mathrm{o}\mathrm{d}q=x,p,q互质，则一定存在则整数s,t，使得a=s\ast p+x,a=t\ast q+x, \\ 所以,s\ast p=t\ast q,则一定存在整数r，使s=r\ast q, \\ 所以，a=r\ast p\ast q+x，得出a\mathrm{m}\mathrm{o}\mathrm{d}p\ast q=x \end{gathered}$$

最大公约数

二进制算法（非递归求GCD）

inline int gcd(int x,int y){
	int i,j;
	if(x==0) return y;
	if(y==0) return x;
	for(i=0;(x&1)==0;i++) x>>=1; //去掉所有的2
	for(j=0;(y&1)==0;j++) y>>=1; //去掉所有的2
	if(j<i) i=j;
	while(1){
		if(x<y) x^=y^=x^=y;//若x<y，交换x,y
		if((x-=y)==0) return y<<i;
		//若x==y，gcd==x==y（就是在辗转减，while(1)控制）
		while((x&1)==0) x>>=1; //去掉所有的2 
	}
}

扩展欧几里得

int ex_gcd(int a,int b,int &x,int &y){
	if(!b){
		x=1;y=0;
		return a;
	}
	int d=ex_gcd(b,a%b,y,x);
	tmp=x;
	x=y;
	y=tmp-(a/b)*y;
	return d;
}

求解线性同余方程

//用扩展欧几里得算法解线性方程：ax+by=c 
bool linearEquation(int a,int b,int c,int &x,int &y){
	int d=ex_gcd(a,b,x,y);
	if(c%d) return false; //如果gcd(a,b)|c才有解 
	int k=c/d;
	x*=k; //+ t*b;
	y*=k; //-t*a;
	//求的只是其中一个解
	return true; 
}

快速幂

int fpow(int a,int p,int mod){ //计算a的p次方在模mod下的值
    int res=1;
    while(p){ //将a^p分治成 a^(p/2)*a^(p/2)，当p为奇数时乘上a
        if(p&1) res=res*a%mod;
        a=a*a%mod;
        p>>=1;
    }
    return res;
}

逆元

$a\ast b\equiv 1(\mathrm{m}\mathrm{o}\mathrm{d}b),a,b互质，则称x为a的逆元，记为a^{-1}$

求逆元的四种方法

①拓展欧几里得

②快速幂+费马小定理

③递推

inv[1]=1;
for(ull i=2;i<=n;i++){
	inv[i]=(ull)(p-p/i)*inv[p%i]%p;
}

唯一分解定理

任意正整数都有且只有一种方式写出其素因子的乘积表达式。

$A=p_1^{k_1}\ast p_2^{k_2}\ast ...\ast p_n^{k_n}，其中p_i均为素数$

约数和方程

对于已经分解的整数$A=p_1^{k_1}\ast p_2^{k_2}\ast ...\ast p_n^{k_n}$

有A的所有因子之和为$S=(1+p_1+p_1^2+...+p_1^{k_1})\ast (1+p_2+p_2^2+...+p_2^{k_2})\ast ...\ast (1+p_n+p_n^2+...+p_n^{k_n})$

中国剩余定理(CRT)

“三人同行七十稀，五树梅花廿一枝，七子团圆月正半，除百零五便得知。”

首先引入同余方程
$$\{\begin{array}{l}x\equiv 2\mathrm{m}\mathrm{o}\mathrm{d}3\\ x\equiv 3\mathrm{m}\mathrm{o}\mathrm{d}5\\ x\equiv 2\mathrm{m}\mathrm{o}\mathrm{d}7\end{array}$$
小模数分别为5*7=35，3*7=21，3*5=15，找乘法逆元
$$\{\begin{array}{l}35a\equiv 1\mathrm{m}\mathrm{o}\mathrm{d}3\\ 21b\equiv 1\mathrm{m}\mathrm{o}\mathrm{d}5\\ 15c\equiv 1\mathrm{m}\mathrm{o}\mathrm{d}7\end{array}$$

我们可以解得逆元a=2,b=1,c=1，于是有下列式子
$$\{\begin{array}{l}x\equiv 70\mathrm{m}\mathrm{o}\mathrm{d}105\\ x\equiv 21\mathrm{m}\mathrm{o}\mathrm{d}105\\ x\equiv 15\mathrm{m}\mathrm{o}\mathrm{d}105\end{array}$$

这里70=2*5*7，21=3*7，15=3*5（都是不同方程里面的素因子）

我们可以求得$2\ast 70+3\ast 21+2\ast 15\mathrm{m}\mathrm{o}\mathrm{d}105=23$为x的解

定义

$$\begin{gathered} 设m_1,m_2,...,m_r两两互素，并记N=m_1\ast m_2\ast ...\ast m_r，则同余方程 \\ \{\begin{array}{l}x\equiv b_1\mathrm{m}\mathrm{o}\mathrm{d}{m}_1\\ x\equiv b_2\mathrm{m}\mathrm{o}\mathrm{d}{m}_2\\ ...x\equiv b_r\mathrm{m}\mathrm{o}\mathrm{d}{m}_r\end{array} \\ 在模N同余的意义下有唯一解。 \end{gathered}$$

$$由于m_i两两互质，令x=(N/m_i)\ast y，方程组等同于解同余方程(N/m_i)y\equiv 1\mathrm{m}\mathrm{o}\mathrm{d}{m}_i，得到特解y_i，则方程组的解为：x_0=b_1{x}_1+b_2{x}_2+...+b_r{x}_r\mathrm{m}\mathrm{o}\mathrm{d}N,在模N意义下唯一。$$

模板

#include<bits/stdc++.h>
using namespace std;
typedef long long ll;
//P1495 【模板】中国剩余定理(CRT)/曹冲养猪

long long n,mul=1,m[16],a[16],mi[16],x,ans=0;

void exgcd(ll a,ll b,ll &x,ll &y){
	if(!b){
		x=1;y=0;
		return;
	}
	exgcd(b,a%b,x,y);
	ll tmp=x;
	x=y;
	y=tmp-a/b*y;
}

int main(){
	scanf("%lld",&n); //同余方程个数
	for(int i=1;i<=n;i++){
		scanf("%lld%lld",&x,&a[i]); //余数和模数
		m[i]=x; 
		mul*=x; //记录大模数N
	}
	for(int i=1;i<=n;i++){
		mi[i]=mul/m[i]; //
		ll x,y;
		exgcd(mi[i],m[i],x,y);
		ans+=a[i]*mi[i]*(x<0?x+m[i]:x);
	}
	printf("%lld",ans%mul);
	return (0^0);
}

总结步骤

求大模数M

对于每个小模数p=M/m,求模m意义下的逆元i，那么p*i*a就是满足方程的最小数

每个方程的最小数相加模M，就是方程组最小解。

其他

**完全数（完美数）：**全部因数之和等于他本身

**盈数：**全部因数之和大于他本身

**亏数：**全部因数之和小于他本身

**亲和数：**一个属的真因子的和等于另一个数

参考代码

luoguP3518 [POI2011]SEJ-Strongbox:https://www.luogu.com.cn/record/54555918

UVA374 Big Mod：https://www.luogu.com.cn/record/54557937

UVA11105 H-半素数 Semi-prime H-numbers：https://www.luogu.com.cn/record/54559621

UVA756 Biorhythms：https://blog.nowcoder.net/n/037e60b2ee124f6286e2b105b7e4c9bf

参考资料

信息学奥赛之数学一本通 （林厚从）

https://www.bilibili.com/video/BV1o5411T7Np?from=search&seid=3644460168410828006