题目描述
输入M、N,1 < M < N < 1000000,求区间[M,N]内的所有素数的个数。素数定义:除了1以外,只能被1和自己整除的自然数称为素数
输入描述:
两个整数M,N
输出描述:
区间内素数的个数
使用埃拉托斯特尼筛法
基本思想:素数的倍数一定不是素数
实现方法:用一个长度为N+1的数组保存信息(0表示素数,1表示非素数),先假设所有的数都是素数(初始化为0),
从第一个素数2开始,把2的倍数都标记为非素数(置为1),一直到大于N;
然后进行下一趟,找到2后面的下一个素数3,进行同样的处理,直到最后,数组中依然为0的数即为素数。
如果n*n > 范围最大值就跳出。
m,n=list(map(int,input().split()))
def prime(n):
original = [0 for i in range(n)] # 列表中的n个元素都为false
for i in range(2,n):
if(not original[i]):
for j in range(i*i,n,i): #优化start从i*i开始而不是i*2,例如5,只要标记5*5,5*6,5*7…为合数,因为5*2,5*3,5*4…已经被之前出现的数的倍数标记过了
original[j]=1
return n-sum(original)
print(prime(n)-prime(m))
埃拉托斯特尼筛法
# python 原生实现
def primes(n):
P = []
f = []
for i in range(n+1):
if i > 2 and i%2 == 0:
f.append(1)
else:
f.append(0)
i = 3
while i*i <= n:
if f[i] == 0:
j = i*i
while j <= n:
f[j] = 1
j += i+i
i += 2
P.append(2)
for x in range(3,n+1,2):
if f[x] == 0:
P.append(x)
return P
n = 100 #100以内的素数
P = primes(n)
print (P)