数论基础

质数

质数这东西让我联想到了原子,一样不可拆分,所以我们说质数具有原子性。
我们设 π ( x ) \pi(x) π(x)为不大于x的质数个数,可得 π ( x ) = O ( x ln ⁡ ( x ) ) \pi(x)=O(\frac x{\ln(x)}) π(x)=O(ln(x)x)

唯一分解定理

任何一个 N ( N > 1 ) N(N>1) N(N>1)能被质数分解,且只有一种分解方式,即
N = p 1 c 1 × p 2 c 2 × … … × p m c m ( N > 1 ) N=p_1^{c_1}\times p_2^{c_2}\times……\times p_m^{c_m}(N>1) N=p1c1×p2c2××pmcm(N>1)
其中 p i p_i pi 为素数, c i c_i ci 为正整数,
满足 p 1 < p 2 < … … < p m p_1<p_2<……<p_m p1<p2<<pm

筛质数

可得:对于 n n n,最多只有一个质因子大于 n \sqrt n n
证明:若两个质因子 a > n , b > n a>\sqrt n,b>\sqrt n a>n b>n
           ~~~~~~~~~~            a × b > n a\times b>n a×b>n(不成立),即原命题正确。

暴力试除法

根据上述结论我们对于一个数,只需要枚举1~ n \sqrt n n 即可,
这里就不给代码了,时间复杂度O( n n n\sqrt n nn )
埃式筛法
枚举到一个质数,我们只需要将它的倍数全部标记即可,没被标记的就是质数,这里给出核心代码。

int prime[100005],flag[100005],tot;
void find(int n)
{
	for(int i=2;i<=n;i++)
		if(!flag[i])
		{
			flag[i]=1;
			prime[++tot]=i;
			for(int j=2;j<=n/i;j++)
				flag[i*j]=1;
		}
	for(int i=1;i<=tot;i++)
		printf("%d ",prime[i]);
}

欧拉筛法(线性筛法)
埃式筛法已经接近 O ( n ) O(n) O(n)的复杂度了,但是优化的余地。
对于一个数,可能被多次标记,我们可以减少它们的次数,来达到优化的目的。
例如:30,它可以被分解成:
30 = 2 × 15 = 3 × 10 = 5 × 6 30=2\times15=3\times10=5\times6 30=2×15=3×10=5×6
根据上述的唯一分解定理,我们只需要用最小的质因子来分解。
例如:30分解成 30 = 2 × 15 30=2\times15 30=2×15

int prime[100005],flag[100005],tot;
void find(int n)
{
	for(int i=2;i<=n;i++)
	{
		if(!flag[i]) prime[++tot]=i;
		for(int j=1;j<=tot&&prime[j]*i<=n;j++)
		{
			flag[i*prime[j]]=1;
			if(i%prime[j]==0) break;
		}
	}
	for(int i=1;i<=tot;i++)
		printf("%d ",prime[i]);
}

整除

a ∣ b a∣b ab表示 a a a整除 b b b,其中 b b b a a a的倍数, a a a b b b的约数。

一些性质:
1、若 d ∣ a d∣a da,有 d ∣ k a d∣ka dka
2、对称性:若 a ∣ b , b ∣ a a∣b,b∣a abba,有 a = b a = b a=b
3、传递性:若 c ∣ b , b ∣ a c∣b,b∣a cbba,则 c ∣ a c∣a ca
4、若 c ∣ a , d ∣ b c∣a,d∣b cadb,则 c d ∣ a b cd∣ab cdab
5、若 c ∣ a , c ∣ b c∣a,c∣b cacb,则 c ∣ ( m a + n b ) c∣(ma+nb) c(ma+nb)

约数:
N = p 1 c 1 × p 2 c 2 × … … × p m c m ( N > 1 ) N=p_1^{c_1}\times p_2^{c_2}\times……\times p_m^{c_m}(N>1) N=p1c1×p2c2××pmcm(N>1)
其中 p i p_i pi为素数, c i c_i ci为正整数,
满足 p 1 < p 2 < … … < p m p_1<p_2<……<p_m p1<p2<<pm
则N 的正约数集合 { p 1 b 1 , p 2 b 2 , … … , p m b m } \{p_1^{b_1},p_2^{b_2},……,p_m^{b_m}\} {p1b1p2b2pmbm}
其中 0 ≤ b i ≤ c i 0\le b_i\le c_i 0bici
由此可推出如下两条公式

  1. N的约数个数: ∏ i = 1 m c i + 1 \prod^m_{i=1}c_i+1 i=1mci+1
  2. N的约数的和: ∏ i = 1 m ( ∑ j = 0 c j ( p i ) j ) \prod_{i=1}^m(\sum^{c_j}_{j=0}(p_i)^j) i=1m(j=0cj(pi)j)

