C/C++素数筛

判断一个数是不是素数的代码：

int isprime(long long n)
{
    int flag=0;
    if(n== 0 || n==1)  flag=1;
    int i;
    for(i=2; i<=sqrt(n); i++)
    {
        if(n%i==0)
        {
            flag=1;
            break;
        }
    }
    return flag;
}

**返回0是素数，返回1不是素数。**

素数筛题目：
1.http://acm.nefu.edu.cn/problemShow.php?problem_id=586
重点是如何求出每一个去掉最高位之后的每一个数，方法是将这个数倒序存放在一个数组里面，然后用循环的方法计算出去掉最高位之后的每一个数。

#include <bits/stdc++.h>
using namespace std;
const int N=1e6;
int n,cnt,prime[N+10];
bool vis[N+10],flag[N+10];
void get_prime()
{
    memset(vis,1,sizeof(vis));
    vis[1]=0;
    for(int i=2;i<=N;i++)
    {
        if(vis[i])prime[++cnt]=i;
        for(int j=1;j<=cnt&&i*prime[j]<=N;j++)
        {
            vis[i*prime[j]]=0;
            if(i%prime[j]==0)break;
        }
    }
}
bool judge(int x)
{
    int a[8],len=0;
    memset(a,0,sizeof(a));
    while(x)
    {
        a[++len]=x%10;
        x=x/10;
    }
    for(int i=1;i<=len;i++)
    {
        int s=0;
        for(int j=i;j>=1;j--)
            s=s*10+a[j];
        if(!vis[s])return 0;
    }
    return 1;
}
int main()
{
    ios::sync_with_stdio(false);
    get_prime();
    while(cin>>n)
    {
        if(judge(n))printf("YES\n");
        else printf("NO\n");
    }
    return 0;
}

2.http://acm.nefu.edu.cn/problemShow.php?problem_id=585
这个问题的关键是要在打出素数表的同时用一个数组来记录素数的个数，假设这个数组名是vis,那么vis[y]-vis[x]代表着（x,y]的素数个数！！！

#include <bits/stdc++.h>
using namespace std;
const int N=1e6+10;
int n,cnt,prime[N],pre[N];//flag标记合数为0，pre[i]记录[1,i]中有多少个素数
bool flag[N];
void get_prime()
{
    memset(flag,1,sizeof(flag));//如果要用memset函数把数组初始化为1，必须开bool数组，不能开int数组
    flag[1]=cnt=0;//标记1为合数
    for(int i=2;i<=N;i++)
    {
        if(flag[i])prime[++cnt]=i;
        for(int j=1;j<=cnt&&prime[j]*i<=N;j++)
        {
            flag[prime[j]*i]=0;
            if(i%prime[j]==0)break;//保证每个合数只会被它的最小质因数筛去，因此每个数只会被标记一次
        }
        pre[i]=cnt;
    }
   /* for(int i=1;i<=N;i++)
    {
        if(flag[i])pre[i]=pre[i-1]+1;
        else pre[i]=pre[i-1];
    }*/

}
int main()
{
    get_prime();
    while(cin>>n)
    {
        if(n==1){printf("0\n");continue;}
        for(int i=n;i>=2;i--)//注意是i>=2，不能写i>=sqrt(n)或者i>=n/2
        {
            if(n%i==0&&flag[i])
            {printf("%d\n",pre[i]);break;}
        }
    }
    return 0;
}

3 http://acm.nefu.edu.cn/problemShow.php?problem_id=585

#include <bits/stdc++.h>
using namespace std;
const int N=1e6;
int n,cnt,prime[N+10],pre[N+10];
bool vis[N+10];
void get_prime()
{
    memset(vis,1,sizeof(vis));
    pre[1]=0;vis[1]=0;
    for(int i=2;i<=N;i++)
    {
        if(vis[i])prime[++cnt]=i;
        for(int j=1;j<=cnt&&i*prime[j]<=N;j++)
        {
            vis[i*prime[j]]=0;
            if(i%prime[j]==0)break;
        }
        pre[i]=cnt;//pre[i]记录【1，i】有多少个素数
    }
}
int main()
{
    ios::sync_with_stdio(false);
    get_prime();
    while(cin>>n)
    {
        if(n==1){printf("0\n");continue;}
        for(int i=n;i>=2;i--)
            if(n%i==0&&vis[i]){printf("%d\n",pre[i]);break;}
    }
    return 0;
}

4.http://acm.nefu.edu.cn/forumDiscuss.php?problem_id=1669

#include <bits/stdc++.h>
using namespace std;
long long x,y,z,num=0,i;
int main()
{
    scanf("%lld%lld",&x,&y);
        z=__gcd(x,y);
    for(i=1;i<sqrt(z);i++)
    {
     if(z%i==0)
     num=num+2;
 }
 if(i==sqrt(z))//注意理解为什么不可以写成z%i==0 ！！！
 num+=1;
 printf("%lld",num);
    return 0;
}

4.http://acm.nefu.edu.cn/problemShow.php?problem_id=2113

#include <bits/stdc++.h>
using namespace std;
const int N=1e8;
int n,k,t,q,cnt,prime[N+10];
bool vis[N+10];
void get_prime()
{
    memset(vis,1,sizeof(vis));
    vis[1]=0;
    for(int i=2;i<=n;i++)
    {
        if(vis[i])prime[++cnt]=i;
        for(int j=1;j<=cnt&&i*prime[j]<=N;j++)
        {
            vis[i*prime[j]]=0;
            if(i%prime[j]==0)break;
        }
    }
}
int main()
{
    ios::sync_with_stdio(false);
    cin>>n>>q;
    get_prime();
    while(q--)
    {
        cin>>k;
        printf("%d\n",prime[k]);
    }
    return 0;
}

5.http://acm.nefu.edu.cn/problemShow.php?problem_id=1262
求一个大数的因子和

#include <bits/stdc++.h>
using namespace std;
typedef long long ll;
const int N=1e6+3;
//小于1e12的最大素数 999999999989
//大于1e6的最小素数 1e6+3  注意打表素数筛的N>=1e6+3 否则会被999999999989这个数据卡掉
ll T,n,cnt,prime[N+10];
bool vis[N+10];
void get_prime()
{
    memset(vis,1,sizeof(vis));
    vis[1]=0;
    for(int i=2;i<=N;i++)
    {
        if(vis[i])prime[++cnt]=i;
        for(int j=1;j<=cnt&&i*prime[j]<=N;j++)
        {
            vis[i*prime[j]]=0;
            if(i%prime[j]==0)break;
        }
    }
}
ll f(ll x)//求x的因子和
{
    ll s=1;
    for(ll i=1;prime[i]*prime[i]<=x;i++)//找x的素因子
    {
        if(x%prime[i]==0)
        {
            ll t=1;
            while(x%prime[i]==0)
                t=t*prime[i],x=x/prime[i];
            s*=(t*prime[i]-1)/(prime[i]-1);
        }
    }
    if(x>1)s*=(1+x);//若x未被除尽到1，说明最后的x还剩下一个素数没有乘到答案里面，把它乘进去
    return s;
}
int main()
{
    ios::sync_with_stdio(false);
    get_prime();
    cin>>T;
    while(T--)
    {
        cin>>n;
        printf("%lld\n",f(n));
    }
    return 0;
}