素数对猜想中运行超时问题的解决

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

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

输入格式:

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

输出格式:

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

输入样例:20

输出样例:4

最开始使用的代码如下：

#include <stdio.h>

int judgesushu(int number)
{
    for(int i=2;i<number;i++)
    {
        if(number%i==0)
            return 0;
    }
    return 1;
}
int main()
{
    int n,x,p1,p2;
    int sum=0;
    scanf("%d",&n);
    for(x=2;x<=n-2;x++)
    {
        p1=judgesushu(x);
        p2=judgesushu(x+2);
        if(p1==1&&p2==1)
            sum++;
    }
    printf("%d",sum);
    return 0;
}

大部分测试用例是通过了的，但是有那么两个死活过不了。原因是运行超时。

原来是因为普通的素数判断是遍历2到该数前一个数，当这个数较大时，就会导致运行超时，非常低效。

然后发现了埃拉托斯特尼筛法（Sieve of Eratosthenes）和Miller-Rabin素数测试这两种比较高效的方法。

埃拉托斯特尼筛法（Sieve of Eratosthenes）

是一种用于筛选素数的算法。它的基本思想是从小到大遍历所有的数，如果当前数是素数，则将其所有的倍数标记为非素数。通过这种方式，最终剩下的未被标记的数就是素数。

具体的步骤如下：

1. 创建一个长度为n+1的布尔数组a，初始时将所有元素都标记为素数（即设为true）。
2. 从2开始遍历到n，如果a[i]为素数（即true），则将其所有的倍数（大于i且小于等于n）标记为非素数（即设为false）。
3. 遍历完所有的数后，剩下的未被标记的数就是素数。

埃拉托斯特尼筛法的时间复杂度为O(nloglogn)，效率较高。但是，它只适用于筛选较小范围内的素数。

Miller-Rabin素数测试

是一种用于判断一个数是否为素数的概率性算法。它基于一个数论的定理，即对于一个大于2的奇数n，如果存在一个整数a，使得a(n-1) ≡ 1 (mod n)且a((n-1)/2) ≡ -1 (mod n)，则n可能是一个素数。

具体的步骤如下：

1. 将待判断的数n分解为n-1 = 2^s * d，其中d为奇数。
2. 选择一个随机数a，满足1 < a < n-1。
3. 计算x = a^d mod n。
4. 如果x ≡ 1 (mod n)或x ≡ -1 (mod n)，则n可能是一个素数。
5. 重复步骤3和步骤4，共进行s次。
6. 如果在s次迭代中，都满足x ≡ 1 (mod n)或x ≡ -1 (mod n)，则n可能是一个素数。
7. 如果在s次迭代中存在x不满足x ≡ 1 (mod n)且x ≡ -1 (mod n)，则n不是一个素数。

Miller-Rabin素数测试的时间复杂度取决于迭代次数s。通常情况下，迭代次数越多，判断结果越准确，但也会增加运算时间。因此，可以根据需求来选择合适的迭代次数。

所以接下来选择了埃拉托斯特尼筛法，以下是改进后的代码：

#include <stdio.h>
#include <math.h>
int main()
{
    int a[100000]={0};
    int n,sum=0;
    scanf("%d",&n);
    for(int i=2;i<=n;i++)
    {
        if(a[i]==0)
        {
            for(int j=i*2;j<=n;j+=i)
                a[j]=1;
        }
    };
    a[0]=1;
    a[1]=1;
    for(int i=1;i<=n-2;i++)
    {
        if(a[i]==0&&a[i+2]==0)
            sum++;
    }
    printf("%d",sum);
}