最大公约数和最小公倍数

我们设 g c d ( a , b ) gcd(a,b) gcd(a,b) a 、 b a、b ab的最大公约数。
           l c m ( a , b ) ~~~~~~~~~~lcm(a,b)           lcm(a,b) a 、 b a、b ab的最下公倍数。

l c m ( a , b ) lcm(a,b) lcm(a,b) g c d ( a , b )      i n t e r e s t i n g gcd(a,b)\;\;interesting gcd(a,b)interesting的方法:
将a,b分解质因数,得:
a = p 1 a 1 × p 2 a 2 × … … × p m a m a=p_1^{a_1}\times p_2^{a_2}\times……\times p_m^{a_m} a=p1a1×p2a2××pmam
b = p 1 b 1 × p 2 b 2 × … … × p m b m b=p_1^{b_1}\times p_2^{b_2}\times……\times p_m^{b_m} b=p1b1×p2b2××pmbm
那么 g c d ( a , b ) = p 1 m i n ( a 1 , b 1 ) × p 2 m i n ( a 2 , b 2 ) × … … × p m m i n ( a m , b m ) gcd(a,b)=p_1^{min(a_1,b_1)}\times p_2^{min(a_2,b_2)}\times……\times p_m^{min(a_m,b_m)} gcd(a,b)=p1min(a1,b1)×p2min(a2,b2)××pmmin(am,bm)
       l c m ( a , b ) = p 1 m a x ( a 1 , b 1 ) × p 2 m a x ( a 2 , b 2 ) × … … × p m m a x ( a m , b m ) ~~~~~~lcm(a,b)=p_1^{max(a_1,b_1)}\times p_2^{max(a_2,b_2)}\times……\times p_m^{max(a_m,b_m)}       lcm(a,b)=p1max(a1,b1)×p2max(a2,b2)××pmmax(am,bm)

欧几里得算法
g c d ( a , b ) = g c d ( b , a   m o d   b ) gcd(a,b)=gcd(b,a\bmod b) gcd(a,b)=gcd(b,amodb)
我们可以根据这个求出 g c d ( a , b ) gcd(a,b) gcd(a,b)

更相减损法
g c d ( a , b ) = g c d ( b , a − b ) = g c d ( a , a − b ) gcd(a,b)=gcd(b,a-b)=gcd(a,a-b) gcd(a,b)=gcd(b,ab)=gcd(a,ab)
g c d ( 2 a , 2 b ) = 2 g c d ( a , b ) gcd(2a,2b)=2gcd(a,b) gcd(2a,2b)=2gcd(a,b)
根据上述可得 g c d ( a , b ) gcd(a,b) gcd(a,b)

互质
g c d ( a , b ) = 1 : a , b gcd(a,b)=1:a,b gcd(a,b)=1ab互质。
g c d ( a , b , c ) = 1 : a , b , c gcd(a,b,c)=1:a,b,c gcd(a,b,c)=1abc互质。
g c d ( a , b ) = g c d ( b , c ) = g c d ( a , c ) = 1 gcd(a,b)=gcd(b,c)=gcd(a,c)=1 gcd(a,b)=gcd(b,c)=gcd(a,c)=1:三者两两互质。

欧拉函数
欧拉函数 φ ( n ) \varphi(n) φ(n),表示 1   n 1~n 1 n中与 n n n互质的个数。
N = p 1 c 1 × p 2 c 2 × … … × p m c m N=p_1^{c_1}\times p_2^{c_2}\times……\times p_m^{c_m} N=p1c1×p2c2××pmcm
其中 p i p_i pi为素数, c i c_i ci 为正整数,
满足 p 1 < p 2 < … … < p m p_1<p_2<……<p_m p1<p2<<pm
φ ( n ) = N × p 1 − 1 p 1 × p 2 − 1 p 2 × … … × p m − 1 p m = N × ∏ ( 1 − 1 p ) \varphi(n)=N\times \frac{p_1-1}{p_1} \times \frac{p_2-1}{p_2}\times……\times\frac{p_m-1}{p_m}=N\times \prod(1-\frac1p) φ(n)=N×p1p11×p2p21××pmpm1=N×(1p1)

