素数相关(结合回文数,合数)线性筛素数(欧拉筛法)Euler【算法模板笔记】

news2024/11/29 22:42:13
一、朴素筛法(埃拉托斯特尼筛法)
  • Eratosthenes 筛法(埃拉托斯特尼筛法,简称埃氏筛法)

  • 时间复杂度是O(nloglogn)

不常用,被欧拉筛代替,略

二、 线性筛素数(欧拉筛法)
  • 简介

  • 线性筛法 也称为 Euler 筛法(欧拉筛法)

  • 对比

  • 普通的筛法就是1到n的倍数来筛,线性筛就是 用1到n中的素数的倍数来筛

  • 时间复杂度 O(n)

  • buff

  • 筛法求素数的同时也得到了每个数的最小质因子。

  • 原理

  • 中心思想:每个数只能被自己的最小质因子筛掉一次

  • 算法原理解释

  • 原理:对于任意合数,必定可以有最小质因子乘以最大因子的分解方式。因此,对于每个合数,只要用最大因子筛一遍,枚举时只要枚举最小质因子即可。由于每个合数都只被标记一次,达到了线性

  • // 任意合数,必定可以有最小质因子乘以最大因子的分解方式。

  • // 所以我们只要保证每个合数都由最小质因子筛掉就行了

  • //用两层循环枚举最小质因子和最大因子。i是最大因子,prime[j]是最小质因子

  • break整除中断原因

  • 解释1

  • // 这里为什么要break

  • // 因为如果不break

  • // i * p[j + 1] 的最小质因子就不是p[j + 1] 了

  • // 而是 i 的最小质因子 p[j]。

  • 解释2

  • //这个break发生时,这个primes[j]的值,是i的最小质因子

  • //保证合数只被最小质因子划掉一次

  • //如果i是质数,则最多枚举到自身中断

  • //如果i是合数,则最多枚举到自身的最小质数中断

  • 文章

  • 【算法/数论】欧拉筛法详解:过程详述、正确性证明、复杂度证明_CSDN博客_欧拉筛

  • 例题

  • P3383 【模板】线性筛素数 - 洛谷 | 计算机科学教育新生态 (luogu.com.cn)

  • Sherlock and his girlfriend - 洛谷 | 计算机科学教育新生态 (luogu.com.cn)

三、与回文数结合
  • 回文数

  • 判断回文数

  • 构造回文数

  • 回文质数

  • 偶数肯定不是质数。这样至少排除一半多的数据量

  • 知识点: “偶数长度的回文数”中只有11是素数,其他的都可以被11整除。

  • 偶数位数回文数(除11)必定不是质数

  • 例题

  • 回文素数 - 回文素数 - 力扣(LeetCode)

  • P1217 [USACO1.5]回文质数 Prime Palindromes - 洛谷 | 计算机科学教育新生态 (luogu.com.cn)

四、与合数结合
  • 方法

  • 用欧拉筛反向筛出合数

  • 合数相关理论

  • 例题

  • P1835 素数密度 - 洛谷 | 计算机科学教育新生态 (luogu.com.cn)

模板
/*prime_nums_Euler*/
#include<bits/stdc++.h>
#define MAXN 100000005
using namespace std;

//P3383 【模板】线性筛素数

bool isprime[MAXN]; // isprime[i]表示i是不是素数
int prime[MAXN]; // 现在已经筛出的素数列表
int nums = 0; // 已经筛出的素数个数

void euler_old(int n){
    memset(isprime, true, sizeof(isprime)); // 先全部标记为素数
    isprime[1] = false; // 1不是素数
    for(int i = 2; i <= n; ++i){ // i从2循环到n(外层循环)
        if(isprime[i]) prime[++nums] = i;// 如果i没有被前面的数筛掉,则i是素数
        for(int j = 1; j <= nums && prime[j] <= n/i; ++j){//枚举已经记录的质数(内层循环)
        //i * prime[j] <= n:由于题目只求小于n的数是否是质数,所以大于n的就不管了【越界中断】
            isprime[i * prime[j]] = false;// 倍数标记为合数,也就是i用prime[j]把i * prime[j]筛掉了
            if(i % prime[j] == 0) break;//【整除中断】
        }
    }
}

//以下是acwing的模板

int primes[MAXN], cnt;  // primes[]存储所有素数
bool st[MAXN];  // st[x]存储x是否被筛掉:true表示要被筛掉(没有处理1的特判)

