pd(n)={
  my(f=factor(n),v=f[,1]~);
  for(i=1,#v,
    while(f[i,2]--,
      v=concat(v,f[i,1])
    )
  );
  vecsort(v)
};