让我们定义dn为:dn =pn+1−pn,其中pi是第i个素数。显然有d1=1,且对于n>1有dn是偶数。“素数对猜想”认为“存在无穷多对相邻且差为2的素数”。
现给定任意正整数N(<105 ),请计算不超过N的满足猜想的素数对的个数。
输入格式:
输入在一行给出正整数N。
输出格式:
在一行中输出不超过N的满足猜想的素数对的个数。
输入样例:
20
输出样例:
4
PTA运行超时,我写的:
num = int(input())
count = 0
for i in range(3, num, 2):
for j in range(2, i):
if i % j == 0:
break
else:
count += 1
print((count+1)//2)
请教了一个大神,得到下面的算法,运行几次可以得到通过,PTA神奇了…
大神写的:
num = int(input())
count = 0
P = 3
for i in range(3, num+1, 2):
for j in range(3, int(i**0.5)+1, 2):
if i % j == 0:
break
else:
if i - P == 2:
count += 1
P = i
print(count)
CDN大牛写的:
import math
def is_prime(number, primes):
up_data = int(math.sqrt(number))
for prime in primes:
if prime > up_data:
break
if number % prime == 0:
return False
return True
up_limit = int(input())
last_prime, next_prime = 2, 3
count = 0
primes = [2, 3]
if up_limit > 4:
for number in range(5, up_limit+1):
if is_prime(number, primes):
primes.append(number)
last_prime = next_prime
next_prime = number
if (next_prime - last_prime) == 2:
count += 1
print(str(count))
求优化:请不吝赐教。