软工日报5.28

8-7 【Python0008】筛法求素数

分数 10

全屏浏览

作者 doublebest

单位 石家庄铁道大学

【题目描述】用户输入整数n和m（1<n<m<1000），应用筛法求[n,m]范围内的所有素数。
【练习要求】请给出源代码程序和运行测试结果，源代码程序要求添加必要的注释。
【输入格式】在一行中输入整数n和m，数据间以空格为间隔。
【输出格式】按照一行5个的格式输出[n,m]范围内的所有素数，其中每个素数的显示格式为：占5位，右对齐。
【输入样例】50 100
【输出样例】 53 59 61 67 71
73 79 83 89 97

 

def sieve_of_eratosthenes(n, m):

 

"""

 

使用筛法求解 \[n,m\] 范围内的所有素数

 

 

 

参数:

 

n (int): 范围的左端点

 

m (int): 范围的右端点

 

 

 

返回:

 

list: 素数列表

 

"""

 

# 初始化标记数组

 

is_prime = \[True\] \* (m + 1)

 

is_prime\[0\] = is_prime\[1\] = False

 

 

 

# 标记非素数

 

for i in range(2, int(m \*\* 0.5) + 1):

 

if is_prime\[i\]:

 

# 将 i 的倍数都标记为非素数

 

for j in range(i \* i, m + 1, i):

 

is_prime\[j\] = False

 

 

 

# 返回素数列表

 

primes = \[i for i in range(n, m + 1) if is_prime\[i\]\]

 

return primes

 

\# 获取用户输入

 

n, m = map(int, input("请输入 n 和 m (1 < n < m < 1000): ").split())

 

\# 使用筛法求解素数并输出结果

 

primes = sieve_of_eratosthenes(n, m)

 

for i, prime in enumerate(primes):

 

# 使用占 5 位、右对齐的格式输出每个素数

 

print(f"{prime:5}", end="")

 

# 每输出 5 个素数就换行

 

if (i + 1) % 5 == 0:

 

print()

 

\# 如果最后一行没有输出 5 个素数，就换行

 

if len(primes) % 5 != 0:

 

print()

 代码量 50行

时间 0.5小时