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

char *name = "From Warwick Tucker's Validated Numerics";

double f_if() {
        float r8824 = 333.75;
        float r8825 = 33096;
        float r8826 = 6;
        float r8827 = pow(r8825, r8826);
        float r8828 = r8824 * r8827;
        float r8829 = 77617;
        float r8830 = r8829 * r8829;
        float r8831 = 11;
        float r8832 = r8831 * r8830;
        float r8833 = r8825 * r8825;
        float r8834 = r8832 * r8833;
        float r8835 = -r8827;
        float r8836 = r8834 + r8835;
        float r8837 = -121;
        float r8838 = 4;
        float r8839 = pow(r8825, r8838);
        float r8840 = r8837 * r8839;
        float r8841 = r8836 + r8840;
        float r8842 = -2;
        float r8843 = r8841 + r8842;
        float r8844 = r8830 * r8843;
        float r8845 = r8828 + r8844;
        float r8846 = 5.5;
        float r8847 = 8;
        float r8848 = pow(r8825, r8847);
        float r8849 = r8846 * r8848;
        float r8850 = r8845 + r8849;
        float r8851 = 2;
        float r8852 = r8851 * r8825;
        float r8853 = r8829 / r8852;
        float r8854 = r8850 + r8853;
        return r8854;
}

double f_id() {
        double r8855 = 333.75;
        double r8856 = 33096;
        double r8857 = 6;
        double r8858 = pow(r8856, r8857);
        double r8859 = r8855 * r8858;
        double r8860 = 77617;
        double r8861 = r8860 * r8860;
        double r8862 = 11;
        double r8863 = r8862 * r8861;
        double r8864 = r8856 * r8856;
        double r8865 = r8863 * r8864;
        double r8866 = -r8858;
        double r8867 = r8865 + r8866;
        double r8868 = -121;
        double r8869 = 4;
        double r8870 = pow(r8856, r8869);
        double r8871 = r8868 * r8870;
        double r8872 = r8867 + r8871;
        double r8873 = -2;
        double r8874 = r8872 + r8873;
        double r8875 = r8861 * r8874;
        double r8876 = r8859 + r8875;
        double r8877 = 5.5;
        double r8878 = 8;
        double r8879 = pow(r8856, r8878);
        double r8880 = r8877 * r8879;
        double r8881 = r8876 + r8880;
        double r8882 = 2;
        double r8883 = r8882 * r8856;
        double r8884 = r8860 / r8883;
        double r8885 = r8881 + r8884;
        return r8885;
}


double f_of() {
        float r8886 = 5.5;
        float r8887 = 33096;
        float r8888 = 8;
        float r8889 = pow(r8887, r8888);
        float r8890 = r8886 * r8889;
        float r8891 = 3;
        float r8892 = pow(r8890, r8891);
        float r8893 = 77617;
        float r8894 = r8893 * r8893;
        float r8895 = 11;
        float r8896 = r8887 * r8895;
        float r8897 = r8896 * r8894;
        float r8898 = r8887 * r8897;
        float r8899 = 6;
        float r8900 = pow(r8887, r8899);
        float r8901 = -2;
        float r8902 = r8900 - r8901;
        float r8903 = -121;
        float r8904 = 4;
        float r8905 = pow(r8887, r8904);
        float r8906 = r8903 * r8905;
        float r8907 = r8902 - r8906;
        float r8908 = r8898 - r8907;
        float r8909 = r8894 * r8908;
        float r8910 = 333.75;
        float r8911 = r8900 * r8910;
        float r8912 = r8909 + r8911;
        float r8913 = r8912 * r8912;
        float r8914 = r8913 * r8912;
        float r8915 = r8892 + r8914;
        float r8916 = r8889 * r8886;
        float r8917 = r8916 * r8916;
        float r8918 = r8887 * r8896;
        float r8919 = r8918 * r8894;
        float r8920 = r8919 - r8907;
        float r8921 = r8894 * r8920;
        float r8922 = r8921 + r8911;
        float r8923 = pow(r8916, r8891);
        float r8924 = cbrt(r8923);
        float r8925 = r8924 - r8911;
        float r8926 = r8925 - r8921;
        float r8927 = r8922 * r8926;
        float r8928 = r8917 - r8927;
        float r8929 = r8915 / r8928;
        float r8930 = 2;
        float r8931 = r8930 * r8887;
        float r8932 = r8893 / r8931;
        float r8933 = r8929 + r8932;
        return r8933;
}

