PAT乙级 1007 素数对猜想(孪生素数对) (20分) ---（C语言 + 详细注释）

让我们定义d​n​​为：d​n​​=p​n+1​​−p​n​​，其中p​i​​是第i个素数。显然有d​1​​=1，且对于n>1有d​n​​是偶数。“素数对猜想”认为“存在无穷多对相邻且差为2的素数”。

现给定任意正整数N(<10​5​​)，请计算不超过N的满足猜想的素数对的个数。

输入格式:

输入在一行给出正整数N。

输出格式:

在一行中输出不超过N的满足猜想的素数对的个数。

输入样例:

20

输出样例:

4

//大体思路：先找到一个素数p1，判断 p1 + 2是否是素数，如果是，就是一个素数对，计数器+1

#include<stdio.h>
#include<string.h>
#define MAX 100001
int prime[MAX];
int checkprime(int n);

int main() {
    int i, k, twin, n;
    scanf("%d", &n);              //读入上界n
    checkprime(n);                //埃拉托色尼筛法（详见https://blog.csdn.net/qq_45472866/article/details/104051475）
    twin = 2;                          //第一个素数是2
    k = 0;                              //计数器
    for (i = 2; i <= n; i++)
        if (prime[i]) {                 //如果是找到了一个素数
            if (i - twin == 2)        //判断这个素数和上一个素数是否相差2
                k++;
            twin = i;                    //把twin更新为目前的素数
        }
    printf("%d", k);
}

int checkprime(int n) {
    int i, j;
    for (i = 0; i <= n; i++)
        prime[i] = 1;
    for (i = 2; i * i <= n; i++)
        if (prime[i])
            for (j = i; j * i <= n; j++)
                prime[j * i] = 0;

}