Problem 838

f(N) = smallest positive integer not coprime to any n<=N whose LSD is 3. Prints ln(f(10^6)) rounded to 6 digits. Ported from native C to pure Flow.

Answer250591.442792
Output250591.442792
StatusPASS
Native helperno
Runtime0 ms
Peak memory4864 KB
Time complexityO(n^2) (estimated)
Space complexityO(n^2) (estimated)

Performance comparison

MetricOur solutionBest known
Time complexityO(n^2)?
Space complexityO(n^2)?
ApproachFlow solutionNot curated
VerdictUnknown

Flow source

# Project Euler 838: Not Coprime.
# f(N) = smallest positive integer not coprime to any n<=N whose LSD is 3.
# Prints ln(f(10^6)) rounded to 6 digits.
# Ported from native C to pure Flow.

extern {
    function calloc(n: i64, size: i64) -> ptr<void>
    function free(p: ptr<void>) -> void
    function memset(p: ptr<void>, c: i32, n: i64) -> ptr<void>
    function log(x: f64) -> f64
}

const N: i64 = 1000000
const SCALE: i64 = 1000000000000

# ---- Dinic max-flow globals (linked-list adjacency) ----
let mut g_head: ptr<i64> = null
let mut e_to: ptr<i64> = null
let mut e_cap: ptr<i64> = null
let mut e_rev: ptr<i64> = null
let mut e_next: ptr<i64> = null
let mut g_level: ptr<i64> = null
let mut g_it: ptr<i64> = null
let mut g_q: ptr<i64> = null
let mut g_n: i64 = 0
let mut g_edge_cnt: i64 = 0

function add_edge(u: i64, v: i64, cap: i64) -> void {
    let e1: i64 = g_edge_cnt
    g_edge_cnt = g_edge_cnt + 1
    let e2: i64 = g_edge_cnt
    g_edge_cnt = g_edge_cnt + 1
    e_to[e1] = v
    e_cap[e1] = cap
    e_rev[e1] = e2
    e_next[e1] = g_head[u]
    g_head[u] = e1
    e_to[e2] = u
    e_cap[e2] = 0
    e_rev[e2] = e1
    e_next[e2] = g_head[v]
    g_head[v] = e2
}

function bfs(s: i64, t: i64) -> bool {
    let mut i: i64 = 0
    while i < g_n {
        g_level[i] = -1
        i = i + 1
    }
    let mut qh: i64 = 0
    let mut qt: i64 = 0
    g_level[s] = 0
    g_q[qt] = s
    qt = qt + 1
    while qh < qt {
        let u: i64 = g_q[qh]
        qh = qh + 1
        let mut e: i64 = g_head[u]
        while e != -1 {
            let v: i64 = e_to[e]
            if e_cap[e] > 0 && g_level[v] < 0 {
                g_level[v] = g_level[u] + 1
                g_q[qt] = v
                qt = qt + 1
            }
            e = e_next[e]
        }
    }
    return g_level[t] >= 0
}

function dfs(u: i64, t: i64, f: i64) -> i64 {
    if u == t { return f }
    let mut e: i64 = g_it[u]
    while e != -1 {
        let v: i64 = e_to[e]
        if e_cap[e] > 0 && g_level[v] == g_level[u] + 1 {
            let mut push: i64 = f
            if e_cap[e] < push { push = e_cap[e] }
            let pushed: i64 = dfs(v, t, push)
            if pushed > 0 {
                e_cap[e] = e_cap[e] - pushed
                e_cap[e_rev[e]] = e_cap[e_rev[e]] + pushed
                g_it[u] = e
                return pushed
            }
        }
        e = e_next[e]
        g_it[u] = e
    }
    return 0
}

function max_flow(s: i64, t: i64) -> i64 {
    let mut flow: i64 = 0
    let INF: i64 = 1000000000000000000
    while bfs(s, t) {
        let mut i: i64 = 0
        while i < g_n {
            g_it[i] = g_head[i]
            i = i + 1
        }
        let mut pushed: i64 = dfs(s, t, INF)
        while pushed > 0 {
            flow = flow + pushed
            pushed = dfs(s, t, INF)
        }
    }
    return flow
}

function reachable(s: i64, vis: ptr<i8>) -> void {
    let mut i: i64 = 0
    while i < g_n {
        vis[i] = 0
        i = i + 1
    }
    let mut qh: i64 = 0
    let mut qt: i64 = 0
    vis[s] = 1
    g_q[qt] = s
    qt = qt + 1
    while qh < qt {
        let u: i64 = g_q[qh]
        qh = qh + 1
        let mut e: i64 = g_head[u]
        while e != -1 {
            let v: i64 = e_to[e]
            if e_cap[e] > 0 && vis[v] == 0 {
                vis[v] = 1
                g_q[qt] = v
                qt = qt + 1
            }
            e = e_next[e]
        }
    }
}

function iroot3(n: i64) -> i64 {
    let mut r: i64 = 0
    let mut cont: bool = true
    while cont {
        let nr: i64 = r + 1
        let c: i128 = (nr as i128) * (nr as i128) * (nr as i128)
        if c <= (n as i128) {
            r = nr
        } else {
            cont = false
        }
    }
    return r
}

function main() -> i32 {
    # Sieve primes up to N
    let is_prime: ptr<i8> = calloc(N + 1, 1) as ptr<i8>
    memset(is_prime, 1, N + 1)
    is_prime[0] = 0
    is_prime[1] = 0
    let mut p: i64 = 2
    while p * p <= N {
        if is_prime[p] != 0 {
            let mut j: i64 = p * p
            while j <= N {
                is_prime[j] = 0
                j = j + p
            }
        }
        p = p + 1
    }

    # Count and collect primes
    let mut pc: i64 = 0
    let mut i: i64 = 2
    while i <= N {
        if is_prime[i] != 0 { pc = pc + 1 }
        i = i + 1
    }
    let primes: ptr<i64> = calloc(pc, 8) as ptr<i64>
    pc = 0
    i = 2
    while i <= N {
        if is_prime[i] != 0 {
            primes[pc] = i
            pc = pc + 1
        }
        i = i + 1
    }

    let cbrt: i64 = iroot3(N)

    # Build left list: p%10==7, p*19<=N, p>cbrt (not forced)
    let left: ptr<i64> = calloc(pc, 8) as ptr<i64>
    let mut nL: i64 = 0
    i = 0
    while i < pc {
        let pp: i64 = primes[i]
        if pp % 10 == 7 && pp * 19 <= N && pp > cbrt {
            left[nL] = pp
            nL = nL + 1
        }
        i = i + 1
    }

    let mut ln_total: f64 = 0.0

    # Add forced primes' log contributions
    i = 0
    while i < pc {
        let pp: i64 = primes[i]
        if pp % 10 == 3 {
            ln_total = ln_total + log(pp as f64)
        } else {
            if pp % 10 == 7 && pp <= cbrt {
                ln_total = ln_total + log(pp as f64)
            }
        }
        i = i + 1
    }

    if nL == 0 {
        printf("%.6f\n", ln_total)
        free(is_prime)
        free(primes)
        free(left)
        return 0
    }

    let pmin: i64 = left[0]
    # Build right list: q%10==9, q <= N/pmin
    let right: ptr<i64> = calloc(pc, 8) as ptr<i64>
    let mut nR: i64 = 0
    let rlim: i64 = N / pmin
    i = 0
    while i < pc {
        let qq: i64 = primes[i]
        if qq % 10 == 9 && qq <= rlim {
            right[nR] = qq
            nR = nR + 1
        }
        i = i + 1
    }

    if nR == 0 {
        printf("%.6f\n", ln_total)
        free(is_prime)
        free(primes)
        free(left)
        free(right)
        return 0
    }

    # prefix lengths: for each left[i], count right[j] <= N/left[i]
    let pref_len: ptr<i64> = calloc(nL, 8) as ptr<i64>
    let mut total_pref: i64 = 0
    i = 0
    while i < nL {
        let lim: i64 = N / left[i]
        let mut lo: i64 = 0
        let mut hi: i64 = nR
        while lo < hi {
            let mid: i64 = (lo + hi) / 2
            if right[mid] <= lim {
                lo = mid + 1
            } else {
                hi = mid
            }
        }
        pref_len[i] = lo
        total_pref = total_pref + lo
        i = i + 1
    }

    # Weights scaled by 1e12
    let wL: ptr<i64> = calloc(nL, 8) as ptr<i64>
    let wR: ptr<i64> = calloc(nR, 8) as ptr<i64>
    let mut total_w: i64 = 0
    i = 0
    while i < nL {
        wL[i] = (log(left[i] as f64) * (SCALE as f64) + 0.5) as i64
        total_w = total_w + wL[i]
        i = i + 1
    }
    i = 0
    while i < nR {
        wR[i] = (log(right[i] as f64) * (SCALE as f64) + 0.5) as i64
        total_w = total_w + wR[i]
        i = i + 1
    }
    let INF: i64 = total_w + 1

    # Build flow network: S=0, offL=1, offR=1+nL, T=1+nL+nR
    let nNodes: i64 = 2 + nL + nR
    let S: i64 = 0
    let offL: i64 = 1
    let offR: i64 = 1 + nL
    let T: i64 = 1 + nL + nR

    let total_edges: i64 = 2 * (nL + nR + total_pref)
    g_n = nNodes
    g_edge_cnt = 0
    g_head = calloc(nNodes, 8) as ptr<i64>
    i = 0
    while i < nNodes {
        g_head[i] = -1
        i = i + 1
    }
    e_to = calloc(total_edges, 8) as ptr<i64>
    e_cap = calloc(total_edges, 8) as ptr<i64>
    e_rev = calloc(total_edges, 8) as ptr<i64>
    e_next = calloc(total_edges, 8) as ptr<i64>
    g_level = calloc(nNodes, 8) as ptr<i64>
    g_it = calloc(nNodes, 8) as ptr<i64>
    g_q = calloc(nNodes, 8) as ptr<i64>

    i = 0
    while i < nL {
        add_edge(S, offL + i, wL[i])
        i = i + 1
    }
    i = 0
    while i < nR {
        add_edge(offR + i, T, wR[i])
        i = i + 1
    }
    i = 0
    while i < nL {
        let u: i64 = offL + i
        let k: i64 = pref_len[i]
        let mut j: i64 = 0
        while j < k {
            add_edge(u, offR + j, INF)
            j = j + 1
        }
        i = i + 1
    }

    max_flow(S, T)

    let vis: ptr<i8> = calloc(nNodes, 1) as ptr<i8>
    reachable(S, vis)

    # min vertex cover = (Left not reachable) U (Right reachable)
    i = 0
    while i < nL {
        if vis[offL + i] == 0 {
            ln_total = ln_total + log(left[i] as f64)
        }
        i = i + 1
    }
    i = 0
    while i < nR {
        if vis[offR + i] != 0 {
            ln_total = ln_total + log(right[i] as f64)
        }
        i = i + 1
    }

    printf("%.6f\n", ln_total)

    free(vis)
    free(g_head)
    free(e_to)
    free(e_cap)
    free(e_rev)
    free(e_next)
    free(g_level)
    free(g_it)
    free(g_q)
    free(wL)
    free(wR)
    free(pref_len)
    free(left)
    free(right)
    free(primes)
    free(is_prime)
    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; }

void add_edge_i64_i64_i64(int64_t u, int64_t v, int64_t cap);
bool bfs_i64_i64(int64_t s, int64_t t);
int64_t dfs_i64_i64_i64(int64_t u, int64_t t, int64_t f);
int64_t max_flow_i64_i64(int64_t s, int64_t t);
void reachable_i64_ptr_i8(int64_t s, int8_t* vis);
int64_t iroot3_i64(int64_t n);
int32_t main(void);

static const int64_t N = 1000000;
static const int64_t SCALE = 1000000000000;

/* Module statics */
static int64_t* g_head = NULL;
static int64_t* e_to = NULL;
static int64_t* e_cap = NULL;
static int64_t* e_rev = NULL;
static int64_t* e_next = NULL;
static int64_t* g_level = NULL;
static int64_t* g_it = NULL;
static int64_t* g_q = NULL;
static int64_t g_n = 0;
static int64_t g_edge_cnt = 0;





void add_edge_i64_i64_i64(int64_t u, int64_t v, int64_t cap) {
    int64_t e1 = g_edge_cnt;
    g_edge_cnt = (g_edge_cnt + 1);
    int64_t e2 = g_edge_cnt;
    g_edge_cnt = (g_edge_cnt + 1);
    e_to[e1] = v;
    e_cap[e1] = cap;
    e_rev[e1] = e2;
    e_next[e1] = g_head[u];
    g_head[u] = e1;
    e_to[e2] = u;
    e_cap[e2] = 0;
    e_rev[e2] = e1;
    e_next[e2] = g_head[v];
    g_head[v] = e2;
}

bool bfs_i64_i64(int64_t s, int64_t t) {
    int64_t i = 0;
    while (i < g_n) {
        g_level[i] = (-1);
        i = (i + 1);
    }
    int64_t qh = 0;
    int64_t qt = 0;
    g_level[s] = 0;
    g_q[qt] = s;
    qt = (qt + 1);
    while (qh < qt) {
        int64_t u = g_q[qh];
        qh = (qh + 1);
        int64_t e = g_head[u];
        while (e != (-1)) {
            int64_t v = e_to[e];
            if ((e_cap[e] > 0 && g_level[v] < 0)) {
                g_level[v] = (g_level[u] + 1);
                g_q[qt] = v;
                qt = (qt + 1);
            }
            e = e_next[e];
        }
    }
    return g_level[t] >= 0;
}

int64_t dfs_i64_i64_i64(int64_t u, int64_t t, int64_t f) {
    if (u == t) {
        return f;
    }
    int64_t e = g_it[u];
    while (e != (-1)) {
        int64_t v = e_to[e];
        if ((e_cap[e] > 0 && g_level[v] == (g_level[u] + 1))) {
            int64_t push = f;
            if (e_cap[e] < push) {
                push = e_cap[e];
            }
            int64_t pushed = dfs_i64_i64_i64(v, t, push);
            if (pushed > 0) {
                e_cap[e] = (e_cap[e] - pushed);
                e_cap[e_rev[e]] = (e_cap[e_rev[e]] + pushed);
                g_it[u] = e;
                return pushed;
            }
        }
        e = e_next[e];
        g_it[u] = e;
    }
    return 0;
}