void get_primes(int n){
    for(int i = 2; i <= n; i++){
        if(!st[i]) primes[cnt++] = i;  //如果没有被筛掉,直接录入,并且cnt++
        for(int j = 0; primes[j] <= n/i; j++){  //从小到大枚举所有质数
            st[primes[j] * i] = 1;  //筛掉质数的倍数的数
            if (i % primes[j] == 0) break;  //【精髓:整除中断】
                //这个break发生时,这个primes[j]的值,是i的最小质因子
        }
    }
}

int main(){
    int n; // 上限,即筛出<=n的素数
    scanf("%d", &n);
    get_primes(n);
    int q,k;//q次询问第 k 小的素数
    scanf("%d", &q);
    while(q--){
        scanf("%d", &k);
        printf("%d\n", primes[k-1]);
    }
    // for(int i = 2; i<=n; i++) if(isprime[i]) printf("%d ", i);
    // for(int i = 2; i<=n; i++) if(!st[i]) printf("%d ", i);

    return 0;
}

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

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

相关文章

C++007-C++循环结构

文章目录C007-C循环结构for循环for循环举例for循环格式题目描述 输出十次手机号题目描述 打印区间内的整数题目描述 打印字符之间的所有字符题目描述 打印区间内符合条件的整数数数量作业在线练习&#xff1a;总结C007-C循环结构 在线练习&#xff1a; http://noi.openjudge.cn…

HappyAI 算法任务调度平台 - 开篇

HappyAI 算法任务调度平台 - 开篇 HappyAI 算法任务调度平台 - 接入 HappyAI 算法任务调度平台 - 开发 1. 支持不同算法侧接入即算法集群&#xff08;如&#xff1a;paddle算法平台&#xff1b;mmdetection算法平台&#xff09; 2. 支持不同相机取流&#xff08;如&#xff…

启动您的 MySQL 查询:ChatGPT 如何帮助您检索 MySQL 数据

本文涉及到的相关软件&#xff0c;dbForge Studio for MySQL和ChatGPT,在您尝操作前&#xff0c;请务必先确保两个程序都正常运作哦&#xff0c;废话不多说&#xff0c;开始上干货&#xff01; Devart 提供包括Oracle、SQL Server、MySQL、PostgreSQL、InterBase以及Firebird在…

msvcr120.dll下载官网,如何处理msvcr120.dll丢失

msvcr120.dll下载官网&#xff1f;目前msvcr120.dll是没有官网的&#xff01;MSVCR120.dll是Microsoft Visual C Redistributable安装程序的一部分&#xff0c;通常在运行依赖于Visual C库的应用程序时需要该文件。所以不要相信网上那些写着所谓的官网&#xff0c;下面给大家分…

Sentinel源码阅读

基础介绍 Sentinel 的使用可以分为两个部分: 核心库&#xff08;Java 客户端&#xff09;&#xff1a;不依赖任何框架/库&#xff0c;能够运行于 Java 8 及以上的版本的运行时环境&#xff0c;同时对 Dubbo / Spring Cloud 等框架也有较好的支持&#xff08;见 主流框架适配&…

电子纸墨水屏的现实应用场景

电子纸挺好个东西&#xff0c;大家都把注意力集中在商超场景 其实还有更多有趣的场景方案可用&#xff0c;价值也不小&#xff0c;比如&#xff1a; 一、仓库场景 通过亮灯拣选&#xff0c;提高仓库作业效率 二、仓库循环使用标签 做NFC类发卡式应用&#xff0c;替代传统纸…

Python中ArcPy将多张栅格图像分割为多个网格小块的方法

本文介绍基于Python中的ArcPy模块&#xff0c;依据渔网的矢量数据文件或通过手动划分小方格的方法&#xff0c;批量将大量栅格图像分割为多个矩形小栅格的方法。 首先明确一下我们的需求。现有同一区域内的多张栅格遥感影像&#xff0c;我们希望对于每一景栅格遥感影像而言&…

QML- 信号和事件系统

QML- 信号和事件系统一、概述二、信号事件处理系统使用1. 用信号处理器接收信号2. 信号处理程序3. 使用Connections 类型4. 附加信号处理程序三、向自定义QML类型添加信号四、将信号连接到方法和信号1. 信号与信号的连接一、概述 应用程序和用户界面组件需要相互通信。例如&am…

数据迁移——技术选型

日常我们在开发中&#xff0c;随着业务需求的变更&#xff0c;重构系统是很常见的事情。重构系统常见的一个场景是变更底层数据模型与存储结构。这种情况下就要对数据进行迁移&#xff0c;从而使业务能正常支行。 背景如下&#xff1a;老系统中使用了mongo数据库&#xff0c;由…

