Add a function to decompose a number as a product of powers of primes
Importing solution for p5