int64_t max_flow_i64_i64(int64_t s, int64_t t) {
    int64_t flow = 0;
    int64_t INF = 1000000000000000000;
    while (bfs_i64_i64(s, t)) {
        int64_t i = 0;
        while (i < g_n) {
            g_it[i] = g_head[i];
            i = (i + 1);
        }
        int64_t pushed = dfs_i64_i64_i64(s, t, INF);
        while (pushed > 0) {
            flow = (flow + pushed);
            pushed = dfs_i64_i64_i64(s, t, INF);
        }
    }
    return flow;
}

void reachable_i64_ptr_i8(int64_t s, int8_t* vis) {
    int64_t i = 0;
    while (i < g_n) {
        vis[i] = 0;
        i = (i + 1);
    }
    int64_t qh = 0;
    int64_t qt = 0;
    vis[s] = 1;
    g_q[qt] = s;
    qt = (qt + 1);
    while (qh < qt) {
        int64_t u = g_q[qh];
        qh = (qh + 1);
        int64_t e = g_head[u];
        while (e != (-1)) {
            int64_t v = e_to[e];
            if ((e_cap[e] > 0 && vis[v] == 0)) {
                vis[v] = 1;
                g_q[qt] = v;
                qt = (qt + 1);
            }
            e = e_next[e];
        }
    }
}

int64_t iroot3_i64(int64_t n) {
    int64_t r = 0;
    bool cont = 1;
    while (cont) {
        int64_t nr = (r + 1);
        __int128 c = ((((__int128)(nr)) * ((__int128)(nr))) * ((__int128)(nr)));
        if (c <= ((__int128)(n))) {
            r = nr;
        } else {
            cont = 0;
        }
    }
    return r;
}

int32_t main(void) {
    int8_t* is_prime = (int8_t*)(((int8_t*)(calloc((N + 1), 1))));
    memset(is_prime, 1, (N + 1));
    is_prime[0] = 0;
    is_prime[1] = 0;
    int64_t p = 2;
    while ((p * p) <= N) {
        if (is_prime[p] != 0) {
            int64_t j = (p * p);
            while (j <= N) {
                is_prime[j] = 0;
                j = (j + p);
            }
        }
        p = (p + 1);
    }
    int64_t pc = 0;
    int64_t i = 2;
    while (i <= N) {
        if (is_prime[i] != 0) {
            pc = (pc + 1);
        }
        i = (i + 1);
    }
    int64_t* primes = (int64_t*)(((int64_t*)(calloc(pc, 8))));
    pc = 0;
    i = 2;
    while (i <= N) {
        if (is_prime[i] != 0) {
            primes[pc] = i;
            pc = (pc + 1);
        }
        i = (i + 1);
    }
    int64_t cbrt = iroot3_i64(N);
    int64_t* left = (int64_t*)(((int64_t*)(calloc(pc, 8))));
    int64_t nL = 0;
    i = 0;
    while (i < pc) {
        int64_t pp = primes[i];
        if (((FLOW_CHECKED_MOD((pp), (10)) == 7 && (pp * 19) <= N) && pp > cbrt)) {
            left[nL] = pp;
            nL = (nL + 1);
        }
        i = (i + 1);
    }
    double ln_total = 0.0;
    i = 0;
    while (i < pc) {
        int64_t pp = primes[i];
        if (FLOW_CHECKED_MOD((pp), (10)) == 3) {
            ln_total = (ln_total + log(((double)(pp))));
        } else {
            if ((FLOW_CHECKED_MOD((pp), (10)) == 7 && pp <= cbrt)) {
                ln_total = (ln_total + log(((double)(pp))));
            }
        }
        i = (i + 1);
    }
    if (nL == 0) {
        printf("%.6f\n", ln_total);
        free(is_prime);
        free(primes);
        free(left);
        return 0;
    }
    int64_t pmin = left[0];
    int64_t* right = (int64_t*)(((int64_t*)(calloc(pc, 8))));
    int64_t nR = 0;
    int64_t rlim = FLOW_CHECKED_DIV((N), (pmin));
    i = 0;
    while (i < pc) {
        int64_t qq = primes[i];
        if ((FLOW_CHECKED_MOD((qq), (10)) == 9 && qq <= rlim)) {
            right[nR] = qq;
            nR = (nR + 1);
        }
        i = (i + 1);
    }
    if (nR == 0) {
        printf("%.6f\n", ln_total);
        free(is_prime);
        free(primes);
        free(left);
        free(right);
        return 0;
    }
    int64_t* pref_len = (int64_t*)(((int64_t*)(calloc(nL, 8))));
    int64_t total_pref = 0;
    i = 0;
    while (i < nL) {
        int64_t lim = FLOW_CHECKED_DIV((N), (left[i]));
        int64_t lo = 0;
        int64_t hi = nR;
        while (lo < hi) {
            int64_t mid = FLOW_CHECKED_DIV(((lo + hi)), (2));
            if (right[mid] <= lim) {
                lo = (mid + 1);
            } else {
                hi = mid;
            }
        }
        pref_len[i] = lo;
        total_pref = (total_pref + lo);
        i = (i + 1);
    }
    int64_t* wL = (int64_t*)(((int64_t*)(calloc(nL, 8))));
    int64_t* wR = (int64_t*)(((int64_t*)(calloc(nR, 8))));
    int64_t total_w = 0;
    i = 0;
    while (i < nL) {
        wL[i] = ((int64_t)(((log(((double)(left[i]))) * ((double)(SCALE))) + 0.5)));
        total_w = (total_w + wL[i]);
        i = (i + 1);
    }
    i = 0;
    while (i < nR) {
        wR[i] = ((int64_t)(((log(((double)(right[i]))) * ((double)(SCALE))) + 0.5)));
        total_w = (total_w + wR[i]);
        i = (i + 1);
    }
    int64_t INF = (total_w + 1);
    int64_t nNodes = ((2 + nL) + nR);
    int64_t S = 0;
    int64_t offL = 1;
    int64_t offR = (1 + nL);
    int64_t T = ((1 + nL) + nR);
    int64_t total_edges = (2 * ((nL + nR) + total_pref));
    g_n = nNodes;
    g_edge_cnt = 0;
    g_head = ((int64_t*)(calloc(nNodes, 8)));
    i = 0;
    while (i < nNodes) {
        g_head[i] = (-1);
        i = (i + 1);
    }
    e_to = ((int64_t*)(calloc(total_edges, 8)));
    e_cap = ((int64_t*)(calloc(total_edges, 8)));
    e_rev = ((int64_t*)(calloc(total_edges, 8)));
    e_next = ((int64_t*)(calloc(total_edges, 8)));
    g_level = ((int64_t*)(calloc(nNodes, 8)));
    g_it = ((int64_t*)(calloc(nNodes, 8)));
    g_q = ((int64_t*)(calloc(nNodes, 8)));
    i = 0;
    while (i < nL) {
        add_edge_i64_i64_i64(S, (offL + i), wL[i]);
        i = (i + 1);
    }
    i = 0;
    while (i < nR) {
        add_edge_i64_i64_i64((offR + i), T, wR[i]);
        i = (i + 1);
    }
    i = 0;
    while (i < nL) {
        int64_t u = (offL + i);
        int64_t k = pref_len[i];
        int64_t j = 0;
        while (j < k) {
            add_edge_i64_i64_i64(u, (offR + j), INF);
            j = (j + 1);
        }
        i = (i + 1);
    }
    max_flow_i64_i64(S, T);
    int8_t* vis = (int8_t*)(((int8_t*)(calloc(nNodes, 1))));
    reachable_i64_ptr_i8(S, vis);
    i = 0;
    while (i < nL) {
        if (vis[(offL + i)] == 0) {
            ln_total = (ln_total + log(((double)(left[i]))));
        }
        i = (i + 1);
    }
    i = 0;
    while (i < nR) {
        if (vis[(offR + i)] != 0) {
            ln_total = (ln_total + log(((double)(right[i]))));
        }
        i = (i + 1);
    }
    printf("%.6f\n", ln_total);
    free(vis);
    free(g_head);
    free(e_to);
    free(e_cap);
    free(e_rev);
    free(e_next);
    free(g_level);
    free(g_it);
    free(g_q);
    free(wL);
    free(wR);
    free(pref_len);
    free(left);
    free(right);
    free(primes);
    free(is_prime);
    return 0;
}

Generated MLIR

