python求小于n的所有素数_Python素数检测的方法

本文实例讲述了Python素数检测的方法。分享给大家供大家参考。具体如下:

因子检测:

检测因子,时间复杂度O(n^(1/2))

1

2

3

4

5

6

7

def is_prime(n):

if n <2:

return False

for iin xrange(2,int(n**0.5+1)):

if n%i== 0:

return False

return True

费马小定理:

如果n是一个素数,a是小于n的任意正整数,那么a的n次方与a模n同余

实现方法:

选择一个底数(例如2),对于大整数p,如果2^(p-1)与1不是模p同余数,则p一定不是素数;否则,则p很可能是一个素数

2**(n-1)%n 不是一个容易计算的数字

模运算规则:

1

2

(a^b)% p= ((a% p)^b)% p

(a* b)% p= (a% p* b% p)% p

计算X^N(% P)

可以

如果N是偶数,那么X^N =(X*X)^[N/2];

如果N是奇数,那么X^N = X*X^(N-1) = X *(X*X)^[N/2];

1

2

3

4

5

6

7

def xn_mod_p(x, n, p):

if n== 0:

return 1

res= xn_mod_p((x*x)%p, n>>1, p)

if n&1 != 0:

res= (res*x)%p

return res

也可以归纳为下面的算法 两个函数是一样的

1

2

3

4

5

6

7

8

def xn_mod_p2(x, n, p):

res= 1

n_bin= bin(n)[2:]

for iin range(0,len(n_bin)):

res= res**2 % p

if n_bin[i]== '1':

res= res* x% p

return res

有了模幂运算快速处理就可以实现费马检测

费马测试当给出否定结论时,是准确的,但是肯定结论有可能是错误的,对于大整数的效率很高,并且误判率随着整数的增大而降低

1

2

3

4

5

6

7

def fermat_test_prime(n):

if n== 1:

return False

if n== 2:

return True

res= xn_mod_p(2, n-1, n)

return res== 1

MILLER-RABIN检测

Miller-Rabin检测是目前应用比较广泛的一种

二次探测定理:如果p是一个素数,且0

费马小定理:a^(p-1) ≡ 1(mod p)

这就是Miller-Rabin素性测试的方法。不断地提取指数n-1中的因子2,把n-1表示成d*2^r(其中d是一个奇数)。那么我们需要计算的东西就变成了a的d*2^r次方除以n的余数。于是,a^(d * 2^(r-1))要么等于1,要么等于n-1。如果a^(d * 2^(r-1))等于1,定理继续适用于a^(d * 2^(r-2)),这样不断开方开下去,直到对于某个i满足a^(d * 2^i) mod n = n-1或者最后指数中的2用完了得到的a^d mod n=1或n-1。这样,Fermat小定理加强为如下形式:

尽可能提取因子2,把n-1表示成d*2^r,如果n是一个素数,那么或者a^d mod n=1,或者存在某个i使得a^(d*2^i) mod n=n-1 ( 0<=i

定理:若n是素数,a是小于n的正整数,则n对以a为基的Miller测试,结果为真.

Miller测试进行k次,将合数当成素数处理的错误概率最多不会超过4^(-k)

1

2

3

4

5

6

7

8

9

10

11

12

13

14

15

16

17

18

19

20

21

22

23

24

25

26

27

28

29

30

def miller_rabin_witness(a, p):

if p== 1:

return False

if p== 2:

return True

#p-1 = u*2^t 求解 u, t

n= p- 1

t= int(math.floor(math.log(n,2)))

u= 1

while t >0:

u= n/ 2**t

if n% 2**t== 0 and u% 2 == 1:

break

t= t- 1

b1= b2= xn_mod_p2(a, u, p)

for iin range(1, t+ 1):

b2= b1**2 % p

if b2== 1 and b1 != 1 and b1 != (p- 1):

return False

b1= b2

if b1 != 1:

return False

return True

def prime_test_miller_rabin(p, k):

while k >0:

a= randint(1, p- 1)

if not miller_rabin_witness(a, p):

return False

k= k- 1

return True

希望本文所述对大家的Python程序设计有所帮助。

  • 0
    点赞
  • 1
    收藏
    觉得还不错? 一键收藏
  • 0
    评论
大学生参加学科竞赛有着诸多好处,不仅有助于个人综合素质的提升,还能为未来职业发展奠定良好基础。以下是一些分析: 首先,学科竞赛是提高专业知识和技能水平的有效途径。通过参与竞赛,学生不仅能够深入学习相关专业知识,还能够接触到最新的科研成果和技术发展趋势。这有助于拓展学生的学科视野,使其对专业领域有更深刻的理解。在竞赛过程中,学生通常需要解决实际问题,这锻炼了他们独立思考和解决问题的能力。 其次,学科竞赛培养了学生的团队合作精神。许多竞赛项目需要团队协作来完成,这促使学生学会有效地与他人合作、协调分工。在团队合作中,学生们能够学到如何有效沟通、共同制定目标和分工合作,这对于日后进入职场具有重要意义。 此外,学科竞赛是提高学生综合能力的一种途径。竞赛项目通常会涉及到理论知识、实际操作和创新思维等多个方面,要参赛者具备全面的素质。在竞赛过程中,学生不仅需要展现自己的专业知识,还需要具备创新意识和解决问题的能力。这种全面的综合能力培养对于未来从事各类职业都具有积极作用。 此外,学科竞赛可以为学生提供展示自我、树立信心的机会。通过比赛的舞台,学生有机会展现自己在专业领域的优势,得到他人的认可和赞誉。这对于培养学生的自信心和自我价值感非常重要,有助于他们更加积极主动地投入学习和未来的职业生涯。 最后,学科竞赛对于个人职业发展具有积极的助推作用。在竞赛中脱颖而出的学生通常能够引起企业、研究机构等用人单位的关注。获得竞赛奖项不仅可以作为个人履历的亮点,还可以为进入理想的工作岗位提供有力的支持。
评论
添加红包

请填写红包祝福语或标题

红包个数最小为10个

红包金额最低5元

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

抵扣说明:

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

余额充值