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 r3582263 = 1.0;
        double r3582264 = s;
        double r3582265 = -r3582264;
        double r3582266 = exp(r3582265);
        double r3582267 = r3582263 + r3582266;
        double r3582268 = r3582263 / r3582267;
        double r3582269 = c_p;
        double r3582270 = pow(r3582268, r3582269);
        double r3582271 = r3582263 - r3582268;
        double r3582272 = c_n;
        double r3582273 = pow(r3582271, r3582272);
        double r3582274 = r3582270 * r3582273;
        double r3582275 = t;
        double r3582276 = -r3582275;
        double r3582277 = exp(r3582276);
        double r3582278 = r3582263 + r3582277;
        double r3582279 = r3582263 / r3582278;
        double r3582280 = pow(r3582279, r3582269);
        double r3582281 = r3582263 - r3582279;
        double r3582282 = pow(r3582281, r3582272);
        double r3582283 = r3582280 * r3582282;
        double r3582284 = r3582274 / r3582283;
        return r3582284;
}