质数（素数）、约数

最新推荐文章于 2023-02-04 14:00:09 发布

_Shmily

最新推荐文章于 2023-02-04 14:00:09 发布

阅读量1.3k

点赞数 2

分类专栏：数论基本算法

本文链接：https://blog.csdn.net/zhaoxinxin1234/article/details/89528636

版权

本文介绍了质数（素数）的概念，包括它们的定义、数量特征和获取素数表的方法，如Eratosthenes筛法和欧拉线性筛。接着讨论了质因数分解的策略，以及如何快速找到正整数的约数。此外，还涉及了最大公约数的计算，如欧几里得算法。

摘要由CSDN通过智能技术生成

一、质数：
若一个正整数无法被除了1和它自身之外的任何自然数整除，则称该数为质数（或素数），
否则称该正整数为合数。
1不是素数，2是最小的素数。
在整个自然数集合中，质数的数量不多，分布比较稀疏，
对于一个最够大的整数 N，不超过 N 的质数大约有 $\frac{N}{lnN}$ 个，即每lnN个数大约有1个质数。

二、素数表的获取——Eratosthenes筛法：
小于x*x的x的倍数在扫描更小的数时就已经被标记过了。
对于某一个合数，其倍数一定在某个比其小的质数（其质因数）筛掉了。
时间复杂度为O（N loglogN）接近线性。

#include<iostream>
#include<cstring>
#include<cstdio>
using namespace std;
const int maxn=1000005;
bool v[maxn];
void init(int n)
{
   
    memset(v,0,sizeof(v));
    v[1]=true;
    for(int i=2;i<=n;i++)
    {
   
        if(v[i]) continue;
        for(int j=i;j<=n/i;j++)
            v[i*j]=true;
    }
    return ;
}

int main(void)
{
   
    int n;
    scanf("%d",&n);
    init(n);
    for(int i=1;i<=n;i++)
    {
   
        if(!v[i])
            printf("%d ",i);
    }
    return 0;
}

三、素数表的获取——欧拉线性筛：
1.

#include<iostream>
#include<cstring>
#include<cstdio>
using namespace std;
const int maxn=1000005;
int v[maxn],prime[maxn];
int cnt=0;
void init(int n)
{
   
    memset(v,0,sizeof(v));
    v[1]=1;
    for(int i=2;i<=n;i++)
    {
   
        if(v[i]==0)
        {
   
            v[i]=i;
            prime[cnt++]=i;
        }
        for(int j=0;j<cnt;j++)
        {
   
            if(prime[j]>v[i]||prime[j]>n/i) break;
            v[i*prime[j]]=prime[j];
        }
    }
    return ;
}

int main(void)
{
   
    int n;
    scanf("%d",&n);
    init(n);
    for(int i=0;i<cnt;i++)
        printf("%d ",prime[i]);

    return 0;
}

#include<iostream>
#include<cstring>
#include<cstdio>
using namespace std;
const int maxn=1000005;
int prime[maxn];
bool ha[maxn];
int cnt=0;
void Prime(int n)
{
   
    memset(ha,0,sizeof(ha));
    ha[1]=true;
    for(int i=2;i<=n;i++)
    {
   
        if(!ha[i])
        prime[cnt++]=i;
        for(int j=0;j<cnt&&prime[j]<=n/i;j++)