acw数学知识

数论基础知识

1. 质数

质数：在大于1的整数中，只包含1和它本身这两个约数的数。
$n$中最多只包含一个大于$sqrt(n)$的质因子。

1.1 质数的判定

试除法 $O(sqrt(n))$

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

1.2 分解质因数

试除法 $O(log(n))$~$O(sqrt(n))$

void divide(int n)
{
	for(int i=2;i<=n/i;i++)
		if(n%i==0)//i一定是质数
		{
			int s=0;
			while(n%i==0)
			{
				n/=i;
				s++;
			}
			printf("%d^%d",i,s);
			if(n!=1) printf(" * ");
		}
	if(n>1) printf("%d^%d\n",n,1);
}

1.3 筛质数

1.3.1 朴素筛法

int prime[N],pNum;
bool p[N];
void find_prime(int n)
{
	for(int i=2;i<=n;i++)
	{
		if(!p[i])
		{
			prime[pNum++]=i;
			for(int j=i+i;j<=n;j+=i)
				p[j]=true;
		}
	}
}

1.3.2 线性筛（欧拉筛）

int prime[N],pNum;
bool p[N];
void find_prime(int n)
{
	for(int i=2;i<=n;i++)
	{
		if(!p[i])
			prime[pNum++]=i;
		for(int j=0;prime[j]<=n/i;j++)//从小到大枚举最小质数
		{
			p[prime[j]*i]=true;
			if(i%prime[j]==0) break;//prime[j]一定是i的最小质因子
		}
	}
}

2. 约数

2.1 试除法求所有约数

vector<int> get_divisors(int n)
{
	vector<int> res;
	for(int i=1;i<=n/i;i++)
	{
		if(n%i==0)
		{
			res.push_back(i);
			if(i!=n/i) res.push_back(n/i);
		}
	}
	sort(res.begin(),res.end());
	return res;
}
int main()
{
	int n;cin>>n;
	auto res=get_divisors(n);
	for(auto t:res) printf("%d ",t);
}

2.2 约数个数

对$N$质因数分解

$N$=$p_1^{{\alpha }_1}$$p_2^{{\alpha }_2}$$p_3^{{\alpha }_3}$…$p_n^{{\alpha }_n}$

其中$p$为质数，1~$N$对$p_i$的倍数有$p_i^0$,$p_i^1$…$p_i^{{\alpha }_i}$

约数个数=(${\alpha }_1$+1)(${\alpha }_2$+1)(${\alpha }_3$+1)…(${\alpha }_n$+1)

#include<bits/stdc++.h>
using namespace std;
typedef long long ll;
const int mod=1e9+7;
int main()
{
	int n;cin>>n;
	unordered_map<int,int> primes;
	for(int i=2;i<=n/i;i++)
	{
		while(n%i==0) 
		{
			n/=i;primes[i]++;
		}
	}
	if(n>1) primes[n]++;
	ll res=1;
	for(auto prime:primes) 
		res=res*(prime.second+1)%mod;
	printf("%lld\n",res);
	return 0;
}

2.3 约数之和

