while n > 1:
    if not (n % i) and IsPrime(i):
        if i in b:
            continue
        else:
            b[i] = 1
            n /= i
            while not (n % i):
                b[i] += 1
                n /= i
    i += 1