输出1000以内所有素数

输出1000以内所有的素数，即只有两个正因数（1和本身）的自然数，可以使用一种常见的算法，称为埃拉托斯特尼筛法（Sieve of Eratosthenes）。以下是Python的一个简单实现：

def find_primes(n):
    primes = [True] * (n+1)  # 初始化一个布尔数组，所有元素都设为True表示都是素数
    p = 2  # 首个质数从2开始
    while p * p <= n:
        if primes[p]:  # 如果p是素数
            for i in range(p * p, n + 1, p):  # 将p的所有倍数标记为非素数
                primes[i] = False
        p += 1
    # 打印出1000以内的素数
    for p in range(2, n):
        if primes[p]:
            print(p)

find_primes(1000)

运行这段代码会打印出1到1000之间的所有素数。