#include <tgmath.h>
#include <gmp.h>
#include <mpfr.h>
#include <stdio.h>
#include <stdbool.h>

char *name = "2nthrt (problem 3.4.6)";

double f_if(float x, float n) {
        float r21074 = x;
        float r21075 = 1;
        float r21076 = r21074 + r21075;
        float r21077 = n;
        float r21078 = r21075 / r21077;
        float r21079 = pow(r21076, r21078);
        float r21080 = pow(r21074, r21078);
        float r21081 = r21079 - r21080;
        return r21081;
}

double f_id(double x, double n) {
        double r21082 = x;
        double r21083 = 1;
        double r21084 = r21082 + r21083;
        double r21085 = n;
        double r21086 = r21083 / r21085;
        double r21087 = pow(r21084, r21086);
        double r21088 = pow(r21082, r21086);
        double r21089 = r21087 - r21088;
        return r21089;
}


double f_of(float x, float n) {
        float r21090 = 1;
        float r21091 = n;
        float r21092 = r21090 / r21091;
        float r21093 = -3.6227047746916703e-09;
        bool r21094 = r21092 <= r21093;
        float r21095 = x;
        float r21096 = r21090 + r21095;
        float r21097 = pow(r21096, r21092);
        float r21098 = pow(r21095, r21092);
        float r21099 = r21097 - r21098;
        float r21100 = exp(r21099);
        float r21101 = log(r21100);
        float r21102 = 6.875959904884827e-20;
        bool r21103 = r21092 <= r21102;
        float r21104 = r21095 * r21091;
        float r21105 = r21090 / r21104;
        float r21106 = 1/2;
        float r21107 = r21106 / r21095;
        float r21108 = r21107 / r21104;
        float r21109 = r21105 - r21108;
        float r21110 = log(r21095);
        float r21111 = r21110 / r21091;
        float r21112 = r21111 / r21104;
        float r21113 = r21109 - r21112;
        float r21114 = r21095 + r21090;
        float r21115 = pow(r21114, r21092);
        float r21116 = r21115 * r21115;
        float r21117 = r21098 * r21098;
        float r21118 = r21116 - r21117;
        float r21119 = r21115 + r21098;
        float r21120 = r21118 / r21119;
        float r21121 = r21103 ? r21113 : r21120;
        float r21122 = r21094 ? r21101 : r21121;
        return r21122;
}

double f_od(double x, double n) {
        double r21123 = 1;
        double r21124 = n;
        double r21125 = r21123 / r21124;
        double r21126 = -3.6227047746916703e-09;
        bool r21127 = r21125 <= r21126;
        double r21128 = x;
        double r21129 = r21123 + r21128;
        double r21130 = pow(r21129, r21125);
        double r21131 = pow(r21128, r21125);
        double r21132 = r21130 - r21131;
        double r21133 = exp(r21132);
        double r21134 = log(r21133);
        double r21135 = 6.875959904884827e-20;
        bool r21136 = r21125 <= r21135;
        double r21137 = r21128 * r21124;
        double r21138 = r21123 / r21137;
        double r21139 = 1/2;
        double r21140 = r21139 / r21128;
        double r21141 = r21140 / r21137;
        double r21142 = r21138 - r21141;
        double r21143 = log(r21128);
        double r21144 = r21143 / r21124;
        double r21145 = r21144 / r21137;
        double r21146 = r21142 - r21145;
        double r21147 = r21128 + r21123;
        double r21148 = pow(r21147, r21125);
        double r21149 = r21148 * r21148;
        double r21150 = r21131 * r21131;
        double r21151 = r21149 - r21150;
        double r21152 = r21148 + r21131;
        double r21153 = r21151 / r21152;
        double r21154 = r21136 ? r21146 : r21153;
        double r21155 = r21127 ? r21134 : r21154;
        return r21155;
}

void mpfr_fmod2(mpfr_t r, mpfr_t n, mpfr_t d, mpfr_rnd_t rmd) {
        mpfr_fmod(r, n, d, rmd);
        if (mpfr_cmp_ui(r, 0) < 0) mpfr_add(r, r, d, rmd);
}


static mpfr_t r21156, r21157, r21158, r21159, r21160, r21161, r21162, r21163;

void setup_mpfr_f_im() {
        mpfr_set_default_prec(1424);
        mpfr_init(r21156);
        mpfr_init_set_str(r21157, "1", 10, MPFR_RNDN);
        mpfr_init(r21158);
        mpfr_init(r21159);
        mpfr_init(r21160);
        mpfr_init(r21161);
        mpfr_init(r21162);
        mpfr_init(r21163);
}

