2021-01-27

N以内孪生素数个数（想了三个小时）

#include<stdio.h>
#include<math.h>
#define N 100001
int vit[100001]= {99,99},prime[40000],twinprime[100001]= {0},count=0;
void creat_prime()
{
    int m=floor(sqrt(N)+0.5);
    for(int i=2; i<=m; i++)
    {
        if(vit[i]==0)
            for(int j=i*i; j<=N; j=j+i)
                vit[j]=1;
    }
    for(int i=2; i<=N; i++)
        if(vit[i]==0)
            prime[count++]=i;
}
void statistics_twinprime()
{
    int k=0,flag=3;
    for(int i=1; i<count; i++)
    {
        if(prime[i]-prime[i-1]==2)
        {
            k++;
            twinprime[prime[i]]=k;
            for(int j=flag+1; j<prime[i]; j++)
                twinprime[j]=k-1;
                flag=prime[i];
        }
    }
    for(int i=flag+1;i<100000;i++)
        twinprime[i]=k;
}
int main()
{
    int n;
    creat_prime();
    statistics_twinprime();
    while(scanf("%d",&n)&&n>=0)
    {
        printf("%d\n",twinprime[n]);
    }
    return 0;
}