第二十三章 数论——质数(1)(超级详细的推导)

news2024/12/23 2:16:03

第二十三章 数论——质数

  • 一、什么是质数
  • 二、质数的判断
    • 1、试除法(朴素版)
    • 2、试除法(优化版)
  • 三、分解质因数

一、什么是质数

质数也叫素数。如果一个大于等于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) (p1)的数字整除。

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=ab,那么必定满足不等式: a ≤ p ≤ b a\leq \sqrt p \leq b ap 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 ab<p

如果 a > p a > \sqrt p a>p 同时 b > p b > \sqrt p b>p ,那么 a ∗ b > p a*b > p ab>p

因此,为了让 a ∗ b = p a*b=p ab=p,只能是: a ≤ p ≤ b a\leq \sqrt p \leq b ap b


由于,我们只是想判断一个数是不是质数,所以我们只需要找到一个 2 2 2 ~ ( p − 1 ) (p-1) (p1)的数能够整除 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=ab

由于 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=p1kp2mp3x...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=bp1kp2mp3x...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=p1kp2mp3x...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=p1rp2tp3y...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=p2tp3y...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=p1qp2w

同时 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=p2wp3y...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);
    }
    
}

本文来自互联网用户投稿,该文观点仅代表作者本人,不代表本站立场。本站仅提供信息存储空间服务,不拥有所有权,不承担相关法律责任。如若转载,请注明出处:http://www.coloradmin.cn/o/113947.html

如若内容造成侵权/违法违规/事实不符,请联系多彩编程网进行投诉反馈,一经查实,立即删除!

相关文章

RepNAS: 基于NAS的结构重参数化技术

1. 介绍 在过去几年里&#xff0c;NAS技术取得了长足进展。然而&#xff0c;由于搜索约束与实际推理之间的差异导致高效网络搜索仍极具挑战性。为搜索一个具有高性能、低推理延迟的模型&#xff0c;已有方案往往在算法中添加计算复杂度约束。然而&#xff0c;推理速度会受多种…

【强化学习笔记】马尔可夫过程、马尔可夫奖励过程

文章目录1.马尔可夫过程1.1.随机过程1.2.马尔可夫性质1.3.马尔可夫过程2. 马尔可夫奖励过程2.1.回报2.2.价值函数3.马尔可夫决策过程1.马尔可夫过程 马尔可夫过程&#xff08;Markov process&#xff09; 指具有 马尔可夫性质 的 随机过程 &#xff0c;也被称为马尔可夫链&…

C++GUI之wxWidgets(4)-编写应用涉及的类和方法(2)-wxDialog,wxCloseEvent

目录wxDialog包含类继承具体描述模态和无模态支持样式此类发出的事件wxWindow:&#xff1a;Close()wxCloseEvent具体描述使用此类的事件wxDialog 包含 #include <wx/dialog.h>类继承 描述主 具体描述 对话框是一个带有标题栏的窗口&#xff0c;有时还有一个系统菜单…

python-多线程、网络编程、正则表达式

目录 闭包 多线程 主线程 线程阻塞 同步锁 网络编程 正则表达式 re.match函数 re.search方法 re.match与re.search的区别 re.findall()方法 正则表达式的特殊规则 闭包 account0 def atm(num,flag):global accountif flag:accountnumaccountprint(account)else:acco…

免费开源的高精度OCR文本提取,支持 100 多种语言、自动文本定位和脚本检测,几行代码即可实现离线使用(附源码)

免费开源的高精度OCR文本提取,支持 100 多种语言、自动文本定位和脚本检测,几行代码即可实现离线使用(附源码)。 要从图像、照片中提取文本吗?是否刚刚拍了讲义的照片并想将其转换为文本?那么您将需要一个可以通过 OCR(光学字符识别)识别文本的应用程序。 图片文字识…

html圣诞树代码

一、前言 想做一个圣诞树&#xff0c;通过html实现了下 二、效果展示 三、代码 <!DOCTYPE html> <html> <head> <meta http-equiv"Content-Type" content"text/html; charsetutf-8" /> <meta name"viewport" cont…

245. 你能回答这些问题吗——线段树

给定长度为 N 的数列 A&#xff0c;以及 M 条指令&#xff0c;每条指令可能是以下两种之一&#xff1a; 1 x y&#xff0c;查询区间 [x,y] 中的最大连续子段和&#xff0c; 2 x y&#xff0c;把 A[x] 改成 y。 对于每个查询指令&#xff0c;输出一个整数表示答案。 输入格式…

RabbitMQ 第二天 高级 9 RabbitMQ 集群搭建 9.3 集群管理 9.5 负载均衡-HAProxy

RabbitMQ 【黑马程序员RabbitMQ全套教程&#xff0c;rabbitmq消息中间件到实战】 文章目录RabbitMQ第二天 高级9 RabbitMQ 集群搭建9.3 集群管理9.5 负载均衡-HAProxy9.5.1 安装HAProxy9.5.2 配置HAProxy第二天 高级 9 RabbitMQ 集群搭建 9.3 集群管理 rabbitmqctl join_cl…

