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 r17872801 = 1.0;
        double r17872802 = s;
        double r17872803 = -r17872802;
        double r17872804 = exp(r17872803);
        double r17872805 = r17872801 + r17872804;
        double r17872806 = r17872801 / r17872805;
        double r17872807 = c_p;
        double r17872808 = pow(r17872806, r17872807);
        double r17872809 = r17872801 - r17872806;
        double r17872810 = c_n;
        double r17872811 = pow(r17872809, r17872810);
        double r17872812 = r17872808 * r17872811;
        double r17872813 = t;
        double r17872814 = -r17872813;
        double r17872815 = exp(r17872814);
        double r17872816 = r17872801 + r17872815;
        double r17872817 = r17872801 / r17872816;
        double r17872818 = pow(r17872817, r17872807);
        double r17872819 = r17872801 - r17872817;
        double r17872820 = pow(r17872819, r17872810);
        double r17872821 = r17872818 * r17872820;
        double r17872822 = r17872812 / r17872821;
        return r17872822;
}