Problem 180

Sum of distinct s=x+y+z for golden triples of order 35; answer = num+den of total.

Answer285196020571078987
Output285196020571078987
StatusPASS
Native helperno
Runtime20 ms
Peak memory2704 KB
Time complexityO(n^3) (estimated)
Space complexityO(n^2) (estimated)

Performance comparison

MetricOur solutionBest known
Time complexityO(n^3)O(n^2)
Space complexityO(n^2)O(n)
ApproachFlow solutionTriple enumeration
VerdictSuboptimal

Flow source

# Project Euler 180
# Sum of distinct s=x+y+z for golden triples of order 35; answer = num+den of total.

import euler.nt { isqrt }

extern {
    function calloc(n: i64, size: i64) -> ptr<void>
    function free(p: ptr<void>) -> void
}

function gcd_abs(a: i64, b: i64) -> i64 {
    let mut x: i64 = a
    if x < 0 { x = 0 - x }
    let mut y: i64 = b
    if y < 0 { y = 0 - y }
    while y != 0 {
        let t: i64 = x % y
        x = y
        y = t
    }
    return x
}

# Store fractions as num/den in arrays; hash set of sums as reduced fractions
# Max rationals for k=35: sum_{b=2}^{35} phi(b) ≈ 385

function main() -> i32 {
    let K: i64 = 35
    let nums: ptr<i64> = calloc(400, 8)
    let dens: ptr<i64> = calloc(400, 8)
    if nums == null || dens == null { return 1 }
    let mut qc: i64 = 0
    let mut b: i64 = 2
    while b <= K {
        let mut a: i64 = 1
        while a < b {
            if gcd_abs(a, b) == 1 {
                nums[qc] = a
                dens[qc] = b
                qc = qc + 1
            }
            a = a + 1
        }
        b = b + 1
    }

    # set of sums: sn[i]/sd[i]
    let MAXS: i64 = 200000
    let sn: ptr<i64> = calloc(MAXS, 8)
    let sd: ptr<i64> = calloc(MAXS, 8)
    if sn == null || sd == null { return 1 }
    let mut sc: i64 = 0

    let mut i: i64 = 0
    while i < qc {
        let mut j: i64 = i
        while j < qc {
            let xn: i64 = nums[i]
            let xd: i64 = dens[i]
            let yn: i64 = nums[j]
            let yd: i64 = dens[j]

            # helper to check if zn/zd is in Q and add sum
            # n=1: z = x+y = (xn*yd + yn*xd)/(xd*yd)
            let mut zn: i64 = xn * yd + yn * xd
            let mut zd: i64 = xd * yd
            let g: i64 = gcd_abs(zn, zd)
            zn = zn / g
            zd = zd / g
            # check in Q: 0 < zn < zd <= K
            if zn > 0 && zn < zd && zd <= K {
                let mut snum: i64 = xn * yd * zd + yn * xd * zd + zn * xd * yd
                let mut sden: i64 = xd * yd * zd
                let g2: i64 = gcd_abs(snum, sden)
                snum = snum / g2
                sden = sden / g2
                # insert if new
                let mut found: bool = false
                let mut k: i64 = 0
                while k < sc {
                    if sn[k] == snum && sd[k] == sden { found = true; break }
                    k = k + 1
                }
                if !found {
                    sn[sc] = snum
                    sd[sc] = sden
                    sc = sc + 1
                }
            }

            # n=-1: z = xy/(x+y)
            let mut den_sum: i64 = xn * yd + yn * xd
            zn = xn * yn
            zd = den_sum
            # z = (xn/xd)*(yn/yd) / ((xn*yd+yn*xd)/(xd*yd)) = xn*yn*xd*yd / (xd*yd*den_sum) wait
            # z = (xn*yn)/(xd*yd) * (xd*yd)/den_sum = (xn*yn)/den_sum
            zn = xn * yn
            zd = den_sum
            let g3: i64 = gcd_abs(zn, zd)
            zn = zn / g3
            zd = zd / g3
            if zn > 0 && zn < zd && zd <= K {
                let mut snum: i64 = xn * yd * zd + yn * xd * zd + zn * xd * yd
                let mut sden: i64 = xd * yd * zd
                let g2: i64 = gcd_abs(snum, sden)
                snum = snum / g2
                sden = sden / g2
                let mut found: bool = false
                let mut k: i64 = 0
                while k < sc {
                    if sn[k] == snum && sd[k] == sden { found = true; break }
                    k = k + 1
                }
                if !found {
                    sn[sc] = snum
                    sd[sc] = sden
                    sc = sc + 1
                }
            }

            # n=2: z^2 = x^2+y^2
            let z2n: i64 = xn * xn * yd * yd + yn * yn * xd * xd
            let z2d: i64 = xd * xd * yd * yd
            let g4: i64 = gcd_abs(z2n, z2d)
            let z2nr: i64 = z2n / g4
            let z2dr: i64 = z2d / g4
            let sn2: i64 = isqrt(z2nr)
            let sd2: i64 = isqrt(z2dr)
            if sn2 * sn2 == z2nr && sd2 * sd2 == z2dr {
                zn = sn2
                zd = sd2
                let g5: i64 = gcd_abs(zn, zd)
                zn = zn / g5
                zd = zd / g5
                if zn > 0 && zn < zd && zd <= K {
                    let mut snum: i64 = xn * yd * zd + yn * xd * zd + zn * xd * yd
                    let mut sden: i64 = xd * yd * zd
                    let g2: i64 = gcd_abs(snum, sden)
                    snum = snum / g2
                    sden = sden / g2
                    let mut found: bool = false
                    let mut k: i64 = 0
                    while k < sc {
                        if sn[k] == snum && sd[k] == sden { found = true; break }
                        k = k + 1
                    }
                    if !found {
                        sn[sc] = snum
                        sd[sc] = sden
                        sc = sc + 1
                    }
                }
            }

            # n=-2: z^2 = (xy)^2 / (x^2+y^2)
            # z2 = (xn*yn)^2 / (z2n/z2d wait) = (xn*yn)^2 * z2d / z2n but use reduced
            # x^2+y^2 = z2nr/z2dr already computed above before reduce... recompute
            let x2y2n: i64 = xn * xn * yd * yd + yn * yn * xd * xd
            let x2y2d: i64 = xd * xd * yd * yd
            # z2 = (xn*yn/(xd*yd))^2 / (x2y2n/x2y2d) = (xn*yn)^2 * x2y2d / ((xd*yd)^2 * x2y2n)
            let mut zz_n: i64 = xn * yn * xn * yn * x2y2d
            let mut zz_d: i64 = xd * yd * xd * yd * x2y2n
            let g6: i64 = gcd_abs(zz_n, zz_d)
            zz_n = zz_n / g6
            zz_d = zz_d / g6
            let sn3: i64 = isqrt(zz_n)
            let sd3: i64 = isqrt(zz_d)
            if sn3 * sn3 == zz_n && sd3 * sd3 == zz_d {
                zn = sn3
                zd = sd3
                let g5: i64 = gcd_abs(zn, zd)
                zn = zn / g5
                zd = zd / g5
                if zn > 0 && zn < zd && zd <= K {
                    let mut snum: i64 = xn * yd * zd + yn * xd * zd + zn * xd * yd
                    let mut sden: i64 = xd * yd * zd
                    let g2: i64 = gcd_abs(snum, sden)
                    snum = snum / g2
                    sden = sden / g2
                    let mut found: bool = false
                    let mut k: i64 = 0
                    while k < sc {
                        if sn[k] == snum && sd[k] == sden { found = true; break }
                        k = k + 1
                    }
                    if !found {
                        sn[sc] = snum
                        sd[sc] = sden
                        sc = sc + 1
                    }
                }
            }

            j = j + 1
        }
        i = i + 1
    }

    # sum all fractions
    let mut tnum: i128 = 0
    let mut tden: i128 = 1
    i = 0
    while i < sc {
        # tnum/tden + sn/sd
        let a: i128 = sn[i] as i128
        let b: i128 = sd[i] as i128
        tnum = tnum * b + tden * a
        tden = tden * b
        # reduce
        let mut x: i128 = tnum
        if x < (0 as i128) { x = 0 as i128 - x }
        let mut y: i128 = tden
        while y != (0 as i128) {
            let t: i128 = x % y
            x = y
            y = t
        }
        let g: i128 = x
        tnum = tnum / g
        tden = tden / g
        i = i + 1
    }
    let ans: i64 = (tnum + tden) as i64
    printf("%lld\n", ans)
    free(nums); free(dens); free(sn); free(sd)
    return 0
}

Generated C

#include <stdint.h>
#include <stdbool.h>
#include <stdio.h>
#include <stdlib.h>
#include <string.h>

/* Flow runtime helpers */
typedef struct flow_temp_node { struct flow_temp_node* next; } flow_temp_node;
static flow_temp_node* flow_temp_head = NULL;
static int flow_temp_atexit_set = 0;
__attribute__((unused)) static void flow_temp_free_all(void) {
    while (flow_temp_head) {
        flow_temp_node* n = flow_temp_head;
        flow_temp_head = n->next;
        free(n);
    }
}
__attribute__((unused)) static void* flow_temp_alloc(size_t nbytes) {
    flow_temp_node* node = (flow_temp_node*)malloc(sizeof(flow_temp_node) + nbytes);
    if (!node) return NULL;
    node->next = flow_temp_head;
    flow_temp_head = node;
    if (!flow_temp_atexit_set) {
        flow_temp_atexit_set = 1;
        atexit(flow_temp_free_all);
    }
    return (void*)(node + 1);
}
#ifndef FLOW_DIAG
#define FLOW_DIAG(msg) fprintf(stderr, "%s", (msg))
#endif
#ifndef FLOW_LOG
#define FLOW_LOG(fmt, ...) printf(fmt, __VA_ARGS__)
#endif
#ifndef FLOW_LOG_EMPTY
#define FLOW_LOG_EMPTY(fmt) printf(fmt)
#endif
static char* flow_strcat(const char* a, const char* b) {
    size_t la = strlen(a ? a : ""), lb = strlen(b ? b : "");
    char* r = (char*)flow_temp_alloc(la + lb + 1);
    if (!r) return NULL;
    if (la) memcpy(r, a, la);
    if (lb) memcpy(r + la, b, lb);
    r[la + lb] = '\0';
    return r;
}

#define __flow_in_arr(arr, val) __extension__ ({ \
    int _found = 0; \
    size_t _n = sizeof(arr)/sizeof((arr)[0]); \
    for (size_t _i = 0; _i < _n; _i++) { \
        if ((arr)[_i] == (val)) { _found = 1; break; } \
    } _found; })

/* Unified fault handler (MISRA #279) — override with -DFLOW_FAULT_HANDLER=fn */
#ifndef FLOW_FAULT_HANDLER
__attribute__((unused)) static inline void flow_fault_handler(const char* msg) {
    fprintf(stderr, "flow: %s\n", msg ? msg : "fault");
    abort();
#if defined(__GNUC__) || defined(__clang__)
    __builtin_unreachable();
#endif
}
#else
#define flow_fault_handler FLOW_FAULT_HANDLER
#endif
#define flow_div_by_zero_handler() flow_fault_handler("division by zero")
#define flow_shift_ub_handler() flow_fault_handler("invalid shift (amount out of range or left-shift of negative)")