MariaDB上市:MySQL之父奋斗13年终敲钟 要写代码写到100岁

雷递网 雷建平 12月24日云数据库公司MariaDB日前与特殊目的公司Angel Pond Holdings完成合并&#xff0c;并在纽交所上市&#xff0c;新公司更名为MariaDB。MariaDB是2022年初与Angel Pond Holdings达成合并协议&#xff0c;对新公司的作价为6.72亿美元。MariaDB是MySQL之父Mic…

【技术应用】java基于UNIX域套接字(unix domain socket)连接redis

【技术应用】java基于UNIX域套接字unix domain socket连接redis一、前言二、实现思路三、代码实现1、java socket基于redis.sock连接redis2、Lettuce框架基于redis.sock连接redis一、前言 在公司工作中经常涉及到一些中小型项目&#xff0c;这些项目都会涉及使用redis数据库&a…

Redis5.0+——持久化——RDBAOF

Redis持久化-RDB 1.实现目标&#xff1a; 在redis持久化时&#xff0c;持久化dump.rdb文件放入到redis解压目录下的data目录下的6379目录下 2.前期准备 1.在redis-5.0.3解压目录下新建data数据目录 2.编辑前面配置的/etc/redis.conf配置文件 修改持久化文件位置 (1) 进入安…

MySQL热备之PXB备份与恢复

&#x1f4e2;&#x1f4e2;&#x1f4e2;&#x1f4e3;&#x1f4e3;&#x1f4e3; 哈喽&#xff01;大家好&#xff0c;我是【IT邦德】&#xff0c;江湖人称jeames007&#xff0c;10余年DBA工作经验 一位上进心十足的【大数据领域博主】&#xff01;&#x1f61c;&#x1f61…

【语音处理】使用块反射器的基于DFT的系统中用于旁瓣抑制的正交预编码(Matlab代码实现)

&#x1f468;‍&#x1f393;个人主页&#xff1a;研学社的博客 &#x1f4a5;&#x1f4a5;&#x1f49e;&#x1f49e;欢迎来到本博客❤️❤️&#x1f4a5;&#x1f4a5; &#x1f3c6;博主优势&#xff1a;&#x1f31e;&#x1f31e;&#x1f31e;博客内容尽量做到思维缜…

Java中的多线程(下)

作者&#xff1a;~小明学编程 文章专栏&#xff1a;JavaEE 格言&#xff1a;热爱编程的&#xff0c;终将被编程所厚爱。 目录 多线程案例 单例模式 饿汉模式 懒汉模式 阻塞式队列 为什么要引入阻塞队列 Java中的阻塞队列 模拟实现阻塞队列 定时器 标准库中的定时器 …

docker的虚悬镜像是什么?

虚悬镜像是什么? 答:仓库名、标签都是<none>的镜像,俗称:dangling image 我们使用Dockerfile写一个: 1:编写 from ubuntu CMD echo action is success2:构建 docker build . 注意没有 -t 产生原因: 1:构建时候因为编写错误导致 2:删除的时候 对于这样…

数据溢出的二进制原理

char 类型的数据占一个字节&#xff0c;一个字节有 8 位&#xff0c;最高位为符号位&#xff0c;1表示负数&#xff0c;0表示正数。在计算机中&#xff0c;数据用补码表示&#xff0c;正数的补码是它本身&#xff0c;负数的补码为 “符号位不变&#xff0c;其他位取反后再加1”…

Spring Cloud 系列之OpenFeign:(4)集成OpenFeign

目录 传送门 服务间调用 集成OpenFeign 说明文档 添加pom依赖 启用OpenFeign 声明OpenFeign接口 改造远程调用 定义OpenFeign接口 测试OpenFeign调用 传送门 Spring Cloud Alibaba系列之nacos&#xff1a;(1)安装 Spring Cloud Alibaba系列之nacos&#xff1a;(2)单…

离散数学数理逻辑部分【2】

文章目录命题逻辑等值演算公式的使用【重点】析取范式和合取范式【重点】范式存在定义【了解】求公式A的范式的步骤&#xff1a;【重点】极大项和极小项【重点】主合取范式和主析取范式【重点】等式演算求主析取范式【重点】真值表求主析取范式【了解】主范式的应用【重点】推理…

Python : 使用python实现教务管理系统(GUI界面+数据库)

一、设计目的 1.熟悉Python和相关软件的操作。 2.基于本学期所学Python知识&#xff0c;熟练应用掌握&#xff0c;制作符合要求的教务管理系统。 3.会对程序运行中的错误代码进行分析&#xff0c;找出合理的解决方案。 4.掌握tkinter开发流程&#xff0c;布局方法和主要组件&a…

C语言位域

如果程序的结构中包含多个开关量&#xff0c;只有 TRUE/FALSE 变量&#xff0c;如下&#xff1a; struct {unsigned int widthValidated;unsigned int heightValidated; } status; 这种结构需要 8 字节的内存空间&#xff0c;但在实际上&#xff0c;在每个变量中&#xff0c;…