将三位数看成由三个数组成，并判断其各位数字立方和是否等于该数本身:

#!/usr/bin/python
# -*- coding: UTF-8 -*-

for x in range(1,10):
    for y in range(0,10):
        for z in range(0,10):
            s1=x*100+y*10+z
            s2=pow(x,3)+pow(y,3)+pow(z,3)
            if s1==s2:
                print "水仙花数有：%7ld" %(s1)

水仙花数及其个数

#!/usr/bin/env python
# -*- coding: UTF-8 -*-

L = []
for n in range(100,1000):
	i = n / 100
	j = n /10 % 10
	k = n % 10
	if n == i ** 3 + j ** 3 + k ** 3:
		if n == pow(i,3) + pow(j,3) +pow(k,3):
			L.append(n)
print(L)
print len(L)