factorint_as_string()

["Ruud H.G. van Tol" <rvtol@isolution.nl>]

sf() = my(v=Vec(factorint(n)~)); if(!#v, v=Vec([1;1])); for(i=1, #v, 
v[i]=if(v[i][2]>1, strjoin(v[i],"^"), v[i][1])); strjoin(["(", 
strjoin(v," * "), ")"])

Example:

? sf(3087)
%527 = "(3^2 * 7^3)"

-- Ruud