一些性质

  1. φ ( 1 ) = 1 \varphi(1)=1 φ(1)=1

  2. p是质数, φ ( p ) = p − 1 \varphi(p)=p-1 φ(p)=p1

  3. p是质数,n = p k p^k pk φ ( n ) = p k − p k − 1 = ( p − 1 ) p k − 1 \varphi(n)=p^k-p^{k-1}=(p-1)p^{k-1} φ(n)=pkpk1=(p1)pk1

  4. g c d ( a , b ) = 1 , φ ( a b ) = φ ( a ) × φ ( b ) gcd(a,b)=1,\varphi(ab)=\varphi(a)\times\varphi(b) gcd(a,b)=1φ(ab)=φ(a)×φ(b)(积性函数)。

  5. 对于质数p,有
          n % p = 0 , φ ( n p ) = φ ( n ) × p ~~~~~n\%p=0,\varphi(np)=\varphi(n)\times p      n%p=0φ(np)=φ(n)×p
          n % p ≠ 0 , φ ( n p ) = φ ( n ) × ( p − 1 ) = φ ( n ) × φ ( p ) ~~~~~n\%p\neq0,\varphi(np)=\varphi(n)\times(p-1)=\varphi(n)\times\varphi(p)      n%p=0φ(np)=φ(n)×(p1)=φ(n)×φ(p)

  6. 对于n>1, 1 − n 1-n 1n中与n互质的数的和为 n × φ ( n ) 2 n\times\frac{\varphi(n)}2 n×2φ(n)

  7. ∑ d ∣ n φ ( d ) = n \sum_{d|n}\varphi(d)=n dnφ(d)=n

欧拉筛求欧拉函数
根据上述性质:
每次枚举出一个质数i, p h i i = i − 1 phi_i=i-1 phii=i1
根据第5个性质,即可推出。

int prime[100005],flag[100005],phi[100005],tot;
void find(int n)
{
	phi[1]=1;
	for(int i=2;i<=n;i++)
	{
		if(!flag[i]) prime[++tot]=i,phi[i]=i-1;
		for(int j=1;j<=tot&&prime[j]*i<=n;j++)
		{
			flag[i*prime[j]]=1;
			if(i%prime[j]==0)
			{
				phi[i*prime[j]]=phi[i]*prime[j];
				break;
			}
			else phi[i*prime[j]]=phi[i]*phi[prime[j]];
		}
	}
	for(int i=1;i<=n;i++)
		printf("%d ",phi[i]);
}

同余

用m除整数a,b余数相同,称a和b对模m同余,记作 a ≡ b ( m o d m ) a\equiv b\pmod m ab(modm)

两个概念:

  1. 标准剩余系,每个同余系的任意两个整数模m同余。
  2. 既约剩余系,每个同余系的任意整数与m互质。

一些性质1:

  1. a ≡ b ( m o d m ) a\equiv b\pmod m ab(modm),当且仅当 m ∣ ( a − b ) m|(a-b) m(ab)
  2. a ≡ b ( m o d m ) a\equiv b\pmod m ab(modm),当且仅当最在整数k,使得 a = b + k m a=b+km a=b+km

一些性质2:
有整数a,b,c,自然数m,n

  1. 自反性: a ≡ a ( m o d m ) a\equiv a\pmod m aa(modm)
  2. 对称性:若 a ≡ b ( m o d m ) a\equiv b\pmod m ab(modm),则 b ≡ a ( m o d m ) b\equiv a\pmod m ba(modm)
  3. 传递性: a ≡ b ( m o d m ) , b ≡ c ( m o d m ) , a ≡ c ( m o d m ) a\equiv b\pmod m,b\equiv c\pmod m,a\equiv c\pmod m ab(modm)bc(modm)ac(modm)
  4. 同加性: a ≡ b ( m o d m ) , a + c ≡ b + c ( m o d m ) a\equiv b\pmod m,a+c\equiv b+c\pmod m ab(modm)a+cb+c(modm)
  5. 同乘性: a ≡ b ( m o d m ) , a c ≡ b c ( m o d m ) a\equiv b\pmod m,ac\equiv bc\pmod m ab(modm)acbc(modm)
  6. 同幂性: a ≡ b ( m o d m ) , a n ≡ b n ( m o d m ) a\equiv b\pmod m,a^n\equiv b^n\pmod m ab(modm)anbn(modm)
  7. a   m o d   p = x , a   m o d   q = x , a   m o d   ( p q ) = x a \bmod p=x,a\bmod q=x,a\bmod (pq)=x amodp=x,amodq=x,amod(pq)=x p , q p,q pq互质。

