求区间素数的埃拉托斯特尼筛法

题目描述
输入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)