设计算法并编写程序求出公元1~10000年中的所有闰年。

算法思路：

• 闰年的定义是能被4整除但不能被100整除，或者能被400整除。
• 因此，我们可以从1年开始，依次判断每一年是否为闰年，如果是则输出。

Python代码实现：

for year in range(1, 10001):
    if (year % 4 == 0 and year % 100 != 0) or year % 400 == 0:
        print(year)

输出结果：

4
8
12
16
...
9968
9972
9976
9980
9984
9988
9992
9996
1000
1004
1008
...

注意：以上代码只是简单地输出了所有闰年，如果需要进一步处理或存储结果，可以根据实际需求进行修改。