#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 r20225 = x;
        float r20226 = 1;
        float r20227 = r20225 + r20226;
        float r20228 = n;
        float r20229 = r20226 / r20228;
        float r20230 = pow(r20227, r20229);
        float r20231 = pow(r20225, r20229);
        float r20232 = r20230 - r20231;
        return r20232;
}

double f_id(double x, double n) {
        double r20233 = x;
        double r20234 = 1;
        double r20235 = r20233 + r20234;
        double r20236 = n;
        double r20237 = r20234 / r20236;
        double r20238 = pow(r20235, r20237);
        double r20239 = pow(r20233, r20237);
        double r20240 = r20238 - r20239;
        return r20240;
}


double f_of(float x, float n) {
        float r20241 = n;
        float r20242 = -78372191.85757636;
        bool r20243 = r20241 <= r20242;
        float r20244 = 1;
        float r20245 = x;
        float r20246 = r20245 * r20241;
        float r20247 = r20244 / r20246;
        float r20248 = 1/2;
        float r20249 = r20248 / r20245;
        float r20250 = r20249 / r20246;
        float r20251 = r20247 - r20250;
        float r20252 = log(r20245);
        float r20253 = r20252 / r20241;
        float r20254 = r20253 / r20246;
        float r20255 = r20251 - r20254;
        float r20256 = 27520717438324964.0;
        bool r20257 = r20241 <= r20256;
        float r20258 = r20245 + r20244;
        float r20259 = sqrt(r20258);
        float r20260 = r20244 / r20241;
        float r20261 = pow(r20259, r20260);
        float r20262 = pow(r20245, r20260);
        float r20263 = sqrt(r20262);
        float r20264 = r20261 + r20263;
        float r20265 = r20244 + r20245;
        float r20266 = sqrt(r20265);
        float r20267 = pow(r20266, r20260);
        float r20268 = r20267 - r20263;
        float r20269 = exp(r20268);
        float r20270 = log(r20269);
        float r20271 = r20264 * r20270;
        float r20272 = r20257 ? r20271 : r20255;
        float r20273 = r20243 ? r20255 : r20272;
        return r20273;
}

double f_od(double x, double n) {
        double r20274 = n;
        double r20275 = -78372191.85757636;
        bool r20276 = r20274 <= r20275;
        double r20277 = 1;
        double r20278 = x;
        double r20279 = r20278 * r20274;
        double r20280 = r20277 / r20279;
        double r20281 = 1/2;
        double r20282 = r20281 / r20278;
        double r20283 = r20282 / r20279;
        double r20284 = r20280 - r20283;
        double r20285 = log(r20278);
        double r20286 = r20285 / r20274;
        double r20287 = r20286 / r20279;
        double r20288 = r20284 - r20287;
        double r20289 = 27520717438324964.0;
        bool r20290 = r20274 <= r20289;
        double r20291 = r20278 + r20277;
        double r20292 = sqrt(r20291);
        double r20293 = r20277 / r20274;
        double r20294 = pow(r20292, r20293);
        double r20295 = pow(r20278, r20293);
        double r20296 = sqrt(r20295);
        double r20297 = r20294 + r20296;
        double r20298 = r20277 + r20278;
        double r20299 = sqrt(r20298);
        double r20300 = pow(r20299, r20293);
        double r20301 = r20300 - r20296;
        double r20302 = exp(r20301);
        double r20303 = log(r20302);
        double r20304 = r20297 * r20303;
        double r20305 = r20290 ? r20304 : r20288;
        double r20306 = r20276 ? r20288 : r20305;
        return r20306;
}

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 r20307, r20308, r20309, r20310, r20311, r20312, r20313, r20314;

void setup_mpfr_f_im() {
        mpfr_set_default_prec(1424);
        mpfr_init(r20307);
        mpfr_init_set_str(r20308, "1", 10, MPFR_RNDN);
        mpfr_init(r20309);
        mpfr_init(r20310);
        mpfr_init(r20311);
        mpfr_init(r20312);
        mpfr_init(r20313);
        mpfr_init(r20314);
}

double f_im(double x, double n) {
        mpfr_set_d(r20307, x, MPFR_RNDN);
        ;
        mpfr_add(r20309, r20307, r20308, MPFR_RNDN);
        mpfr_set_d(r20310, n, MPFR_RNDN);
        mpfr_div(r20311, r20308, r20310, MPFR_RNDN);
        mpfr_pow(r20312, r20309, r20311, MPFR_RNDN);
        mpfr_pow(r20313, r20307, r20311, MPFR_RNDN);
        mpfr_sub(r20314, r20312, r20313, MPFR_RNDN);
        return mpfr_get_d(r20314, MPFR_RNDN);
}

