void prodprim(unsigned int n, unsigned int &p)
{
    p = 1;
    int div = 2, exp = 0;
    while(n > 1)
    {
        while(n % div == 0)
        {
            n = n / div;
            exp++;
        }
        if(exp > 0)
        {
            p = p * div;
            cout << div << " ^ " << exp << "\n";
        }
        div = div + 1;
        exp = 0;
    }
}