webpack 的热更新是如何做到的?原理是什么?

Hot Module Replacement&#xff0c;简称 HMR&#xff0c;在不需要刷新整个页面的同时更新模块&#xff0c;能够提升开发的效率和体验。热更新时只会局部刷新页面上发生了变化的模块&#xff0c;同时可以保留当前页面的状态&#xff0c;比如复选框的选中状态等。 在 webpack 中…

Maxwell工作流程详解

要介绍maxwell的工作原理&#xff0c;首先需要讲一下mysql主从复制的原理 mysql主从复制原理&#xff1a; 如上图&#xff0c;左边是master主节点&#xff0c;右边是slave从节点 工作流程&#xff1a; 1.往主节点mysql的数据库中写入数据&#xff0c;产生数据变化&#xff0c…

让ChatGPT介绍一下ChatGPT

申请新必应内测通过了&#xff0c;我在New Bing中使用下ChatGPT&#xff0c;让ChatGPT介绍一下ChatGPT 问题1&#xff1a;帮我生成一篇介绍chatGPT的文章&#xff0c;不少于2000字 回答&#xff1a; chatGPT是什么&#xff1f;它有什么特点和用途&#xff1f; chatGPT是一种…

IO流概述

&#x1f3e1;个人主页 &#xff1a; 守夜人st &#x1f680;系列专栏&#xff1a;Java …持续更新中敬请关注… &#x1f649;博主简介&#xff1a;软件工程专业&#xff0c;在校学生&#xff0c;写博客是为了总结回顾一些所学知识点 目录IO流概述IO 流的分类总结流的四大类字…

TCP传输文件

传输文件和传输信息的区别&#xff1a; 传输信息&#xff0c;只是一条数据&#xff0c;传输文件是多条数据传输信息传输过去一般都会显示&#xff0c;传输文件一般不会显示&#xff0c;一般只是存放在文件中传输文件需要传输&#xff0c;文件大小和文件名称&#xff08;不然不知…

C语言(字符串函数)

这章的内容记得引用<string.h>头文件 目录 1.strlen&#xff08;&#xff09; 2.strcat() 3.strncat() 4.strcmp() 5.strncmp() 6.strcpy() 7.strncpy() 8.sprintf() 8.strchr() 9.strpbrk() 10.strrchr() 11.strstr() 1.strlen&#xff08;&#xff09; 用于统计字符串的…

CentOS Stream 9尝鲜安装教程

作者&#xff1a;IT圈黎俊杰 一、下载CentOS Stream 9安装介质 在CentOS官网可以下载到CentOS Stream 9的安装介质&#xff0c;正面列出ISO介质的下载链接地址&#xff1a; https://download.cf.centos.org/9-stream/BaseOS/x86_64/iso/CentOS-Stream-9-20221019.0-x86_64-dv…

157、【动态规划】leetcode ——377. 组合总和 Ⅳ(C++版本)

题目描述 原题链接&#xff1a;377. 组合总和 Ⅳ 解题思路 一维滚动数组思路与 518. 零钱兑换 II&#xff08;二维数组一维滚动数组&#xff09; 大致相同&#xff0c;区别之处在于&#xff0c;零钱兑换解决的是组合问题&#xff08;元素在不同位置都代表同一个含义&#xff…

Java基础之IO流详解及使用方式(建议收藏)

目录前言一.初识IO流二.流的分类三.基本字节流3.1 FileOutputStream3.2 FileInputStream四.文件拷贝五.字符集六.乱码原因七.基本字符流7.1 FileReader7.2 FileWriter经典案例Demo八.高级流8.1 缓冲流8.1.1 字节缓冲流8.1.2 字符缓冲流8.2 转换流8.3 序列化流8.4 打印流8.4.1 字…

【算法】前缀和

作者&#xff1a;指针不指南吗 专栏&#xff1a;算法篇 &#x1f43e;要学会在纸上打草稿&#xff0c;这个很重要&#x1f43e; 文章目录1.什么是前缀和&#xff1f;2.怎么求前缀和&#xff1f;3.前缀和有什么用&#xff1f;4.进阶二维:矩阵和前缀和 主打一个记公式 1.什么是前…

四数之和-力扣18-java排序+双指针

一、题目描述给你一个由 n 个整数组成的数组 nums &#xff0c;和一个目标值 target 。请你找出并返回满足下述全部条件且不重复的四元组 [nums[a], nums[b], nums[c], nums[d]] &#xff08;若两个四元组元素一一对应&#xff0c;则认为两个四元组重复&#xff09;&#xff1a…