Cannot sample enough valid points. (more)

\[0 \lt c_p \land 0 \lt c_n\]
\[\frac{{\left(\frac{1}{1 + e^{-s}}\right)}^{c_p} \cdot {\left(1 - \frac{1}{1 + e^{-s}}\right)}^{c_n}}{{\left(\frac{1}{1 + e^{-t}}\right)}^{c_p} \cdot {\left(1 - \frac{1}{1 + e^{-t}}\right)}^{c_n}}\]
\frac{{\left(\frac{1}{1 + e^{-s}}\right)}^{c_p} \cdot {\left(1 - \frac{1}{1 + e^{-s}}\right)}^{c_n}}{{\left(\frac{1}{1 + e^{-t}}\right)}^{c_p} \cdot {\left(1 - \frac{1}{1 + e^{-t}}\right)}^{c_n}}
double f(double c_p, double c_n, double t, double s) {
        double r4793567 = 1.0;
        double r4793568 = s;
        double r4793569 = -r4793568;
        double r4793570 = exp(r4793569);
        double r4793571 = r4793567 + r4793570;
        double r4793572 = r4793567 / r4793571;
        double r4793573 = c_p;
        double r4793574 = pow(r4793572, r4793573);
        double r4793575 = r4793567 - r4793572;
        double r4793576 = c_n;
        double r4793577 = pow(r4793575, r4793576);
        double r4793578 = r4793574 * r4793577;
        double r4793579 = t;
        double r4793580 = -r4793579;
        double r4793581 = exp(r4793580);
        double r4793582 = r4793567 + r4793581;
        double r4793583 = r4793567 / r4793582;
        double r4793584 = pow(r4793583, r4793573);
        double r4793585 = r4793567 - r4793583;
        double r4793586 = pow(r4793585, r4793576);
        double r4793587 = r4793584 * r4793586;
        double r4793588 = r4793578 / r4793587;
        return r4793588;
}