因數分解

for n in [1500,3500,6000]:

    nb = n

    a = []

    for i in range(2,int(n**0.5)+1):

        while n%i==0:

            a.append(i)   

            n = n /i

    

    out = []

    s = set(a)

    for i in s:

        t = f'{i}^{a.count(i)}'

        out.append(t)

    out = '*'.join(out)

    out = str(nb) + '=' + out

    out = out.replace('^1','')

    print(out)