double f_od() {
        double r8934 = 5.5;
        double r8935 = 33096;
        double r8936 = 8;
        double r8937 = pow(r8935, r8936);
        double r8938 = r8934 * r8937;
        double r8939 = 3;
        double r8940 = pow(r8938, r8939);
        double r8941 = 77617;
        double r8942 = r8941 * r8941;
        double r8943 = 11;
        double r8944 = r8935 * r8943;
        double r8945 = r8944 * r8942;
        double r8946 = r8935 * r8945;
        double r8947 = 6;
        double r8948 = pow(r8935, r8947);
        double r8949 = -2;
        double r8950 = r8948 - r8949;
        double r8951 = -121;
        double r8952 = 4;
        double r8953 = pow(r8935, r8952);
        double r8954 = r8951 * r8953;
        double r8955 = r8950 - r8954;
        double r8956 = r8946 - r8955;
        double r8957 = r8942 * r8956;
        double r8958 = 333.75;
        double r8959 = r8948 * r8958;
        double r8960 = r8957 + r8959;
        double r8961 = r8960 * r8960;
        double r8962 = r8961 * r8960;
        double r8963 = r8940 + r8962;
        double r8964 = r8937 * r8934;
        double r8965 = r8964 * r8964;
        double r8966 = r8935 * r8944;
        double r8967 = r8966 * r8942;
        double r8968 = r8967 - r8955;
        double r8969 = r8942 * r8968;
        double r8970 = r8969 + r8959;
        double r8971 = pow(r8964, r8939);
        double r8972 = cbrt(r8971);
        double r8973 = r8972 - r8959;
        double r8974 = r8973 - r8969;
        double r8975 = r8970 * r8974;
        double r8976 = r8965 - r8975;
        double r8977 = r8963 / r8976;
        double r8978 = 2;
        double r8979 = r8978 * r8935;
        double r8980 = r8941 / r8979;
        double r8981 = r8977 + r8980;
        return r8981;
}

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 r8982, r8983, r8984, r8985, r8986, r8987, r8988, r8989, r8990, r8991, r8992, r8993, r8994, r8995, r8996, r8997, r8998, r8999, r9000, r9001, r9002, r9003, r9004, r9005, r9006, r9007, r9008, r9009, r9010, r9011, r9012;

void setup_mpfr_f_im() {
        mpfr_set_default_prec(592);
        mpfr_init_set_str(r8982, "333.75", 10, MPFR_RNDN);
        mpfr_init_set_str(r8983, "33096", 10, MPFR_RNDN);
        mpfr_init_set_str(r8984, "6", 10, MPFR_RNDN);
        mpfr_init(r8985);
        mpfr_init(r8986);
        mpfr_init_set_str(r8987, "77617", 10, MPFR_RNDN);
        mpfr_init(r8988);
        mpfr_init_set_str(r8989, "11", 10, MPFR_RNDN);
        mpfr_init(r8990);
        mpfr_init(r8991);
        mpfr_init(r8992);
        mpfr_init(r8993);
        mpfr_init(r8994);
        mpfr_init_set_str(r8995, "-121", 10, MPFR_RNDN);
        mpfr_init_set_str(r8996, "4", 10, MPFR_RNDN);
        mpfr_init(r8997);
        mpfr_init(r8998);
        mpfr_init(r8999);
        mpfr_init_set_str(r9000, "-2", 10, MPFR_RNDN);
        mpfr_init(r9001);
        mpfr_init(r9002);
        mpfr_init(r9003);
        mpfr_init_set_str(r9004, "5.5", 10, MPFR_RNDN);
        mpfr_init_set_str(r9005, "8", 10, MPFR_RNDN);
        mpfr_init(r9006);
        mpfr_init(r9007);
        mpfr_init(r9008);
        mpfr_init_set_str(r9009, "2", 10, MPFR_RNDN);
        mpfr_init(r9010);
        mpfr_init(r9011);
        mpfr_init(r9012);
}