double f_im(double x, double n) {
        mpfr_set_d(r21156, x, MPFR_RNDN);
        ;
        mpfr_add(r21158, r21156, r21157, MPFR_RNDN);
        mpfr_set_d(r21159, n, MPFR_RNDN);
        mpfr_div(r21160, r21157, r21159, MPFR_RNDN);
        mpfr_pow(r21161, r21158, r21160, MPFR_RNDN);
        mpfr_pow(r21162, r21156, r21160, MPFR_RNDN);
        mpfr_sub(r21163, r21161, r21162, MPFR_RNDN);
        return mpfr_get_d(r21163, MPFR_RNDN);
}

static mpfr_t r21164, r21165, r21166, r21167, r21168, r21169, r21170, r21171, r21172, r21173, r21174, r21175, r21176, r21177, r21178, r21179, r21180, r21181, r21182, r21183, r21184, r21185, r21186, r21187, r21188, r21189, r21190, r21191, r21192, r21193, r21194, r21195, r21196;

void setup_mpfr_f_fm() {
        mpfr_set_default_prec(1424);
        mpfr_init_set_str(r21164, "1", 10, MPFR_RNDN);
        mpfr_init(r21165);
        mpfr_init(r21166);
        mpfr_init_set_str(r21167, "-3.6227047746916703e-09", 10, MPFR_RNDN);
        mpfr_init(r21168);
        mpfr_init(r21169);
        mpfr_init(r21170);
        mpfr_init(r21171);
        mpfr_init(r21172);
        mpfr_init(r21173);
        mpfr_init(r21174);
        mpfr_init(r21175);
        mpfr_init_set_str(r21176, "6.875959904884827e-20", 10, MPFR_RNDN);
        mpfr_init(r21177);
        mpfr_init(r21178);
        mpfr_init(r21179);
        mpfr_init_set_str(r21180, "1/2", 10, MPFR_RNDN);
        mpfr_init(r21181);
        mpfr_init(r21182);
        mpfr_init(r21183);
        mpfr_init(r21184);
        mpfr_init(r21185);
        mpfr_init(r21186);
        mpfr_init(r21187);
        mpfr_init(r21188);
        mpfr_init(r21189);
        mpfr_init(r21190);
        mpfr_init(r21191);
        mpfr_init(r21192);
        mpfr_init(r21193);
        mpfr_init(r21194);
        mpfr_init(r21195);
        mpfr_init(r21196);
}

double f_fm(double x, double n) {
        ;
        mpfr_set_d(r21165, n, MPFR_RNDN);
        mpfr_div(r21166, r21164, r21165, MPFR_RNDN);
        ;
        mpfr_set_si(r21168, mpfr_cmp(r21166, r21167) <= 0, MPFR_RNDN);
        mpfr_set_d(r21169, x, MPFR_RNDN);
        mpfr_add(r21170, r21164, r21169, MPFR_RNDN);
        mpfr_pow(r21171, r21170, r21166, MPFR_RNDN);
        mpfr_pow(r21172, r21169, r21166, MPFR_RNDN);
        mpfr_sub(r21173, r21171, r21172, MPFR_RNDN);
        mpfr_exp(r21174, r21173, MPFR_RNDN);
        mpfr_log(r21175, r21174, MPFR_RNDN);
        ;
        mpfr_set_si(r21177, mpfr_cmp(r21166, r21176) <= 0, MPFR_RNDN);
        mpfr_mul(r21178, r21169, r21165, MPFR_RNDN);
        mpfr_div(r21179, r21164, r21178, MPFR_RNDN);
        ;
        mpfr_div(r21181, r21180, r21169, MPFR_RNDN);
        mpfr_div(r21182, r21181, r21178, MPFR_RNDN);
        mpfr_sub(r21183, r21179, r21182, MPFR_RNDN);
        mpfr_log(r21184, r21169, MPFR_RNDN);
        mpfr_div(r21185, r21184, r21165, MPFR_RNDN);
        mpfr_div(r21186, r21185, r21178, MPFR_RNDN);
        mpfr_sub(r21187, r21183, r21186, MPFR_RNDN);
        mpfr_add(r21188, r21169, r21164, MPFR_RNDN);
        mpfr_pow(r21189, r21188, r21166, MPFR_RNDN);
        mpfr_mul(r21190, r21189, r21189, MPFR_RNDN);
        mpfr_mul(r21191, r21172, r21172, MPFR_RNDN);
        mpfr_sub(r21192, r21190, r21191, MPFR_RNDN);
        mpfr_add(r21193, r21189, r21172, MPFR_RNDN);
        mpfr_div(r21194, r21192, r21193, MPFR_RNDN);
        if (mpfr_get_si(r21177, MPFR_RNDN)) { mpfr_set(r21195, r21187, MPFR_RNDN); } else { mpfr_set(r21195, r21194, MPFR_RNDN); };
        if (mpfr_get_si(r21168, MPFR_RNDN)) { mpfr_set(r21196, r21175, MPFR_RNDN); } else { mpfr_set(r21196, r21195, MPFR_RNDN); };
        return mpfr_get_d(r21196, MPFR_RNDN);
}