约数之和= ($p_1^0$+$p_1^1$+$p_1^2$…$p_1^{{\alpha }_1}$)($p_2^0$+$p_2^1$+$p_2^2$…$p_2^{{\alpha }_1}$) …($p_n^0$+$p_n^1$+$p_n^2$…$p_n^{{\alpha }_n}$）

$t$=1

$t$=$p_1$+1
…
$t$=$p_1$*( $p_1^{{\alpha }_1-2}$+…+ $p_1^0$ )+1

$t$=$p_1$* ( $p_1$*( $p_1^{{\alpha }_1-2}$+…+ $p_1^0$ )+1 ) +1

$t$=$p_1^{{\alpha }_1}$+…+ $p_1^2$+$p_1^1$+$p_1^0$

#include<bits/stdc++.h>
using namespace std;
typedef long long ll;
const int mod=1e9+7;
int main()
{
	int n;cin>>n;
	unordered_map<int,int> primes;
	for(int i=2;i<=n/i;i++)
	{
		while(n%i==0) 
		{
			n/=i;primes[i]++;
		}
	}
	if(n>1) primes[n]++;
	ll res=1;
	for(auto prime:primes) 
	{
		int p=prime.first,a=prime.second;
		ll t=1;
		while(a--) t=(t*p+1)%mod;
		res=res*t%mod;
	}	 
	printf("%lld\n",res);
	return 0;
}

2.4 最大公约数（欧几里得算法）

$d$可以整数$a$,$d$可以整数$b$,那么$d$可以整数$a\ast x+b\ast y$

0除以任何一个数都等于0，所以任何数都可以整除0

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

2.5 最小公倍数

int lcm(int a,int b)
{
	return a/gcd(a,b)*b;
}

3. 欧拉函数

3.1 定义

定义 1~$N$中与$N$互质的数的个数被称为欧拉函数
N= $p_1^{{\alpha }_1}$ $p_2^{{\alpha }_2}$ ··· $p_k^{{\alpha }_k}$

$\Phi (N)$=N(1-$\frac{1}{p_1}$)(1-$\frac{1}{p_2}$)···(1-$\frac{1}{p_k}$）

3.2 证明

3.2.1 证明

根据容斥原理
$\Phi (N)$=$N$
-$\frac{N}{p_1}$ -$\frac{N}{p_2}$ -···- $\frac{N}{p_k}$
+$\frac{N}{p_1\ast p_2}$+$\frac{N}{p_1\ast p_3}$ +···+$\frac{N}{p_{k-1}\ast p_k}$
-$\frac{N}{p_1\ast p_2\ast p_3}$ -$\frac{N}{p_1\ast p_2\ast p_4}$ -···-$\frac{N}{p_{k-2}\ast p_{k-1}\ast p_k}$
+$\frac{N}{p_1\ast p_2\ast p_3\ast p_4}$+$\frac{N}{p_1\ast p_2\ast p_3\ast p_5}$+··· $\frac{N}{p_{k-3}\ast p_{k-2}\ast p_{k-1}\ast p_k}$

3.2.2代码

#include<bits/stdc++.h>
using namespace std;
int main()
{
	int n;
	cin>>n;
	int res=n;
	for(int i=2;i<=n/i;i++)
	{
		if(n%i==0)
		{
			res=res/i*(i-1);//(a-1)/a 
			while(n%i==0) n/=i;
		}
		if(n>1) res=res/n*(n-1);
	}
	cout<<res<<endl;
	return 0;
}

3.3筛法求欧拉函数

3.3.1推导
（1）$i$为质数，$\varphi (i)$= $i$ - 1
（2）$i$ $mod$ $p_j$ $==0$ 时,$\varphi (p_j\ast i)$ = $p_j$ * $\varphi (i)$
（3）$i$ $mod$ $p_j$ $!=0$ 时，$\varphi (p_j\ast i)$= $\varphi (p_j)$ * $\varphi (i)$ = $\varphi (i)$ * ($p_j$-1)

3.3.2代码

#include<bits/stdc++.h>
using namespace std;
typedef long long ll;
const int N=1e6+10;
int primes[N],cnt;
int phi[N];
bool p[N];

ll get_eulers(int n)
{
	phi[1]=1;
	for(int i=2;i<=n;i++)
	{
		if(!p[i])
		{
			primes[cnt++]=i;
			phi[i]=i-1;
		}	
		for(int j=0;primes[j]<=n/i;j++)
		{
			p[primes[j]*i]=true;
			if(i%primes[j]==0) 
			{
				phi[primes[j]*i]=phi[i]*primes[j];
				break;	
			}
			phi[primes[j]*i]=phi[i]*(primes[j]-1);
		}
	}
	ll res=0;
	for(int i=1;i<=n;i++) res+=phi[i];
	return res;
}
int main()
{
	int n;cin>>n;
	cout<<get_eulers(n)<<endl;
}

4.欧拉定理

4.1 定理

欧拉定理（Euler Theorem，也称费马-欧拉定理或欧拉函数定理）：
设$a,m$ $\in$ $N^{+}$，且 $gcd(a,m)==1$，我们有
$a^{\varphi (m)}$ $\equiv 1$ $(mod$ $m$)

4.2 证明

$a,n$互质
$1-n$中，有$a_1$,$a_2$···,$a_{\varphi (n)}$这$\varphi (n)$个数与n互质
那么，a*$a_1$ , a*$a_2$ ,···, a*$a_{\varphi (n)}$也与n互质

( a$a_1$ $mod$ $n$ ) (a$a_2$ $mod$ $n$)··· (a$a_{\varphi (n)}$ $mod$ $n$)
=( a$a_1$ * a$a_2$ ··· a$a_{\varphi (n)}$) $mod$ $n$
=$a^{\varphi (n)}$ ($a_1$$a_2$···$a_{\varphi (n)}$) $mod$ $n$

即 $a^{\varphi (n)}$ ($a_1$$a_2$···$a_{\varphi (n)}$) ≡ $a_1$$a_2$···$a_{\varphi (n)}$ （$mod$ $n$ ）
所以，$a^{\varphi (n)}$ ≡ 1 ($mod$ $n$ )

4.3 简单推论

费马小定理： 当n是质数时，$a^{\varphi (p)}$ ≡ 1 (mod p) ，即$a^{p-1}$ $\equiv 1$ $(mod$ $p)$

5. 快速幂

5.1 推导

${\alpha }^k$
=${\alpha }^{2^{x_1}}$ * ${\alpha }^{2^{x_2}}$ ··· ${\alpha }^{2^{x_k}}$
=${\alpha }^{2^{x_1}+2^{x_2}+\cdot \cdot \cdot +2^{x_t}}$
即k=$2^{x_1}$+$2^{x_2}$+···+$2^{x_t}$

${\alpha }^1$
${\alpha }^{2^1}$=${\alpha }^{2^0}$ * ${\alpha }^{2^0}$
${\alpha }^{2^2}$=${\alpha }^{2^1}$ * ${\alpha }^{2^1}$
···
${\alpha }^{2^{lo{g}_{2}k}}$=${\alpha }^{2^{lo{g}_{2}k-1}}$ * ${\alpha }^{2^{lo{g}_{2}k-1}}$

例：$4^5$=$4^{(101{)}_2}$=$4^0$ * $4^2$
$a=4,k=5,p=10,res=1$
$res=res\ast a$ $mod$ $p;$
$k>>=1;$
$a=a\ast a$ $mod$ $p$;

5.2 代码

typedef long long ll;
//返回a^k mod p
int qmi(int a,int k,int p)
{
	int res=1;
	while(k)
	{
		if(k&1) res=(ll)res*a%p; 
		k>>=1;
		a=(ll)a*a%p;
	}
	return res;
}

6. 逆元

6.1 乘法逆元

定义 若整数$b,m$互质，并且$b|a$,则存在一个整数$x$，使得$a/b$ $\equiv a\ast x$ $(mod$ $m)$，则称$x$为$b$的模$m$乘法逆元，记为$b^{-1}$ ($mod$ $m$)。
$b$ 存在乘法逆元 $<==>$ $b$与$m$互质
根据费马小定理，当$m$为质数时，$b^{m-2}$即为$b$的乘法逆元。
$b\ast x$ $\equiv 1$ $(mod$ $m)$
$b^{m-1}$ $\equiv 1$ $(mod$ $m)$
$b\ast b^{m-2}$ $\equiv 1$ $(mod$ $m)$
$-->$ $x$ $=$ $b^{m-2}$

7. 扩展欧几里得

7.1 裴蜀定理

定理 若$a,b$是整数，且 $gcd(a,b)=d$，那么对于任意的整数$x,y$，$a\ast x+b\ast y$ 都一定是$d$的倍数。特别地，一定存在整数$x,y$，使$a\ast x+b\ast y=d$ 成立。
证明
因为$d|a,d|b,$ 所以$d|(ax+by)$

7.2 扩展欧几里得

证明
$a\ast x_1+b\ast y_1=gcd(a,b)$
$b\ast x_2+(a$ $mod$ $b)\ast y_2$ $=gcd(b,a$ $mod$ $b)$
因为 $gcd(a,b)$ $=$ $gcd(b,$ $a$ $mod$ $b)$
其中，$a$ $mod$ $b$ $=a-$ $\lfloor a/b\rfloor$$\ast b$
那么，$b\ast x_2+$(a $mod$ $b)\ast y_2$
$=b\ast x_2+$(a-$\lfloor \frac{a}{b}\rfloor$$\ast b)\ast y_2$
$=a\ast y_2$$+$$b\ast (x_2$ -$\lfloor \frac{a}{b}\rfloor$*$y_2$$)$

推得
$x_1=y_2$
$y_1=$ $x_2$ $-$ $\lfloor \frac{a}{b}\rfloor$*$y_2$

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; 
}
int main()
{
	int a,b,x,y;
	scanf("%d%d",&a,&b);
	exgcd(a,b,x,y);
	printf("%d %d\n",x,y);
	return 0;
}

8. 线性同余方程

给定$a,b,m$，求$x$，使得$a\ast x\equiv b$ $(mod$ $m)$
解
$\exists$ $y$ $\in$ $Z$ ，$st.$ $a\ast x\equiv$ $m\ast y+b$
$a\ast x-m\ast y=b$
令$y^{\prime }=-y$
根据裴蜀定理可知
$a\ast x+m\ast y^{\prime }=b$ 当且仅当$gcd(a,m)|b$时有解
解为 $x=x\ast$ $\frac{b}{gcd(a,m)}$ $mod$ $m$

9.中国剩余定理（孙子定理）

9.1 定理说明

$$\{\begin{array}{l}x\equiv a_1(mod{m}_1)\\ x\equiv a_2(mod{m}_2)\\ \cdot \cdot \cdot \\ x\equiv a_n(mod{m}_n)\end{array}$$

假定整数$m_1,m_2\cdot \cdot \cdot m_n$两两互质，则对任意的整数：$a_1,a_2\cdot \cdot \cdot a_n$,方程组$S$有解。
方程的解为：
$$x=\sum _{i=1}^n{a}_i\ast t_i\ast M_i$$
其中，
$M={\prod }_{i=1}^n{m}_i$
$M_i=\frac{M}{m_i}$
$t_i=m_i^{-1}(mod{M}_i)$