质数判定
对于一个数是否是质数,只用看它能否被 2 → N 2\to \sqrt N 2→N中任意的整数整除即可。
inline bool is_prime(int n){
if(n<2)return false;
for(int i=2;i<=sqrt(n);i++)
if(!(n%i))return false;
return true;
}
质因数分解
引入:算术基本定理
任何一个大于1的正整数都能唯一分解为有限个质数的乘积,可以写作:
N
=
p
1
c
1
p
2
c
2
p
3
c
3
.
.
.
p
m
c
m
N=p_1^{c_1}p_2^{c_2}p_3^{c_3}...p_m^{c_m}
N=p1c1p2c2p3c3...pmcm
其中
c
i
c_i
ci都是正整数,
p
i
p_i
pi都是质数,且满足
p
1
<
p
2
<
.
.
.
<
p
m
p_1\lt p_2\lt ...\lt p_m
p1<p2<...<pm。
方法:试除法
结合质数判定法与埃氏筛法,我们可以扫描
2
→
N
2\to \sqrt N
2→N的每一个数
d
d
d,如果
d
d
d能整除
N
N
N,则从
N
N
N中除掉所有的因子
d
d
d,同时累积除去
d
d
d的个数。时间复杂度为
O
(
s
q
r
t
(
N
)
)
O(sqrt(N))
O(sqrt(N))。
特别的,若
N
N
N没有被任何
2
→
N
2\to \sqrt N
2→N的数整除,则
N
N
N为质数,无需分解。
inline void divide(int n){
int m=0;
for(int i=2;i<=sqrt(n);i++){
if(n%i==0){
p[++m]=i;c[m]=0;
while(n%i==0)n/=i,c[m]++;
}
}
if(n>1)p[++m]=n;c[m]=1;
for(int i=1;i<=m;i++)cout<<p[i]<<"^"<<c[i]<<'\n';
}