Euler 27 二次“素数生成”多项式

最新推荐文章于 2022-05-10 20:19:35 发布

银-TaMa

最新推荐文章于 2022-05-10 20:19:35 发布

阅读量791

点赞数

分类专栏： Euler ｃ学习

本文链接：https://blog.csdn.net/weixin_39722329/article/details/86093530

版权

博客介绍了如何解决Euler 27问题，涉及关键点包括使用米勒测试快速判断素数以及优化剪枝策略。通过米勒测试根据费马小定理判断素数，确保函数输入范围限定，并讨论了剪枝策略，如当i=0时b必须为素数，i=1时a+b+1也需为素数。最后，提出了一个素数筛的变形来找到合数的最小素因子。

摘要由CSDN通过智能技术生成

题目链接

关于这道题目

需要解决的点主要是

判断一个数是不是素数
如何有效的减枝

判断一个数是不是素数

判断一个数是不是素数有很多种做法，比如试除法

今天呢使用一个更快的办法米勒测试

主要的原理是根据费马小定理

定理描述

当且仅当　Ｐ　为素数时：

ａ^p-1 mod P 为　１

１　<= a <= p - 1

那么我们只需要选取若干个ａ，代入以上公式

求得结果，若均为１，说明　Ｐ　在很大概率上是个素数

代码如下

/*************************************************************************
	> File Name: 1.c
	> Author:Gin.TaMa 
	> Mail:1137554811@qq.com 
	> Created Time: 2019年01月08日 星期二 11时52分13秒
 ************************************************************************/

#include<stdio.h>
#include<stdlib.h>
#include<time.h>
#define test_round 30

int r_m_test(int x){
   
    if(x <= 1)return 0;
    //注意对函数的输入范围进行限制，这才是一个良好的函数定义方法
    int a , m, ans ;
    for(int i = 0;i < test_round;i ++){
   
        a = rand()%(x - 1) + 1;
        m = x - 1;
        ans = 1;
        while(m){
   
            if(m & 1)ans = ans * a % x;
            a = a * a % x;
            m >>= 1;
        }
        if(ans != 1

最低0.47元/天解锁文章