#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 r20232 = x;
        float r20233 = 1;
        float r20234 = r20232 + r20233;
        float r20235 = n;
        float r20236 = r20233 / r20235;
        float r20237 = pow(r20234, r20236);
        float r20238 = pow(r20232, r20236);
        float r20239 = r20237 - r20238;
        return r20239;
}

double f_id(double x, double n) {
        double r20240 = x;
        double r20241 = 1;
        double r20242 = r20240 + r20241;
        double r20243 = n;
        double r20244 = r20241 / r20243;
        double r20245 = pow(r20242, r20244);
        double r20246 = pow(r20240, r20244);
        double r20247 = r20245 - r20246;
        return r20247;
}


double f_of(float x, float n) {
        float r20248 = n;
        float r20249 = -5450781.891258214;
        bool r20250 = r20248 <= r20249;
        float r20251 = 1;
        float r20252 = x;
        float r20253 = r20252 * r20248;
        float r20254 = r20251 / r20253;
        float r20255 = 1/2;
        float r20256 = r20255 / r20252;
        float r20257 = r20256 / r20253;
        float r20258 = r20254 - r20257;
        float r20259 = log(r20252);
        float r20260 = r20259 / r20248;
        float r20261 = r20260 / r20253;
        float r20262 = r20258 - r20261;
        float r20263 = 5.83627980787353e+16;
        bool r20264 = r20248 <= r20263;
        float r20265 = r20252 + r20251;
        float r20266 = r20251 / r20248;
        float r20267 = pow(r20265, r20266);
        float r20268 = pow(r20252, r20266);
        float r20269 = r20267 - r20268;
        float r20270 = cbrt(r20269);
        float r20271 = r20270 * r20270;
        float r20272 = 3;
        float r20273 = pow(r20268, r20272);
        float r20274 = cbrt(r20273);
        float r20275 = r20267 - r20274;
        float r20276 = cbrt(r20275);
        float r20277 = r20271 * r20276;
        float r20278 = r20264 ? r20277 : r20262;
        float r20279 = r20250 ? r20262 : r20278;
        return r20279;
}

double f_od(double x, double n) {
        double r20280 = n;
        double r20281 = -5450781.891258214;
        bool r20282 = r20280 <= r20281;
        double r20283 = 1;
        double r20284 = x;
        double r20285 = r20284 * r20280;
        double r20286 = r20283 / r20285;
        double r20287 = 1/2;
        double r20288 = r20287 / r20284;
        double r20289 = r20288 / r20285;
        double r20290 = r20286 - r20289;
        double r20291 = log(r20284);
        double r20292 = r20291 / r20280;
        double r20293 = r20292 / r20285;
        double r20294 = r20290 - r20293;
        double r20295 = 5.83627980787353e+16;
        bool r20296 = r20280 <= r20295;
        double r20297 = r20284 + r20283;
        double r20298 = r20283 / r20280;
        double r20299 = pow(r20297, r20298);
        double r20300 = pow(r20284, r20298);
        double r20301 = r20299 - r20300;
        double r20302 = cbrt(r20301);
        double r20303 = r20302 * r20302;
        double r20304 = 3;
        double r20305 = pow(r20300, r20304);
        double r20306 = cbrt(r20305);
        double r20307 = r20299 - r20306;
        double r20308 = cbrt(r20307);
        double r20309 = r20303 * r20308;
        double r20310 = r20296 ? r20309 : r20294;
        double r20311 = r20282 ? r20294 : r20310;
        return r20311;
}

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 r20312, r20313, r20314, r20315, r20316, r20317, r20318, r20319;

void setup_mpfr_f_im() {
        mpfr_set_default_prec(1424);
        mpfr_init(r20312);
        mpfr_init_set_str(r20313, "1", 10, MPFR_RNDN);
        mpfr_init(r20314);
        mpfr_init(r20315);
        mpfr_init(r20316);
        mpfr_init(r20317);
        mpfr_init(r20318);
        mpfr_init(r20319);
}

double f_im(double x, double n) {
        mpfr_set_d(r20312, x, MPFR_RNDN);
        ;
        mpfr_add(r20314, r20312, r20313, MPFR_RNDN);
        mpfr_set_d(r20315, n, MPFR_RNDN);
        mpfr_div(r20316, r20313, r20315, MPFR_RNDN);
        mpfr_pow(r20317, r20314, r20316, MPFR_RNDN);
        mpfr_pow(r20318, r20312, r20316, MPFR_RNDN);
        mpfr_sub(r20319, r20317, r20318, MPFR_RNDN);
        return mpfr_get_d(r20319, MPFR_RNDN);
}

static mpfr_t 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, r20348, r20349, r20350, r20351;

void setup_mpfr_f_fm() {
        mpfr_set_default_prec(1424);
        mpfr_init(r20320);
        mpfr_init_set_str(r20321, "-5450781.891258214", 10, MPFR_RNDN);
        mpfr_init(r20322);
        mpfr_init_set_str(r20323, "1", 10, MPFR_RNDN);
        mpfr_init(r20324);
        mpfr_init(r20325);
        mpfr_init(r20326);
        mpfr_init_set_str(r20327, "1/2", 10, MPFR_RNDN);
        mpfr_init(r20328);
        mpfr_init(r20329);
        mpfr_init(r20330);
        mpfr_init(r20331);
        mpfr_init(r20332);
        mpfr_init(r20333);
        mpfr_init(r20334);
        mpfr_init_set_str(r20335, "5.83627980787353e+16", 10, MPFR_RNDN);
        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_set_str(r20344, "3", 10, MPFR_RNDN);
        mpfr_init(r20345);
        mpfr_init(r20346);
        mpfr_init(r20347);
        mpfr_init(r20348);
        mpfr_init(r20349);
        mpfr_init(r20350);
        mpfr_init(r20351);
}

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

static mpfr_t 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, r20381, r20382, r20383;

void setup_mpfr_f_dm() {
        mpfr_set_default_prec(1424);
        mpfr_init(r20352);
        mpfr_init_set_str(r20353, "-5450781.891258214", 10, MPFR_RNDN);
        mpfr_init(r20354);
        mpfr_init_set_str(r20355, "1", 10, MPFR_RNDN);
        mpfr_init(r20356);
        mpfr_init(r20357);
        mpfr_init(r20358);
        mpfr_init_set_str(r20359, "1/2", 10, MPFR_RNDN);
        mpfr_init(r20360);
        mpfr_init(r20361);
        mpfr_init(r20362);
        mpfr_init(r20363);
        mpfr_init(r20364);
        mpfr_init(r20365);
        mpfr_init(r20366);
        mpfr_init_set_str(r20367, "5.83627980787353e+16", 10, MPFR_RNDN);
        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_set_str(r20376, "3", 10, MPFR_RNDN);
        mpfr_init(r20377);
        mpfr_init(r20378);
        mpfr_init(r20379);
        mpfr_init(r20380);
        mpfr_init(r20381);
        mpfr_init(r20382);
        mpfr_init(r20383);
}

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