Simulates from the generalized gamma distribution with density
g(x) = exp(-b*x^p)*x^(a-1)/K_3, x>0,
where K_3 is a normalizing constant. The mathodology is explained in Section 4 of Grabchak (2021). This distribution is needed to simulate p-RDTS random variables with negative alpha values.