Project Eurler - Problem 7

By listing the first six prime numbers: 2, 3, 5, 7, 11, and 13, we can see that the 6th prime is 13. What is the 10001st prime number?

这里采用简单的办法求素数,不过还是要讲究一点效率:

1、2是唯一的偶数素数
2、所有比3大的素数可以表示为6k+1或者6k-1
3、判定一个数是否为素数的较快办法为从2到n的平方根,是否存在能整除n的因子,不存在则为素数。

def isPrime(n):
    if n == 1:
        return False
    elif n < 4:
        return True          # 2,3 都是素数
    elif n % 2 == 0:
        return False         # 除2外的偶数不是素数
    elif n < 9 == 0:
        return True          # 到此只可能是7了,是素数。
    elif n % 3 == 0:
        return False         # 使剩下的部分表示为6k加减1
    else:
        r = math.floor(math.sqrt(n))
        f = 5                  # 第一个6k-1
        while f <= r:
            if n % f == 0:
                return False
            if n % (f+2) == 0:  # 6k-1+2成为6k+1
                return False
            f += 6            # 步伐6大跃进
        return True

def the10001stPrime():
    p = 1
    count = 1
    while count < 10001:
        p += 2
        if isPrime(p):
            count += 1
    return p

  • 0
    点赞
  • 0
    收藏
    觉得还不错? 一键收藏
  • 0
    评论

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

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

请填写红包祝福语或标题

红包个数最小为10个

红包金额最低5元

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

抵扣说明:

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

余额充值