刷题记录-自测-2 素数对猜想

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

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

输入格式:

输入在一行给出正整数N

输出格式:

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

输入样例:

20

输出样例:

4

 

#include <stdio.h>
#include <math.h>

int prime(int num);

int main()
{
    int N,i,past=2,cnt=0;
    //输入范围设置
    scanf("%d",&N);
    //不断查询奇数是否为素数
    for(i=3;i<=N;i+=2)
    {
        if(prime(i))
        {
            //若差为2则奇数加1
            if((i-past)==2)
                cnt++;
            past=i;
        }
    }
    printf("%d",cnt);
    return 0;
}

int prime(int num)
{
    int i;
    for(i=2;i<=sqrt(num);i++)
    {
        if(num%i==0)
            return 0;
    }
    return 1;
}
  • 0
    点赞
  • 0
    收藏
    觉得还不错? 一键收藏
  • 0
    评论
好的,我可以为您提供一个C++程序,用三种方法对奇素数的素性进行判断,分别是Fermat测试、Solovay-Strassen测试和Miller-Rabin测试。 Fermat测试的原理是利用费马小定理,对于素数p和任意整数a,有$a^{p-1} \equiv 1 \pmod{p}$,因此如果对于某个整数a,有$a^{p-1} \not\equiv 1 \pmod{p}$,则p不是素数。Fermat测试的缺点是存在卡迈克尔数(Carmichael number),即有些合数也会通过Fermat测试。 Solovay-Strassen测试的原理是利用欧拉判别准则(Euler's criterion),对于素数p和任意整数a,有$a^{\frac{p-1}{2}} \equiv Legendre(a,p) \pmod{p}$,其中Legendre符号是二次剩余的符号。因此,如果对于某个整数a,有$a^{\frac{p-1}{2}} \not\equiv Legendre(a,p) \pmod{p}$,则p不是素数。 Miller-Rabin测试的原理是利用二次探测(quadratic probing)和费马小定理,对于奇素数p和任意整数a,可以将p-1表示为$2^s \times d$的形式,其中d是奇数。然后对于任意整数r,如果存在某个整数t满足$a^{2^tr} \equiv -1 \pmod{p}$,则p可能是素数;否则,如果对于所有的r,有$a^{2^tr} \not\equiv -1 \pmod{p}$,则p不是素数。 下面是完整的C++代码实现: ```c++ #include <iostream> #include <cmath> using namespace std; // 计算幂次方 int power(int a, int b, int p) { int ans = 1 % p; a %= p; while (b) { if (b & 1) ans = (long long)ans * a % p; a = (long long)a * a % p; b >>= 1; } return ans; } // 计算Legendre符号 int legendre(int a, int p) { return power(a, (p - 1) / 2, p); } // Fermat测试 bool fermatTest(int n, int k) { if (n == 2) return true; if (n <= 1 || n % 2 == 0) return false; for (int i = 0; i < k; i++) { int a = rand() % (n - 2) + 2; if (power(a, n - 1, n) != 1) return false; } return true; } // Solovay-Strassen测试 bool solovayStrassenTest(int n, int k) { if (n == 2) return true; if (n <= 1 || n % 2 == 0) return false; for (int i = 0; i < k; i++) { int a = rand() % (n - 2) + 2; int x = power(a, (n - 1) / 2, n); if (x == 0 || x != legendre(a, n)) return false; } return true; } // Miller-Rabin测试 bool millerRabinTest(int n, int k) { if (n == 2) return true; if (n <= 1 || n % 2 == 0) return false; int s = 0, d = n - 1; while (d % 2 == 0) { s++; d /= 2; } for (int i = 0; i < k; i++) { int a = rand() % (n - 2) + 2; int x = power(a, d, n); if (x == 1 || x == n - 1) continue; for (int j = 0; j < s - 1; j++) { x = (long long)x * x % n; if (x == n - 1) break; } if (x != n - 1) return false; } return true; } int main() { srand(time(NULL)); // 随机数种子 int n; cout << "输入一个奇素数(n>=3):" << endl; cin >> n; if (fermatTest(n, 10)) { cout << n << "是素数(Fermat测试)" << endl; } else { cout << n << "不是素数(Fermat测试)" << endl; } if (solovayStrassenTest(n, 10)) { cout << n << "是素数(Solovay-Strassen测试)" << endl; } else { cout << n << "不是素数(Solovay-Strassen测试)" << endl; } if (millerRabinTest(n, 10)) { cout << n << "是素数(Miller-Rabin测试)" << endl; } else { cout << n << "不是素数(Miller-Rabin测试)" << endl; } return 0; } ``` 这段代码可以接受一个奇素数作为输入,并分别使用Fermat测试、Solovay-Strassen测试和Miller-Rabin测试判断是否为素数。您可以根据需要进行修改和扩展。注意,此代码仅用于演示目的,不适用于实际应用。

“相关推荐”对你有帮助么?

  • 非常没帮助
  • 没帮助
  • 一般
  • 有帮助
  • 非常有帮助
提交
评论
添加红包

请填写红包祝福语或标题

红包个数最小为10个

红包金额最低5元

当前余额3.43前往充值 >
需支付:10.00
成就一亿技术人!
领取后你会自动成为博主和红包主的粉丝 规则
hope_wisdom
发出的红包
实付
使用余额支付
点击重新获取
扫码支付
钱包余额 0

抵扣说明:

1.余额是钱包充值的虚拟货币,按照1:1的比例进行支付金额的抵扣。
2.余额无法直接购买下载,可以购买VIP、付费专栏及课程。

余额充值