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 r49542010 = 1.0;
        double r49542011 = s;
        double r49542012 = -r49542011;
        double r49542013 = exp(r49542012);
        double r49542014 = r49542010 + r49542013;
        double r49542015 = r49542010 / r49542014;
        double r49542016 = c_p;
        double r49542017 = pow(r49542015, r49542016);
        double r49542018 = r49542010 - r49542015;
        double r49542019 = c_n;
        double r49542020 = pow(r49542018, r49542019);
        double r49542021 = r49542017 * r49542020;
        double r49542022 = t;
        double r49542023 = -r49542022;
        double r49542024 = exp(r49542023);
        double r49542025 = r49542010 + r49542024;
        double r49542026 = r49542010 / r49542025;
        double r49542027 = pow(r49542026, r49542016);
        double r49542028 = r49542010 - r49542026;
        double r49542029 = pow(r49542028, r49542019);
        double r49542030 = r49542027 * r49542029;
        double r49542031 = r49542021 / r49542030;
        return r49542031;
}