让我们定义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)
是一种用于筛选素数的算法。它的基本思想是从小到大遍历所有的数,如果当前数是素数,则将其所有的倍数标记为非素数。通过这种方式,最终剩下的未被标记的数就是素数。
具体的步骤如下:
- 创建一个长度为n+1的布尔数组a,初始时将所有元素都标记为素数(即设为true)。
- 从2开始遍历到n,如果a[i]为素数(即true),则将其所有的倍数(大于i且小于等于n)标记为非素数(即设为false)。
- 遍历完所有的数后,剩下的未被标记的数就是素数。
埃拉托斯特尼筛法的时间复杂度为O(nloglogn),效率较高。但是,它只适用于筛选较小范围内的素数。
Miller-Rabin素数测试
是一种用于判断一个数是否为素数的概率性算法。它基于一个数论的定理,即对于一个大于2的奇数n,如果存在一个整数a,使得a(n-1) ≡ 1 (mod n)且a((n-1)/2) ≡ -1 (mod n),则n可能是一个素数。
具体的步骤如下:
- 将待判断的数n分解为n-1 = 2^s * d,其中d为奇数。
- 选择一个随机数a,满足1 < a < n-1。
- 计算x = a^d mod n。
- 如果x ≡ 1 (mod n)或x ≡ -1 (mod n),则n可能是一个素数。
- 重复步骤3和步骤4,共进行s次。
- 如果在s次迭代中,都满足x ≡ 1 (mod n)或x ≡ -1 (mod n),则n可能是一个素数。
- 如果在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);
}