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 r75723747 = 1.0;
        double r75723748 = s;
        double r75723749 = -r75723748;
        double r75723750 = exp(r75723749);
        double r75723751 = r75723747 + r75723750;
        double r75723752 = r75723747 / r75723751;
        double r75723753 = c_p;
        double r75723754 = pow(r75723752, r75723753);
        double r75723755 = r75723747 - r75723752;
        double r75723756 = c_n;
        double r75723757 = pow(r75723755, r75723756);
        double r75723758 = r75723754 * r75723757;
        double r75723759 = t;
        double r75723760 = -r75723759;
        double r75723761 = exp(r75723760);
        double r75723762 = r75723747 + r75723761;
        double r75723763 = r75723747 / r75723762;
        double r75723764 = pow(r75723763, r75723753);
        double r75723765 = r75723747 - r75723763;
        double r75723766 = pow(r75723765, r75723756);
        double r75723767 = r75723764 * r75723766;
        double r75723768 = r75723758 / r75723767;
        return r75723768;
}