def pri_num(val): #求出該數字的所有素數
for i in range(2,val+1):
if i<=4:
j=i
else:
j=i//2
for k in range(2,j):
if i%k==0:
break
else:
num.append(i)
def sol_factor(val,k): #求因數
for i in k: #遍歷素數
if val%i==0:
return i
break
else:
return 1
for input_num in range(2,1000):
sum_num=1
num=[] #該數字范圍內的所有素數
fac_num=[] #該數字的因數
pri_num(input_num) #求出該數字的所有素數
decom_num=input_num #要分解的數
while True:
res=sol_factor(decom_num,num)
if res==decom_num: #遍歷結束
sum_num +=res
fac_num.append(str(res))
break
else:
sum_num +=res
fac_num.append(str(res))
decom_num= decom_num/res
if sum_num==input_num:
print(input_num,'=1+','+'.join(fac_num)) #分解后的表達式