d=[]

for i in range(1,200001):

 b=str(i)

 c=0

 for j in b:

  a=1

  for l in range(1,int(j)+1):

   a=l*a

  c=c+a

 if i==c:

  d.append(i)

print(d)

d=[1,2,145,40585]

，，，，，，，，，，，，，，，

def fun(x):

if x==1 or x==0:

return 1

else:

return x*fun(x-1)

d=[]

for i in range(1,200001):

b=str(i)

c=0

for j in b:

k=int(j)

c=fun(k)+c

if i==c:

d.append(i)

print(d)