do(n)= my(k,z); if (k=ispower(n,,&z), print([k,z]));
{v = [
3^101*2^20,
3^101*2^101,
3^101*2^101*5,
2^20*3^15*5^5,
2^3*103^6,
2^5*103^10,
2^7*103^10,
2^15*103^10,
2^21*103^42,
2^35*103^70,
2^105*103^210,
103^100003
]; }

for (i=1, #v, do(v[i]))
for (i=1, 10^6, do(i))
\\#1259
do(-16)
do(-64)
do(-8)
do(-8 * 103^6)
do(-1/64)
