大一寒假集训（10）---素数筛

大一寒假集训（10）—素数筛

知识点

1.

筛法的思想是去除要求范围内所有的合数，剩下的就是素数 了，而任何合数都可以表示为素数的乘积，因此如果已知一 个数为素数，则它的倍数都为合数。

long long su[MAX],cnt;
bool isprime[MAX];
void prime()
{
    cnt=1;
    memset(isprime,1,sizeof(isprime));//初始化，将所有数认为为素数并标记为1
    isprime[0]=isprime[1]=0;//0和1不是素数
    for(long long i=2;i<=MAX;i++)
    {
        if(isprime[i])//保存素数
        {
            su[cnt++];
        }
        for (long long j=i*2;j<=MAX;j+=i)//素数的倍数都为合数
        {
            isprime[j]=0;
        }
    }
}

2.埃氏筛

合数的倍数一定会在筛素数倍数时候被筛掉，所以只筛素数就好，只把质数的素 数筛掉 就是找到一个质数，把它的倍数全部标记为合数。但是你会发现有的数 字会被标记多次，比如 12 被 2，3，都标记，这样会浪费时间。
⭐️ 最常用的方法

const int N=1e7+1;
int prime[N];
int b[N];
int cnt=0,max1=1e7;
int init()
{
    memset(b,1,sizeof(b));
    b[0]=b[1]=0;
    for (int i=2;i<=max1;i++)
    {
        if (b[i])
        {
            prime[++cnt]=i;
            for (int j=2;j*i<=max1;j++)
            b[i*j]=0;
        }
    }
}

3.线性筛

素数筛可以优化，普通的线性筛法虽然大大缩短了求素数的时间， 但是实际上还是做了许多重复运算，比如2*3=6，在素数2的时候筛选了一遍， 在素数为3时又筛选了一遍。如果只筛选小于等于素数i的素数与i的乘积， 既不会造成重复筛选，又不会遗漏。时间复杂度几乎是线性的

const int N=1e7+1;
int prime[N];
int b[N];
int cnt=0,max1=1e7;
int init()
{
    memset(b,1,sizeof(b));
    b[0]=b[1]=0;
    for (int i=2;i<=max1;i++)
    {
        if (b[i])
        prime[++cnt]=i;
        for(int j=1;j<=cnt&&prime[j]*i<=max1;j++)
        {
            b[prime[j]*i]=0;
            if (i%prime[j]==0)
                break;
        }
    }
    return 0;
}

1.知否知否，应是绿肥红瘦 nefu 1704

#include <bits/stdc++.h>
using namespace std;
typedef long long LL;
const int N=1e7+1;
int prime[N];
int b[N];
int cnt=0;
int max1=1e7;
int isprime()
{
    memset(b,1,sizeof(b));
    b[0]=b[1]=0;
    for (int i=2;i<=max1;i++)
    {
        if(b[i])
            prime[++cnt]=i;
        for (int j=1;j<=cnt&&prime[j]*i<=max1;j++)
        {
            b[prime[j]*i]=0;
            if (i%prime[j]==0)
                break;
        }
    }
}
int su(LL n)
{
    int flag=0;
    for (int i=1;prime[i]<=sqrt(n*1.0);i++)
        if (n%prime[i]==0)
       {
            flag=1;
            break;
       }
    if (n==1)
        flag=1;
    return flag;
}
int main()
{
    int n;
    LL x,y,z;
    isprime();
    cin>>n;
    while (n--)
    {
        cin>>x>>y>>z;
        if (su(x+y-z)==0)
            cout<<"yes"<<endl;
        else
            cout<<"no"<<endl;
    }
    return 0;
}

2.函数版素数判定 nefu 825

本题注意用scanf和printf，且将b[N]的类型改为bool,否则会TLE

#include <bits/stdc++.h>
using namespace std;
const int N=1e7;
typedef long long LL;
int cnt=0;
bool b[N];
int n;
int prime[N];
void isprime()
{
    memset(b,1,sizeof(b));
    b[0]=b[1]=0;
    for(int i=2;i<=N;i++)
    {
        if (b[i])
            prime[++cnt]=i;
        for (int j=1;j<=cnt&&prime[j]*i<=N;j++)
        {
            b[prime[j]*i]=0;
            if (i%prime[j]==0)
                break;
        }
    }
}
int main()
{
    isprime();
    while(scanf("%d",&n)!=EOF)
    {
        if(b[n]==0)
            printf("NO\n");
        else
            printf("YES\n");
    }
    return 0;
}

3.高木同学的因子 nefu 1669

本题求因子的时候做过
此题为当时做的时候的题解

#include <bits/stdc++.h>

using namespace std;
int main()
{
    ios::sync_with_stdio(false);
    long long x,y;
    long long n;
    int i;
    cin>>x>>y;
    int ans=0;
    n=__gcd(x,y);//c++中可直接使用  ————gcd(x,y)
    for (i=1;i*i<n;i++)
    {
        if (n%i==0)
            ans+=2;
    }
    if (i*i==n) ans++;
    cout<<ans<<endl;
    return 0;
}

4.半素数 nefu 587

本题只需改一下judge函数即可

#include <bits/stdc++.h>
using namespace std;
const int N=1e6;
bool a[N];
int prime[N];
int n;
int cnt;
void isprime()
{
    memset(a,1,sizeof(a));
    a[1]=0;
    for (int i=2;i<=N;i++)
    {
        if(a[i])
            prime[++cnt]=i;
        for(int j=1;j<=cnt&&i*prime[j]<=N;j++)
        {
            a[i*prime[j]]=0;
            if(i%prime[j]==0)
                break;
        }
    }

}
bool judge(int n)
{
    for(int i=2;i*i<=n;i++)
    {
        if (n%i==0&&a[i]==1&&a[n/i]==1)
            return 1;
    }
    return 0;
}
int main()
{
    isprime();
    while(scanf("%d",&n)!=EOF)
    {
        if(judge(n)==1)
            printf("YES\n");
        else
            printf("NO\n");
    }
    return 0;
}

5.纯素数 nefu 586

#include <bits/stdc++.h>
using namespace std;
const int N=1e6;
bool b[N+1];
int prime[N+1];
int n,cnt;
void isprime()
{
    memset(b,1,sizeof(b));
    b[1]=0;
    for (int i=2;i<=N;i++)
    {
        if(b[i])
            prime[++cnt]=i;
        for (int j=1;j<=cnt&&i*prime[j]<=N;j++)
        {
            b[i*prime[j]]=0;
            if(i%prime[j]==0)
                break;
        }
    }
}
int x;
bool judge(int x)
{
    int r;
    r=0;
    int c[10];
    memset(c,0,sizeof(c));
    //int s=0;
    while(x)
    {
        c[++r]=x%10;
        x=x/10;
    }
    for (int i=1;i<=r;i++)
    {
        int s=0;
        for (int j=i;j>=1;j--)
            s=s*10+c[j];
        if(!b[s])
            return 0;
    }
    return 1;
}
int main()
{
    isprime();
    while(scanf("%d",&n)!=EOF)
    {
        if(judge(n))
            printf("YES\n");
        else
            printf("NO\n");
    }
    return 0;
}