#ifndef FLOW_CHECKED_DIV
#define FLOW_CHECKED_DIV(L, R) (((R) != 0) ? ((L) / (R)) : (flow_div_by_zero_handler(), (L) * 0))
#endif
#ifndef FLOW_CHECKED_MOD
#define FLOW_CHECKED_MOD(L, R) (((R) != 0) ? ((L) % (R)) : (flow_div_by_zero_handler(), (L) * 0))
#endif
#ifndef FLOW_CHECKED_SHL
#define FLOW_CHECKED_SHL(L, R) ((((R) >= 0) && ((unsigned long long)(R) < (sizeof(L) * 8ull)) && ((L) >= 0)) ? ((L) << (R)) : (flow_shift_ub_handler(), (L) * 0))
#endif
#ifndef FLOW_CHECKED_SHR
#define FLOW_CHECKED_SHR(L, R) ((((R) >= 0) && ((unsigned long long)(R) < (sizeof(L) * 8ull))) ? ((L) >> (R)) : (flow_shift_ub_handler(), (L) * 0))
#endif

#include <math.h>

void* _ui_state = NULL;

static inline float i32_to_f32(int32_t v) { return (float)v; }

/* Host stub for @gpu kernels (device codegen replaces this). */
static inline int32_t gpu_thread_id(void) { return 0; }

int64_t gcd_i64_i64(int64_t a0, int64_t b0);
int64_t lcm_i64_i64(int64_t a, int64_t b);
int64_t isqrt_i64(int64_t n);
int64_t mulmod_i64_i64_i64(int64_t a0, int64_t b0, int64_t mod);
int64_t mod_pow_i64_i64_i64(int64_t base, int64_t exp, int64_t mod);
bool is_prime_i64(int64_t n);
int64_t gcd_abs_i64_i64(int64_t a, int64_t b);
int32_t main(void);

int64_t gcd_i64_i64(int64_t a0, int64_t b0) {
    int64_t a = a0;
    int64_t b = b0;
    while (b != 0) {
        int64_t t = FLOW_CHECKED_MOD((a), (b));
        a = b;
        b = t;
    }
    return a;
}

int64_t lcm_i64_i64(int64_t a, int64_t b) {
    if ((a == 0 || b == 0)) {
        return 0;
    }
    return (FLOW_CHECKED_DIV((a), (gcd_i64_i64(a, b))) * b);
}

int64_t isqrt_i64(int64_t n) {
    if (n < 2) {
        return n;
    }
    int64_t x = n;
    int64_t y = FLOW_CHECKED_DIV(((x + 1)), (2));
    while (y < x) {
        x = y;
        y = FLOW_CHECKED_DIV(((x + FLOW_CHECKED_DIV((n), (x)))), (2));
    }
    return x;
}

int64_t mulmod_i64_i64_i64(int64_t a0, int64_t b0, int64_t mod) {
    int64_t a = FLOW_CHECKED_MOD((a0), (mod));
    int64_t b = FLOW_CHECKED_MOD((b0), (mod));
    int64_t result = 0;
    while (b > 0) {
        if (FLOW_CHECKED_MOD((b), (2)) == 1) {
            result = FLOW_CHECKED_MOD(((result + a)), (mod));
        }
        a = FLOW_CHECKED_MOD(((a * 2)), (mod));
        b = FLOW_CHECKED_DIV((b), (2));
    }
    return result;
}

int64_t mod_pow_i64_i64_i64(int64_t base, int64_t exp, int64_t mod) {
    if (mod == 1) {
        return 0;
    }
    int64_t result = 1;
    int64_t b = FLOW_CHECKED_MOD((base), (mod));
    int64_t e = exp;
    while (e > 0) {
        if (FLOW_CHECKED_MOD((e), (2)) == 1) {
            result = mulmod_i64_i64_i64(result, b, mod);
        }
        b = mulmod_i64_i64_i64(b, b, mod);
        e = FLOW_CHECKED_DIV((e), (2));
    }
    return result;
}

bool is_prime_i64(int64_t n) {
    if (n < 2) {
        return 0;
    }
    if (n < 4) {
        return 1;
    }
    if ((FLOW_CHECKED_MOD((n), (2)) == 0 || FLOW_CHECKED_MOD((n), (3)) == 0)) {
        return 0;
    }
    int64_t i = 5;
    while ((i * i) <= n) {
        if ((FLOW_CHECKED_MOD((n), (i)) == 0 || FLOW_CHECKED_MOD((n), ((i + 2))) == 0)) {
            return 0;
        }
        i = (i + 6);
    }
    return 1;
}



int64_t gcd_abs_i64_i64(int64_t a, int64_t b) {
    int64_t x = a;
    if (x < 0) {
        x = (0 - x);
    }
    int64_t y = b;
    if (y < 0) {
        y = (0 - y);
    }
    while (y != 0) {
        int64_t t = FLOW_CHECKED_MOD((x), (y));
        x = y;
        y = t;
    }
    return x;
}

int32_t main(void) {
    int64_t K = 35;
    int64_t* nums = (int64_t*)(calloc(400, 8));
    int64_t* dens = (int64_t*)(calloc(400, 8));
    if ((nums == NULL || dens == NULL)) {
        return 1;
    }
    int64_t qc = 0;
    int64_t b = 2;
    while (b <= K) {
        int64_t a = 1;
        while (a < b) {
            if (gcd_abs_i64_i64(a, b) == 1) {
                nums[qc] = a;
                dens[qc] = b;
                qc = (qc + 1);
            }
            a = (a + 1);
        }
        b = (b + 1);
    }
    int64_t MAXS = 200000;
    int64_t* sn = (int64_t*)(calloc(MAXS, 8));
    int64_t* sd = (int64_t*)(calloc(MAXS, 8));
    if ((sn == NULL || sd == NULL)) {
        return 1;
    }
    int64_t sc = 0;
    int64_t i = 0;
    while (i < qc) {
        int64_t j = i;
        while (j < qc) {
            int64_t xn = nums[i];
            int64_t xd = dens[i];
            int64_t yn = nums[j];
            int64_t yd = dens[j];
            int64_t zn = ((xn * yd) + (yn * xd));
            int64_t zd = (xd * yd);
            int64_t g = gcd_abs_i64_i64(zn, zd);
            zn = FLOW_CHECKED_DIV((zn), (g));
            zd = FLOW_CHECKED_DIV((zd), (g));
            if (((zn > 0 && zn < zd) && zd <= K)) {
                int64_t snum = ((((xn * yd) * zd) + ((yn * xd) * zd)) + ((zn * xd) * yd));
                int64_t sden = ((xd * yd) * zd);
                int64_t g2 = gcd_abs_i64_i64(snum, sden);
                snum = FLOW_CHECKED_DIV((snum), (g2));
                sden = FLOW_CHECKED_DIV((sden), (g2));
                bool found = 0;
                int64_t k = 0;
                while (k < sc) {
                    if ((sn[k] == snum && sd[k] == sden)) {
                        found = 1;
                        break;
                    }
                    k = (k + 1);
                }
                if ((!(found))) {
                    sn[sc] = snum;
                    sd[sc] = sden;
                    sc = (sc + 1);
                }
            }
            int64_t den_sum = ((xn * yd) + (yn * xd));
            zn = (xn * yn);
            zd = den_sum;
            zn = (xn * yn);
            zd = den_sum;
            int64_t g3 = gcd_abs_i64_i64(zn, zd);
            zn = FLOW_CHECKED_DIV((zn), (g3));
            zd = FLOW_CHECKED_DIV((zd), (g3));
            if (((zn > 0 && zn < zd) && zd <= K)) {
                int64_t snum = ((((xn * yd) * zd) + ((yn * xd) * zd)) + ((zn * xd) * yd));
                int64_t sden = ((xd * yd) * zd);
                int64_t g2 = gcd_abs_i64_i64(snum, sden);
                snum = FLOW_CHECKED_DIV((snum), (g2));
                sden = FLOW_CHECKED_DIV((sden), (g2));
                bool found = 0;
                int64_t k = 0;
                while (k < sc) {
                    if ((sn[k] == snum && sd[k] == sden)) {
                        found = 1;
                        break;
                    }
                    k = (k + 1);
                }
                if ((!(found))) {
                    sn[sc] = snum;
                    sd[sc] = sden;
                    sc = (sc + 1);
                }
            }
            int64_t z2n = ((((xn * xn) * yd) * yd) + (((yn * yn) * xd) * xd));
            int64_t z2d = (((xd * xd) * yd) * yd);
            int64_t g4 = gcd_abs_i64_i64(z2n, z2d);
            int64_t z2nr = FLOW_CHECKED_DIV((z2n), (g4));
            int64_t z2dr = FLOW_CHECKED_DIV((z2d), (g4));
            int64_t sn2 = isqrt_i64(z2nr);
            int64_t sd2 = isqrt_i64(z2dr);
            if (((sn2 * sn2) == z2nr && (sd2 * sd2) == z2dr)) {
                zn = sn2;
                zd = sd2;
                int64_t g5 = gcd_abs_i64_i64(zn, zd);
                zn = FLOW_CHECKED_DIV((zn), (g5));
                zd = FLOW_CHECKED_DIV((zd), (g5));
                if (((zn > 0 && zn < zd) && zd <= K)) {
                    int64_t snum = ((((xn * yd) * zd) + ((yn * xd) * zd)) + ((zn * xd) * yd));
                    int64_t sden = ((xd * yd) * zd);
                    int64_t g2 = gcd_abs_i64_i64(snum, sden);
                    snum = FLOW_CHECKED_DIV((snum), (g2));
                    sden = FLOW_CHECKED_DIV((sden), (g2));
                    bool found = 0;
                    int64_t k = 0;
                    while (k < sc) {
                        if ((sn[k] == snum && sd[k] == sden)) {
                            found = 1;
                            break;
                        }
                        k = (k + 1);
                    }
                    if ((!(found))) {
                        sn[sc] = snum;
                        sd[sc] = sden;
                        sc = (sc + 1);
                    }
                }
            }
            int64_t x2y2n = ((((xn * xn) * yd) * yd) + (((yn * yn) * xd) * xd));
            int64_t x2y2d = (((xd * xd) * yd) * yd);
            int64_t zz_n = ((((xn * yn) * xn) * yn) * x2y2d);
            int64_t zz_d = ((((xd * yd) * xd) * yd) * x2y2n);
            int64_t g6 = gcd_abs_i64_i64(zz_n, zz_d);
            zz_n = FLOW_CHECKED_DIV((zz_n), (g6));
            zz_d = FLOW_CHECKED_DIV((zz_d), (g6));
            int64_t sn3 = isqrt_i64(zz_n);
            int64_t sd3 = isqrt_i64(zz_d);
            if (((sn3 * sn3) == zz_n && (sd3 * sd3) == zz_d)) {
                zn = sn3;
                zd = sd3;
                int64_t g5 = gcd_abs_i64_i64(zn, zd);
                zn = FLOW_CHECKED_DIV((zn), (g5));
                zd = FLOW_CHECKED_DIV((zd), (g5));
                if (((zn > 0 && zn < zd) && zd <= K)) {
                    int64_t snum = ((((xn * yd) * zd) + ((yn * xd) * zd)) + ((zn * xd) * yd));
                    int64_t sden = ((xd * yd) * zd);
                    int64_t g2 = gcd_abs_i64_i64(snum, sden);
                    snum = FLOW_CHECKED_DIV((snum), (g2));
                    sden = FLOW_CHECKED_DIV((sden), (g2));
                    bool found = 0;
                    int64_t k = 0;
                    while (k < sc) {
                        if ((sn[k] == snum && sd[k] == sden)) {
                            found = 1;
                            break;
                        }
                        k = (k + 1);
                    }
                    if ((!(found))) {
                        sn[sc] = snum;
                        sd[sc] = sden;
                        sc = (sc + 1);
                    }
                }
            }
            j = (j + 1);
        }
        i = (i + 1);
    }
    __int128 tnum = 0;
    __int128 tden = 1;
    i = 0;
    while (i < sc) {
        __int128 a = ((__int128)(sn[i]));
        __int128 b = ((__int128)(sd[i]));
        tnum = ((tnum * b) + (tden * a));
        tden = (tden * b);
        __int128 x = tnum;
        if (x < ((__int128)(0))) {
            x = (((__int128)(0)) - x);
        }
        __int128 y = tden;
        while (y != ((__int128)(0))) {
            __int128 t = FLOW_CHECKED_MOD((x), (y));
            x = y;
            y = t;
        }
        __int128 g = x;
        tnum = FLOW_CHECKED_DIV((tnum), (g));
        tden = FLOW_CHECKED_DIV((tden), (g));
        i = (i + 1);
    }
    int64_t ans = ((int64_t)((tnum + tden)));
    printf("%lld\n", ans);
    free(nums);
    free(dens);
    free(sn);
    free(sd);
    return 0;
}