static mpfr_t r20315, r20316, r20317, r20318, r20319, r20320, r20321, r20322, r20323, r20324, r20325, r20326, r20327, r20328, r20329, r20330, r20331, r20332, r20333, r20334, r20335, r20336, r20337, r20338, r20339, r20340, r20341, r20342, r20343, r20344, r20345, r20346, r20347;

void setup_mpfr_f_fm() {
        mpfr_set_default_prec(1424);
        mpfr_init(r20315);
        mpfr_init_set_str(r20316, "-78372191.85757636", 10, MPFR_RNDN);
        mpfr_init(r20317);
        mpfr_init_set_str(r20318, "1", 10, MPFR_RNDN);
        mpfr_init(r20319);
        mpfr_init(r20320);
        mpfr_init(r20321);
        mpfr_init_set_str(r20322, "1/2", 10, MPFR_RNDN);
        mpfr_init(r20323);
        mpfr_init(r20324);
        mpfr_init(r20325);
        mpfr_init(r20326);
        mpfr_init(r20327);
        mpfr_init(r20328);
        mpfr_init(r20329);
        mpfr_init_set_str(r20330, "27520717438324964.0", 10, MPFR_RNDN);
        mpfr_init(r20331);
        mpfr_init(r20332);
        mpfr_init(r20333);
        mpfr_init(r20334);
        mpfr_init(r20335);
        mpfr_init(r20336);
        mpfr_init(r20337);
        mpfr_init(r20338);
        mpfr_init(r20339);
        mpfr_init(r20340);
        mpfr_init(r20341);
        mpfr_init(r20342);
        mpfr_init(r20343);
        mpfr_init(r20344);
        mpfr_init(r20345);
        mpfr_init(r20346);
        mpfr_init(r20347);
}

double f_fm(double x, double n) {
        mpfr_set_d(r20315, n, MPFR_RNDN);
        ;
        mpfr_set_si(r20317, mpfr_cmp(r20315, r20316) <= 0, MPFR_RNDN);
        ;
        mpfr_set_d(r20319, x, MPFR_RNDN);
        mpfr_mul(r20320, r20319, r20315, MPFR_RNDN);
        mpfr_div(r20321, r20318, r20320, MPFR_RNDN);
        ;
        mpfr_div(r20323, r20322, r20319, MPFR_RNDN);
        mpfr_div(r20324, r20323, r20320, MPFR_RNDN);
        mpfr_sub(r20325, r20321, r20324, MPFR_RNDN);
        mpfr_log(r20326, r20319, MPFR_RNDN);
        mpfr_div(r20327, r20326, r20315, MPFR_RNDN);
        mpfr_div(r20328, r20327, r20320, MPFR_RNDN);
        mpfr_sub(r20329, r20325, r20328, MPFR_RNDN);
        ;
        mpfr_set_si(r20331, mpfr_cmp(r20315, r20330) <= 0, MPFR_RNDN);
        mpfr_add(r20332, r20319, r20318, MPFR_RNDN);
        mpfr_sqrt(r20333, r20332, MPFR_RNDN);
        mpfr_div(r20334, r20318, r20315, MPFR_RNDN);
        mpfr_pow(r20335, r20333, r20334, MPFR_RNDN);
        mpfr_pow(r20336, r20319, r20334, MPFR_RNDN);
        mpfr_sqrt(r20337, r20336, MPFR_RNDN);
        mpfr_add(r20338, r20335, r20337, MPFR_RNDN);
        mpfr_add(r20339, r20318, r20319, MPFR_RNDN);
        mpfr_sqrt(r20340, r20339, MPFR_RNDN);
        mpfr_pow(r20341, r20340, r20334, MPFR_RNDN);
        mpfr_sub(r20342, r20341, r20337, MPFR_RNDN);
        mpfr_exp(r20343, r20342, MPFR_RNDN);
        mpfr_log(r20344, r20343, MPFR_RNDN);
        mpfr_mul(r20345, r20338, r20344, MPFR_RNDN);
        if (mpfr_get_si(r20331, MPFR_RNDN)) { mpfr_set(r20346, r20345, MPFR_RNDN); } else { mpfr_set(r20346, r20329, MPFR_RNDN); };
        if (mpfr_get_si(r20317, MPFR_RNDN)) { mpfr_set(r20347, r20329, MPFR_RNDN); } else { mpfr_set(r20347, r20346, MPFR_RNDN); };
        return mpfr_get_d(r20347, MPFR_RNDN);
}