double f_im() {
        ;
        ;
        ;
        mpfr_pow(r8985, r8983, r8984, MPFR_RNDN);
        mpfr_mul(r8986, r8982, r8985, MPFR_RNDN);
        ;
        mpfr_mul(r8988, r8987, r8987, MPFR_RNDN);
        ;
        mpfr_mul(r8990, r8989, r8988, MPFR_RNDN);
        mpfr_mul(r8991, r8983, r8983, MPFR_RNDN);
        mpfr_mul(r8992, r8990, r8991, MPFR_RNDN);
        mpfr_neg(r8993, r8985, MPFR_RNDN);
        mpfr_add(r8994, r8992, r8993, MPFR_RNDN);
        ;
        ;
        mpfr_pow(r8997, r8983, r8996, MPFR_RNDN);
        mpfr_mul(r8998, r8995, r8997, MPFR_RNDN);
        mpfr_add(r8999, r8994, r8998, MPFR_RNDN);
        ;
        mpfr_add(r9001, r8999, r9000, MPFR_RNDN);
        mpfr_mul(r9002, r8988, r9001, MPFR_RNDN);
        mpfr_add(r9003, r8986, r9002, MPFR_RNDN);
        ;
        ;
        mpfr_pow(r9006, r8983, r9005, MPFR_RNDN);
        mpfr_mul(r9007, r9004, r9006, MPFR_RNDN);
        mpfr_add(r9008, r9003, r9007, MPFR_RNDN);
        ;
        mpfr_mul(r9010, r9009, r8983, MPFR_RNDN);
        mpfr_div(r9011, r8987, r9010, MPFR_RNDN);
        mpfr_add(r9012, r9008, r9011, MPFR_RNDN);
        return mpfr_get_d(r9012, MPFR_RNDN);
}

static mpfr_t r9013, r9014, r9015, r9016, r9017, r9018, r9019, r9020, r9021, r9022, r9023, r9024, r9025, r9026, r9027, r9028, r9029, r9030, r9031, r9032, r9033, r9034, r9035, r9036, r9037, r9038, r9039, r9040, r9041, r9042, r9043, r9044, r9045, r9046, r9047, r9048, r9049, r9050, r9051, r9052, r9053, r9054, r9055, r9056, r9057, r9058, r9059, r9060;

void setup_mpfr_f_fm() {
        mpfr_set_default_prec(592);
        mpfr_init_set_str(r9013, "5.5", 10, MPFR_RNDN);
        mpfr_init_set_str(r9014, "33096", 10, MPFR_RNDN);
        mpfr_init_set_str(r9015, "8", 10, MPFR_RNDN);
        mpfr_init(r9016);
        mpfr_init(r9017);
        mpfr_init_set_str(r9018, "3", 10, MPFR_RNDN);
        mpfr_init(r9019);
        mpfr_init_set_str(r9020, "77617", 10, MPFR_RNDN);
        mpfr_init(r9021);
        mpfr_init_set_str(r9022, "11", 10, MPFR_RNDN);
        mpfr_init(r9023);
        mpfr_init(r9024);
        mpfr_init(r9025);
        mpfr_init_set_str(r9026, "6", 10, MPFR_RNDN);
        mpfr_init(r9027);
        mpfr_init_set_str(r9028, "-2", 10, MPFR_RNDN);
        mpfr_init(r9029);
        mpfr_init_set_str(r9030, "-121", 10, MPFR_RNDN);
        mpfr_init_set_str(r9031, "4", 10, MPFR_RNDN);
        mpfr_init(r9032);
        mpfr_init(r9033);
        mpfr_init(r9034);
        mpfr_init(r9035);
        mpfr_init(r9036);
        mpfr_init_set_str(r9037, "333.75", 10, MPFR_RNDN);
        mpfr_init(r9038);
        mpfr_init(r9039);
        mpfr_init(r9040);
        mpfr_init(r9041);
        mpfr_init(r9042);
        mpfr_init(r9043);
        mpfr_init(r9044);
        mpfr_init(r9045);
        mpfr_init(r9046);
        mpfr_init(r9047);
        mpfr_init(r9048);
        mpfr_init(r9049);
        mpfr_init(r9050);
        mpfr_init(r9051);
        mpfr_init(r9052);
        mpfr_init(r9053);
        mpfr_init(r9054);
        mpfr_init(r9055);
        mpfr_init(r9056);
        mpfr_init_set_str(r9057, "2", 10, MPFR_RNDN);
        mpfr_init(r9058);
        mpfr_init(r9059);
        mpfr_init(r9060);
}