static mpfr_t r21197, r21198, r21199, r21200, r21201, r21202, r21203, r21204, r21205, r21206, r21207, r21208, r21209, r21210, r21211, r21212, r21213, r21214, r21215, r21216, r21217, r21218, r21219, r21220, r21221, r21222, r21223, r21224, r21225, r21226, r21227, r21228, r21229;

void setup_mpfr_f_dm() {
        mpfr_set_default_prec(1424);
        mpfr_init_set_str(r21197, "1", 10, MPFR_RNDN);
        mpfr_init(r21198);
        mpfr_init(r21199);
        mpfr_init_set_str(r21200, "-3.6227047746916703e-09", 10, MPFR_RNDN);
        mpfr_init(r21201);
        mpfr_init(r21202);
        mpfr_init(r21203);
        mpfr_init(r21204);
        mpfr_init(r21205);
        mpfr_init(r21206);
        mpfr_init(r21207);
        mpfr_init(r21208);
        mpfr_init_set_str(r21209, "6.875959904884827e-20", 10, MPFR_RNDN);
        mpfr_init(r21210);
        mpfr_init(r21211);
        mpfr_init(r21212);
        mpfr_init_set_str(r21213, "1/2", 10, MPFR_RNDN);
        mpfr_init(r21214);
        mpfr_init(r21215);
        mpfr_init(r21216);
        mpfr_init(r21217);
        mpfr_init(r21218);
        mpfr_init(r21219);
        mpfr_init(r21220);
        mpfr_init(r21221);
        mpfr_init(r21222);
        mpfr_init(r21223);
        mpfr_init(r21224);
        mpfr_init(r21225);
        mpfr_init(r21226);
        mpfr_init(r21227);
        mpfr_init(r21228);
        mpfr_init(r21229);
}

double f_dm(double x, double n) {
        ;
        mpfr_set_d(r21198, n, MPFR_RNDN);
        mpfr_div(r21199, r21197, r21198, MPFR_RNDN);
        ;
        mpfr_set_si(r21201, mpfr_cmp(r21199, r21200) <= 0, MPFR_RNDN);
        mpfr_set_d(r21202, x, MPFR_RNDN);
        mpfr_add(r21203, r21197, r21202, MPFR_RNDN);
        mpfr_pow(r21204, r21203, r21199, MPFR_RNDN);
        mpfr_pow(r21205, r21202, r21199, MPFR_RNDN);
        mpfr_sub(r21206, r21204, r21205, MPFR_RNDN);
        mpfr_exp(r21207, r21206, MPFR_RNDN);
        mpfr_log(r21208, r21207, MPFR_RNDN);
        ;
        mpfr_set_si(r21210, mpfr_cmp(r21199, r21209) <= 0, MPFR_RNDN);
        mpfr_mul(r21211, r21202, r21198, MPFR_RNDN);
        mpfr_div(r21212, r21197, r21211, MPFR_RNDN);
        ;
        mpfr_div(r21214, r21213, r21202, MPFR_RNDN);
        mpfr_div(r21215, r21214, r21211, MPFR_RNDN);
        mpfr_sub(r21216, r21212, r21215, MPFR_RNDN);
        mpfr_log(r21217, r21202, MPFR_RNDN);
        mpfr_div(r21218, r21217, r21198, MPFR_RNDN);
        mpfr_div(r21219, r21218, r21211, MPFR_RNDN);
        mpfr_sub(r21220, r21216, r21219, MPFR_RNDN);
        mpfr_add(r21221, r21202, r21197, MPFR_RNDN);
        mpfr_pow(r21222, r21221, r21199, MPFR_RNDN);
        mpfr_mul(r21223, r21222, r21222, MPFR_RNDN);
        mpfr_mul(r21224, r21205, r21205, MPFR_RNDN);
        mpfr_sub(r21225, r21223, r21224, MPFR_RNDN);
        mpfr_add(r21226, r21222, r21205, MPFR_RNDN);
        mpfr_div(r21227, r21225, r21226, MPFR_RNDN);
        if (mpfr_get_si(r21210, MPFR_RNDN)) { mpfr_set(r21228, r21220, MPFR_RNDN); } else { mpfr_set(r21228, r21227, MPFR_RNDN); };
        if (mpfr_get_si(r21201, MPFR_RNDN)) { mpfr_set(r21229, r21208, MPFR_RNDN); } else { mpfr_set(r21229, r21228, MPFR_RNDN); };
        return mpfr_get_d(r21229, MPFR_RNDN);
}