static mpfr_t r20348, r20349, r20350, r20351, r20352, r20353, r20354, r20355, r20356, r20357, r20358, r20359, r20360, r20361, r20362, r20363, r20364, r20365, r20366, r20367, r20368, r20369, r20370, r20371, r20372, r20373, r20374, r20375, r20376, r20377, r20378, r20379, r20380;

void setup_mpfr_f_dm() {
        mpfr_set_default_prec(1424);
        mpfr_init(r20348);
        mpfr_init_set_str(r20349, "-78372191.85757636", 10, MPFR_RNDN);
        mpfr_init(r20350);
        mpfr_init_set_str(r20351, "1", 10, MPFR_RNDN);
        mpfr_init(r20352);
        mpfr_init(r20353);
        mpfr_init(r20354);
        mpfr_init_set_str(r20355, "1/2", 10, MPFR_RNDN);
        mpfr_init(r20356);
        mpfr_init(r20357);
        mpfr_init(r20358);
        mpfr_init(r20359);
        mpfr_init(r20360);
        mpfr_init(r20361);
        mpfr_init(r20362);
        mpfr_init_set_str(r20363, "27520717438324964.0", 10, MPFR_RNDN);
        mpfr_init(r20364);
        mpfr_init(r20365);
        mpfr_init(r20366);
        mpfr_init(r20367);
        mpfr_init(r20368);
        mpfr_init(r20369);
        mpfr_init(r20370);
        mpfr_init(r20371);
        mpfr_init(r20372);
        mpfr_init(r20373);
        mpfr_init(r20374);
        mpfr_init(r20375);
        mpfr_init(r20376);
        mpfr_init(r20377);
        mpfr_init(r20378);
        mpfr_init(r20379);
        mpfr_init(r20380);
}

double f_dm(double x, double n) {
        mpfr_set_d(r20348, n, MPFR_RNDN);
        ;
        mpfr_set_si(r20350, mpfr_cmp(r20348, r20349) <= 0, MPFR_RNDN);
        ;
        mpfr_set_d(r20352, x, MPFR_RNDN);
        mpfr_mul(r20353, r20352, r20348, MPFR_RNDN);
        mpfr_div(r20354, r20351, r20353, MPFR_RNDN);
        ;
        mpfr_div(r20356, r20355, r20352, MPFR_RNDN);
        mpfr_div(r20357, r20356, r20353, MPFR_RNDN);
        mpfr_sub(r20358, r20354, r20357, MPFR_RNDN);
        mpfr_log(r20359, r20352, MPFR_RNDN);
        mpfr_div(r20360, r20359, r20348, MPFR_RNDN);
        mpfr_div(r20361, r20360, r20353, MPFR_RNDN);
        mpfr_sub(r20362, r20358, r20361, MPFR_RNDN);
        ;
        mpfr_set_si(r20364, mpfr_cmp(r20348, r20363) <= 0, MPFR_RNDN);
        mpfr_add(r20365, r20352, r20351, MPFR_RNDN);
        mpfr_sqrt(r20366, r20365, MPFR_RNDN);
        mpfr_div(r20367, r20351, r20348, MPFR_RNDN);
        mpfr_pow(r20368, r20366, r20367, MPFR_RNDN);
        mpfr_pow(r20369, r20352, r20367, MPFR_RNDN);
        mpfr_sqrt(r20370, r20369, MPFR_RNDN);
        mpfr_add(r20371, r20368, r20370, MPFR_RNDN);
        mpfr_add(r20372, r20351, r20352, MPFR_RNDN);
        mpfr_sqrt(r20373, r20372, MPFR_RNDN);
        mpfr_pow(r20374, r20373, r20367, MPFR_RNDN);
        mpfr_sub(r20375, r20374, r20370, MPFR_RNDN);
        mpfr_exp(r20376, r20375, MPFR_RNDN);
        mpfr_log(r20377, r20376, MPFR_RNDN);
        mpfr_mul(r20378, r20371, r20377, MPFR_RNDN);
        if (mpfr_get_si(r20364, MPFR_RNDN)) { mpfr_set(r20379, r20378, MPFR_RNDN); } else { mpfr_set(r20379, r20362, MPFR_RNDN); };
        if (mpfr_get_si(r20350, MPFR_RNDN)) { mpfr_set(r20380, r20362, MPFR_RNDN); } else { mpfr_set(r20380, r20379, MPFR_RNDN); };
        return mpfr_get_d(r20380, MPFR_RNDN);
}