module {
  llvm.func @printf(!llvm.ptr, ...) -> i32
  llvm.mlir.global internal constant @str_0("%.6f\n\00") {addr_space = 0 : i32} : !llvm.array<6 x i8>
  func.func private @calloc(i64, i64) -> !llvm.ptr
  func.func private @free(!llvm.ptr) -> ()
  func.func private @memset(!llvm.ptr, i32, i64) -> !llvm.ptr
  func.func private @log(f64) -> f64
  // Constant: N
  llvm.mlir.global internal constant @N(1000000 : i64) : i64
  // Constant: SCALE
  llvm.mlir.global internal constant @SCALE(1000000000000 : i64) : i64
  // Module static: g_head
  llvm.mlir.global internal @g_head() {addr_space = 0 : i32} : !llvm.ptr {
    %0 = llvm.mlir.zero : !llvm.ptr
    llvm.return %0 : !llvm.ptr
  }
  // Module static: e_to
  llvm.mlir.global internal @e_to() {addr_space = 0 : i32} : !llvm.ptr {
    %1 = llvm.mlir.zero : !llvm.ptr
    llvm.return %1 : !llvm.ptr
  }
  // Module static: e_cap
  llvm.mlir.global internal @e_cap() {addr_space = 0 : i32} : !llvm.ptr {
    %2 = llvm.mlir.zero : !llvm.ptr
    llvm.return %2 : !llvm.ptr
  }
  // Module static: e_rev
  llvm.mlir.global internal @e_rev() {addr_space = 0 : i32} : !llvm.ptr {
    %3 = llvm.mlir.zero : !llvm.ptr
    llvm.return %3 : !llvm.ptr
  }
  // Module static: e_next
  llvm.mlir.global internal @e_next() {addr_space = 0 : i32} : !llvm.ptr {
    %4 = llvm.mlir.zero : !llvm.ptr
    llvm.return %4 : !llvm.ptr
  }
  // Module static: g_level
  llvm.mlir.global internal @g_level() {addr_space = 0 : i32} : !llvm.ptr {
    %5 = llvm.mlir.zero : !llvm.ptr
    llvm.return %5 : !llvm.ptr
  }
  // Module static: g_it
  llvm.mlir.global internal @g_it() {addr_space = 0 : i32} : !llvm.ptr {
    %6 = llvm.mlir.zero : !llvm.ptr
    llvm.return %6 : !llvm.ptr
  }
  // Module static: g_q
  llvm.mlir.global internal @g_q() {addr_space = 0 : i32} : !llvm.ptr {
    %7 = llvm.mlir.zero : !llvm.ptr
    llvm.return %7 : !llvm.ptr
  }
  // Module static: g_n
  llvm.mlir.global internal @g_n(0 : i64) : i64
  // Module static: g_edge_cnt
  llvm.mlir.global internal @g_edge_cnt(0 : i64) : i64
  func.func @add_edge(%arg0: i64, %arg1: i64, %arg2: i64) -> () {
    %8 = llvm.mlir.addressof @g_edge_cnt : !llvm.ptr
    %9 = llvm.load %8 : !llvm.ptr -> i64
    %10 = llvm.mlir.addressof @g_edge_cnt : !llvm.ptr
    %11 = llvm.load %10 : !llvm.ptr -> i64
    %12 = arith.constant 1 : i32
    %14 = arith.extsi %12 : i32 to i64
    %13 = arith.addi %11, %14 : i64
    %15 = llvm.mlir.addressof @g_edge_cnt : !llvm.ptr
    llvm.store %13, %15 : i64, !llvm.ptr
    %16 = llvm.mlir.addressof @g_edge_cnt : !llvm.ptr
    %17 = llvm.load %16 : !llvm.ptr -> i64
    %18 = llvm.mlir.addressof @g_edge_cnt : !llvm.ptr
    %19 = llvm.load %18 : !llvm.ptr -> i64
    %20 = arith.constant 1 : i32
    %22 = arith.extsi %20 : i32 to i64
    %21 = arith.addi %19, %22 : i64
    %23 = llvm.mlir.addressof @g_edge_cnt : !llvm.ptr
    llvm.store %21, %23 : i64, !llvm.ptr
    %24 = llvm.mlir.addressof @e_to : !llvm.ptr
    %25 = llvm.load %24 : !llvm.ptr -> !llvm.ptr
    %26 = llvm.getelementptr %25[%9] : (!llvm.ptr, i64) -> !llvm.ptr, i64
    llvm.store %arg1, %26 : i64, !llvm.ptr
    %27 = llvm.mlir.addressof @e_cap : !llvm.ptr
    %28 = llvm.load %27 : !llvm.ptr -> !llvm.ptr
    %29 = llvm.getelementptr %28[%9] : (!llvm.ptr, i64) -> !llvm.ptr, i64
    llvm.store %arg2, %29 : i64, !llvm.ptr
    %30 = llvm.mlir.addressof @e_rev : !llvm.ptr
    %31 = llvm.load %30 : !llvm.ptr -> !llvm.ptr
    %32 = llvm.getelementptr %31[%9] : (!llvm.ptr, i64) -> !llvm.ptr, i64
    llvm.store %17, %32 : i64, !llvm.ptr
    %34 = llvm.mlir.addressof @g_head : !llvm.ptr
    %35 = llvm.load %34 : !llvm.ptr -> !llvm.ptr
    %36 = llvm.getelementptr %35[%arg0] : (!llvm.ptr, i64) -> !llvm.ptr, i64
    %33 = llvm.load %36 : !llvm.ptr -> i64
    %37 = llvm.mlir.addressof @e_next : !llvm.ptr
    %38 = llvm.load %37 : !llvm.ptr -> !llvm.ptr
    %39 = llvm.getelementptr %38[%9] : (!llvm.ptr, i64) -> !llvm.ptr, i64
    llvm.store %33, %39 : i64, !llvm.ptr
    %40 = llvm.mlir.addressof @g_head : !llvm.ptr
    %41 = llvm.load %40 : !llvm.ptr -> !llvm.ptr
    %42 = llvm.getelementptr %41[%arg0] : (!llvm.ptr, i64) -> !llvm.ptr, i64
    llvm.store %9, %42 : i64, !llvm.ptr
    %43 = llvm.mlir.addressof @e_to : !llvm.ptr
    %44 = llvm.load %43 : !llvm.ptr -> !llvm.ptr
    %45 = llvm.getelementptr %44[%17] : (!llvm.ptr, i64) -> !llvm.ptr, i64
    llvm.store %arg0, %45 : i64, !llvm.ptr
    %46 = arith.constant 0 : i32
    %47 = llvm.mlir.addressof @e_cap : !llvm.ptr
    %48 = llvm.load %47 : !llvm.ptr -> !llvm.ptr
    %49 = arith.extsi %46 : i32 to i64
    %50 = llvm.getelementptr %48[%17] : (!llvm.ptr, i64) -> !llvm.ptr, i64
    llvm.store %49, %50 : i64, !llvm.ptr
    %51 = llvm.mlir.addressof @e_rev : !llvm.ptr
    %52 = llvm.load %51 : !llvm.ptr -> !llvm.ptr
    %53 = llvm.getelementptr %52[%17] : (!llvm.ptr, i64) -> !llvm.ptr, i64
    llvm.store %9, %53 : i64, !llvm.ptr
    %55 = llvm.mlir.addressof @g_head : !llvm.ptr
    %56 = llvm.load %55 : !llvm.ptr -> !llvm.ptr
    %57 = llvm.getelementptr %56[%arg1] : (!llvm.ptr, i64) -> !llvm.ptr, i64
    %54 = llvm.load %57 : !llvm.ptr -> i64
    %58 = llvm.mlir.addressof @e_next : !llvm.ptr
    %59 = llvm.load %58 : !llvm.ptr -> !llvm.ptr
    %60 = llvm.getelementptr %59[%17] : (!llvm.ptr, i64) -> !llvm.ptr, i64
    llvm.store %54, %60 : i64, !llvm.ptr
    %61 = llvm.mlir.addressof @g_head : !llvm.ptr
    %62 = llvm.load %61 : !llvm.ptr -> !llvm.ptr
    %63 = llvm.getelementptr %62[%arg1] : (!llvm.ptr, i64) -> !llvm.ptr, i64
    llvm.store %17, %63 : i64, !llvm.ptr
    func.return
  }
  func.func @bfs(%arg0: i64, %arg1: i64) -> i1 {
    %64 = arith.constant 0 : i32
    %65 = arith.extsi %64 : i32 to i64
    %66 = llvm.mlir.constant(1 : i64) : i64
    %67 = llvm.alloca %66 x i64 : (i64) -> !llvm.ptr
    llvm.store %65, %67 : i64, !llvm.ptr
    cf.br ^bb0
    ^bb0:
    %68 = llvm.load %67 : !llvm.ptr -> i64
    %69 = llvm.mlir.addressof @g_n : !llvm.ptr
    %70 = llvm.load %69 : !llvm.ptr -> i64
    %71 = arith.cmpi slt, %68, %70 : i64
    cf.cond_br %71, ^bb1, ^bb2
    ^bb1:
      %72 = arith.constant 1 : i32
      %74 = arith.constant 0 : i32
      %73 = arith.subi %74, %72 : i32
      %75 = llvm.mlir.addressof @g_level : !llvm.ptr
      %76 = llvm.load %75 : !llvm.ptr -> !llvm.ptr
      %77 = llvm.load %67 : !llvm.ptr -> i64
      %78 = arith.extsi %73 : i32 to i64
      %79 = llvm.getelementptr %76[%77] : (!llvm.ptr, i64) -> !llvm.ptr, i64
      llvm.store %78, %79 : i64, !llvm.ptr
      %80 = llvm.load %67 : !llvm.ptr -> i64
      %81 = arith.constant 1 : i32
      %83 = arith.extsi %81 : i32 to i64
      %82 = arith.addi %80, %83 : i64
      llvm.store %82, %67 : i64, !llvm.ptr
      cf.br ^bb0
    ^bb2:
    %84 = arith.constant 0 : i32
    %85 = arith.extsi %84 : i32 to i64
    %86 = llvm.mlir.constant(1 : i64) : i64
    %87 = llvm.alloca %86 x i64 : (i64) -> !llvm.ptr
    llvm.store %85, %87 : i64, !llvm.ptr
    %88 = arith.constant 0 : i32
    %89 = arith.extsi %88 : i32 to i64
    %90 = llvm.mlir.constant(1 : i64) : i64
    %91 = llvm.alloca %90 x i64 : (i64) -> !llvm.ptr
    llvm.store %89, %91 : i64, !llvm.ptr
    %92 = arith.constant 0 : i32
    %93 = llvm.mlir.addressof @g_level : !llvm.ptr
    %94 = llvm.load %93 : !llvm.ptr -> !llvm.ptr
    %95 = arith.extsi %92 : i32 to i64
    %96 = llvm.getelementptr %94[%arg0] : (!llvm.ptr, i64) -> !llvm.ptr, i64
    llvm.store %95, %96 : i64, !llvm.ptr
    %97 = llvm.mlir.addressof @g_q : !llvm.ptr
    %98 = llvm.load %97 : !llvm.ptr -> !llvm.ptr
    %99 = llvm.load %91 : !llvm.ptr -> i64
    %100 = llvm.getelementptr %98[%99] : (!llvm.ptr, i64) -> !llvm.ptr, i64
    llvm.store %arg0, %100 : i64, !llvm.ptr
    %101 = llvm.load %91 : !llvm.ptr -> i64
    %102 = arith.constant 1 : i32
    %104 = arith.extsi %102 : i32 to i64
    %103 = arith.addi %101, %104 : i64
    llvm.store %103, %91 : i64, !llvm.ptr
    cf.br ^bb3
    ^bb3:
    %105 = llvm.load %87 : !llvm.ptr -> i64
    %106 = llvm.load %91 : !llvm.ptr -> i64
    %107 = arith.cmpi slt, %105, %106 : i64
    cf.cond_br %107, ^bb4, ^bb5
    ^bb4:
      %109 = llvm.mlir.addressof @g_q : !llvm.ptr
      %110 = llvm.load %109 : !llvm.ptr -> !llvm.ptr
      %111 = llvm.load %87 : !llvm.ptr -> i64
      %112 = llvm.getelementptr %110[%111] : (!llvm.ptr, i64) -> !llvm.ptr, i64
      %108 = llvm.load %112 : !llvm.ptr -> i64
      %113 = llvm.load %87 : !llvm.ptr -> i64
      %114 = arith.constant 1 : i32
      %116 = arith.extsi %114 : i32 to i64
      %115 = arith.addi %113, %116 : i64
      llvm.store %115, %87 : i64, !llvm.ptr
      %118 = llvm.mlir.addressof @g_head : !llvm.ptr
      %119 = llvm.load %118 : !llvm.ptr -> !llvm.ptr
      %120 = llvm.getelementptr %119[%108] : (!llvm.ptr, i64) -> !llvm.ptr, i64
      %117 = llvm.load %120 : !llvm.ptr -> i64
      %121 = llvm.mlir.constant(1 : i64) : i64
      %122 = llvm.alloca %121 x i64 : (i64) -> !llvm.ptr
      llvm.store %117, %122 : i64, !llvm.ptr
      cf.br ^bb6
      ^bb6:
      %123 = llvm.load %122 : !llvm.ptr -> i64
      %124 = arith.constant 1 : i32
      %126 = arith.constant 0 : i32
      %125 = arith.subi %126, %124 : i32
      %128 = arith.extsi %125 : i32 to i64
      %127 = arith.cmpi ne, %123, %128 : i64
      cf.cond_br %127, ^bb7, ^bb8
      ^bb7:
        %130 = llvm.mlir.addressof @e_to : !llvm.ptr
        %131 = llvm.load %130 : !llvm.ptr -> !llvm.ptr
        %132 = llvm.load %122 : !llvm.ptr -> i64
        %133 = llvm.getelementptr %131[%132] : (!llvm.ptr, i64) -> !llvm.ptr, i64
        %129 = llvm.load %133 : !llvm.ptr -> i64
        %135 = llvm.mlir.addressof @e_cap : !llvm.ptr
        %136 = llvm.load %135 : !llvm.ptr -> !llvm.ptr
        %137 = llvm.load %122 : !llvm.ptr -> i64
        %138 = llvm.getelementptr %136[%137] : (!llvm.ptr, i64) -> !llvm.ptr, i64
        %134 = llvm.load %138 : !llvm.ptr -> i64
        %139 = arith.constant 0 : i32
        %141 = arith.extsi %139 : i32 to i64
        %140 = arith.cmpi sgt, %134, %141 : i64
        %142 = scf.if %140 -> (i1) {
          %144 = llvm.mlir.addressof @g_level : !llvm.ptr
          %145 = llvm.load %144 : !llvm.ptr -> !llvm.ptr
          %146 = llvm.getelementptr %145[%129] : (!llvm.ptr, i64) -> !llvm.ptr, i64
          %143 = llvm.load %146 : !llvm.ptr -> i64
          %147 = arith.constant 0 : i32
          %149 = arith.extsi %147 : i32 to i64
          %148 = arith.cmpi slt, %143, %149 : i64
          scf.yield %148 : i1
        } else {
          %150 = arith.constant false
          scf.yield %150 : i1
        }
        cf.cond_br %142, ^bb9, ^bb10
        ^bb9:
          %152 = llvm.mlir.addressof @g_level : !llvm.ptr
          %153 = llvm.load %152 : !llvm.ptr -> !llvm.ptr
          %154 = llvm.getelementptr %153[%108] : (!llvm.ptr, i64) -> !llvm.ptr, i64
          %151 = llvm.load %154 : !llvm.ptr -> i64
          %155 = arith.constant 1 : i32
          %157 = arith.extsi %155 : i32 to i64
          %156 = arith.addi %151, %157 : i64
          %158 = llvm.mlir.addressof @g_level : !llvm.ptr
          %159 = llvm.load %158 : !llvm.ptr -> !llvm.ptr
          %160 = llvm.getelementptr %159[%129] : (!llvm.ptr, i64) -> !llvm.ptr, i64
          llvm.store %156, %160 : i64, !llvm.ptr
          %161 = llvm.mlir.addressof @g_q : !llvm.ptr
          %162 = llvm.load %161 : !llvm.ptr -> !llvm.ptr
          %163 = llvm.load %91 : !llvm.ptr -> i64
          %164 = llvm.getelementptr %162[%163] : (!llvm.ptr, i64) -> !llvm.ptr, i64
          llvm.store %129, %164 : i64, !llvm.ptr
          %165 = llvm.load %91 : !llvm.ptr -> i64
          %166 = arith.constant 1 : i32
          %168 = arith.extsi %166 : i32 to i64
          %167 = arith.addi %165, %168 : i64
          llvm.store %167, %91 : i64, !llvm.ptr
          cf.br ^bb11
        ^bb10:
          cf.br ^bb11
        ^bb11:
        %170 = llvm.mlir.addressof @e_next : !llvm.ptr
        %171 = llvm.load %170 : !llvm.ptr -> !llvm.ptr
        %172 = llvm.load %122 : !llvm.ptr -> i64
        %173 = llvm.getelementptr %171[%172] : (!llvm.ptr, i64) -> !llvm.ptr, i64
        %169 = llvm.load %173 : !llvm.ptr -> i64
        llvm.store %169, %122 : i64, !llvm.ptr
        cf.br ^bb6
      ^bb8:
      cf.br ^bb3
    ^bb5:
    %175 = llvm.mlir.addressof @g_level : !llvm.ptr
    %176 = llvm.load %175 : !llvm.ptr -> !llvm.ptr
    %177 = llvm.getelementptr %176[%arg1] : (!llvm.ptr, i64) -> !llvm.ptr, i64
    %174 = llvm.load %177 : !llvm.ptr -> i64
    %178 = arith.constant 0 : i32
    %180 = arith.extsi %178 : i32 to i64
    %179 = arith.cmpi sge, %174, %180 : i64
    func.return %179 : i1
  }
  func.func @dfs(%arg0: i64, %arg1: i64, %arg2: i64) -> i64 {
    %181 = arith.cmpi eq, %arg0, %arg1 : i64
    cf.cond_br %181, ^bb12, ^bb13
    ^bb12:
      func.return %arg2 : i64
    ^bb13:
      cf.br ^bb14
    ^bb14:
    %183 = llvm.mlir.addressof @g_it : !llvm.ptr
    %184 = llvm.load %183 : !llvm.ptr -> !llvm.ptr
    %185 = llvm.getelementptr %184[%arg0] : (!llvm.ptr, i64) -> !llvm.ptr, i64
    %182 = llvm.load %185 : !llvm.ptr -> i64
    %186 = llvm.mlir.constant(1 : i64) : i64
    %187 = llvm.alloca %186 x i64 : (i64) -> !llvm.ptr
    llvm.store %182, %187 : i64, !llvm.ptr
    cf.br ^bb15
    ^bb15:
    %188 = llvm.load %187 : !llvm.ptr -> i64
    %189 = arith.constant 1 : i32
    %191 = arith.constant 0 : i32
    %190 = arith.subi %191, %189 : i32
    %193 = arith.extsi %190 : i32 to i64
    %192 = arith.cmpi ne, %188, %193 : i64
    cf.cond_br %192, ^bb16, ^bb17
    ^bb16:
      %195 = llvm.mlir.addressof @e_to : !llvm.ptr
      %196 = llvm.load %195 : !llvm.ptr -> !llvm.ptr
      %197 = llvm.load %187 : !llvm.ptr -> i64
      %198 = llvm.getelementptr %196[%197] : (!llvm.ptr, i64) -> !llvm.ptr, i64
      %194 = llvm.load %198 : !llvm.ptr -> i64
      %200 = llvm.mlir.addressof @e_cap : !llvm.ptr
      %201 = llvm.load %200 : !llvm.ptr -> !llvm.ptr
      %202 = llvm.load %187 : !llvm.ptr -> i64
      %203 = llvm.getelementptr %201[%202] : (!llvm.ptr, i64) -> !llvm.ptr, i64
      %199 = llvm.load %203 : !llvm.ptr -> i64
      %204 = arith.constant 0 : i32
      %206 = arith.extsi %204 : i32 to i64
      %205 = arith.cmpi sgt, %199, %206 : i64
      %207 = scf.if %205 -> (i1) {
        %209 = llvm.mlir.addressof @g_level : !llvm.ptr
        %210 = llvm.load %209 : !llvm.ptr -> !llvm.ptr
        %211 = llvm.getelementptr %210[%194] : (!llvm.ptr, i64) -> !llvm.ptr, i64
        %208 = llvm.load %211 : !llvm.ptr -> i64
        %213 = llvm.mlir.addressof @g_level : !llvm.ptr
        %214 = llvm.load %213 : !llvm.ptr -> !llvm.ptr
        %215 = llvm.getelementptr %214[%arg0] : (!llvm.ptr, i64) -> !llvm.ptr, i64
        %212 = llvm.load %215 : !llvm.ptr -> i64
        %216 = arith.constant 1 : i32
        %218 = arith.extsi %216 : i32 to i64
        %217 = arith.addi %212, %218 : i64
        %219 = arith.cmpi eq, %208, %217 : i64
        scf.yield %219 : i1
      } else {
        %220 = arith.constant false
        scf.yield %220 : i1
      }
      cf.cond_br %207, ^bb18, ^bb19
      ^bb18:
        %221 = llvm.mlir.constant(1 : i64) : i64
        %222 = llvm.alloca %221 x i64 : (i64) -> !llvm.ptr
        llvm.store %arg2, %222 : i64, !llvm.ptr
        %224 = llvm.mlir.addressof @e_cap : !llvm.ptr
        %225 = llvm.load %224 : !llvm.ptr -> !llvm.ptr
        %226 = llvm.load %187 : !llvm.ptr -> i64
        %227 = llvm.getelementptr %225[%226] : (!llvm.ptr, i64) -> !llvm.ptr, i64
        %223 = llvm.load %227 : !llvm.ptr -> i64
        %228 = llvm.load %222 : !llvm.ptr -> i64
        %229 = arith.cmpi slt, %223, %228 : i64
        cf.cond_br %229, ^bb21, ^bb22
        ^bb21:
          %231 = llvm.mlir.addressof @e_cap : !llvm.ptr
          %232 = llvm.load %231 : !llvm.ptr -> !llvm.ptr
          %233 = llvm.load %187 : !llvm.ptr -> i64
          %234 = llvm.getelementptr %232[%233] : (!llvm.ptr, i64) -> !llvm.ptr, i64
          %230 = llvm.load %234 : !llvm.ptr -> i64
          llvm.store %230, %222 : i64, !llvm.ptr
          cf.br ^bb23
        ^bb22:
          cf.br ^bb23
        ^bb23:
        %236 = llvm.load %222 : !llvm.ptr -> i64
        %235 = func.call @dfs(%194, %arg1, %236) : (i64, i64, i64) -> i64
        %237 = arith.constant 0 : i32
        %239 = arith.extsi %237 : i32 to i64
        %238 = arith.cmpi sgt, %235, %239 : i64
        cf.cond_br %238, ^bb24, ^bb25
        ^bb24:
          %241 = llvm.mlir.addressof @e_cap : !llvm.ptr
          %242 = llvm.load %241 : !llvm.ptr -> !llvm.ptr
          %243 = llvm.load %187 : !llvm.ptr -> i64
          %244 = llvm.getelementptr %242[%243] : (!llvm.ptr, i64) -> !llvm.ptr, i64
          %240 = llvm.load %244 : !llvm.ptr -> i64
          %245 = arith.subi %240, %235 : i64
          %246 = llvm.mlir.addressof @e_cap : !llvm.ptr
          %247 = llvm.load %246 : !llvm.ptr -> !llvm.ptr
          %248 = llvm.load %187 : !llvm.ptr -> i64
          %249 = llvm.getelementptr %247[%248] : (!llvm.ptr, i64) -> !llvm.ptr, i64
          llvm.store %245, %249 : i64, !llvm.ptr
          %251 = llvm.mlir.addressof @e_cap : !llvm.ptr
          %252 = llvm.load %251 : !llvm.ptr -> !llvm.ptr
          %254 = llvm.mlir.addressof @e_rev : !llvm.ptr
          %255 = llvm.load %254 : !llvm.ptr -> !llvm.ptr
          %256 = llvm.load %187 : !llvm.ptr -> i64
          %257 = llvm.getelementptr %255[%256] : (!llvm.ptr, i64) -> !llvm.ptr, i64
          %253 = llvm.load %257 : !llvm.ptr -> i64
          %258 = llvm.getelementptr %252[%253] : (!llvm.ptr, i64) -> !llvm.ptr, i64
          %250 = llvm.load %258 : !llvm.ptr -> i64
          %259 = arith.addi %250, %235 : i64
          %260 = llvm.mlir.addressof @e_cap : !llvm.ptr
          %261 = llvm.load %260 : !llvm.ptr -> !llvm.ptr
          %263 = llvm.mlir.addressof @e_rev : !llvm.ptr
          %264 = llvm.load %263 : !llvm.ptr -> !llvm.ptr
          %265 = llvm.load %187 : !llvm.ptr -> i64
          %266 = llvm.getelementptr %264[%265] : (!llvm.ptr, i64) -> !llvm.ptr, i64
          %262 = llvm.load %266 : !llvm.ptr -> i64
          %267 = llvm.getelementptr %261[%262] : (!llvm.ptr, i64) -> !llvm.ptr, i64
          llvm.store %259, %267 : i64, !llvm.ptr
          %268 = llvm.load %187 : !llvm.ptr -> i64
          %269 = llvm.mlir.addressof @g_it : !llvm.ptr
          %270 = llvm.load %269 : !llvm.ptr -> !llvm.ptr
          %271 = llvm.getelementptr %270[%arg0] : (!llvm.ptr, i64) -> !llvm.ptr, i64
          llvm.store %268, %271 : i64, !llvm.ptr
          func.return %235 : i64
        ^bb25:
          cf.br ^bb26
        ^bb26:
        cf.br ^bb20
      ^bb19:
        cf.br ^bb20
      ^bb20:
      %273 = llvm.mlir.addressof @e_next : !llvm.ptr
      %274 = llvm.load %273 : !llvm.ptr -> !llvm.ptr
      %275 = llvm.load %187 : !llvm.ptr -> i64
      %276 = llvm.getelementptr %274[%275] : (!llvm.ptr, i64) -> !llvm.ptr, i64
      %272 = llvm.load %276 : !llvm.ptr -> i64
      llvm.store %272, %187 : i64, !llvm.ptr
      %277 = llvm.load %187 : !llvm.ptr -> i64
      %278 = llvm.mlir.addressof @g_it : !llvm.ptr
      %279 = llvm.load %278 : !llvm.ptr -> !llvm.ptr
      %280 = llvm.getelementptr %279[%arg0] : (!llvm.ptr, i64) -> !llvm.ptr, i64
      llvm.store %277, %280 : i64, !llvm.ptr
      cf.br ^bb15
    ^bb17:
    %281 = arith.constant 0 : i32
    %282 = arith.extsi %281 : i32 to i64
    func.return %282 : i64
  }
  func.func @max_flow(%arg0: i64, %arg1: i64) -> i64 {
    %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
    %287 = arith.constant 999999995705032704 : i32
    %288 = arith.extsi %287 : i32 to i64
    cf.br ^bb27
    ^bb27:
    %289 = func.call @bfs(%arg0, %arg1) : (i64, i64) -> i1
    cf.cond_br %289, ^bb28, ^bb29
    ^bb28:
      %290 = arith.constant 0 : i32
      %291 = arith.extsi %290 : i32 to i64
      %292 = llvm.mlir.constant(1 : i64) : i64
      %293 = llvm.alloca %292 x i64 : (i64) -> !llvm.ptr
      llvm.store %291, %293 : i64, !llvm.ptr
      cf.br ^bb30
      ^bb30:
      %294 = llvm.load %293 : !llvm.ptr -> i64
      %295 = llvm.mlir.addressof @g_n : !llvm.ptr
      %296 = llvm.load %295 : !llvm.ptr -> i64
      %297 = arith.cmpi slt, %294, %296 : i64
      cf.cond_br %297, ^bb31, ^bb32
      ^bb31:
        %299 = llvm.mlir.addressof @g_head : !llvm.ptr
        %300 = llvm.load %299 : !llvm.ptr -> !llvm.ptr
        %301 = llvm.load %293 : !llvm.ptr -> i64
        %302 = llvm.getelementptr %300[%301] : (!llvm.ptr, i64) -> !llvm.ptr, i64
        %298 = llvm.load %302 : !llvm.ptr -> i64
        %303 = llvm.mlir.addressof @g_it : !llvm.ptr
        %304 = llvm.load %303 : !llvm.ptr -> !llvm.ptr
        %305 = llvm.load %293 : !llvm.ptr -> i64
        %306 = llvm.getelementptr %304[%305] : (!llvm.ptr, i64) -> !llvm.ptr, i64
        llvm.store %298, %306 : i64, !llvm.ptr
        %307 = llvm.load %293 : !llvm.ptr -> i64
        %308 = arith.constant 1 : i32
        %310 = arith.extsi %308 : i32 to i64
        %309 = arith.addi %307, %310 : i64
        llvm.store %309, %293 : i64, !llvm.ptr
        cf.br ^bb30
      ^bb32:
      %311 = func.call @dfs(%arg0, %arg1, %288) : (i64, i64, i64) -> i64
      %312 = llvm.mlir.constant(1 : i64) : i64
      %313 = llvm.alloca %312 x i64 : (i64) -> !llvm.ptr
      llvm.store %311, %313 : i64, !llvm.ptr
      cf.br ^bb33
      ^bb33:
      %314 = llvm.load %313 : !llvm.ptr -> i64
      %315 = arith.constant 0 : i32
      %317 = arith.extsi %315 : i32 to i64
      %316 = arith.cmpi sgt, %314, %317 : i64
      cf.cond_br %316, ^bb34, ^bb35
      ^bb34:
        %318 = llvm.load %286 : !llvm.ptr -> i64
        %319 = llvm.load %313 : !llvm.ptr -> i64
        %320 = arith.addi %318, %319 : i64
        llvm.store %320, %286 : i64, !llvm.ptr
        %321 = func.call @dfs(%arg0, %arg1, %288) : (i64, i64, i64) -> i64
        llvm.store %321, %313 : i64, !llvm.ptr
        cf.br ^bb33
      ^bb35:
      cf.br ^bb27
    ^bb29:
    %322 = llvm.load %286 : !llvm.ptr -> i64
    func.return %322 : i64
  }
  func.func @reachable(%arg0: i64, %arg1: !llvm.ptr) -> () {
    %323 = arith.constant 0 : i32
    %324 = arith.extsi %323 : i32 to i64
    %325 = llvm.mlir.constant(1 : i64) : i64
    %326 = llvm.alloca %325 x i64 : (i64) -> !llvm.ptr
    llvm.store %324, %326 : i64, !llvm.ptr
    cf.br ^bb36
    ^bb36:
    %327 = llvm.load %326 : !llvm.ptr -> i64
    %328 = llvm.mlir.addressof @g_n : !llvm.ptr
    %329 = llvm.load %328 : !llvm.ptr -> i64
    %330 = arith.cmpi slt, %327, %329 : i64
    cf.cond_br %330, ^bb37, ^bb38
    ^bb37:
      %331 = arith.constant 0 : i32
      %332 = llvm.load %326 : !llvm.ptr -> i64
      %333 = arith.trunci %331 : i32 to i8
      %334 = llvm.getelementptr %arg1[%332] : (!llvm.ptr, i64) -> !llvm.ptr, i8
      llvm.store %333, %334 : i8, !llvm.ptr
      %335 = llvm.load %326 : !llvm.ptr -> i64
      %336 = arith.constant 1 : i32
      %338 = arith.extsi %336 : i32 to i64
      %337 = arith.addi %335, %338 : i64
      llvm.store %337, %326 : i64, !llvm.ptr
      cf.br ^bb36
    ^bb38:
    %339 = arith.constant 0 : i32
    %340 = arith.extsi %339 : i32 to i64
    %341 = llvm.mlir.constant(1 : i64) : i64
    %342 = llvm.alloca %341 x i64 : (i64) -> !llvm.ptr
    llvm.store %340, %342 : i64, !llvm.ptr
    %343 = arith.constant 0 : i32
    %344 = arith.extsi %343 : i32 to i64
    %345 = llvm.mlir.constant(1 : i64) : i64
    %346 = llvm.alloca %345 x i64 : (i64) -> !llvm.ptr
    llvm.store %344, %346 : i64, !llvm.ptr
    %347 = arith.constant 1 : i32
    %348 = arith.trunci %347 : i32 to i8
    %349 = llvm.getelementptr %arg1[%arg0] : (!llvm.ptr, i64) -> !llvm.ptr, i8
    llvm.store %348, %349 : i8, !llvm.ptr
    %350 = llvm.mlir.addressof @g_q : !llvm.ptr
    %351 = llvm.load %350 : !llvm.ptr -> !llvm.ptr
    %352 = llvm.load %346 : !llvm.ptr -> i64
    %353 = llvm.getelementptr %351[%352] : (!llvm.ptr, i64) -> !llvm.ptr, i64
    llvm.store %arg0, %353 : i64, !llvm.ptr
    %354 = llvm.load %346 : !llvm.ptr -> i64
    %355 = arith.constant 1 : i32
    %357 = arith.extsi %355 : i32 to i64
    %356 = arith.addi %354, %357 : i64
    llvm.store %356, %346 : i64, !llvm.ptr
    cf.br ^bb39
    ^bb39:
    %358 = llvm.load %342 : !llvm.ptr -> i64
    %359 = llvm.load %346 : !llvm.ptr -> i64
    %360 = arith.cmpi slt, %358, %359 : i64
    cf.cond_br %360, ^bb40, ^bb41
    ^bb40:
      %362 = llvm.mlir.addressof @g_q : !llvm.ptr
      %363 = llvm.load %362 : !llvm.ptr -> !llvm.ptr
      %364 = llvm.load %342 : !llvm.ptr -> i64
      %365 = llvm.getelementptr %363[%364] : (!llvm.ptr, i64) -> !llvm.ptr, i64
      %361 = llvm.load %365 : !llvm.ptr -> i64
      %366 = llvm.load %342 : !llvm.ptr -> i64
      %367 = arith.constant 1 : i32
      %369 = arith.extsi %367 : i32 to i64
      %368 = arith.addi %366, %369 : i64
      llvm.store %368, %342 : i64, !llvm.ptr
      %371 = llvm.mlir.addressof @g_head : !llvm.ptr
      %372 = llvm.load %371 : !llvm.ptr -> !llvm.ptr
      %373 = llvm.getelementptr %372[%361] : (!llvm.ptr, i64) -> !llvm.ptr, i64
      %370 = llvm.load %373 : !llvm.ptr -> i64
      %374 = llvm.mlir.constant(1 : i64) : i64
      %375 = llvm.alloca %374 x i64 : (i64) -> !llvm.ptr
      llvm.store %370, %375 : i64, !llvm.ptr
      cf.br ^bb42
      ^bb42:
      %376 = llvm.load %375 : !llvm.ptr -> i64
      %377 = arith.constant 1 : i32
      %379 = arith.constant 0 : i32
      %378 = arith.subi %379, %377 : i32
      %381 = arith.extsi %378 : i32 to i64
      %380 = arith.cmpi ne, %376, %381 : i64
      cf.cond_br %380, ^bb43, ^bb44
      ^bb43:
        %383 = llvm.mlir.addressof @e_to : !llvm.ptr
        %384 = llvm.load %383 : !llvm.ptr -> !llvm.ptr
        %385 = llvm.load %375 : !llvm.ptr -> i64
        %386 = llvm.getelementptr %384[%385] : (!llvm.ptr, i64) -> !llvm.ptr, i64
        %382 = llvm.load %386 : !llvm.ptr -> i64
        %388 = llvm.mlir.addressof @e_cap : !llvm.ptr
        %389 = llvm.load %388 : !llvm.ptr -> !llvm.ptr
        %390 = llvm.load %375 : !llvm.ptr -> i64
        %391 = llvm.getelementptr %389[%390] : (!llvm.ptr, i64) -> !llvm.ptr, i64
        %387 = llvm.load %391 : !llvm.ptr -> i64
        %392 = arith.constant 0 : i32
        %394 = arith.extsi %392 : i32 to i64
        %393 = arith.cmpi sgt, %387, %394 : i64
        %395 = scf.if %393 -> (i1) {
          %397 = llvm.getelementptr %arg1[%382] : (!llvm.ptr, i64) -> !llvm.ptr, i8
          %396 = llvm.load %397 : !llvm.ptr -> i8
          %398 = arith.constant 0 : i32
          %400 = arith.extsi %396 : i8 to i32
          %399 = arith.cmpi eq, %400, %398 : i32
          scf.yield %399 : i1
        } else {
          %401 = arith.constant false
          scf.yield %401 : i1
        }
        cf.cond_br %395, ^bb45, ^bb46
        ^bb45:
          %402 = arith.constant 1 : i32
          %403 = arith.trunci %402 : i32 to i8
          %404 = llvm.getelementptr %arg1[%382] : (!llvm.ptr, i64) -> !llvm.ptr, i8
          llvm.store %403, %404 : i8, !llvm.ptr
          %405 = llvm.mlir.addressof @g_q : !llvm.ptr
          %406 = llvm.load %405 : !llvm.ptr -> !llvm.ptr
          %407 = llvm.load %346 : !llvm.ptr -> i64
          %408 = llvm.getelementptr %406[%407] : (!llvm.ptr, i64) -> !llvm.ptr, i64
          llvm.store %382, %408 : i64, !llvm.ptr
          %409 = llvm.load %346 : !llvm.ptr -> i64
          %410 = arith.constant 1 : i32
          %412 = arith.extsi %410 : i32 to i64
          %411 = arith.addi %409, %412 : i64
          llvm.store %411, %346 : i64, !llvm.ptr
          cf.br ^bb47
        ^bb46:
          cf.br ^bb47
        ^bb47:
        %414 = llvm.mlir.addressof @e_next : !llvm.ptr
        %415 = llvm.load %414 : !llvm.ptr -> !llvm.ptr
        %416 = llvm.load %375 : !llvm.ptr -> i64
        %417 = llvm.getelementptr %415[%416] : (!llvm.ptr, i64) -> !llvm.ptr, i64
        %413 = llvm.load %417 : !llvm.ptr -> i64
        llvm.store %413, %375 : i64, !llvm.ptr
        cf.br ^bb42
      ^bb44:
      cf.br ^bb39
    ^bb41:
    func.return
  }
  func.func @iroot3(%arg0: i64) -> i64 {
    %418 = arith.constant 0 : i32
    %419 = arith.extsi %418 : i32 to i64
    %420 = llvm.mlir.constant(1 : i64) : i64
    %421 = llvm.alloca %420 x i64 : (i64) -> !llvm.ptr
    llvm.store %419, %421 : i64, !llvm.ptr
    %422 = arith.constant 1 : i1
    %423 = llvm.mlir.constant(1 : i64) : i64
    %424 = llvm.alloca %423 x i1 : (i64) -> !llvm.ptr
    llvm.store %422, %424 : i1, !llvm.ptr
    cf.br ^bb48
    ^bb48:
    %425 = llvm.load %424 : !llvm.ptr -> i1
    cf.cond_br %425, ^bb49, ^bb50
    ^bb49:
      %426 = llvm.load %421 : !llvm.ptr -> i64
      %427 = arith.constant 1 : i32
      %429 = arith.extsi %427 : i32 to i64
      %428 = arith.addi %426, %429 : i64
      %430 = arith.extsi %428 : i64 to i128
      %431 = arith.extsi %428 : i64 to i128
      %433 = arith.trunci %430 : i128 to i64
      %434 = arith.trunci %431 : i128 to i64
      %432 = arith.muli %433, %434 : i64
      %435 = arith.extsi %428 : i64 to i128
      %437 = arith.trunci %435 : i128 to i64
      %436 = arith.muli %432, %437 : i64
      %438 = arith.extsi %436 : i64 to i128
      %439 = arith.extsi %arg0 : i64 to i128
      %441 = arith.trunci %438 : i128 to i64
      %442 = arith.trunci %439 : i128 to i64
      %440 = arith.cmpi sle, %441, %442 : i64
      cf.cond_br %440, ^bb51, ^bb52
      ^bb51:
        llvm.store %428, %421 : i64, !llvm.ptr
        cf.br ^bb53
      ^bb52:
        %443 = arith.constant 0 : i1
        llvm.store %443, %424 : i1, !llvm.ptr
        cf.br ^bb53
      ^bb53:
      cf.br ^bb48
    ^bb50:
    %444 = llvm.load %421 : !llvm.ptr -> i64
    func.return %444 : i64
  }
  func.func @main() -> i32 {
    %446 = llvm.mlir.addressof @N : !llvm.ptr
    %447 = llvm.load %446 : !llvm.ptr -> i64
    %448 = arith.constant 1 : i32
    %450 = arith.extsi %448 : i32 to i64
    %449 = arith.addi %447, %450 : i64
    %451 = arith.constant 1 : i32
    %452 = arith.extsi %451 : i32 to i64
    %445 = func.call @calloc(%449, %452) : (i64, i64) -> !llvm.ptr
    %454 = arith.constant 1 : i32
    %455 = llvm.mlir.addressof @N : !llvm.ptr
    %456 = llvm.load %455 : !llvm.ptr -> i64
    %457 = arith.constant 1 : i32
    %459 = arith.extsi %457 : i32 to i64
    %458 = arith.addi %456, %459 : i64
    %453 = func.call @memset(%445, %454, %458) : (!llvm.ptr, i32, i64) -> !llvm.ptr
    %460 = arith.constant 0 : i32
    %461 = arith.constant 0 : i32
    %462 = arith.trunci %460 : i32 to i8
    %463 = arith.extsi %461 : i32 to i64
    %464 = llvm.getelementptr %445[%463] : (!llvm.ptr, i64) -> !llvm.ptr, i8
    llvm.store %462, %464 : i8, !llvm.ptr
    %465 = arith.constant 0 : i32
    %466 = arith.constant 1 : i32
    %467 = arith.trunci %465 : i32 to i8
    %468 = arith.extsi %466 : i32 to i64
    %469 = llvm.getelementptr %445[%468] : (!llvm.ptr, i64) -> !llvm.ptr, i8
    llvm.store %467, %469 : i8, !llvm.ptr
    %470 = arith.constant 2 : i32
    %471 = arith.extsi %470 : i32 to i64
    %472 = llvm.mlir.constant(1 : i64) : i64
    %473 = llvm.alloca %472 x i64 : (i64) -> !llvm.ptr
    llvm.store %471, %473 : i64, !llvm.ptr
    cf.br ^bb54
    ^bb54:
    %474 = llvm.load %473 : !llvm.ptr -> i64
    %475 = llvm.load %473 : !llvm.ptr -> i64
    %476 = arith.muli %474, %475 : i64
    %477 = llvm.mlir.addressof @N : !llvm.ptr
    %478 = llvm.load %477 : !llvm.ptr -> i64
    %479 = arith.cmpi sle, %476, %478 : i64
    cf.cond_br %479, ^bb55, ^bb56
    ^bb55:
      %481 = llvm.load %473 : !llvm.ptr -> i64
      %482 = llvm.getelementptr %445[%481] : (!llvm.ptr, i64) -> !llvm.ptr, i8
      %480 = llvm.load %482 : !llvm.ptr -> i8
      %483 = arith.constant 0 : i32
      %485 = arith.extsi %480 : i8 to i32
      %484 = arith.cmpi ne, %485, %483 : i32
      cf.cond_br %484, ^bb57, ^bb58
      ^bb57:
        %486 = llvm.load %473 : !llvm.ptr -> i64
        %487 = llvm.load %473 : !llvm.ptr -> i64
        %488 = arith.muli %486, %487 : i64
        %489 = llvm.mlir.constant(1 : i64) : i64
        %490 = llvm.alloca %489 x i64 : (i64) -> !llvm.ptr
        llvm.store %488, %490 : i64, !llvm.ptr
        cf.br ^bb60
        ^bb60:
        %491 = llvm.load %490 : !llvm.ptr -> i64
        %492 = llvm.mlir.addressof @N : !llvm.ptr
        %493 = llvm.load %492 : !llvm.ptr -> i64
        %494 = arith.cmpi sle, %491, %493 : i64
        cf.cond_br %494, ^bb61, ^bb62
        ^bb61:
          %495 = arith.constant 0 : i32
          %496 = llvm.load %490 : !llvm.ptr -> i64
          %497 = arith.trunci %495 : i32 to i8
          %498 = llvm.getelementptr %445[%496] : (!llvm.ptr, i64) -> !llvm.ptr, i8
          llvm.store %497, %498 : i8, !llvm.ptr
          %499 = llvm.load %490 : !llvm.ptr -> i64
          %500 = llvm.load %473 : !llvm.ptr -> i64
          %501 = arith.addi %499, %500 : i64
          llvm.store %501, %490 : i64, !llvm.ptr
          cf.br ^bb60
        ^bb62:
        cf.br ^bb59
      ^bb58:
        cf.br ^bb59
      ^bb59:
      %502 = llvm.load %473 : !llvm.ptr -> i64
      %503 = arith.constant 1 : i32
      %505 = arith.extsi %503 : i32 to i64
      %504 = arith.addi %502, %505 : i64
      llvm.store %504, %473 : i64, !llvm.ptr
      cf.br ^bb54
    ^bb56:
    %506 = arith.constant 0 : i32
    %507 = arith.extsi %506 : i32 to i64
    %508 = llvm.mlir.constant(1 : i64) : i64
    %509 = llvm.alloca %508 x i64 : (i64) -> !llvm.ptr
    llvm.store %507, %509 : i64, !llvm.ptr
    %510 = arith.constant 2 : i32
    %511 = arith.extsi %510 : i32 to i64
    %512 = llvm.mlir.constant(1 : i64) : i64
    %513 = llvm.alloca %512 x i64 : (i64) -> !llvm.ptr
    llvm.store %511, %513 : i64, !llvm.ptr
    cf.br ^bb63
    ^bb63:
    %514 = llvm.load %513 : !llvm.ptr -> i64
    %515 = llvm.mlir.addressof @N : !llvm.ptr
    %516 = llvm.load %515 : !llvm.ptr -> i64
    %517 = arith.cmpi sle, %514, %516 : i64
    cf.cond_br %517, ^bb64, ^bb65
    ^bb64:
      %519 = llvm.load %513 : !llvm.ptr -> i64
      %520 = llvm.getelementptr %445[%519] : (!llvm.ptr, i64) -> !llvm.ptr, i8
      %518 = llvm.load %520 : !llvm.ptr -> i8
      %521 = arith.constant 0 : i32
      %523 = arith.extsi %518 : i8 to i32
      %522 = arith.cmpi ne, %523, %521 : i32
      cf.cond_br %522, ^bb66, ^bb67
      ^bb66:
        %524 = llvm.load %509 : !llvm.ptr -> i64
        %525 = arith.constant 1 : i32
        %527 = arith.extsi %525 : i32 to i64
        %526 = arith.addi %524, %527 : i64
        llvm.store %526, %509 : i64, !llvm.ptr
        cf.br ^bb68
      ^bb67:
        cf.br ^bb68
      ^bb68:
      %528 = llvm.load %513 : !llvm.ptr -> i64
      %529 = arith.constant 1 : i32
      %531 = arith.extsi %529 : i32 to i64
      %530 = arith.addi %528, %531 : i64
      llvm.store %530, %513 : i64, !llvm.ptr
      cf.br ^bb63
    ^bb65:
    %533 = llvm.load %509 : !llvm.ptr -> i64
    %534 = arith.constant 8 : i32
    %535 = arith.extsi %534 : i32 to i64
    %532 = func.call @calloc(%533, %535) : (i64, i64) -> !llvm.ptr
    %536 = arith.constant 0 : i32
    %537 = arith.extsi %536 : i32 to i64
    llvm.store %537, %509 : i64, !llvm.ptr
    %538 = arith.constant 2 : i32
    %539 = arith.extsi %538 : i32 to i64
    llvm.store %539, %513 : i64, !llvm.ptr
    cf.br ^bb69
    ^bb69:
    %540 = llvm.load %513 : !llvm.ptr -> i64
    %541 = llvm.mlir.addressof @N : !llvm.ptr
    %542 = llvm.load %541 : !llvm.ptr -> i64
    %543 = arith.cmpi sle, %540, %542 : i64
    cf.cond_br %543, ^bb70, ^bb71
    ^bb70:
      %545 = llvm.load %513 : !llvm.ptr -> i64
      %546 = llvm.getelementptr %445[%545] : (!llvm.ptr, i64) -> !llvm.ptr, i8
      %544 = llvm.load %546 : !llvm.ptr -> i8
      %547 = arith.constant 0 : i32
      %549 = arith.extsi %544 : i8 to i32
      %548 = arith.cmpi ne, %549, %547 : i32
      cf.cond_br %548, ^bb72, ^bb73
      ^bb72:
        %550 = llvm.load %513 : !llvm.ptr -> i64
        %551 = llvm.load %509 : !llvm.ptr -> i64
        %552 = llvm.getelementptr %532[%551] : (!llvm.ptr, i64) -> !llvm.ptr, i64
        llvm.store %550, %552 : i64, !llvm.ptr
        %553 = llvm.load %509 : !llvm.ptr -> i64
        %554 = arith.constant 1 : i32
        %556 = arith.extsi %554 : i32 to i64
        %555 = arith.addi %553, %556 : i64
        llvm.store %555, %509 : i64, !llvm.ptr
        cf.br ^bb74
      ^bb73:
        cf.br ^bb74
      ^bb74:
      %557 = llvm.load %513 : !llvm.ptr -> i64
      %558 = arith.constant 1 : i32
      %560 = arith.extsi %558 : i32 to i64
      %559 = arith.addi %557, %560 : i64
      llvm.store %559, %513 : i64, !llvm.ptr
      cf.br ^bb69
    ^bb71:
    %562 = llvm.mlir.addressof @N : !llvm.ptr
    %563 = llvm.load %562 : !llvm.ptr -> i64
    %561 = func.call @iroot3(%563) : (i64) -> i64
    %565 = llvm.load %509 : !llvm.ptr -> i64
    %566 = arith.constant 8 : i32
    %567 = arith.extsi %566 : i32 to i64
    %564 = func.call @calloc(%565, %567) : (i64, i64) -> !llvm.ptr
    %568 = arith.constant 0 : i32
    %569 = arith.extsi %568 : i32 to i64
    %570 = llvm.mlir.constant(1 : i64) : i64
    %571 = llvm.alloca %570 x i64 : (i64) -> !llvm.ptr
    llvm.store %569, %571 : i64, !llvm.ptr
    %572 = arith.constant 0 : i32
    %573 = arith.extsi %572 : i32 to i64
    llvm.store %573, %513 : i64, !llvm.ptr
    cf.br ^bb75
    ^bb75:
    %574 = llvm.load %513 : !llvm.ptr -> i64
    %575 = llvm.load %509 : !llvm.ptr -> i64
    %576 = arith.cmpi slt, %574, %575 : i64
    cf.cond_br %576, ^bb76, ^bb77
    ^bb76:
      %578 = llvm.load %513 : !llvm.ptr -> i64
      %579 = llvm.getelementptr %532[%578] : (!llvm.ptr, i64) -> !llvm.ptr, i64
      %577 = llvm.load %579 : !llvm.ptr -> i64
      %580 = arith.constant 10 : i32
      %582 = arith.extsi %580 : i32 to i64
      %581 = arith.remsi %577, %582 : i64
      %583 = arith.constant 7 : i32
      %585 = arith.extsi %583 : i32 to i64
      %584 = arith.cmpi eq, %581, %585 : i64
      %586 = scf.if %584 -> (i1) {
        %587 = arith.constant 19 : i32
        %589 = arith.extsi %587 : i32 to i64
        %588 = arith.muli %577, %589 : i64
        %590 = llvm.mlir.addressof @N : !llvm.ptr
        %591 = llvm.load %590 : !llvm.ptr -> i64
        %592 = arith.cmpi sle, %588, %591 : i64
        scf.yield %592 : i1
      } else {
        %593 = arith.constant false
        scf.yield %593 : i1
      }
      %594 = scf.if %586 -> (i1) {
        %595 = arith.cmpi sgt, %577, %561 : i64
        scf.yield %595 : i1
      } else {
        %596 = arith.constant false
        scf.yield %596 : i1
      }
      cf.cond_br %594, ^bb78, ^bb79
      ^bb78:
        %597 = llvm.load %571 : !llvm.ptr -> i64
        %598 = llvm.getelementptr %564[%597] : (!llvm.ptr, i64) -> !llvm.ptr, i64
        llvm.store %577, %598 : i64, !llvm.ptr
        %599 = llvm.load %571 : !llvm.ptr -> i64
        %600 = arith.constant 1 : i32
        %602 = arith.extsi %600 : i32 to i64
        %601 = arith.addi %599, %602 : i64
        llvm.store %601, %571 : i64, !llvm.ptr
        cf.br ^bb80
      ^bb79:
        cf.br ^bb80
      ^bb80:
      %603 = llvm.load %513 : !llvm.ptr -> i64
      %604 = arith.constant 1 : i32
      %606 = arith.extsi %604 : i32 to i64
      %605 = arith.addi %603, %606 : i64
      llvm.store %605, %513 : i64, !llvm.ptr
      cf.br ^bb75
    ^bb77:
    %607 = arith.constant 0.0 : f32
    %608 = arith.extf %607 : f32 to f64
    %609 = llvm.mlir.constant(1 : i64) : i64
    %610 = llvm.alloca %609 x f64 : (i64) -> !llvm.ptr
    llvm.store %608, %610 : f64, !llvm.ptr
    %611 = arith.constant 0 : i32
    %612 = arith.extsi %611 : i32 to i64
    llvm.store %612, %513 : i64, !llvm.ptr
    cf.br ^bb81
    ^bb81:
    %613 = llvm.load %513 : !llvm.ptr -> i64
    %614 = llvm.load %509 : !llvm.ptr -> i64
    %615 = arith.cmpi slt, %613, %614 : i64
    cf.cond_br %615, ^bb82, ^bb83
    ^bb82:
      %617 = llvm.load %513 : !llvm.ptr -> i64
      %618 = llvm.getelementptr %532[%617] : (!llvm.ptr, i64) -> !llvm.ptr, i64
      %616 = llvm.load %618 : !llvm.ptr -> i64
      %619 = arith.constant 10 : i32
      %621 = arith.extsi %619 : i32 to i64
      %620 = arith.remsi %616, %621 : i64
      %622 = arith.constant 3 : i32
      %624 = arith.extsi %622 : i32 to i64
      %623 = arith.cmpi eq, %620, %624 : i64
      cf.cond_br %623, ^bb84, ^bb85
      ^bb84:
        %625 = llvm.load %610 : !llvm.ptr -> f64
        %626 = arith.sitofp %616 : i64 to f64
        %627 = math.log %626 : f64
        %628 = arith.addf %625, %627 : f64
        llvm.store %628, %610 : f64, !llvm.ptr
        cf.br ^bb86
      ^bb85:
        %629 = arith.constant 10 : i32
        %631 = arith.extsi %629 : i32 to i64
        %630 = arith.remsi %616, %631 : i64
        %632 = arith.constant 7 : i32
        %634 = arith.extsi %632 : i32 to i64
        %633 = arith.cmpi eq, %630, %634 : i64
        %635 = scf.if %633 -> (i1) {
          %636 = arith.cmpi sle, %616, %561 : i64
          scf.yield %636 : i1
        } else {
          %637 = arith.constant false
          scf.yield %637 : i1
        }
        cf.cond_br %635, ^bb87, ^bb88
        ^bb87:
          %638 = llvm.load %610 : !llvm.ptr -> f64
          %639 = arith.sitofp %616 : i64 to f64
          %640 = math.log %639 : f64
          %641 = arith.addf %638, %640 : f64
          llvm.store %641, %610 : f64, !llvm.ptr
          cf.br ^bb89
        ^bb88:
          cf.br ^bb89
        ^bb89:
        cf.br ^bb86
      ^bb86:
      %642 = llvm.load %513 : !llvm.ptr -> i64
      %643 = arith.constant 1 : i32
      %645 = arith.extsi %643 : i32 to i64
      %644 = arith.addi %642, %645 : i64
      llvm.store %644, %513 : i64, !llvm.ptr
      cf.br ^bb81
    ^bb83:
    %646 = llvm.load %571 : !llvm.ptr -> i64
    %647 = arith.constant 0 : i32
    %649 = arith.extsi %647 : i32 to i64
    %648 = arith.cmpi eq, %646, %649 : i64
    cf.cond_br %648, ^bb90, ^bb91
    ^bb90:
      %650 = llvm.mlir.addressof @str_0 : !llvm.ptr
      %651 = llvm.load %610 : !llvm.ptr -> f64
      %652 = llvm.call @printf(%650, %651) vararg(!llvm.func<i32 (ptr, ...)>) : (!llvm.ptr, f64) -> i32
      func.call @free(%445) : (!llvm.ptr) -> ()
      func.call @free(%532) : (!llvm.ptr) -> ()
      func.call @free(%564) : (!llvm.ptr) -> ()
      %656 = arith.constant 0 : i32
      func.return %656 : i32
    ^bb91:
      cf.br ^bb92
    ^bb92:
    %658 = arith.constant 0 : i32
    %659 = arith.extsi %658 : i32 to i64
    %660 = llvm.getelementptr %564[%659] : (!llvm.ptr, i64) -> !llvm.ptr, i64
    %657 = llvm.load %660 : !llvm.ptr -> i64
    %662 = llvm.load %509 : !llvm.ptr -> i64
    %663 = arith.constant 8 : i32
    %664 = arith.extsi %663 : i32 to i64
    %661 = func.call @calloc(%662, %664) : (i64, i64) -> !llvm.ptr
    %665 = arith.constant 0 : i32
    %666 = arith.extsi %665 : i32 to i64
    %667 = llvm.mlir.constant(1 : i64) : i64
    %668 = llvm.alloca %667 x i64 : (i64) -> !llvm.ptr
    llvm.store %666, %668 : i64, !llvm.ptr
    %669 = llvm.mlir.addressof @N : !llvm.ptr
    %670 = llvm.load %669 : !llvm.ptr -> i64
    %671 = arith.divsi %670, %657 : i64
    %672 = arith.constant 0 : i32
    %673 = arith.extsi %672 : i32 to i64
    llvm.store %673, %513 : i64, !llvm.ptr
    cf.br ^bb93
    ^bb93:
    %674 = llvm.load %513 : !llvm.ptr -> i64
    %675 = llvm.load %509 : !llvm.ptr -> i64
    %676 = arith.cmpi slt, %674, %675 : i64
    cf.cond_br %676, ^bb94, ^bb95
    ^bb94:
      %678 = llvm.load %513 : !llvm.ptr -> i64
      %679 = llvm.getelementptr %532[%678] : (!llvm.ptr, i64) -> !llvm.ptr, i64
      %677 = llvm.load %679 : !llvm.ptr -> i64
      %680 = arith.constant 10 : i32
      %682 = arith.extsi %680 : i32 to i64
      %681 = arith.remsi %677, %682 : i64
      %683 = arith.constant 9 : i32
      %685 = arith.extsi %683 : i32 to i64
      %684 = arith.cmpi eq, %681, %685 : i64
      %686 = scf.if %684 -> (i1) {
        %687 = arith.cmpi sle, %677, %671 : i64
        scf.yield %687 : i1
      } else {
        %688 = arith.constant false
        scf.yield %688 : i1
      }
      cf.cond_br %686, ^bb96, ^bb97
      ^bb96:
        %689 = llvm.load %668 : !llvm.ptr -> i64
        %690 = llvm.getelementptr %661[%689] : (!llvm.ptr, i64) -> !llvm.ptr, i64
        llvm.store %677, %690 : i64, !llvm.ptr
        %691 = llvm.load %668 : !llvm.ptr -> i64
        %692 = arith.constant 1 : i32
        %694 = arith.extsi %692 : i32 to i64
        %693 = arith.addi %691, %694 : i64
        llvm.store %693, %668 : i64, !llvm.ptr
        cf.br ^bb98
      ^bb97:
        cf.br ^bb98
      ^bb98:
      %695 = llvm.load %513 : !llvm.ptr -> i64
      %696 = arith.constant 1 : i32
      %698 = arith.extsi %696 : i32 to i64
      %697 = arith.addi %695, %698 : i64
      llvm.store %697, %513 : i64, !llvm.ptr
      cf.br ^bb93
    ^bb95:
    %699 = llvm.load %668 : !llvm.ptr -> i64
    %700 = arith.constant 0 : i32
    %702 = arith.extsi %700 : i32 to i64
    %701 = arith.cmpi eq, %699, %702 : i64
    cf.cond_br %701, ^bb99, ^bb100
    ^bb99:
      %703 = llvm.mlir.addressof @str_0 : !llvm.ptr
      %704 = llvm.load %610 : !llvm.ptr -> f64
      %705 = llvm.call @printf(%703, %704) vararg(!llvm.func<i32 (ptr, ...)>) : (!llvm.ptr, f64) -> i32
      func.call @free(%445) : (!llvm.ptr) -> ()
      func.call @free(%532) : (!llvm.ptr) -> ()
      func.call @free(%564) : (!llvm.ptr) -> ()
      func.call @free(%661) : (!llvm.ptr) -> ()
      %710 = arith.constant 0 : i32
      func.return %710 : i32
    ^bb100:
      cf.br ^bb101
    ^bb101:
    %712 = llvm.load %571 : !llvm.ptr -> i64
    %713 = arith.constant 8 : i32
    %714 = arith.extsi %713 : i32 to i64
    %711 = func.call @calloc(%712, %714) : (i64, i64) -> !llvm.ptr
    %715 = arith.constant 0 : i32
    %716 = arith.extsi %715 : i32 to i64
    %717 = llvm.mlir.constant(1 : i64) : i64
    %718 = llvm.alloca %717 x i64 : (i64) -> !llvm.ptr
    llvm.store %716, %718 : i64, !llvm.ptr
    %719 = arith.constant 0 : i32
    %720 = arith.extsi %719 : i32 to i64
    llvm.store %720, %513 : i64, !llvm.ptr
    cf.br ^bb102
    ^bb102:
    %721 = llvm.load %513 : !llvm.ptr -> i64
    %722 = llvm.load %571 : !llvm.ptr -> i64
    %723 = arith.cmpi slt, %721, %722 : i64
    cf.cond_br %723, ^bb103, ^bb104
    ^bb103:
      %724 = llvm.mlir.addressof @N : !llvm.ptr
      %725 = llvm.load %724 : !llvm.ptr -> i64
      %727 = llvm.load %513 : !llvm.ptr -> i64
      %728 = llvm.getelementptr %564[%727] : (!llvm.ptr, i64) -> !llvm.ptr, i64
      %726 = llvm.load %728 : !llvm.ptr -> i64
      %729 = arith.divsi %725, %726 : i64
      %730 = arith.constant 0 : i32
      %731 = arith.extsi %730 : i32 to i64
      %732 = llvm.mlir.constant(1 : i64) : i64
      %733 = llvm.alloca %732 x i64 : (i64) -> !llvm.ptr
      llvm.store %731, %733 : i64, !llvm.ptr
      %734 = llvm.load %668 : !llvm.ptr -> i64
      %735 = llvm.mlir.constant(1 : i64) : i64
      %736 = llvm.alloca %735 x i64 : (i64) -> !llvm.ptr
      llvm.store %734, %736 : i64, !llvm.ptr
      cf.br ^bb105
      ^bb105:
      %737 = llvm.load %733 : !llvm.ptr -> i64
      %738 = llvm.load %736 : !llvm.ptr -> i64
      %739 = arith.cmpi slt, %737, %738 : i64
      cf.cond_br %739, ^bb106, ^bb107
      ^bb106:
        %740 = llvm.load %733 : !llvm.ptr -> i64
        %741 = llvm.load %736 : !llvm.ptr -> i64
        %742 = arith.addi %740, %741 : i64
        %743 = arith.constant 2 : i32
        %745 = arith.extsi %743 : i32 to i64
        %744 = arith.divsi %742, %745 : i64
        %747 = llvm.getelementptr %661[%744] : (!llvm.ptr, i64) -> !llvm.ptr, i64
        %746 = llvm.load %747 : !llvm.ptr -> i64
        %748 = arith.cmpi sle, %746, %729 : i64
        cf.cond_br %748, ^bb108, ^bb109
        ^bb108:
          %749 = arith.constant 1 : i32
          %751 = arith.extsi %749 : i32 to i64
          %750 = arith.addi %744, %751 : i64
          llvm.store %750, %733 : i64, !llvm.ptr
          cf.br ^bb110
        ^bb109:
          llvm.store %744, %736 : i64, !llvm.ptr
          cf.br ^bb110
        ^bb110:
        cf.br ^bb105
      ^bb107:
      %752 = llvm.load %733 : !llvm.ptr -> i64
      %753 = llvm.load %513 : !llvm.ptr -> i64
      %754 = llvm.getelementptr %711[%753] : (!llvm.ptr, i64) -> !llvm.ptr, i64
      llvm.store %752, %754 : i64, !llvm.ptr
      %755 = llvm.load %718 : !llvm.ptr -> i64
      %756 = llvm.load %733 : !llvm.ptr -> i64
      %757 = arith.addi %755, %756 : i64
      llvm.store %757, %718 : i64, !llvm.ptr
      %758 = llvm.load %513 : !llvm.ptr -> i64
      %759 = arith.constant 1 : i32
      %761 = arith.extsi %759 : i32 to i64
      %760 = arith.addi %758, %761 : i64
      llvm.store %760, %513 : i64, !llvm.ptr
      cf.br ^bb102
    ^bb104:
    %763 = llvm.load %571 : !llvm.ptr -> i64
    %764 = arith.constant 8 : i32
    %765 = arith.extsi %764 : i32 to i64
    %762 = func.call @calloc(%763, %765) : (i64, i64) -> !llvm.ptr
    %767 = llvm.load %668 : !llvm.ptr -> i64
    %768 = arith.constant 8 : i32
    %769 = arith.extsi %768 : i32 to i64
    %766 = func.call @calloc(%767, %769) : (i64, i64) -> !llvm.ptr
    %770 = arith.constant 0 : i32
    %771 = arith.extsi %770 : i32 to i64
    %772 = llvm.mlir.constant(1 : i64) : i64
    %773 = llvm.alloca %772 x i64 : (i64) -> !llvm.ptr
    llvm.store %771, %773 : i64, !llvm.ptr
    %774 = arith.constant 0 : i32
    %775 = arith.extsi %774 : i32 to i64
    llvm.store %775, %513 : i64, !llvm.ptr
    cf.br ^bb111
    ^bb111:
    %776 = llvm.load %513 : !llvm.ptr -> i64
    %777 = llvm.load %571 : !llvm.ptr -> i64
    %778 = arith.cmpi slt, %776, %777 : i64
    cf.cond_br %778, ^bb112, ^bb113
    ^bb112:
      %780 = llvm.load %513 : !llvm.ptr -> i64
      %781 = llvm.getelementptr %564[%780] : (!llvm.ptr, i64) -> !llvm.ptr, i64
      %779 = llvm.load %781 : !llvm.ptr -> i64
      %782 = arith.sitofp %779 : i64 to f64
      %783 = math.log %782 : f64
      %784 = llvm.mlir.addressof @SCALE : !llvm.ptr
      %785 = llvm.load %784 : !llvm.ptr -> i64
      %786 = arith.sitofp %785 : i64 to f64
      %787 = arith.mulf %783, %786 : f64
      %788 = arith.constant 0.5 : f32
      %790 = arith.extf %788 : f32 to f64
      %789 = arith.addf %787, %790 : f64
      %791 = arith.fptosi %789 : f64 to i64
      %792 = llvm.load %513 : !llvm.ptr -> i64
      %793 = llvm.getelementptr %762[%792] : (!llvm.ptr, i64) -> !llvm.ptr, i64
      llvm.store %791, %793 : i64, !llvm.ptr
      %794 = llvm.load %773 : !llvm.ptr -> i64
      %796 = llvm.load %513 : !llvm.ptr -> i64
      %797 = llvm.getelementptr %762[%796] : (!llvm.ptr, i64) -> !llvm.ptr, i64
      %795 = llvm.load %797 : !llvm.ptr -> i64
      %798 = arith.addi %794, %795 : i64
      llvm.store %798, %773 : i64, !llvm.ptr
      %799 = llvm.load %513 : !llvm.ptr -> i64
      %800 = arith.constant 1 : i32
      %802 = arith.extsi %800 : i32 to i64
      %801 = arith.addi %799, %802 : i64
      llvm.store %801, %513 : i64, !llvm.ptr
      cf.br ^bb111
    ^bb113:
    %803 = arith.constant 0 : i32
    %804 = arith.extsi %803 : i32 to i64
    llvm.store %804, %513 : i64, !llvm.ptr
    cf.br ^bb114
    ^bb114:
    %805 = llvm.load %513 : !llvm.ptr -> i64
    %806 = llvm.load %668 : !llvm.ptr -> i64
    %807 = arith.cmpi slt, %805, %806 : i64
    cf.cond_br %807, ^bb115, ^bb116
    ^bb115:
      %809 = llvm.load %513 : !llvm.ptr -> i64
      %810 = llvm.getelementptr %661[%809] : (!llvm.ptr, i64) -> !llvm.ptr, i64
      %808 = llvm.load %810 : !llvm.ptr -> i64
      %811 = arith.sitofp %808 : i64 to f64
      %812 = math.log %811 : f64
      %813 = llvm.mlir.addressof @SCALE : !llvm.ptr
      %814 = llvm.load %813 : !llvm.ptr -> i64
      %815 = arith.sitofp %814 : i64 to f64
      %816 = arith.mulf %812, %815 : f64
      %817 = arith.constant 0.5 : f32
      %819 = arith.extf %817 : f32 to f64
      %818 = arith.addf %816, %819 : f64
      %820 = arith.fptosi %818 : f64 to i64
      %821 = llvm.load %513 : !llvm.ptr -> i64
      %822 = llvm.getelementptr %766[%821] : (!llvm.ptr, i64) -> !llvm.ptr, i64
      llvm.store %820, %822 : i64, !llvm.ptr
      %823 = llvm.load %773 : !llvm.ptr -> i64
      %825 = llvm.load %513 : !llvm.ptr -> i64
      %826 = llvm.getelementptr %766[%825] : (!llvm.ptr, i64) -> !llvm.ptr, i64
      %824 = llvm.load %826 : !llvm.ptr -> i64
      %827 = arith.addi %823, %824 : i64
      llvm.store %827, %773 : i64, !llvm.ptr
      %828 = llvm.load %513 : !llvm.ptr -> i64
      %829 = arith.constant 1 : i32
      %831 = arith.extsi %829 : i32 to i64
      %830 = arith.addi %828, %831 : i64
      llvm.store %830, %513 : i64, !llvm.ptr
      cf.br ^bb114
    ^bb116:
    %832 = llvm.load %773 : !llvm.ptr -> i64
    %833 = arith.constant 1 : i32
    %835 = arith.extsi %833 : i32 to i64
    %834 = arith.addi %832, %835 : i64
    %836 = arith.constant 2 : i32
    %837 = llvm.load %571 : !llvm.ptr -> i64
    %839 = arith.extsi %836 : i32 to i64
    %838 = arith.addi %839, %837 : i64
    %840 = llvm.load %668 : !llvm.ptr -> i64
    %841 = arith.addi %838, %840 : i64
    %842 = arith.constant 0 : i32
    %843 = arith.extsi %842 : i32 to i64
    %844 = arith.constant 1 : i32
    %845 = arith.extsi %844 : i32 to i64
    %846 = arith.constant 1 : i32
    %847 = llvm.load %571 : !llvm.ptr -> i64
    %849 = arith.extsi %846 : i32 to i64
    %848 = arith.addi %849, %847 : i64
    %850 = arith.constant 1 : i32
    %851 = llvm.load %571 : !llvm.ptr -> i64
    %853 = arith.extsi %850 : i32 to i64
    %852 = arith.addi %853, %851 : i64
    %854 = llvm.load %668 : !llvm.ptr -> i64
    %855 = arith.addi %852, %854 : i64
    %856 = arith.constant 2 : i32
    %857 = llvm.load %571 : !llvm.ptr -> i64
    %858 = llvm.load %668 : !llvm.ptr -> i64
    %859 = arith.addi %857, %858 : i64
    %860 = llvm.load %718 : !llvm.ptr -> i64
    %861 = arith.addi %859, %860 : i64
    %863 = arith.extsi %856 : i32 to i64
    %862 = arith.muli %863, %861 : i64
    %864 = llvm.mlir.addressof @g_n : !llvm.ptr
    llvm.store %841, %864 : i64, !llvm.ptr
    %865 = arith.constant 0 : i32
    %866 = arith.extsi %865 : i32 to i64
    %867 = llvm.mlir.addressof @g_edge_cnt : !llvm.ptr
    llvm.store %866, %867 : i64, !llvm.ptr
    %869 = arith.constant 8 : i32
    %870 = arith.extsi %869 : i32 to i64
    %868 = func.call @calloc(%841, %870) : (i64, i64) -> !llvm.ptr
    %871 = llvm.mlir.addressof @g_head : !llvm.ptr
    llvm.store %868, %871 : !llvm.ptr, !llvm.ptr
    %872 = arith.constant 0 : i32
    %873 = arith.extsi %872 : i32 to i64
    llvm.store %873, %513 : i64, !llvm.ptr
    cf.br ^bb117
    ^bb117:
    %874 = llvm.load %513 : !llvm.ptr -> i64
    %875 = arith.cmpi slt, %874, %841 : i64
    cf.cond_br %875, ^bb118, ^bb119
    ^bb118:
      %876 = arith.constant 1 : i32
      %878 = arith.constant 0 : i32
      %877 = arith.subi %878, %876 : i32
      %879 = llvm.mlir.addressof @g_head : !llvm.ptr
      %880 = llvm.load %879 : !llvm.ptr -> !llvm.ptr
      %881 = llvm.load %513 : !llvm.ptr -> i64
      %882 = arith.extsi %877 : i32 to i64
      %883 = llvm.getelementptr %880[%881] : (!llvm.ptr, i64) -> !llvm.ptr, i64
      llvm.store %882, %883 : i64, !llvm.ptr
      %884 = llvm.load %513 : !llvm.ptr -> i64
      %885 = arith.constant 1 : i32
      %887 = arith.extsi %885 : i32 to i64
      %886 = arith.addi %884, %887 : i64
      llvm.store %886, %513 : i64, !llvm.ptr
      cf.br ^bb117
    ^bb119:
    %889 = arith.constant 8 : i32
    %890 = arith.extsi %889 : i32 to i64
    %888 = func.call @calloc(%862, %890) : (i64, i64) -> !llvm.ptr
    %891 = llvm.mlir.addressof @e_to : !llvm.ptr
    llvm.store %888, %891 : !llvm.ptr, !llvm.ptr
    %893 = arith.constant 8 : i32
    %894 = arith.extsi %893 : i32 to i64
    %892 = func.call @calloc(%862, %894) : (i64, i64) -> !llvm.ptr
    %895 = llvm.mlir.addressof @e_cap : !llvm.ptr
    llvm.store %892, %895 : !llvm.ptr, !llvm.ptr
    %897 = arith.constant 8 : i32
    %898 = arith.extsi %897 : i32 to i64
    %896 = func.call @calloc(%862, %898) : (i64, i64) -> !llvm.ptr
    %899 = llvm.mlir.addressof @e_rev : !llvm.ptr
    llvm.store %896, %899 : !llvm.ptr, !llvm.ptr
    %901 = arith.constant 8 : i32
    %902 = arith.extsi %901 : i32 to i64
    %900 = func.call @calloc(%862, %902) : (i64, i64) -> !llvm.ptr
    %903 = llvm.mlir.addressof @e_next : !llvm.ptr
    llvm.store %900, %903 : !llvm.ptr, !llvm.ptr
    %905 = arith.constant 8 : i32
    %906 = arith.extsi %905 : i32 to i64
    %904 = func.call @calloc(%841, %906) : (i64, i64) -> !llvm.ptr
    %907 = llvm.mlir.addressof @g_level : !llvm.ptr
    llvm.store %904, %907 : !llvm.ptr, !llvm.ptr
    %909 = arith.constant 8 : i32
    %910 = arith.extsi %909 : i32 to i64
    %908 = func.call @calloc(%841, %910) : (i64, i64) -> !llvm.ptr
    %911 = llvm.mlir.addressof @g_it : !llvm.ptr
    llvm.store %908, %911 : !llvm.ptr, !llvm.ptr
    %913 = arith.constant 8 : i32
    %914 = arith.extsi %913 : i32 to i64
    %912 = func.call @calloc(%841, %914) : (i64, i64) -> !llvm.ptr
    %915 = llvm.mlir.addressof @g_q : !llvm.ptr
    llvm.store %912, %915 : !llvm.ptr, !llvm.ptr
    %916 = arith.constant 0 : i32
    %917 = arith.extsi %916 : i32 to i64
    llvm.store %917, %513 : i64, !llvm.ptr
    cf.br ^bb120
    ^bb120:
    %918 = llvm.load %513 : !llvm.ptr -> i64
    %919 = llvm.load %571 : !llvm.ptr -> i64
    %920 = arith.cmpi slt, %918, %919 : i64
    cf.cond_br %920, ^bb121, ^bb122
    ^bb121:
      %922 = llvm.load %513 : !llvm.ptr -> i64
      %923 = arith.addi %845, %922 : i64
      %925 = llvm.load %513 : !llvm.ptr -> i64
      %926 = llvm.getelementptr %762[%925] : (!llvm.ptr, i64) -> !llvm.ptr, i64
      %924 = llvm.load %926 : !llvm.ptr -> i64
      func.call @add_edge(%843, %923, %924) : (i64, i64, i64) -> ()
      %927 = llvm.load %513 : !llvm.ptr -> i64
      %928 = arith.constant 1 : i32
      %930 = arith.extsi %928 : i32 to i64
      %929 = arith.addi %927, %930 : i64
      llvm.store %929, %513 : i64, !llvm.ptr
      cf.br ^bb120
    ^bb122:
    %931 = arith.constant 0 : i32
    %932 = arith.extsi %931 : i32 to i64
    llvm.store %932, %513 : i64, !llvm.ptr
    cf.br ^bb123
    ^bb123:
    %933 = llvm.load %513 : !llvm.ptr -> i64
    %934 = llvm.load %668 : !llvm.ptr -> i64
    %935 = arith.cmpi slt, %933, %934 : i64
    cf.cond_br %935, ^bb124, ^bb125
    ^bb124:
      %937 = llvm.load %513 : !llvm.ptr -> i64
      %938 = arith.addi %848, %937 : i64
      %940 = llvm.load %513 : !llvm.ptr -> i64
      %941 = llvm.getelementptr %766[%940] : (!llvm.ptr, i64) -> !llvm.ptr, i64
      %939 = llvm.load %941 : !llvm.ptr -> i64
      func.call @add_edge(%938, %855, %939) : (i64, i64, i64) -> ()
      %942 = llvm.load %513 : !llvm.ptr -> i64
      %943 = arith.constant 1 : i32
      %945 = arith.extsi %943 : i32 to i64
      %944 = arith.addi %942, %945 : i64
      llvm.store %944, %513 : i64, !llvm.ptr
      cf.br ^bb123
    ^bb125:
    %946 = arith.constant 0 : i32
    %947 = arith.extsi %946 : i32 to i64
    llvm.store %947, %513 : i64, !llvm.ptr
    cf.br ^bb126
    ^bb126:
    %948 = llvm.load %513 : !llvm.ptr -> i64
    %949 = llvm.load %571 : !llvm.ptr -> i64
    %950 = arith.cmpi slt, %948, %949 : i64
    cf.cond_br %950, ^bb127, ^bb128
    ^bb127:
      %951 = llvm.load %513 : !llvm.ptr -> i64
      %952 = arith.addi %845, %951 : i64
      %954 = llvm.load %513 : !llvm.ptr -> i64
      %955 = llvm.getelementptr %711[%954] : (!llvm.ptr, i64) -> !llvm.ptr, i64
      %953 = llvm.load %955 : !llvm.ptr -> i64
      %956 = arith.constant 0 : i32
      %957 = arith.extsi %956 : i32 to i64
      %958 = llvm.mlir.constant(1 : i64) : i64
      %959 = llvm.alloca %958 x i64 : (i64) -> !llvm.ptr
      llvm.store %957, %959 : i64, !llvm.ptr
      cf.br ^bb129
      ^bb129:
      %960 = llvm.load %959 : !llvm.ptr -> i64
      %961 = arith.cmpi slt, %960, %953 : i64
      cf.cond_br %961, ^bb130, ^bb131
      ^bb130:
        %963 = llvm.load %959 : !llvm.ptr -> i64
        %964 = arith.addi %848, %963 : i64
        func.call @add_edge(%952, %964, %834) : (i64, i64, i64) -> ()
        %965 = llvm.load %959 : !llvm.ptr -> i64
        %966 = arith.constant 1 : i32
        %968 = arith.extsi %966 : i32 to i64
        %967 = arith.addi %965, %968 : i64
        llvm.store %967, %959 : i64, !llvm.ptr
        cf.br ^bb129
      ^bb131:
      %969 = llvm.load %513 : !llvm.ptr -> i64
      %970 = arith.constant 1 : i32
      %972 = arith.extsi %970 : i32 to i64
      %971 = arith.addi %969, %972 : i64
      llvm.store %971, %513 : i64, !llvm.ptr
      cf.br ^bb126
    ^bb128:
    %973 = func.call @max_flow(%843, %855) : (i64, i64) -> i64
    %975 = arith.constant 1 : i32
    %976 = arith.extsi %975 : i32 to i64
    %974 = func.call @calloc(%841, %976) : (i64, i64) -> !llvm.ptr
    func.call @reachable(%843, %974) : (i64, !llvm.ptr) -> ()
    %978 = arith.constant 0 : i32
    %979 = arith.extsi %978 : i32 to i64
    llvm.store %979, %513 : i64, !llvm.ptr
    cf.br ^bb132
    ^bb132:
    %980 = llvm.load %513 : !llvm.ptr -> i64
    %981 = llvm.load %571 : !llvm.ptr -> i64
    %982 = arith.cmpi slt, %980, %981 : i64
    cf.cond_br %982, ^bb133, ^bb134
    ^bb133:
      %984 = llvm.load %513 : !llvm.ptr -> i64
      %985 = arith.addi %845, %984 : i64
      %986 = llvm.getelementptr %974[%985] : (!llvm.ptr, i64) -> !llvm.ptr, i8
      %983 = llvm.load %986 : !llvm.ptr -> i8
      %987 = arith.constant 0 : i32
      %989 = arith.extsi %983 : i8 to i32
      %988 = arith.cmpi eq, %989, %987 : i32
      cf.cond_br %988, ^bb135, ^bb136
      ^bb135:
        %990 = llvm.load %610 : !llvm.ptr -> f64
        %992 = llvm.load %513 : !llvm.ptr -> i64
        %993 = llvm.getelementptr %564[%992] : (!llvm.ptr, i64) -> !llvm.ptr, i64
        %991 = llvm.load %993 : !llvm.ptr -> i64
        %994 = arith.sitofp %991 : i64 to f64
        %995 = math.log %994 : f64
        %996 = arith.addf %990, %995 : f64
        llvm.store %996, %610 : f64, !llvm.ptr
        cf.br ^bb137
      ^bb136:
        cf.br ^bb137
      ^bb137:
      %997 = llvm.load %513 : !llvm.ptr -> i64
      %998 = arith.constant 1 : i32
      %1000 = arith.extsi %998 : i32 to i64
      %999 = arith.addi %997, %1000 : i64
      llvm.store %999, %513 : i64, !llvm.ptr
      cf.br ^bb132
    ^bb134:
    %1001 = arith.constant 0 : i32
    %1002 = arith.extsi %1001 : i32 to i64
    llvm.store %1002, %513 : i64, !llvm.ptr
    cf.br ^bb138
    ^bb138:
    %1003 = llvm.load %513 : !llvm.ptr -> i64
    %1004 = llvm.load %668 : !llvm.ptr -> i64
    %1005 = arith.cmpi slt, %1003, %1004 : i64
    cf.cond_br %1005, ^bb139, ^bb140
    ^bb139:
      %1007 = llvm.load %513 : !llvm.ptr -> i64
      %1008 = arith.addi %848, %1007 : i64
      %1009 = llvm.getelementptr %974[%1008] : (!llvm.ptr, i64) -> !llvm.ptr, i8
      %1006 = llvm.load %1009 : !llvm.ptr -> i8
      %1010 = arith.constant 0 : i32
      %1012 = arith.extsi %1006 : i8 to i32
      %1011 = arith.cmpi ne, %1012, %1010 : i32
      cf.cond_br %1011, ^bb141, ^bb142
      ^bb141:
        %1013 = llvm.load %610 : !llvm.ptr -> f64
        %1015 = llvm.load %513 : !llvm.ptr -> i64
        %1016 = llvm.getelementptr %661[%1015] : (!llvm.ptr, i64) -> !llvm.ptr, i64
        %1014 = llvm.load %1016 : !llvm.ptr -> i64
        %1017 = arith.sitofp %1014 : i64 to f64
        %1018 = math.log %1017 : f64
        %1019 = arith.addf %1013, %1018 : f64
        llvm.store %1019, %610 : f64, !llvm.ptr
        cf.br ^bb143
      ^bb142:
        cf.br ^bb143
      ^bb143:
      %1020 = llvm.load %513 : !llvm.ptr -> i64
      %1021 = arith.constant 1 : i32
      %1023 = arith.extsi %1021 : i32 to i64
      %1022 = arith.addi %1020, %1023 : i64
      llvm.store %1022, %513 : i64, !llvm.ptr
      cf.br ^bb138
    ^bb140:
    %1024 = llvm.mlir.addressof @str_0 : !llvm.ptr
    %1025 = llvm.load %610 : !llvm.ptr -> f64
    %1026 = llvm.call @printf(%1024, %1025) vararg(!llvm.func<i32 (ptr, ...)>) : (!llvm.ptr, f64) -> i32
    func.call @free(%974) : (!llvm.ptr) -> ()
    %1029 = llvm.mlir.addressof @g_head : !llvm.ptr
    %1030 = llvm.load %1029 : !llvm.ptr -> !llvm.ptr
    func.call @free(%1030) : (!llvm.ptr) -> ()
    %1032 = llvm.mlir.addressof @e_to : !llvm.ptr
    %1033 = llvm.load %1032 : !llvm.ptr -> !llvm.ptr
    func.call @free(%1033) : (!llvm.ptr) -> ()
    %1035 = llvm.mlir.addressof @e_cap : !llvm.ptr
    %1036 = llvm.load %1035 : !llvm.ptr -> !llvm.ptr
    func.call @free(%1036) : (!llvm.ptr) -> ()
    %1038 = llvm.mlir.addressof @e_rev : !llvm.ptr
    %1039 = llvm.load %1038 : !llvm.ptr -> !llvm.ptr
    func.call @free(%1039) : (!llvm.ptr) -> ()
    %1041 = llvm.mlir.addressof @e_next : !llvm.ptr
    %1042 = llvm.load %1041 : !llvm.ptr -> !llvm.ptr
    func.call @free(%1042) : (!llvm.ptr) -> ()
    %1044 = llvm.mlir.addressof @g_level : !llvm.ptr
    %1045 = llvm.load %1044 : !llvm.ptr -> !llvm.ptr
    func.call @free(%1045) : (!llvm.ptr) -> ()
    %1047 = llvm.mlir.addressof @g_it : !llvm.ptr
    %1048 = llvm.load %1047 : !llvm.ptr -> !llvm.ptr
    func.call @free(%1048) : (!llvm.ptr) -> ()
    %1050 = llvm.mlir.addressof @g_q : !llvm.ptr
    %1051 = llvm.load %1050 : !llvm.ptr -> !llvm.ptr
    func.call @free(%1051) : (!llvm.ptr) -> ()
    func.call @free(%762) : (!llvm.ptr) -> ()
    func.call @free(%766) : (!llvm.ptr) -> ()
    func.call @free(%711) : (!llvm.ptr) -> ()
    func.call @free(%564) : (!llvm.ptr) -> ()
    func.call @free(%661) : (!llvm.ptr) -> ()
    func.call @free(%532) : (!llvm.ptr) -> ()
    func.call @free(%445) : (!llvm.ptr) -> ()
    %1059 = arith.constant 0 : i32
    func.return %1059 : i32
  }
}