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 r33972993 = 1.0;
        double r33972994 = s;
        double r33972995 = -r33972994;
        double r33972996 = exp(r33972995);
        double r33972997 = r33972993 + r33972996;
        double r33972998 = r33972993 / r33972997;
        double r33972999 = c_p;
        double r33973000 = pow(r33972998, r33972999);
        double r33973001 = r33972993 - r33972998;
        double r33973002 = c_n;
        double r33973003 = pow(r33973001, r33973002);
        double r33973004 = r33973000 * r33973003;
        double r33973005 = t;
        double r33973006 = -r33973005;
        double r33973007 = exp(r33973006);
        double r33973008 = r33972993 + r33973007;
        double r33973009 = r33972993 / r33973008;
        double r33973010 = pow(r33973009, r33972999);
        double r33973011 = r33972993 - r33973009;
        double r33973012 = pow(r33973011, r33973002);
        double r33973013 = r33973010 * r33973012;
        double r33973014 = r33973004 / r33973013;
        return r33973014;
}