当前位置: 首页 > news >正文

C++数论————质数筛法(单独判断一个数,判断N个数) 埃氏筛法

质数想必大家都不陌生

从小学到大

质数的概念:一个数如果除了1和本身之外没有其他的因子,那么这个数被称为质数

今天要讲两个知识点:

  1. 在C++中如何判断一个数是否为质数

  1. 在C++中如何判断1-N之间哪些数为整数

    • 在C++中如何判断一个数是否为质数

这个知识点较为简单

充分利用上面的概念

我们用一个for循环来解决

因为对于每个数来说,1和n都是他们的因子

所以循环开始为2,终止为n-1

for (int i=2;i<n;i++){
    if (n%i==0){
        return false;
    }
}
return true;

But……

时间复杂度也忒高了吧

大概就是O(n)

那如果n>10^9的话,那直接出1秒了

所以还得优化

咱们来举个例子,n=16:

n的因子有

1 2 4 8 16

可以看出,4就像一条分割线,分成了前后两部分

前后一一对应

没错 优化的点找到了

循环前面的部分,如果前面都没有的话,那么后面肯定也没有

那么循环终止条件是什么呢?

4是16的什么数?

因数?

算术平方根!

也就是C++中的sqrt

恭喜你!优化成功

代码就不写了,就是把循环条件的n改成sqrt(n),然后加个cmath头文件

成功把时间复杂度降到O()

第一个问题成功解决

    • 在C++中如何判断1-N之间哪些数为整数

会不会有一些小伙伴说了

哎哟不就是刚才的刚才的函数调用n编吗?至于单列出来一个标题吗~

NONONONONOONONNOONONONONONO

大错特错

这篇文章的主要就在这里

举个例子,10^9都用不上,咱就说n小于等于10^6,还是很常见的吧

那么时间复杂度就是O(10^6*循环次数)

也就那么亿点点的时间就好了

所以还是那句话,优化优化优化

恭迎今天主角:埃氏筛法!

网友:埃氏筛法 满足你的一切需要!

先讲一下他的思想

思想:

  1. 一个合数,一定可以拆成一个质数与另一个数的乘积

  1. 每个质数的倍数都是合数

基于以上这两条思想,就可以实现埃氏筛法

  1. 用一个数组表示每个数的状态,即是否为质数

  1. 循环每个数,如果是质数,就把他所有的倍数全部标为合数

没错就是这么简单

时间复杂度计算方面,我真的很拉

埃氏筛法大概O(n log n)左右吧

(大佬勿喷)

void eratosthemes_prim(int n){
    for (int i=2;i<=n;i++){
        if_p[i]=true;
    }
    for (int i=2;i<=n;i++){
        if (is_p[i]){
            num++;
            p[num]=i;//p数组存放质数,num代表目前质数个数
            for (int j=i<<1;j<=n;j+=i){
                is_p[j]=false;
            }
        }
    }
}

这样是不是快了很多?

别走!来看道简简单单的例题

Goldbach's Conjecture
题目描述

全屏
哥德巴赫猜想:任何大于4的偶数都可以拆成两个奇素数之和。 比如:

8=3+5

20=3+17=7+13

42=5+37=11+31=13+29=19+23

你的任务是:验证小于106的数满足哥德巴赫猜想。

输入
多组数据,每组数据一个n。

读入以0结束。

输出
对于每组数据,输出形如n = a + b,其中a,b是奇素数,=和+两边都有空格。若有多组满足条件的a,b,输出b-a最大的一组。
若无解,输出Goldbach's conjecture is wrong.。

样例输入1
8
20
42
0

样例输出1
8 = 3 + 5
20 = 3 + 17
42 = 5 + 37

提示/说明
【数据范围与提示】

对于全部数据,6≤n≤106。

思想:运用埃氏筛法筛出所有数的情况

循环遍历情况即可

# include <iostream>
# include <cstdio>
# include <cstring>
#define N 1000005
using namespace std;
int n,vis[N],flag;
int main(){
    memset(vis,0,sizeof(vis));
    vis[0]=vis[1]=1;
    for(int i=2;i<=N/2;i++){
        if(vis[i]==0){
            for(int j=2*i;j<=N;j=j+i){
                vis[j]=1;
            }
        }    
    }
    while(scanf("%d",&n)&&n){
        flag=0;
        for(int i=3;i<=N/2;i++){
            if(vis[i]==0&&vis[n-i]==0){
                printf("%d = %d + %d\n",n,i,n-i);
                flag=1;
                break;
            }
        }    
        if(flag==0){
            printf("Goldbach's conjecture is wrong.\n");
        }
    }
    return 0;
}

有些人也许会问

题目里说的是奇素数,可你的代码中没体现出奇数这个条件啊?

问这个问题的,思考一下,偶数可能是素数吗?所以不需要考虑

还有一个问题:陈景润都没证明出哥德巴赫是错的,你还判断一下有必要吗?

这个是可以删除的!

都看到这里了,不点个赞吗?

今天的质数就讲解到这里

欢迎大家关注我

持续更新——

相关文章:

  • 区区几行代码,就能全面实现 Python 自动探索性数据分析
  • C++ 浅谈之智能指针 shared_ptr 循环引用
  • python之selenium入门教程
  • Power BI散点图
  • 分享|2023年全球市场准入认证咨讯
  • synchronized底层原理
  • Java多线程之CAS中的ABA问题与JUC的常见类
  • 【Android -- 开源库】权限适配 RxPermissions 的基本使用
  • Unity-TCP-网络聊天功能(二): 登录与注册
  • 51单片机最强模块化封装(1)
  • Python学习-----起步2(变量与转义符)
  • Java程序设计实验3 | 面向对象(上)
  • 优秀的代码最终选择if else,还是switch case
  • Openharmony的编译构建--进阶篇1
  • 每天一道大厂SQL题【Day02】电商场景TopK统计
  • EMT4J详细介绍与使用,帮你找到Java版本升级带来的问题,让你在项目jdk升级不在头疼
  • 第2章:使用CSS定义样式
  • 【数据结构】动图详解单向链表
  • MySQL基础篇笔记
  • Vue3现状—必然趋势?
  • “你是做什么的?”——人们能否对工作说不?
  • 2025上海半马鸣枪,多个“首次”冲击一城双白金
  • 亲诚惠容行大道,命运与共开新篇——中共中央政治局委员、外交部长王毅谈习近平主席对越南、马来西亚、柬埔寨进行国事访问
  • 重庆发展投资董事长刘小军将任渝农商行董事长
  • 分析|开门红:一季度GDP增长5.4%超预期,市场活力信心增强
  • 深一度|2005年首败来得恰到好处,巴萨三冠王雄心依旧