PAT-B 1007. 素数对猜想

最新推荐文章于 2024-05-23 12:04:54 发布

不明真相的板蓝根

最新推荐文章于 2024-05-23 12:04:54 发布

阅读量600

点赞数

CC 4.0 BY-SA版权

分类专栏： PAT-B 文章标签： PAT 乙级简洁代码

本文链接：https://blog.youkuaiyun.com/plank_root/article/details/51263742

PAT-B 专栏收录该内容

62 篇文章

订阅专栏

本文介绍了一种算法实现方式，用于解决素数对猜想问题，即找出小于等于给定正整数N的所有素数对，这些素数对的差为2。通过生成素数列表并检查相邻素数之间的差值来实现。

题目内容：

让我们定义 $d_n$ 为： $d_n = p_{n+1} - p_n$ ，其中 $p_i$ 是第 $i$ 个素数。显然有 $d_1=1$ 且对于 $n>1$ 有 $d_n$ 是偶数。“素数对猜想”认为“存在无穷多对相邻且差为2的素数”。

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

输入格式：每个测试输入包含1个测试用例，给出正整数N。

输出格式：每个测试用例的输出占一行，不超过N的满足猜想的素数对的个数。

输入样例：
20
输出样例：
4

思路分析：

本题先根据给定的正整数N，生成包含所有在1–N之间的素数的列表。
然后根据生成的列表，查找是否有差为2的相邻两个素数，并计数。

生成素数表的方法：

根据唯一分解定理：任何合数都可以分解为若干个素数的乘积，所以对于任意的正整数n，只要可以不被所有比他小的素数整除，即可证明n是素数。
合数分解得到的若干素数中，最大值不会超过合数的平方根，所以验证到 list[i] * list[i] <= n 就可以了.
因为所有2得倍数都不是素数，所以验证时 i递增2 ，预设2为素数，验证从3开始。

代码：

#include <stdio.h>

int is_prime(int n, int list[])
{
    for (int i = 0; list[i] * list[i] <= n; i++) // list[i] * list[i] <= n时退出循环
        if (n % list[i] == 0) // 验证能否被素数表中的数整除
            return 0;
    return 1;
}

int generate_prime_list(int n, int list[])
{
    int m = 1;
    for (int i = 3; i <= n; i += 2) // 从3开始验证
        if (is_prime(i, list))
            list[m++] = i; // m为已经生成的素数数量
    return m;
}

int main()
{
    int n, size_of_list, prime_couple = 0;
    int prime_list[10000] = {2}; // 预设2为素数表中的第一个数

    scanf("%d", &n);
    size_of_list = generate_prime_list(n, prime_list);

    for (int i = 1; i < size_of_list; i++)
        if (prime_list[i] - prime_list[i - 1] == 2)
            prime_couple++;
    printf("%d", prime_couple);

    return 0;
}

点这里进入试题网页