欧拉定理

g c d ( a , n ) = 1 gcd(a,n)=1 gcd(a,n)=1,则 a φ ( n ) ≡ 1 ( m o d n ) a^{\varphi(n)}\equiv1\pmod n aφ(n)1(modn)

在这里插入图片描述

费马小定理

若p为质数,a为整数,有 a p ≡ a ( m o d p ) a^p\equiv a\pmod p apa(modp)

证明:

  1. g c d ( a , p ) = 1 → a φ ( p ) ≡ 1 ( m o d p ) gcd(a,p)=1\to a^{\varphi(p)}\equiv1\pmod p gcd(a,p)=1aφ(p)1(modp)
       ~~    a p − 1 ≡ 1 ( m o d p ) → a p ≡ a ( m o d p ) a^{p-1}\equiv1\pmod p\to a^p\equiv a\pmod p ap11(modp)apa(modp)
  2. g c d ( a , p ) ≠ 1 gcd(a,p)\neq1 gcd(a,p)=1,即a是p的倍数。
    a p % p = 0 , a a^p\%p=0,a%p=0 ap%p=0a,即 a p ≡ a ( m o d p ) a^p\equiv a\pmod p apa(modp)

裴蜀定理(贝祖定理)

对于任意整数a,b,存在一组整数x,y,满足 a x + b y = g c d ( a , b ) ax+by=gcd(a,b) ax+by=gcd(a,b)

证明:
d = g c d ( a , b ) d=gcd(a,b) d=gcd(a,b),则 a = d a ′ , b = d b ′ a=da^{'},b=db^{'} a=da,b=db
a x + b y = d a ′ x + d b ′ y = d ( a ′ x + b ′ y ) = d ax+by=da^{'}x+db^{'}y=d(a^{'}x+b^{'}y)=d ax+by=dax+dby=d(ax+by)=d
因为 a ′ , b ′ , x , y a^{'},b^{'},x,y a,b,x,y都是整数,
所以只有 a x + b y = g c d ( a , b ) ax+by=gcd(a,b) ax+by=gcd(a,b)或是 g c d ( a , b ) gcd(a,b) gcd(a,b)得倍数时有解。

在这里插入图片描述

扩展欧几里得的通解

对于ax+by=gcd(a,b),通过扩欧得出一组通解 x 0 , y 0 x_0,y_0 x0,y0
得通解: x = x 0 + b ÷ g c d ( a , b ) × t x=x_0+b÷gcd(a,b)\times t x=x0+b÷gcd(a,b)×t
              y = y 0 − a ÷ g c d ( a , b ) × t ~~~~~~~~~~~~~y=y_0-a÷gcd(a,b)\times t              y=y0a÷gcd(a,b)×t

到这里数论的知识就差不多了吧(尽管我还没写中国剩余定理呢)!

  • 12
    点赞
  • 4
    收藏
    觉得还不错? 一键收藏
  • 7
    评论

“相关推荐”对你有帮助么?

  • 非常没帮助
  • 没帮助
  • 一般
  • 有帮助
  • 非常有帮助
提交
评论 7
添加红包

请填写红包祝福语或标题

红包个数最小为10个

红包金额最低5元

当前余额3.43前往充值 >
需支付:10.00
成就一亿技术人!
领取后你会自动成为博主和红包主的粉丝 规则
hope_wisdom
发出的红包
实付
使用余额支付
点击重新获取
扫码支付
钱包余额 0

抵扣说明:

1.余额是钱包充值的虚拟货币,按照1:1的比例进行支付金额的抵扣。
2.余额无法直接购买下载,可以购买VIP、付费专栏及课程。

余额充值