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 r3512257 = 1.0;
        double r3512258 = s;
        double r3512259 = -r3512258;
        double r3512260 = exp(r3512259);
        double r3512261 = r3512257 + r3512260;
        double r3512262 = r3512257 / r3512261;
        double r3512263 = c_p;
        double r3512264 = pow(r3512262, r3512263);
        double r3512265 = r3512257 - r3512262;
        double r3512266 = c_n;
        double r3512267 = pow(r3512265, r3512266);
        double r3512268 = r3512264 * r3512267;
        double r3512269 = t;
        double r3512270 = -r3512269;
        double r3512271 = exp(r3512270);
        double r3512272 = r3512257 + r3512271;
        double r3512273 = r3512257 / r3512272;
        double r3512274 = pow(r3512273, r3512263);
        double r3512275 = r3512257 - r3512273;
        double r3512276 = pow(r3512275, r3512266);
        double r3512277 = r3512274 * r3512276;
        double r3512278 = r3512268 / r3512277;
        return r3512278;
}