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 r1880709 = 1.0;
        double r1880710 = s;
        double r1880711 = -r1880710;
        double r1880712 = exp(r1880711);
        double r1880713 = r1880709 + r1880712;
        double r1880714 = r1880709 / r1880713;
        double r1880715 = c_p;
        double r1880716 = pow(r1880714, r1880715);
        double r1880717 = r1880709 - r1880714;
        double r1880718 = c_n;
        double r1880719 = pow(r1880717, r1880718);
        double r1880720 = r1880716 * r1880719;
        double r1880721 = t;
        double r1880722 = -r1880721;
        double r1880723 = exp(r1880722);
        double r1880724 = r1880709 + r1880723;
        double r1880725 = r1880709 / r1880724;
        double r1880726 = pow(r1880725, r1880715);
        double r1880727 = r1880709 - r1880725;
        double r1880728 = pow(r1880727, r1880718);
        double r1880729 = r1880726 * r1880728;
        double r1880730 = r1880720 / r1880729;
        return r1880730;
}