第二十三章 数论——质数
- 一、什么是质数
- 二、质数的判断
- 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);
    }
    
}


















