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 r6728262 = 1.0;
        double r6728263 = s;
        double r6728264 = -r6728263;
        double r6728265 = exp(r6728264);
        double r6728266 = r6728262 + r6728265;
        double r6728267 = r6728262 / r6728266;
        double r6728268 = c_p;
        double r6728269 = pow(r6728267, r6728268);
        double r6728270 = r6728262 - r6728267;
        double r6728271 = c_n;
        double r6728272 = pow(r6728270, r6728271);
        double r6728273 = r6728269 * r6728272;
        double r6728274 = t;
        double r6728275 = -r6728274;
        double r6728276 = exp(r6728275);
        double r6728277 = r6728262 + r6728276;
        double r6728278 = r6728262 / r6728277;
        double r6728279 = pow(r6728278, r6728268);
        double r6728280 = r6728262 - r6728278;
        double r6728281 = pow(r6728280, r6728271);
        double r6728282 = r6728279 * r6728281;
        double r6728283 = r6728273 / r6728282;
        return r6728283;
}