double f_fm() {
        ;
        ;
        ;
        mpfr_pow(r9016, r9014, r9015, MPFR_RNDN);
        mpfr_mul(r9017, r9013, r9016, MPFR_RNDN);
        ;
        mpfr_pow(r9019, r9017, r9018, MPFR_RNDN);
        ;
        mpfr_mul(r9021, r9020, r9020, MPFR_RNDN);
        ;
        mpfr_mul(r9023, r9014, r9022, MPFR_RNDN);
        mpfr_mul(r9024, r9023, r9021, MPFR_RNDN);
        mpfr_mul(r9025, r9014, r9024, MPFR_RNDN);
        ;
        mpfr_pow(r9027, r9014, r9026, MPFR_RNDN);
        ;
        mpfr_sub(r9029, r9027, r9028, MPFR_RNDN);
        ;
        ;
        mpfr_pow(r9032, r9014, r9031, MPFR_RNDN);
        mpfr_mul(r9033, r9030, r9032, MPFR_RNDN);
        mpfr_sub(r9034, r9029, r9033, MPFR_RNDN);
        mpfr_sub(r9035, r9025, r9034, MPFR_RNDN);
        mpfr_mul(r9036, r9021, r9035, MPFR_RNDN);
        ;
        mpfr_mul(r9038, r9027, r9037, MPFR_RNDN);
        mpfr_add(r9039, r9036, r9038, MPFR_RNDN);
        mpfr_mul(r9040, r9039, r9039, MPFR_RNDN);
        mpfr_mul(r9041, r9040, r9039, MPFR_RNDN);
        mpfr_add(r9042, r9019, r9041, MPFR_RNDN);
        mpfr_mul(r9043, r9016, r9013, MPFR_RNDN);
        mpfr_mul(r9044, r9043, r9043, MPFR_RNDN);
        mpfr_mul(r9045, r9014, r9023, MPFR_RNDN);
        mpfr_mul(r9046, r9045, r9021, MPFR_RNDN);
        mpfr_sub(r9047, r9046, r9034, MPFR_RNDN);
        mpfr_mul(r9048, r9021, r9047, MPFR_RNDN);
        mpfr_add(r9049, r9048, r9038, MPFR_RNDN);
        mpfr_pow(r9050, r9043, r9018, MPFR_RNDN);
        mpfr_cbrt(r9051, r9050, MPFR_RNDN);
        mpfr_sub(r9052, r9051, r9038, MPFR_RNDN);
        mpfr_sub(r9053, r9052, r9048, MPFR_RNDN);
        mpfr_mul(r9054, r9049, r9053, MPFR_RNDN);
        mpfr_sub(r9055, r9044, r9054, MPFR_RNDN);
        mpfr_div(r9056, r9042, r9055, MPFR_RNDN);
        ;
        mpfr_mul(r9058, r9057, r9014, MPFR_RNDN);
        mpfr_div(r9059, r9020, r9058, MPFR_RNDN);
        mpfr_add(r9060, r9056, r9059, MPFR_RNDN);
        return mpfr_get_d(r9060, MPFR_RNDN);
}

static mpfr_t r9061, r9062, r9063, r9064, r9065, r9066, r9067, r9068, r9069, r9070, r9071, r9072, r9073, r9074, r9075, r9076, r9077, r9078, r9079, r9080, r9081, r9082, r9083, r9084, r9085, r9086, r9087, r9088, r9089, r9090, r9091, r9092, r9093, r9094, r9095, r9096, r9097, r9098, r9099, r9100, r9101, r9102, r9103, r9104, r9105, r9106, r9107, r9108;

