编写程序，通过循环结构计算全部水仙花数。提示：水仙花数是一个三位数，该数正好等于组成该三位数的各位数字的立方和

以下 Python 代码实现：

for num in range(100, 1000):
    # 分别计算百位、十位、个位上的数字的立方和
    digit1 = num // 100
    digit2 = (num // 10) % 10
    digit3 = num % 10
    sum_of_cubes = digit1 ** 3 + digit2 ** 3 + digit3 ** 3
    if num == sum_of_cubes:
        print(num)

输出结果：

153
370
371
407