P3383 【模板】线性筛素数

最新推荐文章于 2025-09-07 20:32:33 发布

weixin_30338743

最新推荐文章于 2025-09-07 20:32:33 发布

阅读量50

点赞数

CC 4.0 BY-SA版权

文章标签：数据结构与算法

原文链接：http://www.cnblogs.com/xiaoyezi-wink/p/11003756.html

本文详细介绍了两种高效的素数筛选算法——线性筛与埃氏筛，通过具体实例展示了如何在C++中实现这两种算法，并对比了它们在处理大规模数据时的效率。适合对算法优化及素数筛选感兴趣的读者。

P3383 【模板】线性筛素数

题解

说好的数组不能开很大的呢！？？！这个题开了1e7+10的数组也没事QWQ

埃氏筛也能水过

代码

1.线性筛

#include<bits/stdc++.h>

using namespace std;

const int maxn=1e7+10;
int not_prime[maxn],prime[maxn];
int m,n,cnt_pri;
int a[maxn];

void xxs(int mn)
{
    memset(not_prime,0,sizeof(not_prime));
    not_prime[1]=1;
    for(int i=2;i<=mn;i++)
    {
        if(!not_prime[i])  prime[++cnt_pri]=i;
        for(int j=1;j<=cnt_pri;j++)
        {
            if(prime[j]*i>mn) break;
            not_prime[prime[j]*i]=1;
            if(i%prime[j]==0) break;
        }
    }
}


int main()
{
    scanf("%d%d",&m,&n);
    for(int i=1;i<=n;i++)
      scanf("%d",&a[i]);
    
    xxs(m);
    
    for(int i=1;i<=n;i++)
      if(not_prime[a[i]]) printf("No\n");
      else printf("Yes\n");
    
}

2.埃氏筛

#include<bits/stdc++.h>

using namespace std;

const int maxn=1e7+10;
int is_prime[maxn];
int m,n;
int a[maxn];

int main()
{
    scanf("%d%d",&m,&n);
    for(int i=1;i<=n;i++)
      scanf("%d",&a[i]);
    
    memset(is_prime,1,sizeof(is_prime));
    is_prime[1]=0;
    for(int i=2;i<=m;i++)
    {
        if(is_prime[i])
          for(int j=i+i;j<=m;j+=i)
          is_prime[j]=0;
    }
    
    for(int i=1;i<=n;i++)
      if(is_prime[a[i]]) printf("Yes\n");
      else printf("No\n");
    
}