void setup_mpfr_f_dm() {
        mpfr_set_default_prec(592);
        mpfr_init_set_str(r9061, "5.5", 10, MPFR_RNDN);
        mpfr_init_set_str(r9062, "33096", 10, MPFR_RNDN);
        mpfr_init_set_str(r9063, "8", 10, MPFR_RNDN);
        mpfr_init(r9064);
        mpfr_init(r9065);
        mpfr_init_set_str(r9066, "3", 10, MPFR_RNDN);
        mpfr_init(r9067);
        mpfr_init_set_str(r9068, "77617", 10, MPFR_RNDN);
        mpfr_init(r9069);
        mpfr_init_set_str(r9070, "11", 10, MPFR_RNDN);
        mpfr_init(r9071);
        mpfr_init(r9072);
        mpfr_init(r9073);
        mpfr_init_set_str(r9074, "6", 10, MPFR_RNDN);
        mpfr_init(r9075);
        mpfr_init_set_str(r9076, "-2", 10, MPFR_RNDN);
        mpfr_init(r9077);
        mpfr_init_set_str(r9078, "-121", 10, MPFR_RNDN);
        mpfr_init_set_str(r9079, "4", 10, MPFR_RNDN);
        mpfr_init(r9080);
        mpfr_init(r9081);
        mpfr_init(r9082);
        mpfr_init(r9083);
        mpfr_init(r9084);
        mpfr_init_set_str(r9085, "333.75", 10, MPFR_RNDN);
        mpfr_init(r9086);
        mpfr_init(r9087);
        mpfr_init(r9088);
        mpfr_init(r9089);
        mpfr_init(r9090);
        mpfr_init(r9091);
        mpfr_init(r9092);
        mpfr_init(r9093);
        mpfr_init(r9094);
        mpfr_init(r9095);
        mpfr_init(r9096);
        mpfr_init(r9097);
        mpfr_init(r9098);
        mpfr_init(r9099);
        mpfr_init(r9100);
        mpfr_init(r9101);
        mpfr_init(r9102);
        mpfr_init(r9103);
        mpfr_init(r9104);
        mpfr_init_set_str(r9105, "2", 10, MPFR_RNDN);
        mpfr_init(r9106);
        mpfr_init(r9107);
        mpfr_init(r9108);
}

double f_dm() {
        ;
        ;
        ;
        mpfr_pow(r9064, r9062, r9063, MPFR_RNDN);
        mpfr_mul(r9065, r9061, r9064, MPFR_RNDN);
        ;
        mpfr_pow(r9067, r9065, r9066, MPFR_RNDN);
        ;
        mpfr_mul(r9069, r9068, r9068, MPFR_RNDN);
        ;
        mpfr_mul(r9071, r9062, r9070, MPFR_RNDN);
        mpfr_mul(r9072, r9071, r9069, MPFR_RNDN);
        mpfr_mul(r9073, r9062, r9072, MPFR_RNDN);
        ;
        mpfr_pow(r9075, r9062, r9074, MPFR_RNDN);
        ;
        mpfr_sub(r9077, r9075, r9076, MPFR_RNDN);
        ;
        ;
        mpfr_pow(r9080, r9062, r9079, MPFR_RNDN);
        mpfr_mul(r9081, r9078, r9080, MPFR_RNDN);
        mpfr_sub(r9082, r9077, r9081, MPFR_RNDN);
        mpfr_sub(r9083, r9073, r9082, MPFR_RNDN);
        mpfr_mul(r9084, r9069, r9083, MPFR_RNDN);
        ;
        mpfr_mul(r9086, r9075, r9085, MPFR_RNDN);
        mpfr_add(r9087, r9084, r9086, MPFR_RNDN);
        mpfr_mul(r9088, r9087, r9087, MPFR_RNDN);
        mpfr_mul(r9089, r9088, r9087, MPFR_RNDN);
        mpfr_add(r9090, r9067, r9089, MPFR_RNDN);
        mpfr_mul(r9091, r9064, r9061, MPFR_RNDN);
        mpfr_mul(r9092, r9091, r9091, MPFR_RNDN);
        mpfr_mul(r9093, r9062, r9071, MPFR_RNDN);
        mpfr_mul(r9094, r9093, r9069, MPFR_RNDN);
        mpfr_sub(r9095, r9094, r9082, MPFR_RNDN);
        mpfr_mul(r9096, r9069, r9095, MPFR_RNDN);
        mpfr_add(r9097, r9096, r9086, MPFR_RNDN);
        mpfr_pow(r9098, r9091, r9066, MPFR_RNDN);
        mpfr_cbrt(r9099, r9098, MPFR_RNDN);
        mpfr_sub(r9100, r9099, r9086, MPFR_RNDN);
        mpfr_sub(r9101, r9100, r9096, MPFR_RNDN);
        mpfr_mul(r9102, r9097, r9101, MPFR_RNDN);
        mpfr_sub(r9103, r9092, r9102, MPFR_RNDN);
        mpfr_div(r9104, r9090, r9103, MPFR_RNDN);
        ;
        mpfr_mul(r9106, r9105, r9062, MPFR_RNDN);
        mpfr_div(r9107, r9068, r9106, MPFR_RNDN);
        mpfr_add(r9108, r9104, r9107, MPFR_RNDN);
        return mpfr_get_d(r9108, MPFR_RNDN);
}