第二十三章 数论——质数
- 一、什么是质数
- 二、质数的判断
- 1、试除法(朴素版)
- 2、试除法(优化版)
- 三、分解质因数
- 1、什么是质因数
- 2、算术基本定理
- 3、分解质因数
- (1)问题
- (2)思路
- (3)代码
- (4)分析
- n不断变化会漏解吗?
- i一定是质因数?
- (5)优化
一、什么是质数
质数也叫素数。如果一个大于等于2的整数,只能够被 1 1 1和它本身整除,那么这个数就叫做质数。反之,这个数就称作合数。
二、质数的判断
1、试除法(朴素版)
试除法其实就是从定义出发,假设一个数 p p p是质数,那么 p p p就只能被1和 p p p整除,那我就借此验证一下 p p p是否能够被 2 2 2 ~ ( p − 1 ) (p-1) (p−1)的数字整除。
bool prime(int u)
{
if(u<2)return false;
for(int i=2;i<u;i++)
{
if(u%i==0)return false;
}
return true;
}
但是这种算法的时间复杂度是 O ( n 2 ) O(n^2) O(n2)。这个时间复杂度是很高的。我们如何优化呢?
2、试除法(优化版)
假设 p p p是一个合数,那么 p = a ∗ b p=a*b p=a∗b,那么必定满足不等式: a ≤ p ≤ b a\leq \sqrt p \leq b a≤p≤b
证明:
如果 0 < a < p 0<a < \sqrt p 0<a<p 同时 0 < b < p 0<b < \sqrt p 0<b<p。那么 a ∗ b < p a*b < p a∗b<p。
如果 a > p a > \sqrt p a>p 同时 b > p b > \sqrt p b>p,那么 a ∗ b > p a*b > p a∗b>p。
因此,为了让 a ∗ b = p a*b=p a∗b=p,只能是: a ≤ p ≤ b a\leq \sqrt p \leq b a≤p≤b
由于,我们只是想判断一个数是不是质数,所以我们只需要找到一个 2 2 2 ~ ( p − 1 ) (p-1) (p−1)的数能够整除 p p p就行。因此我们只需要枚举 2 2 2 ~ p \sqrt p p。
bool prime(int u)
{
if(u<2)return false;
for(int i=2;i<=u/i;i++)
{
if(u%i==0)return false;
}
return true;
}
此时代码的时间复杂度就是 O ( n ) O(\sqrt n ) O(n),时间上进行了极大的优化。
三、分解质因数
1、什么是质因数
如果 a a a % b = 0 b=0 b=0,就说 b b b是 a a a的因数,如果这个因数又恰好是质数,那么就称 b b b是 a a a的质因数。
2、算术基本定理
任何一个大于1的自然数 N,如果N不为质数,那么N可以唯一分解成有限个质数的乘积。
证明(反证法)
假设存在一个合数,其无法被分解为有限个质数的乘积,我们将这些数字放在一个集合中,那么这个集合中必存在一个最小的符合上述条件的合数 n n n。
由于 n n n是合数,所以必定存在 n = a ∗ b n=a*b n=a∗b
由于 n n n无法被分解成有限个质数的乘积,那么 a a a和 b b b必定不同时为质数。
a a a和 b b b中至少存在一个合数,所以我们假设 a a a是合数, b b b是质数。
由于我们的 n n n是无法被分解成有限个质数中最小的一个,同时 a a a一定是小于 n n n的。所以 a a a必定是可以分解为有限个质数的乘积的。即, a = p 1 k ∗ p 2 m ∗ p 3 x . . . p n w a=p_1^k*p_2^m*p_3^x...p_n^w a=p1k∗p2m∗p3x...pnw
那么此时我们的 n = b ∗ p 1 k ∗ p 2 m ∗ p 3 x . . . p n w n=b*p_1^k*p_2^m*p_3^x...p_n^w n=b∗p1k∗p2m∗p3x...pnw
我们发现此时 n n n被分解成了有限个质数,和我们的假设矛盾,所以原命题成立。
3、分解质因数
(1)问题
(2)思路
其实很简单,这道题就是考察我们怎么将一个合数写成算术基本定理的样子。
即
a
=
p
1
k
∗
p
2
m
∗
p
3
x
.
.
.
p
n
w
a=p_1^k*p_2^m*p_3^x...p_n^w
a=p1k∗p2m∗p3x...pnw。
因此,我们需要先找到这个数的质因数 m m m,这样我们就确定了底数,然后我们再不断地让 a / m a/m a/m和 a a a% m m m。这样我们就能确定指数。
因为只要指数大于0,那么一定可以被整除,当指数降低为0以下的时候,此时就无法被整除了。
(3)代码
#include<iostream>
using namespace std;
int n;
void Prime_factor(int x)
{
for(int i=2;i<=x;i++)
{
if(x%i==0)
{
int num=0;
while(x%i==0)
{
num++;
x/=i;
}
cout<<i<<" "<<num<<endl;
}
}
puts("");
}
int main()
{
cin>>n;
while(n--)
{
int a=0;
scanf("%d",&a);
Prime_factor(a);
}
}
(4)分析
n不断变化会漏解吗?
我们上述的代码中,直接把 i i i当成了质因数,首先这个 i i i能够整除 n n n,所以它一定是因数。但是我们怎么确定这个数不是合数而是质数呢?
我们首先发现的第一件事其实应该是n是在不断变化的。每次都是去掉一个: p n k p_n^k pnk
所以我们先来思考一下,这样会不会漏解?
我们的
n
=
p
1
r
∗
p
2
t
∗
p
3
y
.
.
.
∗
p
n
u
n=p_1^r*p_2^t*p_3^y...*p_n^u
n=p1r∗p2t∗p3y...∗pnu,同时我们让
p
1
,
p
2
,
p
3
.
.
.
p
n
p_1,p_2,p_3...p_n
p1,p2,p3...pn依次增大。
并且 p 1 , p 2 , p 3 . . . p n p_1,p_2,p_3...p_n p1,p2,p3...pn都是大于等于2的整数,所以 n n n一定是大于任何一个组成 n n n的质数 p n p_n pn
因为, i i i是从小到大依次枚举的,于是先枚举到的一定是 p 1 p_1 p1。当我们将 p 1 r p_1^r p1r去掉后,此时的 n n n记作 n 1 n_1 n1。 n 1 = p 2 t ∗ p 3 y . . . ∗ p n u n_1=p_2^t*p_3^y...*p_n^u n1=p2t∗p3y...∗pnu
此时 i < p 2 < p 3 < . . . < p n < n 1 i<p_2<p_3<...<p_n<n_1 i<p2<p3<...<pn<n1,那么随着 i i i从当前到 n 1 n_1 n1之间逐渐增大的时候,这些质数还是可以被枚举到的。
不会出现漏解的情况。
i一定是质因数?
假设我们存在一个合数
b
b
b符合条件。由于
b
b
b是合数,所以根据算数基本定理,
不妨再假设
b
=
p
1
q
∗
p
2
w
b=p_1^q*p_2^w
b=p1q∗p2w
同时 p 1 < b p_1<b p1<b,此时我们会先枚举到 p 1 p_1 p1。因为 b b b能够整除 n n n,所以 p 1 p_1 p1一定是能够整除 n n n的。因此, i = p 1 i=p_1 i=p1会优先出现。当出现这个质数的时候,我们会更新 n n n,将 n n n中的 p 1 p_1 p1删掉。
当 p 1 p_1 p1删掉后,此时的 n = p 2 w ∗ p 3 y . . . ∗ p n f n=p_2^w*p_3^y...*p_n^f n=p2w∗p3y...∗pnf
很明显,此时更新后的 n n n没有办法被 b b b整除了(因为 n n n中已经没有 p 1 p_1 p1这一项了),此时的合数 b b b不满足条件,和假设矛盾。所以不存在这样的合数。
因此我们枚举到的 i i i都是质因数。
(5)优化
一个数最多有一个 > n > \sqrt n >n的质因数
证明很简单,还是和上次的证明方式相同。如果有两个,那么这个数已经是 > n >n >n的了,不可能是 = n =n =n的。
所以,我们可以利用 n \sqrt n n来枚举,最后特判一下是否存在这样一个 > n > \sqrt n >n的数。
#include<iostream>
using namespace std;
int n;
void Prime_factor(int x)
{
for(int i=2;i<=x/i;i++)
{
if(x%i==0)
{
int num=0;
while(x%i==0)
{
num++;
x/=i;
}
cout<<i<<" "<<num<<endl;
}
}
if(x>1)cout<<x<<" "<<1<<endl;
puts("");
}
int main()
{
cin>>n;
while(n--)
{
int a=0;
scanf("%d",&a);
Prime_factor(a);
}
}