朝题夕解之数论

简介: 朝题夕解之数论

分解质因数


算术基本定理: 对于任何一个大于1的正整数都能唯一分解为有限个质数的乘积,可写作:

N = P1c1 P2c2 …Pmcm

其中ci都是正整数,Pi都是质数,且满足P1 < P2 < … < Pm


例题描述微信图片_20221018130107.jpg

参考代码(C++版本)

#include <iostream>
#include <algorithm>
using namespace std;
void divide(int n)
{
    //试除法枚举所有的数
    for(int i = 2;i <= n / i;i++)
        if(n % i == 0)
        {
            int s = 0;
            while(n % i == 0)
            {
                n /= i;
                s++;
            }
            printf("%d %d\n",i,s); //要清楚i才是底数
        }
    if(n > 1) printf("%d %d\n",n,1);
    puts("");
}
int main()
{
    //输入
    int n;
    scanf("%d",&n);
    while(n--)
    {
        //调函数
        int x;
        scanf("%d",&x);
        divide(x);
    }
    return 0;
}

算法模板


一、算法实现流程图:微信图片_20221018130203.jpg

二、算法的代码实现:

试除法分解质因数
void divide(int x)
{
    for (int i = 2; i <= x / i; i ++ )
        if (x % i == 0)
        {
            int s = 0;
            while (x % i == 0) x /= i, s ++ ;
            cout << i << ' ' << s << endl;
        }
    if (x > 1) cout << x << ' ' << 1 << endl;
    cout << endl;
}

疑难点剖析

一、明白题目需求

微信图片_20221018130256.png题目要的是将传入的数据分解为底数是/质数,指数是整数,分解结果的乘积等于原数据,输出的时候,要从小到大排列。


例如:

6就是21 x 31,因此就输出 2 1和3 1

48就是24 x 31,因此就输出2 4 和 3 1


二、指数的获取

8 = 23 = 2 x 2 x 2;


那么在代码层面,我们可以采用逆向思维,倒着对8进行运算,统计它进行的除法次数,也就是相乘的次数了。同时,也就是指数了。


质数筛选



例题描述

微信图片_20221018130432.jpg

参考代码(C++版本)

// 埃氏筛选法
#include <iostream>
#include <algorithm>
using namespace std;
const int N = 1e6 + 10;
int primes[N],cnt;
bool st[N];
void get_primes(int n)
{
    for(int i = 2;i <= n;i++)
    {
        if(!st[i])
        {
            primes[cnt++] = n;
            //利用现有的质数,将质数的倍数去掉
            for(int j = i+i; j <= n;j += i) st[j] = true;
        }
    }
}
int main()
{
    //输入
    int n;
    cin >>n;
    //调用函数
    get_primes(n);
    //输出
    cout << cnt <<endl;
    return 0;
}
//线性筛选法  算法原理:n只会被最小质因子筛掉
#include <iostream>
#include <algorithm>
using namespace std;
const int N = 1e6 + 10;
int primes[N],cnt;
bool st[N];
void get_primes(int n)
{
    for(int i = 2;i <= n;i++)
    {
        //如果是质数,就把加到数组中去
        if(!st[i]) primes[cnt ++]  = i;
        //从小到大枚举所有的质数
        for(int j = 0; primes[j] <= n/i;j++)
        {
            //把prims[j] * i筛掉
            st[primes[j] * i] = true;
            if(i % primes[j] == 0) break; //primes[j]一定是i的最小质因子
        }
        /*
        1、i % pj == 0 
            pj一定是i的最小质因子,pj一定是pj*i的最小质因子
        2、i % pj != 0
            pj一定小于i的所有质因子,pj也一定是pj * i的最小质因子
        */
    }
}
int main()
{
    //输入
    int n;
    cin >>n;
    //调用函数
    get_primes(n);
    //输出
    cout << cnt <<endl;
    return 0;
}

算法模板


一、埃氏筛法——用当前已有的质数去消去它们的倍数


首先,将2到n范围的所有整数获取到。其中,最小的数字2是质数,将2的所有倍数划去。表中剩余的数字中,最小的是3,它不能被更小的整除,所有它是质数,再将所有3的倍数划去。以此类推如果表中剩余的最小数字是m时,m就是质数。然后将表中所有m的倍数都划去。像这种反复操作,就能够依次枚举n以内的质数。

举个栗子:

微信图片_20221018130534.png

1.1、 算法实现流程如下:

微信图片_20221018130602.jpg

1.2、算法代码实现:

埃氏筛法求素数
int primes[N], cnt;   // primes[]存储所有素数
bool st[N];     // st[x]存储x是否被筛掉
void get_primes(int n)
{
    for (int i = 2; i <= n; i ++ )
    {
        if (st[i]) continue;
        primes[cnt ++ ] = i;
        for (int j = i; j <= n; j += i)
            st[j] = true;
    }
}

1.3、算法的时间复杂度 = O(NloglogN)


二、线性筛法

线性筛选的核心是—— 传入的整数n只会被最小质因子筛掉

2.1、算法实现流程:

image.jpeg

2.2、算法代码实现:

线性筛法求素数
int primes[N], cnt;   // primes[]存储所有素数
bool st[N];     // st[x]存储x是否被筛掉
void get_primes(int n)
{
    for (int i = 2; i <= n; i ++ )
    {
        if (!st[i]) primes[cnt ++ ] = i;
        for (int j = 0; primes[j] <= n / i; j ++ )
        {
            st[primes[j] * i] = true;
            if (i % primes[j] == 0) break;
        }
    }
}

2.3、算法时间复杂度 = O(N)


疑难点剖析


清楚primes数组和st数组是分别用来维护素数集合和被筛除数据的集合

相关文章
|
7月前
数论——高斯消元
数论——高斯消元
37 0
|
机器学习/深度学习 人工智能
数学知识-质数
数学知识-质数
数学知识-约数
数学知识-约数
|
算法
数学知识:中国剩余定理
复习acwing算法基础课的内容,本篇为讲解数学知识:中国剩余定理,关于时间复杂度:目前博主不太会计算,先鸽了,日后一定补上。
125 0
数学知识:中国剩余定理
|
存储 算法 图计算
数学知识:容斥原理
复习acwing算法基础课的内容,本篇为讲解数学知识:容斥原理,关于时间复杂度:目前博主不太会计算,先鸽了,日后一定补上。
82 0
数学知识:容斥原理
|
算法
数学知识:快速幂
复习acwing算法基础课的内容,本篇为讲解数学知识:快速幂,关于时间复杂度:目前博主不太会计算,先鸽了,日后一定补上。
76 0
数学知识:快速幂
|
算法
数学知识:质数(一)
复习acwing算法基础课的内容,本篇为讲解数学知识:质数,关于时间复杂度:目前博主不太会计算,先鸽了,日后一定补上。
114 0
数学知识:质数(一)
|
存储 算法
数学知识:高斯消元(一)
复习acwing算法基础课的内容,本篇为讲解数学知识:高斯消元,关于时间复杂度:目前博主不太会计算,先鸽了,日后一定补上。
118 0
数学知识:高斯消元(一)