return x - (-log1p(- log_cdf ** (1.0 / a))) ** (1.0 / c) * scale