Generated MLIR

module {
  llvm.func @printf(!llvm.ptr, ...) -> i32
  llvm.mlir.global internal constant @str_0("%lld\n\00") {addr_space = 0 : i32} : !llvm.array<6 x i8>
  func.func @gcd(%arg0: i64, %arg1: i64) -> i64 {
    %0 = llvm.mlir.constant(1 : i64) : i64
    %1 = llvm.alloca %0 x i64 : (i64) -> !llvm.ptr
    llvm.store %arg0, %1 : i64, !llvm.ptr
    %2 = llvm.mlir.constant(1 : i64) : i64
    %3 = llvm.alloca %2 x i64 : (i64) -> !llvm.ptr
    llvm.store %arg1, %3 : i64, !llvm.ptr
    cf.br ^bb0
    ^bb0:
    %4 = llvm.load %3 : !llvm.ptr -> i64
    %5 = arith.constant 0 : i32
    %7 = arith.extsi %5 : i32 to i64
    %6 = arith.cmpi ne, %4, %7 : i64
    cf.cond_br %6, ^bb1, ^bb2
    ^bb1:
      %8 = llvm.load %1 : !llvm.ptr -> i64
      %9 = llvm.load %3 : !llvm.ptr -> i64
      %10 = arith.remsi %8, %9 : i64
      %11 = llvm.load %3 : !llvm.ptr -> i64
      llvm.store %11, %1 : i64, !llvm.ptr
      llvm.store %10, %3 : i64, !llvm.ptr
      cf.br ^bb0
    ^bb2:
    %12 = llvm.load %1 : !llvm.ptr -> i64
    func.return %12 : i64
  }
  func.func @lcm(%arg0: i64, %arg1: i64) -> i64 {
    %13 = arith.constant 0 : i32
    %15 = arith.extsi %13 : i32 to i64
    %14 = arith.cmpi eq, %arg0, %15 : i64
    %16 = scf.if %14 -> (i1) {
      %17 = arith.constant true
      scf.yield %17 : i1
    } else {
      %18 = arith.constant 0 : i32
      %20 = arith.extsi %18 : i32 to i64
      %19 = arith.cmpi eq, %arg1, %20 : i64
      scf.yield %19 : i1
    }
    cf.cond_br %16, ^bb3, ^bb4
    ^bb3:
      %21 = arith.constant 0 : i32
      %22 = arith.extsi %21 : i32 to i64
      func.return %22 : i64
    ^bb4:
      cf.br ^bb5
    ^bb5:
    %23 = func.call @gcd(%arg0, %arg1) : (i64, i64) -> i64
    %24 = arith.divsi %arg0, %23 : i64
    %25 = arith.muli %24, %arg1 : i64
    func.return %25 : i64
  }
  func.func @isqrt(%arg0: i64) -> i64 {
    %26 = arith.constant 2 : i32
    %28 = arith.extsi %26 : i32 to i64
    %27 = arith.cmpi slt, %arg0, %28 : i64
    cf.cond_br %27, ^bb6, ^bb7
    ^bb6:
      func.return %arg0 : i64
    ^bb7:
      cf.br ^bb8
    ^bb8:
    %29 = llvm.mlir.constant(1 : i64) : i64
    %30 = llvm.alloca %29 x i64 : (i64) -> !llvm.ptr
    llvm.store %arg0, %30 : i64, !llvm.ptr
    %31 = llvm.load %30 : !llvm.ptr -> i64
    %32 = arith.constant 1 : i32
    %34 = arith.extsi %32 : i32 to i64
    %33 = arith.addi %31, %34 : i64
    %35 = arith.constant 2 : i32
    %37 = arith.extsi %35 : i32 to i64
    %36 = arith.divsi %33, %37 : i64
    %38 = llvm.mlir.constant(1 : i64) : i64
    %39 = llvm.alloca %38 x i64 : (i64) -> !llvm.ptr
    llvm.store %36, %39 : i64, !llvm.ptr
    cf.br ^bb9
    ^bb9:
    %40 = llvm.load %39 : !llvm.ptr -> i64
    %41 = llvm.load %30 : !llvm.ptr -> i64
    %42 = arith.cmpi slt, %40, %41 : i64
    cf.cond_br %42, ^bb10, ^bb11
    ^bb10:
      %43 = llvm.load %39 : !llvm.ptr -> i64
      llvm.store %43, %30 : i64, !llvm.ptr
      %44 = llvm.load %30 : !llvm.ptr -> i64
      %45 = llvm.load %30 : !llvm.ptr -> i64
      %46 = arith.divsi %arg0, %45 : i64
      %47 = arith.addi %44, %46 : i64
      %48 = arith.constant 2 : i32
      %50 = arith.extsi %48 : i32 to i64
      %49 = arith.divsi %47, %50 : i64
      llvm.store %49, %39 : i64, !llvm.ptr
      cf.br ^bb9
    ^bb11:
    %51 = llvm.load %30 : !llvm.ptr -> i64
    func.return %51 : i64
  }
  func.func @mulmod(%arg0: i64, %arg1: i64, %arg2: i64) -> i64 {
    %52 = arith.remsi %arg0, %arg2 : i64
    %53 = llvm.mlir.constant(1 : i64) : i64
    %54 = llvm.alloca %53 x i64 : (i64) -> !llvm.ptr
    llvm.store %52, %54 : i64, !llvm.ptr
    %55 = arith.remsi %arg1, %arg2 : i64
    %56 = llvm.mlir.constant(1 : i64) : i64
    %57 = llvm.alloca %56 x i64 : (i64) -> !llvm.ptr
    llvm.store %55, %57 : i64, !llvm.ptr
    %58 = arith.constant 0 : i32
    %59 = arith.extsi %58 : i32 to i64
    %60 = llvm.mlir.constant(1 : i64) : i64
    %61 = llvm.alloca %60 x i64 : (i64) -> !llvm.ptr
    llvm.store %59, %61 : i64, !llvm.ptr
    cf.br ^bb12
    ^bb12:
    %62 = llvm.load %57 : !llvm.ptr -> i64
    %63 = arith.constant 0 : i32
    %65 = arith.extsi %63 : i32 to i64
    %64 = arith.cmpi sgt, %62, %65 : i64
    cf.cond_br %64, ^bb13, ^bb14
    ^bb13:
      %66 = llvm.load %57 : !llvm.ptr -> i64
      %67 = arith.constant 2 : i32
      %69 = arith.extsi %67 : i32 to i64
      %68 = arith.remsi %66, %69 : i64
      %70 = arith.constant 1 : i32
      %72 = arith.extsi %70 : i32 to i64
      %71 = arith.cmpi eq, %68, %72 : i64
      cf.cond_br %71, ^bb15, ^bb16
      ^bb15:
        %73 = llvm.load %61 : !llvm.ptr -> i64
        %74 = llvm.load %54 : !llvm.ptr -> i64
        %75 = arith.addi %73, %74 : i64
        %76 = arith.remsi %75, %arg2 : i64
        llvm.store %76, %61 : i64, !llvm.ptr
        cf.br ^bb17
      ^bb16:
        cf.br ^bb17
      ^bb17:
      %77 = llvm.load %54 : !llvm.ptr -> i64
      %78 = arith.constant 2 : i32
      %80 = arith.extsi %78 : i32 to i64
      %79 = arith.muli %77, %80 : i64
      %81 = arith.remsi %79, %arg2 : i64
      llvm.store %81, %54 : i64, !llvm.ptr
      %82 = llvm.load %57 : !llvm.ptr -> i64
      %83 = arith.constant 2 : i32
      %85 = arith.extsi %83 : i32 to i64
      %84 = arith.divsi %82, %85 : i64
      llvm.store %84, %57 : i64, !llvm.ptr
      cf.br ^bb12
    ^bb14:
    %86 = llvm.load %61 : !llvm.ptr -> i64
    func.return %86 : i64
  }
  func.func @mod_pow(%arg0: i64, %arg1: i64, %arg2: i64) -> i64 {
    %87 = arith.constant 1 : i32
    %89 = arith.extsi %87 : i32 to i64
    %88 = arith.cmpi eq, %arg2, %89 : i64
    cf.cond_br %88, ^bb18, ^bb19
    ^bb18:
      %90 = arith.constant 0 : i32
      %91 = arith.extsi %90 : i32 to i64
      func.return %91 : i64
    ^bb19:
      cf.br ^bb20
    ^bb20:
    %92 = arith.constant 1 : i32
    %93 = arith.extsi %92 : i32 to i64
    %94 = llvm.mlir.constant(1 : i64) : i64
    %95 = llvm.alloca %94 x i64 : (i64) -> !llvm.ptr
    llvm.store %93, %95 : i64, !llvm.ptr
    %96 = arith.remsi %arg0, %arg2 : i64
    %97 = llvm.mlir.constant(1 : i64) : i64
    %98 = llvm.alloca %97 x i64 : (i64) -> !llvm.ptr
    llvm.store %96, %98 : i64, !llvm.ptr
    %99 = llvm.mlir.constant(1 : i64) : i64
    %100 = llvm.alloca %99 x i64 : (i64) -> !llvm.ptr
    llvm.store %arg1, %100 : i64, !llvm.ptr
    cf.br ^bb21
    ^bb21:
    %101 = llvm.load %100 : !llvm.ptr -> i64
    %102 = arith.constant 0 : i32
    %104 = arith.extsi %102 : i32 to i64
    %103 = arith.cmpi sgt, %101, %104 : i64
    cf.cond_br %103, ^bb22, ^bb23
    ^bb22:
      %105 = llvm.load %100 : !llvm.ptr -> i64
      %106 = arith.constant 2 : i32
      %108 = arith.extsi %106 : i32 to i64
      %107 = arith.remsi %105, %108 : i64
      %109 = arith.constant 1 : i32
      %111 = arith.extsi %109 : i32 to i64
      %110 = arith.cmpi eq, %107, %111 : i64
      cf.cond_br %110, ^bb24, ^bb25
      ^bb24:
        %113 = llvm.load %95 : !llvm.ptr -> i64
        %114 = llvm.load %98 : !llvm.ptr -> i64
        %112 = func.call @mulmod(%113, %114, %arg2) : (i64, i64, i64) -> i64
        llvm.store %112, %95 : i64, !llvm.ptr
        cf.br ^bb26
      ^bb25:
        cf.br ^bb26
      ^bb26:
      %116 = llvm.load %98 : !llvm.ptr -> i64
      %117 = llvm.load %98 : !llvm.ptr -> i64
      %115 = func.call @mulmod(%116, %117, %arg2) : (i64, i64, i64) -> i64
      llvm.store %115, %98 : i64, !llvm.ptr
      %118 = llvm.load %100 : !llvm.ptr -> i64
      %119 = arith.constant 2 : i32
      %121 = arith.extsi %119 : i32 to i64
      %120 = arith.divsi %118, %121 : i64
      llvm.store %120, %100 : i64, !llvm.ptr
      cf.br ^bb21
    ^bb23:
    %122 = llvm.load %95 : !llvm.ptr -> i64
    func.return %122 : i64
  }
  func.func @is_prime(%arg0: i64) -> i1 {
    %123 = arith.constant 2 : i32
    %125 = arith.extsi %123 : i32 to i64
    %124 = arith.cmpi slt, %arg0, %125 : i64
    cf.cond_br %124, ^bb27, ^bb28
    ^bb27:
      %126 = arith.constant 0 : i1
      func.return %126 : i1
    ^bb28:
      cf.br ^bb29
    ^bb29:
    %127 = arith.constant 4 : i32
    %129 = arith.extsi %127 : i32 to i64
    %128 = arith.cmpi slt, %arg0, %129 : i64
    cf.cond_br %128, ^bb30, ^bb31
    ^bb30:
      %130 = arith.constant 1 : i1
      func.return %130 : i1
    ^bb31:
      cf.br ^bb32
    ^bb32:
    %131 = arith.constant 2 : i32
    %133 = arith.extsi %131 : i32 to i64
    %132 = arith.remsi %arg0, %133 : i64
    %134 = arith.constant 0 : i32
    %136 = arith.extsi %134 : i32 to i64
    %135 = arith.cmpi eq, %132, %136 : i64
    %137 = scf.if %135 -> (i1) {
      %138 = arith.constant true
      scf.yield %138 : i1
    } else {
      %139 = arith.constant 3 : i32
      %141 = arith.extsi %139 : i32 to i64
      %140 = arith.remsi %arg0, %141 : i64
      %142 = arith.constant 0 : i32
      %144 = arith.extsi %142 : i32 to i64
      %143 = arith.cmpi eq, %140, %144 : i64
      scf.yield %143 : i1
    }
    cf.cond_br %137, ^bb33, ^bb34
    ^bb33:
      %145 = arith.constant 0 : i1
      func.return %145 : i1
    ^bb34:
      cf.br ^bb35
    ^bb35:
    %146 = arith.constant 5 : i32
    %147 = arith.extsi %146 : i32 to i64
    %148 = llvm.mlir.constant(1 : i64) : i64
    %149 = llvm.alloca %148 x i64 : (i64) -> !llvm.ptr
    llvm.store %147, %149 : i64, !llvm.ptr
    cf.br ^bb36
    ^bb36:
    %150 = llvm.load %149 : !llvm.ptr -> i64
    %151 = llvm.load %149 : !llvm.ptr -> i64
    %152 = arith.muli %150, %151 : i64
    %153 = arith.cmpi sle, %152, %arg0 : i64
    cf.cond_br %153, ^bb37, ^bb38
    ^bb37:
      %154 = llvm.load %149 : !llvm.ptr -> i64
      %155 = arith.remsi %arg0, %154 : i64
      %156 = arith.constant 0 : i32
      %158 = arith.extsi %156 : i32 to i64
      %157 = arith.cmpi eq, %155, %158 : i64
      %159 = scf.if %157 -> (i1) {
        %160 = arith.constant true
        scf.yield %160 : i1
      } else {
        %161 = llvm.load %149 : !llvm.ptr -> i64
        %162 = arith.constant 2 : i32
        %164 = arith.extsi %162 : i32 to i64
        %163 = arith.addi %161, %164 : i64
        %165 = arith.remsi %arg0, %163 : i64
        %166 = arith.constant 0 : i32
        %168 = arith.extsi %166 : i32 to i64
        %167 = arith.cmpi eq, %165, %168 : i64
        scf.yield %167 : i1
      }
      cf.cond_br %159, ^bb39, ^bb40
      ^bb39:
        %169 = arith.constant 0 : i1
        func.return %169 : i1
      ^bb40:
        cf.br ^bb41
      ^bb41:
      %170 = llvm.load %149 : !llvm.ptr -> i64
      %171 = arith.constant 6 : i32
      %173 = arith.extsi %171 : i32 to i64
      %172 = arith.addi %170, %173 : i64
      llvm.store %172, %149 : i64, !llvm.ptr
      cf.br ^bb36
    ^bb38:
    %174 = arith.constant 1 : i1
    func.return %174 : i1
  }
  func.func private @calloc(i64, i64) -> !llvm.ptr
  func.func private @free(!llvm.ptr) -> ()
  func.func @gcd_abs(%arg0: i64, %arg1: i64) -> i64 {
    %175 = llvm.mlir.constant(1 : i64) : i64
    %176 = llvm.alloca %175 x i64 : (i64) -> !llvm.ptr
    llvm.store %arg0, %176 : i64, !llvm.ptr
    %177 = llvm.load %176 : !llvm.ptr -> i64
    %178 = arith.constant 0 : i32
    %180 = arith.extsi %178 : i32 to i64
    %179 = arith.cmpi slt, %177, %180 : i64
    cf.cond_br %179, ^bb42, ^bb43
    ^bb42:
      %181 = arith.constant 0 : i32
      %182 = llvm.load %176 : !llvm.ptr -> i64
      %184 = arith.extsi %181 : i32 to i64
      %183 = arith.subi %184, %182 : i64
      llvm.store %183, %176 : i64, !llvm.ptr
      cf.br ^bb44
    ^bb43:
      cf.br ^bb44
    ^bb44:
    %185 = llvm.mlir.constant(1 : i64) : i64
    %186 = llvm.alloca %185 x i64 : (i64) -> !llvm.ptr
    llvm.store %arg1, %186 : i64, !llvm.ptr
    %187 = llvm.load %186 : !llvm.ptr -> i64
    %188 = arith.constant 0 : i32
    %190 = arith.extsi %188 : i32 to i64
    %189 = arith.cmpi slt, %187, %190 : i64
    cf.cond_br %189, ^bb45, ^bb46
    ^bb45:
      %191 = arith.constant 0 : i32
      %192 = llvm.load %186 : !llvm.ptr -> i64
      %194 = arith.extsi %191 : i32 to i64
      %193 = arith.subi %194, %192 : i64
      llvm.store %193, %186 : i64, !llvm.ptr
      cf.br ^bb47
    ^bb46:
      cf.br ^bb47
    ^bb47:
    cf.br ^bb48
    ^bb48:
    %195 = llvm.load %186 : !llvm.ptr -> i64
    %196 = arith.constant 0 : i32
    %198 = arith.extsi %196 : i32 to i64
    %197 = arith.cmpi ne, %195, %198 : i64
    cf.cond_br %197, ^bb49, ^bb50
    ^bb49:
      %199 = llvm.load %176 : !llvm.ptr -> i64
      %200 = llvm.load %186 : !llvm.ptr -> i64
      %201 = arith.remsi %199, %200 : i64
      %202 = llvm.load %186 : !llvm.ptr -> i64
      llvm.store %202, %176 : i64, !llvm.ptr
      llvm.store %201, %186 : i64, !llvm.ptr
      cf.br ^bb48
    ^bb50:
    %203 = llvm.load %176 : !llvm.ptr -> i64
    func.return %203 : i64
  }
  func.func @main() -> i32 {
    %204 = arith.constant 35 : i32
    %205 = arith.extsi %204 : i32 to i64
    %207 = arith.constant 400 : i32
    %208 = arith.constant 8 : i32
    %209 = arith.extsi %207 : i32 to i64
    %210 = arith.extsi %208 : i32 to i64
    %206 = func.call @calloc(%209, %210) : (i64, i64) -> !llvm.ptr
    %212 = arith.constant 400 : i32
    %213 = arith.constant 8 : i32
    %214 = arith.extsi %212 : i32 to i64
    %215 = arith.extsi %213 : i32 to i64
    %211 = func.call @calloc(%214, %215) : (i64, i64) -> !llvm.ptr
    %216 = llvm.mlir.zero : !llvm.ptr
    %217 = llvm.icmp "eq" %206, %216 : !llvm.ptr
    %218 = scf.if %217 -> (i1) {
      %219 = arith.constant true
      scf.yield %219 : i1
    } else {
      %220 = llvm.mlir.zero : !llvm.ptr
      %221 = llvm.icmp "eq" %211, %220 : !llvm.ptr
      scf.yield %221 : i1
    }
    cf.cond_br %218, ^bb51, ^bb52
    ^bb51:
      %222 = arith.constant 1 : i32
      func.return %222 : i32
    ^bb52:
      cf.br ^bb53
    ^bb53:
    %223 = arith.constant 0 : i32
    %224 = arith.extsi %223 : i32 to i64
    %225 = llvm.mlir.constant(1 : i64) : i64
    %226 = llvm.alloca %225 x i64 : (i64) -> !llvm.ptr
    llvm.store %224, %226 : i64, !llvm.ptr
    %227 = arith.constant 2 : i32
    %228 = arith.extsi %227 : i32 to i64
    %229 = llvm.mlir.constant(1 : i64) : i64
    %230 = llvm.alloca %229 x i64 : (i64) -> !llvm.ptr
    llvm.store %228, %230 : i64, !llvm.ptr
    cf.br ^bb54
    ^bb54:
    %231 = llvm.load %230 : !llvm.ptr -> i64
    %232 = arith.cmpi sle, %231, %205 : i64
    cf.cond_br %232, ^bb55, ^bb56
    ^bb55:
      %233 = arith.constant 1 : i32
      %234 = arith.extsi %233 : i32 to i64
      %235 = llvm.mlir.constant(1 : i64) : i64
      %236 = llvm.alloca %235 x i64 : (i64) -> !llvm.ptr
      llvm.store %234, %236 : i64, !llvm.ptr
      cf.br ^bb57
      ^bb57:
      %237 = llvm.load %236 : !llvm.ptr -> i64
      %238 = llvm.load %230 : !llvm.ptr -> i64
      %239 = arith.cmpi slt, %237, %238 : i64
      cf.cond_br %239, ^bb58, ^bb59
      ^bb58:
        %241 = llvm.load %236 : !llvm.ptr -> i64
        %242 = llvm.load %230 : !llvm.ptr -> i64
        %240 = func.call @gcd_abs(%241, %242) : (i64, i64) -> i64
        %243 = arith.constant 1 : i32
        %245 = arith.extsi %243 : i32 to i64
        %244 = arith.cmpi eq, %240, %245 : i64
        cf.cond_br %244, ^bb60, ^bb61
        ^bb60:
          %246 = llvm.load %236 : !llvm.ptr -> i64
          %247 = llvm.load %226 : !llvm.ptr -> i64
          %248 = llvm.getelementptr %206[%247] : (!llvm.ptr, i64) -> !llvm.ptr, i64
          llvm.store %246, %248 : i64, !llvm.ptr
          %249 = llvm.load %230 : !llvm.ptr -> i64
          %250 = llvm.load %226 : !llvm.ptr -> i64
          %251 = llvm.getelementptr %211[%250] : (!llvm.ptr, i64) -> !llvm.ptr, i64
          llvm.store %249, %251 : i64, !llvm.ptr
          %252 = llvm.load %226 : !llvm.ptr -> i64
          %253 = arith.constant 1 : i32
          %255 = arith.extsi %253 : i32 to i64
          %254 = arith.addi %252, %255 : i64
          llvm.store %254, %226 : i64, !llvm.ptr
          cf.br ^bb62
        ^bb61:
          cf.br ^bb62
        ^bb62:
        %256 = llvm.load %236 : !llvm.ptr -> i64
        %257 = arith.constant 1 : i32
        %259 = arith.extsi %257 : i32 to i64
        %258 = arith.addi %256, %259 : i64
        llvm.store %258, %236 : i64, !llvm.ptr
        cf.br ^bb57
      ^bb59:
      %260 = llvm.load %230 : !llvm.ptr -> i64
      %261 = arith.constant 1 : i32
      %263 = arith.extsi %261 : i32 to i64
      %262 = arith.addi %260, %263 : i64
      llvm.store %262, %230 : i64, !llvm.ptr
      cf.br ^bb54
    ^bb56:
    %264 = arith.constant 200000 : i32
    %265 = arith.extsi %264 : i32 to i64
    %267 = arith.constant 8 : i32
    %268 = arith.extsi %267 : i32 to i64
    %266 = func.call @calloc(%265, %268) : (i64, i64) -> !llvm.ptr
    %270 = arith.constant 8 : i32
    %271 = arith.extsi %270 : i32 to i64
    %269 = func.call @calloc(%265, %271) : (i64, i64) -> !llvm.ptr
    %272 = llvm.mlir.zero : !llvm.ptr
    %273 = llvm.icmp "eq" %266, %272 : !llvm.ptr
    %274 = scf.if %273 -> (i1) {
      %275 = arith.constant true
      scf.yield %275 : i1
    } else {
      %276 = llvm.mlir.zero : !llvm.ptr
      %277 = llvm.icmp "eq" %269, %276 : !llvm.ptr
      scf.yield %277 : i1
    }
    cf.cond_br %274, ^bb63, ^bb64
    ^bb63:
      %278 = arith.constant 1 : i32
      func.return %278 : i32
    ^bb64:
      cf.br ^bb65
    ^bb65:
    %279 = arith.constant 0 : i32
    %280 = arith.extsi %279 : i32 to i64
    %281 = llvm.mlir.constant(1 : i64) : i64
    %282 = llvm.alloca %281 x i64 : (i64) -> !llvm.ptr
    llvm.store %280, %282 : i64, !llvm.ptr
    %283 = arith.constant 0 : i32
    %284 = arith.extsi %283 : i32 to i64
    %285 = llvm.mlir.constant(1 : i64) : i64
    %286 = llvm.alloca %285 x i64 : (i64) -> !llvm.ptr
    llvm.store %284, %286 : i64, !llvm.ptr
    cf.br ^bb66
    ^bb66:
    %287 = llvm.load %286 : !llvm.ptr -> i64
    %288 = llvm.load %226 : !llvm.ptr -> i64
    %289 = arith.cmpi slt, %287, %288 : i64
    cf.cond_br %289, ^bb67, ^bb68
    ^bb67:
      %290 = llvm.load %286 : !llvm.ptr -> i64
      %291 = llvm.mlir.constant(1 : i64) : i64
      %292 = llvm.alloca %291 x i64 : (i64) -> !llvm.ptr
      llvm.store %290, %292 : i64, !llvm.ptr
      cf.br ^bb69
      ^bb69:
      %293 = llvm.load %292 : !llvm.ptr -> i64
      %294 = llvm.load %226 : !llvm.ptr -> i64
      %295 = arith.cmpi slt, %293, %294 : i64
      cf.cond_br %295, ^bb70, ^bb71
      ^bb70:
        %297 = llvm.load %286 : !llvm.ptr -> i64
        %298 = llvm.getelementptr %206[%297] : (!llvm.ptr, i64) -> !llvm.ptr, i64
        %296 = llvm.load %298 : !llvm.ptr -> i64
        %300 = llvm.load %286 : !llvm.ptr -> i64
        %301 = llvm.getelementptr %211[%300] : (!llvm.ptr, i64) -> !llvm.ptr, i64
        %299 = llvm.load %301 : !llvm.ptr -> i64
        %303 = llvm.load %292 : !llvm.ptr -> i64
        %304 = llvm.getelementptr %206[%303] : (!llvm.ptr, i64) -> !llvm.ptr, i64
        %302 = llvm.load %304 : !llvm.ptr -> i64
        %306 = llvm.load %292 : !llvm.ptr -> i64
        %307 = llvm.getelementptr %211[%306] : (!llvm.ptr, i64) -> !llvm.ptr, i64
        %305 = llvm.load %307 : !llvm.ptr -> i64
        %308 = arith.muli %296, %305 : i64
        %309 = arith.muli %302, %299 : i64
        %310 = arith.addi %308, %309 : i64
        %311 = llvm.mlir.constant(1 : i64) : i64
        %312 = llvm.alloca %311 x i64 : (i64) -> !llvm.ptr
        llvm.store %310, %312 : i64, !llvm.ptr
        %313 = arith.muli %299, %305 : i64
        %314 = llvm.mlir.constant(1 : i64) : i64
        %315 = llvm.alloca %314 x i64 : (i64) -> !llvm.ptr
        llvm.store %313, %315 : i64, !llvm.ptr
        %317 = llvm.load %312 : !llvm.ptr -> i64
        %318 = llvm.load %315 : !llvm.ptr -> i64
        %316 = func.call @gcd_abs(%317, %318) : (i64, i64) -> i64
        %319 = llvm.load %312 : !llvm.ptr -> i64
        %320 = arith.divsi %319, %316 : i64
        llvm.store %320, %312 : i64, !llvm.ptr
        %321 = llvm.load %315 : !llvm.ptr -> i64
        %322 = arith.divsi %321, %316 : i64
        llvm.store %322, %315 : i64, !llvm.ptr
        %323 = llvm.load %312 : !llvm.ptr -> i64
        %324 = arith.constant 0 : i32
        %326 = arith.extsi %324 : i32 to i64
        %325 = arith.cmpi sgt, %323, %326 : i64
        %327 = scf.if %325 -> (i1) {
          %328 = llvm.load %312 : !llvm.ptr -> i64
          %329 = llvm.load %315 : !llvm.ptr -> i64
          %330 = arith.cmpi slt, %328, %329 : i64
          scf.yield %330 : i1
        } else {
          %331 = arith.constant false
          scf.yield %331 : i1
        }
        %332 = scf.if %327 -> (i1) {
          %333 = llvm.load %315 : !llvm.ptr -> i64
          %334 = arith.cmpi sle, %333, %205 : i64
          scf.yield %334 : i1
        } else {
          %335 = arith.constant false
          scf.yield %335 : i1
        }
        cf.cond_br %332, ^bb72, ^bb73
        ^bb72:
          %336 = arith.muli %296, %305 : i64
          %337 = llvm.load %315 : !llvm.ptr -> i64
          %338 = arith.muli %336, %337 : i64
          %339 = arith.muli %302, %299 : i64
          %340 = llvm.load %315 : !llvm.ptr -> i64
          %341 = arith.muli %339, %340 : i64
          %342 = arith.addi %338, %341 : i64
          %343 = llvm.load %312 : !llvm.ptr -> i64
          %344 = arith.muli %343, %299 : i64
          %345 = arith.muli %344, %305 : i64
          %346 = arith.addi %342, %345 : i64
          %347 = llvm.mlir.constant(1 : i64) : i64
          %348 = llvm.alloca %347 x i64 : (i64) -> !llvm.ptr
          llvm.store %346, %348 : i64, !llvm.ptr
          %349 = arith.muli %299, %305 : i64
          %350 = llvm.load %315 : !llvm.ptr -> i64
          %351 = arith.muli %349, %350 : i64
          %352 = llvm.mlir.constant(1 : i64) : i64
          %353 = llvm.alloca %352 x i64 : (i64) -> !llvm.ptr
          llvm.store %351, %353 : i64, !llvm.ptr
          %355 = llvm.load %348 : !llvm.ptr -> i64
          %356 = llvm.load %353 : !llvm.ptr -> i64
          %354 = func.call @gcd_abs(%355, %356) : (i64, i64) -> i64
          %357 = llvm.load %348 : !llvm.ptr -> i64
          %358 = arith.divsi %357, %354 : i64
          llvm.store %358, %348 : i64, !llvm.ptr
          %359 = llvm.load %353 : !llvm.ptr -> i64
          %360 = arith.divsi %359, %354 : i64
          llvm.store %360, %353 : i64, !llvm.ptr
          %361 = arith.constant 0 : i1
          %362 = llvm.mlir.constant(1 : i64) : i64
          %363 = llvm.alloca %362 x i1 : (i64) -> !llvm.ptr
          llvm.store %361, %363 : i1, !llvm.ptr
          %364 = arith.constant 0 : i32
          %365 = arith.extsi %364 : i32 to i64
          %366 = llvm.mlir.constant(1 : i64) : i64
          %367 = llvm.alloca %366 x i64 : (i64) -> !llvm.ptr
          llvm.store %365, %367 : i64, !llvm.ptr
          cf.br ^bb75
          ^bb75:
          %368 = llvm.load %367 : !llvm.ptr -> i64
          %369 = llvm.load %282 : !llvm.ptr -> i64
          %370 = arith.cmpi slt, %368, %369 : i64
          cf.cond_br %370, ^bb76, ^bb77
          ^bb76:
            %372 = llvm.load %367 : !llvm.ptr -> i64
            %373 = llvm.getelementptr %266[%372] : (!llvm.ptr, i64) -> !llvm.ptr, i64
            %371 = llvm.load %373 : !llvm.ptr -> i64
            %374 = llvm.load %348 : !llvm.ptr -> i64
            %375 = arith.cmpi eq, %371, %374 : i64
            %376 = scf.if %375 -> (i1) {
              %378 = llvm.load %367 : !llvm.ptr -> i64
              %379 = llvm.getelementptr %269[%378] : (!llvm.ptr, i64) -> !llvm.ptr, i64
              %377 = llvm.load %379 : !llvm.ptr -> i64
              %380 = llvm.load %353 : !llvm.ptr -> i64
              %381 = arith.cmpi eq, %377, %380 : i64
              scf.yield %381 : i1
            } else {
              %382 = arith.constant false
              scf.yield %382 : i1
            }
            cf.cond_br %376, ^bb78, ^bb79
            ^bb78:
              %383 = arith.constant 1 : i1
              llvm.store %383, %363 : i1, !llvm.ptr
              cf.br ^bb77
            ^bb79:
              cf.br ^bb80
            ^bb80:
            %384 = llvm.load %367 : !llvm.ptr -> i64
            %385 = arith.constant 1 : i32
            %387 = arith.extsi %385 : i32 to i64
            %386 = arith.addi %384, %387 : i64
            llvm.store %386, %367 : i64, !llvm.ptr
            cf.br ^bb75
          ^bb77:
          %388 = llvm.load %363 : !llvm.ptr -> i1
          %390 = arith.constant 1 : i1
          %389 = arith.xori %388, %390 : i1
          cf.cond_br %389, ^bb81, ^bb82
          ^bb81:
            %392 = llvm.load %348 : !llvm.ptr -> i64
            %393 = llvm.load %282 : !llvm.ptr -> i64
            %394 = llvm.getelementptr %266[%393] : (!llvm.ptr, i64) -> !llvm.ptr, i64
            llvm.store %392, %394 : i64, !llvm.ptr
            %395 = llvm.load %353 : !llvm.ptr -> i64
            %396 = llvm.load %282 : !llvm.ptr -> i64
            %397 = llvm.getelementptr %269[%396] : (!llvm.ptr, i64) -> !llvm.ptr, i64
            llvm.store %395, %397 : i64, !llvm.ptr
            %398 = llvm.load %282 : !llvm.ptr -> i64
            %399 = arith.constant 1 : i32
            %401 = arith.extsi %399 : i32 to i64
            %400 = arith.addi %398, %401 : i64
            llvm.store %400, %282 : i64, !llvm.ptr
            cf.br ^bb83
          ^bb82:
            cf.br ^bb83
          ^bb83:
          cf.br ^bb74
        ^bb73:
          cf.br ^bb74
        ^bb74:
        %402 = arith.muli %296, %305 : i64
        %403 = arith.muli %302, %299 : i64
        %404 = arith.addi %402, %403 : i64
        %405 = llvm.mlir.constant(1 : i64) : i64
        %406 = llvm.alloca %405 x i64 : (i64) -> !llvm.ptr
        llvm.store %404, %406 : i64, !llvm.ptr
        %407 = arith.muli %296, %302 : i64
        llvm.store %407, %312 : i64, !llvm.ptr
        %408 = llvm.load %406 : !llvm.ptr -> i64
        llvm.store %408, %315 : i64, !llvm.ptr
        %409 = arith.muli %296, %302 : i64
        llvm.store %409, %312 : i64, !llvm.ptr
        %410 = llvm.load %406 : !llvm.ptr -> i64
        llvm.store %410, %315 : i64, !llvm.ptr
        %412 = llvm.load %312 : !llvm.ptr -> i64
        %413 = llvm.load %315 : !llvm.ptr -> i64
        %411 = func.call @gcd_abs(%412, %413) : (i64, i64) -> i64
        %414 = llvm.load %312 : !llvm.ptr -> i64
        %415 = arith.divsi %414, %411 : i64
        llvm.store %415, %312 : i64, !llvm.ptr
        %416 = llvm.load %315 : !llvm.ptr -> i64
        %417 = arith.divsi %416, %411 : i64
        llvm.store %417, %315 : i64, !llvm.ptr
        %418 = llvm.load %312 : !llvm.ptr -> i64
        %419 = arith.constant 0 : i32
        %421 = arith.extsi %419 : i32 to i64
        %420 = arith.cmpi sgt, %418, %421 : i64
        %422 = scf.if %420 -> (i1) {
          %423 = llvm.load %312 : !llvm.ptr -> i64
          %424 = llvm.load %315 : !llvm.ptr -> i64
          %425 = arith.cmpi slt, %423, %424 : i64
          scf.yield %425 : i1
        } else {
          %426 = arith.constant false
          scf.yield %426 : i1
        }
        %427 = scf.if %422 -> (i1) {
          %428 = llvm.load %315 : !llvm.ptr -> i64
          %429 = arith.cmpi sle, %428, %205 : i64
          scf.yield %429 : i1
        } else {
          %430 = arith.constant false
          scf.yield %430 : i1
        }
        cf.cond_br %427, ^bb84, ^bb85
        ^bb84:
          %431 = arith.muli %296, %305 : i64
          %432 = llvm.load %315 : !llvm.ptr -> i64
          %433 = arith.muli %431, %432 : i64
          %434 = arith.muli %302, %299 : i64
          %435 = llvm.load %315 : !llvm.ptr -> i64
          %436 = arith.muli %434, %435 : i64
          %437 = arith.addi %433, %436 : i64
          %438 = llvm.load %312 : !llvm.ptr -> i64
          %439 = arith.muli %438, %299 : i64
          %440 = arith.muli %439, %305 : i64
          %441 = arith.addi %437, %440 : i64
          %442 = llvm.mlir.constant(1 : i64) : i64
          %443 = llvm.alloca %442 x i64 : (i64) -> !llvm.ptr
          llvm.store %441, %443 : i64, !llvm.ptr
          %444 = arith.muli %299, %305 : i64
          %445 = llvm.load %315 : !llvm.ptr -> i64
          %446 = arith.muli %444, %445 : i64
          %447 = llvm.mlir.constant(1 : i64) : i64
          %448 = llvm.alloca %447 x i64 : (i64) -> !llvm.ptr
          llvm.store %446, %448 : i64, !llvm.ptr
          %450 = llvm.load %443 : !llvm.ptr -> i64
          %451 = llvm.load %448 : !llvm.ptr -> i64
          %449 = func.call @gcd_abs(%450, %451) : (i64, i64) -> i64
          %452 = llvm.load %443 : !llvm.ptr -> i64
          %453 = arith.divsi %452, %449 : i64
          llvm.store %453, %443 : i64, !llvm.ptr
          %454 = llvm.load %448 : !llvm.ptr -> i64
          %455 = arith.divsi %454, %449 : i64
          llvm.store %455, %448 : i64, !llvm.ptr
          %456 = arith.constant 0 : i1
          %457 = llvm.mlir.constant(1 : i64) : i64
          %458 = llvm.alloca %457 x i1 : (i64) -> !llvm.ptr
          llvm.store %456, %458 : i1, !llvm.ptr
          %459 = arith.constant 0 : i32
          %460 = arith.extsi %459 : i32 to i64
          %461 = llvm.mlir.constant(1 : i64) : i64
          %462 = llvm.alloca %461 x i64 : (i64) -> !llvm.ptr
          llvm.store %460, %462 : i64, !llvm.ptr
          cf.br ^bb87
          ^bb87:
          %463 = llvm.load %462 : !llvm.ptr -> i64
          %464 = llvm.load %282 : !llvm.ptr -> i64
          %465 = arith.cmpi slt, %463, %464 : i64
          cf.cond_br %465, ^bb88, ^bb89
          ^bb88:
            %467 = llvm.load %462 : !llvm.ptr -> i64
            %468 = llvm.getelementptr %266[%467] : (!llvm.ptr, i64) -> !llvm.ptr, i64
            %466 = llvm.load %468 : !llvm.ptr -> i64
            %469 = llvm.load %443 : !llvm.ptr -> i64
            %470 = arith.cmpi eq, %466, %469 : i64
            %471 = scf.if %470 -> (i1) {
              %473 = llvm.load %462 : !llvm.ptr -> i64
              %474 = llvm.getelementptr %269[%473] : (!llvm.ptr, i64) -> !llvm.ptr, i64
              %472 = llvm.load %474 : !llvm.ptr -> i64
              %475 = llvm.load %448 : !llvm.ptr -> i64
              %476 = arith.cmpi eq, %472, %475 : i64
              scf.yield %476 : i1
            } else {
              %477 = arith.constant false
              scf.yield %477 : i1
            }
            cf.cond_br %471, ^bb90, ^bb91
            ^bb90:
              %478 = arith.constant 1 : i1
              llvm.store %478, %458 : i1, !llvm.ptr
              cf.br ^bb89
            ^bb91:
              cf.br ^bb92
            ^bb92:
            %479 = llvm.load %462 : !llvm.ptr -> i64
            %480 = arith.constant 1 : i32
            %482 = arith.extsi %480 : i32 to i64
            %481 = arith.addi %479, %482 : i64
            llvm.store %481, %462 : i64, !llvm.ptr
            cf.br ^bb87
          ^bb89:
          %483 = llvm.load %458 : !llvm.ptr -> i1
          %485 = arith.constant 1 : i1
          %484 = arith.xori %483, %485 : i1
          cf.cond_br %484, ^bb93, ^bb94
          ^bb93:
            %487 = llvm.load %443 : !llvm.ptr -> i64
            %488 = llvm.load %282 : !llvm.ptr -> i64
            %489 = llvm.getelementptr %266[%488] : (!llvm.ptr, i64) -> !llvm.ptr, i64
            llvm.store %487, %489 : i64, !llvm.ptr
            %490 = llvm.load %448 : !llvm.ptr -> i64
            %491 = llvm.load %282 : !llvm.ptr -> i64
            %492 = llvm.getelementptr %269[%491] : (!llvm.ptr, i64) -> !llvm.ptr, i64
            llvm.store %490, %492 : i64, !llvm.ptr
            %493 = llvm.load %282 : !llvm.ptr -> i64
            %494 = arith.constant 1 : i32
            %496 = arith.extsi %494 : i32 to i64
            %495 = arith.addi %493, %496 : i64
            llvm.store %495, %282 : i64, !llvm.ptr
            cf.br ^bb95
          ^bb94:
            cf.br ^bb95
          ^bb95:
          cf.br ^bb86
        ^bb85:
          cf.br ^bb86
        ^bb86:
        %497 = arith.muli %296, %296 : i64
        %498 = arith.muli %497, %305 : i64
        %499 = arith.muli %498, %305 : i64
        %500 = arith.muli %302, %302 : i64
        %501 = arith.muli %500, %299 : i64
        %502 = arith.muli %501, %299 : i64
        %503 = arith.addi %499, %502 : i64
        %504 = arith.muli %299, %299 : i64
        %505 = arith.muli %504, %305 : i64
        %506 = arith.muli %505, %305 : i64
        %507 = func.call @gcd_abs(%503, %506) : (i64, i64) -> i64
        %508 = arith.divsi %503, %507 : i64
        %509 = arith.divsi %506, %507 : i64
        %510 = func.call @isqrt(%508) : (i64) -> i64
        %511 = func.call @isqrt(%509) : (i64) -> i64
        %512 = arith.muli %510, %510 : i64
        %513 = arith.cmpi eq, %512, %508 : i64
        %514 = scf.if %513 -> (i1) {
          %515 = arith.muli %511, %511 : i64
          %516 = arith.cmpi eq, %515, %509 : i64
          scf.yield %516 : i1
        } else {
          %517 = arith.constant false
          scf.yield %517 : i1
        }
        cf.cond_br %514, ^bb96, ^bb97
        ^bb96:
          llvm.store %510, %312 : i64, !llvm.ptr
          llvm.store %511, %315 : i64, !llvm.ptr
          %519 = llvm.load %312 : !llvm.ptr -> i64
          %520 = llvm.load %315 : !llvm.ptr -> i64
          %518 = func.call @gcd_abs(%519, %520) : (i64, i64) -> i64
          %521 = llvm.load %312 : !llvm.ptr -> i64
          %522 = arith.divsi %521, %518 : i64
          llvm.store %522, %312 : i64, !llvm.ptr
          %523 = llvm.load %315 : !llvm.ptr -> i64
          %524 = arith.divsi %523, %518 : i64
          llvm.store %524, %315 : i64, !llvm.ptr
          %525 = llvm.load %312 : !llvm.ptr -> i64
          %526 = arith.constant 0 : i32
          %528 = arith.extsi %526 : i32 to i64
          %527 = arith.cmpi sgt, %525, %528 : i64
          %529 = scf.if %527 -> (i1) {
            %530 = llvm.load %312 : !llvm.ptr -> i64
            %531 = llvm.load %315 : !llvm.ptr -> i64
            %532 = arith.cmpi slt, %530, %531 : i64
            scf.yield %532 : i1
          } else {
            %533 = arith.constant false
            scf.yield %533 : i1
          }
          %534 = scf.if %529 -> (i1) {
            %535 = llvm.load %315 : !llvm.ptr -> i64
            %536 = arith.cmpi sle, %535, %205 : i64
            scf.yield %536 : i1
          } else {
            %537 = arith.constant false
            scf.yield %537 : i1
          }
          cf.cond_br %534, ^bb99, ^bb100
          ^bb99:
            %538 = arith.muli %296, %305 : i64
            %539 = llvm.load %315 : !llvm.ptr -> i64
            %540 = arith.muli %538, %539 : i64
            %541 = arith.muli %302, %299 : i64
            %542 = llvm.load %315 : !llvm.ptr -> i64
            %543 = arith.muli %541, %542 : i64
            %544 = arith.addi %540, %543 : i64
            %545 = llvm.load %312 : !llvm.ptr -> i64
            %546 = arith.muli %545, %299 : i64
            %547 = arith.muli %546, %305 : i64
            %548 = arith.addi %544, %547 : i64
            %549 = llvm.mlir.constant(1 : i64) : i64
            %550 = llvm.alloca %549 x i64 : (i64) -> !llvm.ptr
            llvm.store %548, %550 : i64, !llvm.ptr
            %551 = arith.muli %299, %305 : i64
            %552 = llvm.load %315 : !llvm.ptr -> i64
            %553 = arith.muli %551, %552 : i64
            %554 = llvm.mlir.constant(1 : i64) : i64
            %555 = llvm.alloca %554 x i64 : (i64) -> !llvm.ptr
            llvm.store %553, %555 : i64, !llvm.ptr
            %557 = llvm.load %550 : !llvm.ptr -> i64
            %558 = llvm.load %555 : !llvm.ptr -> i64
            %556 = func.call @gcd_abs(%557, %558) : (i64, i64) -> i64
            %559 = llvm.load %550 : !llvm.ptr -> i64
            %560 = arith.divsi %559, %556 : i64
            llvm.store %560, %550 : i64, !llvm.ptr
            %561 = llvm.load %555 : !llvm.ptr -> i64
            %562 = arith.divsi %561, %556 : i64
            llvm.store %562, %555 : i64, !llvm.ptr
            %563 = arith.constant 0 : i1
            %564 = llvm.mlir.constant(1 : i64) : i64
            %565 = llvm.alloca %564 x i1 : (i64) -> !llvm.ptr
            llvm.store %563, %565 : i1, !llvm.ptr
            %566 = arith.constant 0 : i32
            %567 = arith.extsi %566 : i32 to i64
            %568 = llvm.mlir.constant(1 : i64) : i64
            %569 = llvm.alloca %568 x i64 : (i64) -> !llvm.ptr
            llvm.store %567, %569 : i64, !llvm.ptr
            cf.br ^bb102
            ^bb102:
            %570 = llvm.load %569 : !llvm.ptr -> i64
            %571 = llvm.load %282 : !llvm.ptr -> i64
            %572 = arith.cmpi slt, %570, %571 : i64
            cf.cond_br %572, ^bb103, ^bb104
            ^bb103:
              %574 = llvm.load %569 : !llvm.ptr -> i64
              %575 = llvm.getelementptr %266[%574] : (!llvm.ptr, i64) -> !llvm.ptr, i64
              %573 = llvm.load %575 : !llvm.ptr -> i64
              %576 = llvm.load %550 : !llvm.ptr -> i64
              %577 = arith.cmpi eq, %573, %576 : i64
              %578 = scf.if %577 -> (i1) {
                %580 = llvm.load %569 : !llvm.ptr -> i64
                %581 = llvm.getelementptr %269[%580] : (!llvm.ptr, i64) -> !llvm.ptr, i64
                %579 = llvm.load %581 : !llvm.ptr -> i64
                %582 = llvm.load %555 : !llvm.ptr -> i64
                %583 = arith.cmpi eq, %579, %582 : i64
                scf.yield %583 : i1
              } else {
                %584 = arith.constant false
                scf.yield %584 : i1
              }
              cf.cond_br %578, ^bb105, ^bb106
              ^bb105:
                %585 = arith.constant 1 : i1
                llvm.store %585, %565 : i1, !llvm.ptr
                cf.br ^bb104
              ^bb106:
                cf.br ^bb107
              ^bb107:
              %586 = llvm.load %569 : !llvm.ptr -> i64
              %587 = arith.constant 1 : i32
              %589 = arith.extsi %587 : i32 to i64
              %588 = arith.addi %586, %589 : i64
              llvm.store %588, %569 : i64, !llvm.ptr
              cf.br ^bb102
            ^bb104:
            %590 = llvm.load %565 : !llvm.ptr -> i1
            %592 = arith.constant 1 : i1
            %591 = arith.xori %590, %592 : i1
            cf.cond_br %591, ^bb108, ^bb109
            ^bb108:
              %594 = llvm.load %550 : !llvm.ptr -> i64
              %595 = llvm.load %282 : !llvm.ptr -> i64
              %596 = llvm.getelementptr %266[%595] : (!llvm.ptr, i64) -> !llvm.ptr, i64
              llvm.store %594, %596 : i64, !llvm.ptr
              %597 = llvm.load %555 : !llvm.ptr -> i64
              %598 = llvm.load %282 : !llvm.ptr -> i64
              %599 = llvm.getelementptr %269[%598] : (!llvm.ptr, i64) -> !llvm.ptr, i64
              llvm.store %597, %599 : i64, !llvm.ptr
              %600 = llvm.load %282 : !llvm.ptr -> i64
              %601 = arith.constant 1 : i32
              %603 = arith.extsi %601 : i32 to i64
              %602 = arith.addi %600, %603 : i64
              llvm.store %602, %282 : i64, !llvm.ptr
              cf.br ^bb110
            ^bb109:
              cf.br ^bb110
            ^bb110:
            cf.br ^bb101
          ^bb100:
            cf.br ^bb101
          ^bb101:
          cf.br ^bb98
        ^bb97:
          cf.br ^bb98
        ^bb98:
        %604 = arith.muli %296, %296 : i64
        %605 = arith.muli %604, %305 : i64
        %606 = arith.muli %605, %305 : i64
        %607 = arith.muli %302, %302 : i64
        %608 = arith.muli %607, %299 : i64
        %609 = arith.muli %608, %299 : i64
        %610 = arith.addi %606, %609 : i64
        %611 = arith.muli %299, %299 : i64
        %612 = arith.muli %611, %305 : i64
        %613 = arith.muli %612, %305 : i64
        %614 = arith.muli %296, %302 : i64
        %615 = arith.muli %614, %296 : i64
        %616 = arith.muli %615, %302 : i64
        %617 = arith.muli %616, %613 : i64
        %618 = llvm.mlir.constant(1 : i64) : i64
        %619 = llvm.alloca %618 x i64 : (i64) -> !llvm.ptr
        llvm.store %617, %619 : i64, !llvm.ptr
        %620 = arith.muli %299, %305 : i64
        %621 = arith.muli %620, %299 : i64
        %622 = arith.muli %621, %305 : i64
        %623 = arith.muli %622, %610 : i64
        %624 = llvm.mlir.constant(1 : i64) : i64
        %625 = llvm.alloca %624 x i64 : (i64) -> !llvm.ptr
        llvm.store %623, %625 : i64, !llvm.ptr
        %627 = llvm.load %619 : !llvm.ptr -> i64
        %628 = llvm.load %625 : !llvm.ptr -> i64
        %626 = func.call @gcd_abs(%627, %628) : (i64, i64) -> i64
        %629 = llvm.load %619 : !llvm.ptr -> i64
        %630 = arith.divsi %629, %626 : i64
        llvm.store %630, %619 : i64, !llvm.ptr
        %631 = llvm.load %625 : !llvm.ptr -> i64
        %632 = arith.divsi %631, %626 : i64
        llvm.store %632, %625 : i64, !llvm.ptr
        %634 = llvm.load %619 : !llvm.ptr -> i64
        %633 = func.call @isqrt(%634) : (i64) -> i64
        %636 = llvm.load %625 : !llvm.ptr -> i64
        %635 = func.call @isqrt(%636) : (i64) -> i64
        %637 = arith.muli %633, %633 : i64
        %638 = llvm.load %619 : !llvm.ptr -> i64
        %639 = arith.cmpi eq, %637, %638 : i64
        %640 = scf.if %639 -> (i1) {
          %641 = arith.muli %635, %635 : i64
          %642 = llvm.load %625 : !llvm.ptr -> i64
          %643 = arith.cmpi eq, %641, %642 : i64
          scf.yield %643 : i1
        } else {
          %644 = arith.constant false
          scf.yield %644 : i1
        }
        cf.cond_br %640, ^bb111, ^bb112
        ^bb111:
          llvm.store %633, %312 : i64, !llvm.ptr
          llvm.store %635, %315 : i64, !llvm.ptr
          %646 = llvm.load %312 : !llvm.ptr -> i64
          %647 = llvm.load %315 : !llvm.ptr -> i64
          %645 = func.call @gcd_abs(%646, %647) : (i64, i64) -> i64
          %648 = llvm.load %312 : !llvm.ptr -> i64
          %649 = arith.divsi %648, %645 : i64
          llvm.store %649, %312 : i64, !llvm.ptr
          %650 = llvm.load %315 : !llvm.ptr -> i64
          %651 = arith.divsi %650, %645 : i64
          llvm.store %651, %315 : i64, !llvm.ptr
          %652 = llvm.load %312 : !llvm.ptr -> i64
          %653 = arith.constant 0 : i32
          %655 = arith.extsi %653 : i32 to i64
          %654 = arith.cmpi sgt, %652, %655 : i64
          %656 = scf.if %654 -> (i1) {
            %657 = llvm.load %312 : !llvm.ptr -> i64
            %658 = llvm.load %315 : !llvm.ptr -> i64
            %659 = arith.cmpi slt, %657, %658 : i64
            scf.yield %659 : i1
          } else {
            %660 = arith.constant false
            scf.yield %660 : i1
          }
          %661 = scf.if %656 -> (i1) {
            %662 = llvm.load %315 : !llvm.ptr -> i64
            %663 = arith.cmpi sle, %662, %205 : i64
            scf.yield %663 : i1
          } else {
            %664 = arith.constant false
            scf.yield %664 : i1
          }
          cf.cond_br %661, ^bb114, ^bb115
          ^bb114:
            %665 = arith.muli %296, %305 : i64
            %666 = llvm.load %315 : !llvm.ptr -> i64
            %667 = arith.muli %665, %666 : i64
            %668 = arith.muli %302, %299 : i64
            %669 = llvm.load %315 : !llvm.ptr -> i64
            %670 = arith.muli %668, %669 : i64
            %671 = arith.addi %667, %670 : i64
            %672 = llvm.load %312 : !llvm.ptr -> i64
            %673 = arith.muli %672, %299 : i64
            %674 = arith.muli %673, %305 : i64
            %675 = arith.addi %671, %674 : i64
            %676 = llvm.mlir.constant(1 : i64) : i64
            %677 = llvm.alloca %676 x i64 : (i64) -> !llvm.ptr
            llvm.store %675, %677 : i64, !llvm.ptr
            %678 = arith.muli %299, %305 : i64
            %679 = llvm.load %315 : !llvm.ptr -> i64
            %680 = arith.muli %678, %679 : i64
            %681 = llvm.mlir.constant(1 : i64) : i64
            %682 = llvm.alloca %681 x i64 : (i64) -> !llvm.ptr
            llvm.store %680, %682 : i64, !llvm.ptr
            %684 = llvm.load %677 : !llvm.ptr -> i64
            %685 = llvm.load %682 : !llvm.ptr -> i64
            %683 = func.call @gcd_abs(%684, %685) : (i64, i64) -> i64
            %686 = llvm.load %677 : !llvm.ptr -> i64
            %687 = arith.divsi %686, %683 : i64
            llvm.store %687, %677 : i64, !llvm.ptr
            %688 = llvm.load %682 : !llvm.ptr -> i64
            %689 = arith.divsi %688, %683 : i64
            llvm.store %689, %682 : i64, !llvm.ptr
            %690 = arith.constant 0 : i1
            %691 = llvm.mlir.constant(1 : i64) : i64
            %692 = llvm.alloca %691 x i1 : (i64) -> !llvm.ptr
            llvm.store %690, %692 : i1, !llvm.ptr
            %693 = arith.constant 0 : i32
            %694 = arith.extsi %693 : i32 to i64
            %695 = llvm.mlir.constant(1 : i64) : i64
            %696 = llvm.alloca %695 x i64 : (i64) -> !llvm.ptr
            llvm.store %694, %696 : i64, !llvm.ptr
            cf.br ^bb117
            ^bb117:
            %697 = llvm.load %696 : !llvm.ptr -> i64
            %698 = llvm.load %282 : !llvm.ptr -> i64
            %699 = arith.cmpi slt, %697, %698 : i64
            cf.cond_br %699, ^bb118, ^bb119
            ^bb118:
              %701 = llvm.load %696 : !llvm.ptr -> i64
              %702 = llvm.getelementptr %266[%701] : (!llvm.ptr, i64) -> !llvm.ptr, i64
              %700 = llvm.load %702 : !llvm.ptr -> i64
              %703 = llvm.load %677 : !llvm.ptr -> i64
              %704 = arith.cmpi eq, %700, %703 : i64
              %705 = scf.if %704 -> (i1) {
                %707 = llvm.load %696 : !llvm.ptr -> i64
                %708 = llvm.getelementptr %269[%707] : (!llvm.ptr, i64) -> !llvm.ptr, i64
                %706 = llvm.load %708 : !llvm.ptr -> i64
                %709 = llvm.load %682 : !llvm.ptr -> i64
                %710 = arith.cmpi eq, %706, %709 : i64
                scf.yield %710 : i1
              } else {
                %711 = arith.constant false
                scf.yield %711 : i1
              }
              cf.cond_br %705, ^bb120, ^bb121
              ^bb120:
                %712 = arith.constant 1 : i1
                llvm.store %712, %692 : i1, !llvm.ptr
                cf.br ^bb119
              ^bb121:
                cf.br ^bb122
              ^bb122:
              %713 = llvm.load %696 : !llvm.ptr -> i64
              %714 = arith.constant 1 : i32
              %716 = arith.extsi %714 : i32 to i64
              %715 = arith.addi %713, %716 : i64
              llvm.store %715, %696 : i64, !llvm.ptr
              cf.br ^bb117
            ^bb119:
            %717 = llvm.load %692 : !llvm.ptr -> i1
            %719 = arith.constant 1 : i1
            %718 = arith.xori %717, %719 : i1
            cf.cond_br %718, ^bb123, ^bb124
            ^bb123:
              %721 = llvm.load %677 : !llvm.ptr -> i64
              %722 = llvm.load %282 : !llvm.ptr -> i64
              %723 = llvm.getelementptr %266[%722] : (!llvm.ptr, i64) -> !llvm.ptr, i64
              llvm.store %721, %723 : i64, !llvm.ptr
              %724 = llvm.load %682 : !llvm.ptr -> i64
              %725 = llvm.load %282 : !llvm.ptr -> i64
              %726 = llvm.getelementptr %269[%725] : (!llvm.ptr, i64) -> !llvm.ptr, i64
              llvm.store %724, %726 : i64, !llvm.ptr
              %727 = llvm.load %282 : !llvm.ptr -> i64
              %728 = arith.constant 1 : i32
              %730 = arith.extsi %728 : i32 to i64
              %729 = arith.addi %727, %730 : i64
              llvm.store %729, %282 : i64, !llvm.ptr
              cf.br ^bb125
            ^bb124:
              cf.br ^bb125
            ^bb125:
            cf.br ^bb116
          ^bb115:
            cf.br ^bb116
          ^bb116:
          cf.br ^bb113
        ^bb112:
          cf.br ^bb113
        ^bb113:
        %731 = llvm.load %292 : !llvm.ptr -> i64
        %732 = arith.constant 1 : i32
        %734 = arith.extsi %732 : i32 to i64
        %733 = arith.addi %731, %734 : i64
        llvm.store %733, %292 : i64, !llvm.ptr
        cf.br ^bb69
      ^bb71:
      %735 = llvm.load %286 : !llvm.ptr -> i64
      %736 = arith.constant 1 : i32
      %738 = arith.extsi %736 : i32 to i64
      %737 = arith.addi %735, %738 : i64
      llvm.store %737, %286 : i64, !llvm.ptr
      cf.br ^bb66
    ^bb68:
    %739 = arith.constant 0 : i32
    %740 = arith.extsi %739 : i32 to i128
    %741 = llvm.mlir.constant(1 : i64) : i64
    %742 = llvm.alloca %741 x i128 : (i64) -> !llvm.ptr
    llvm.store %740, %742 : i128, !llvm.ptr
    %743 = arith.constant 1 : i32
    %744 = arith.extsi %743 : i32 to i128
    %745 = llvm.mlir.constant(1 : i64) : i64
    %746 = llvm.alloca %745 x i128 : (i64) -> !llvm.ptr
    llvm.store %744, %746 : i128, !llvm.ptr
    %747 = arith.constant 0 : i32
    %748 = arith.extsi %747 : i32 to i64
    llvm.store %748, %286 : i64, !llvm.ptr
    cf.br ^bb126
    ^bb126:
    %749 = llvm.load %286 : !llvm.ptr -> i64
    %750 = llvm.load %282 : !llvm.ptr -> i64
    %751 = arith.cmpi slt, %749, %750 : i64
    cf.cond_br %751, ^bb127, ^bb128
    ^bb127:
      %753 = llvm.load %286 : !llvm.ptr -> i64
      %754 = llvm.getelementptr %266[%753] : (!llvm.ptr, i64) -> !llvm.ptr, i64
      %752 = llvm.load %754 : !llvm.ptr -> i64
      %755 = arith.extsi %752 : i64 to i128
      %757 = llvm.load %286 : !llvm.ptr -> i64
      %758 = llvm.getelementptr %269[%757] : (!llvm.ptr, i64) -> !llvm.ptr, i64
      %756 = llvm.load %758 : !llvm.ptr -> i64
      %759 = arith.extsi %756 : i64 to i128
      %760 = llvm.load %742 : !llvm.ptr -> i128
      %762 = arith.trunci %760 : i128 to i64
      %763 = arith.trunci %759 : i128 to i64
      %761 = arith.muli %762, %763 : i64
      %764 = llvm.load %746 : !llvm.ptr -> i128
      %766 = arith.trunci %764 : i128 to i64
      %767 = arith.trunci %755 : i128 to i64
      %765 = arith.muli %766, %767 : i64
      %768 = arith.addi %761, %765 : i64
      %769 = arith.extsi %768 : i64 to i128
      llvm.store %769, %742 : i128, !llvm.ptr
      %770 = llvm.load %746 : !llvm.ptr -> i128
      %772 = arith.trunci %770 : i128 to i64
      %773 = arith.trunci %759 : i128 to i64
      %771 = arith.muli %772, %773 : i64
      %774 = arith.extsi %771 : i64 to i128
      llvm.store %774, %746 : i128, !llvm.ptr
      %775 = llvm.load %742 : !llvm.ptr -> i128
      %776 = llvm.mlir.constant(1 : i64) : i64
      %777 = llvm.alloca %776 x i128 : (i64) -> !llvm.ptr
      llvm.store %775, %777 : i128, !llvm.ptr
      %778 = llvm.load %777 : !llvm.ptr -> i128
      %779 = arith.constant 0 : i32
      %780 = arith.extsi %779 : i32 to i128
      %782 = arith.trunci %778 : i128 to i64
      %783 = arith.trunci %780 : i128 to i64
      %781 = arith.cmpi slt, %782, %783 : i64
      cf.cond_br %781, ^bb129, ^bb130
      ^bb129:
        %784 = arith.constant 0 : i32
        %785 = arith.extsi %784 : i32 to i128
        %786 = llvm.load %777 : !llvm.ptr -> i128
        %788 = arith.trunci %785 : i128 to i64
        %789 = arith.trunci %786 : i128 to i64
        %787 = arith.subi %788, %789 : i64
        %790 = arith.extsi %787 : i64 to i128
        llvm.store %790, %777 : i128, !llvm.ptr
        cf.br ^bb131
      ^bb130:
        cf.br ^bb131
      ^bb131:
      %791 = llvm.load %746 : !llvm.ptr -> i128
      %792 = llvm.mlir.constant(1 : i64) : i64
      %793 = llvm.alloca %792 x i128 : (i64) -> !llvm.ptr
      llvm.store %791, %793 : i128, !llvm.ptr
      cf.br ^bb132
      ^bb132:
      %794 = llvm.load %793 : !llvm.ptr -> i128
      %795 = arith.constant 0 : i32
      %796 = arith.extsi %795 : i32 to i128
      %798 = arith.trunci %794 : i128 to i64
      %799 = arith.trunci %796 : i128 to i64
      %797 = arith.cmpi ne, %798, %799 : i64
      cf.cond_br %797, ^bb133, ^bb134
      ^bb133:
        %800 = llvm.load %777 : !llvm.ptr -> i128
        %801 = llvm.load %793 : !llvm.ptr -> i128
        %803 = arith.trunci %800 : i128 to i64
        %804 = arith.trunci %801 : i128 to i64
        %802 = arith.remsi %803, %804 : i64
        %805 = arith.extsi %802 : i64 to i128
        %806 = llvm.load %793 : !llvm.ptr -> i128
        llvm.store %806, %777 : i128, !llvm.ptr
        llvm.store %805, %793 : i128, !llvm.ptr
        cf.br ^bb132
      ^bb134:
      %807 = llvm.load %777 : !llvm.ptr -> i128
      %808 = llvm.load %742 : !llvm.ptr -> i128
      %810 = arith.trunci %808 : i128 to i64
      %811 = arith.trunci %807 : i128 to i64
      %809 = arith.divsi %810, %811 : i64
      %812 = arith.extsi %809 : i64 to i128
      llvm.store %812, %742 : i128, !llvm.ptr
      %813 = llvm.load %746 : !llvm.ptr -> i128
      %815 = arith.trunci %813 : i128 to i64
      %816 = arith.trunci %807 : i128 to i64
      %814 = arith.divsi %815, %816 : i64
      %817 = arith.extsi %814 : i64 to i128
      llvm.store %817, %746 : i128, !llvm.ptr
      %818 = llvm.load %286 : !llvm.ptr -> i64
      %819 = arith.constant 1 : i32
      %821 = arith.extsi %819 : i32 to i64
      %820 = arith.addi %818, %821 : i64
      llvm.store %820, %286 : i64, !llvm.ptr
      cf.br ^bb126
    ^bb128:
    %822 = llvm.load %742 : !llvm.ptr -> i128
    %823 = llvm.load %746 : !llvm.ptr -> i128
    %825 = arith.trunci %822 : i128 to i64
    %826 = arith.trunci %823 : i128 to i64
    %824 = arith.addi %825, %826 : i64
    %827 = llvm.mlir.addressof @str_0 : !llvm.ptr
    %828 = llvm.call @printf(%827, %824) vararg(!llvm.func<i32 (ptr, ...)>) : (!llvm.ptr, i64) -> i32
    func.call @free(%206) : (!llvm.ptr) -> ()
    func.call @free(%211) : (!llvm.ptr) -> ()
    func.call @free(%266) : (!llvm.ptr) -> ()
    func.call @free(%269) : (!llvm.ptr) -> ()
    %833 = arith.constant 0 : i32
    func.return %833 : i32
  }
}