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 r50548705 = 1.0;
        double r50548706 = s;
        double r50548707 = -r50548706;
        double r50548708 = exp(r50548707);
        double r50548709 = r50548705 + r50548708;
        double r50548710 = r50548705 / r50548709;
        double r50548711 = c_p;
        double r50548712 = pow(r50548710, r50548711);
        double r50548713 = r50548705 - r50548710;
        double r50548714 = c_n;
        double r50548715 = pow(r50548713, r50548714);
        double r50548716 = r50548712 * r50548715;
        double r50548717 = t;
        double r50548718 = -r50548717;
        double r50548719 = exp(r50548718);
        double r50548720 = r50548705 + r50548719;
        double r50548721 = r50548705 / r50548720;
        double r50548722 = pow(r50548721, r50548711);
        double r50548723 = r50548705 - r50548721;
        double r50548724 = pow(r50548723, r50548714);
        double r50548725 = r50548722 * r50548724;
        double r50548726 = r50548716 / r50548725;
        return r50548726;
}