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 r2373864 = 1.0;
        double r2373865 = s;
        double r2373866 = -r2373865;
        double r2373867 = exp(r2373866);
        double r2373868 = r2373864 + r2373867;
        double r2373869 = r2373864 / r2373868;
        double r2373870 = c_p;
        double r2373871 = pow(r2373869, r2373870);
        double r2373872 = r2373864 - r2373869;
        double r2373873 = c_n;
        double r2373874 = pow(r2373872, r2373873);
        double r2373875 = r2373871 * r2373874;
        double r2373876 = t;
        double r2373877 = -r2373876;
        double r2373878 = exp(r2373877);
        double r2373879 = r2373864 + r2373878;
        double r2373880 = r2373864 / r2373879;
        double r2373881 = pow(r2373880, r2373870);
        double r2373882 = r2373864 - r2373880;
        double r2373883 = pow(r2373882, r2373873);
        double r2373884 = r2373881 * r2373883;
        double r2373885 = r2373875 / r2373884;
        return r2373885;
}