Problem 495

Writing n as Product of k Squarefree Integers W(10000!, 30) mod 1_000_000_007 via cycle-index generating functions.

Answer789107601
Output789107601
StatusPASS
Native helperno
Runtime290 ms
Peak memory40320 KB
Time complexityO(n^2) (estimated)
Space complexityO(n^2) (estimated)

Performance comparison

MetricOur solutionBest known
Time complexityO(n^2)O(n log log n)
Space complexityO(n^2)O(n)
ApproachFlow solutionMobius sieve
VerdictSuboptimal

Flow source

# Project Euler 495
# Writing n as Product of k Squarefree Integers
# W(10000!, 30) mod 1_000_000_007 via cycle-index generating functions.

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

const MOD: i64 = 1000000007
const K: i64 = 30
const NFACT: i64 = 10000

let mut MAX_E: i64 = 0
let mut NEXP: i64 = 0
let mut EXP_LIST: ptr<i64> = null
let mut EXP_CNT: ptr<i64> = null
let mut FREQ1: i64 = 0
let mut BASE: ptr<i64> = null
let mut STRIDE: i64 = 0
let mut R2MAX: i64 = 0
let mut INV_FACT: ptr<i64> = null
let mut INV_POW: ptr<i64> = null
let mut PARTS: ptr<i64> = null
let mut BPART: i64 = 0
let mut TOTAL: i64 = 0
let mut DP_TMP: ptr<i64> = null

function modpow(base0: i64, exp0: i64) -> i64 {
    let mut r: i64 = 1
    let mut b: i64 = base0 % MOD
    let mut e: i64 = exp0
    while e > 0 {
        if e % 2 == 1 {
            r = (r * b) % MOD
        }
        b = (b * b) % MOD
        e = e / 2
    }
    return r
}

function base_off(r1: i64, r2: i64) -> i64 {
    return (r1 * R2MAX + r2) * STRIDE
}

function apply_parts() -> void {
    let mut r1: i64 = 0
    let mut r2: i64 = 0
    let mut i: i64 = BPART - 1
    while i >= 0 {
        if PARTS[i] != 1 { break }
        r1 = r1 + 1
        i = i - 1
    }
    while i >= 0 {
        if PARTS[i] != 2 { break }
        r2 = r2 + 1
        i = i - 1
    }
    if r1 == 0 {
        if FREQ1 != 0 { return }
    }
    let off: i64 = base_off(r1, r2)
    let mut t: i64 = 0
    while t <= MAX_E {
        DP_TMP[t] = BASE[off + t]
        t = t + 1
    }
    let mut j: i64 = 0
    while j <= i {
        let m: i64 = PARTS[j]
        t = m
        while t <= MAX_E {
            let x: i64 = DP_TMP[t] + DP_TMP[t - m]
            if x >= MOD {
                DP_TMP[t] = x - MOD
            } else {
                DP_TMP[t] = x
            }
            t = t + 1
        }
        j = j + 1
    }
    let mut fval: i64 = 1
    j = 0
    while j < NEXP {
        let c: i64 = DP_TMP[EXP_LIST[j]]
        if c == 0 {
            fval = 0
            break
        }
        let fe: i64 = EXP_CNT[j]
        if fe == 1 {
            fval = (fval * c) % MOD
        } else {
            if c != 1 {
                fval = (fval * modpow(c, fe)) % MOD
            }
        }
        j = j + 1
    }
    if fval == 0 { return }
    let mut w: i64 = 1
    j = 0
    while j < BPART {
        let s: i64 = PARTS[j]
        let mut jj: i64 = j + 1
        while jj < BPART {
            if PARTS[jj] != s { break }
            jj = jj + 1
        }
        let cnt: i64 = jj - j
        w = (w * INV_POW[s * (K + 1) + cnt]) % MOD
        w = (w * INV_FACT[cnt]) % MOD
        j = jj
    }
    if ((K - BPART) & 1) != 0 {
        w = MOD - w
    }
    TOTAL = (TOTAL + w * fval) % MOD
}

function gen(rem: i64, maxp0: i64) -> void {
    if rem == 0 {
        apply_parts()
        return
    }
    let mut maxp: i64 = maxp0
    if maxp > rem { maxp = rem }
    let mut first: i64 = maxp
    while first >= 1 {
        PARTS[BPART] = first
        BPART = BPART + 1
        gen(rem - first, first)
        BPART = BPART - 1
        first = first - 1
    }
}

function main() -> i32 {
    let freq_arr: ptr<i64> = calloc(NFACT + 1, 8)
    let sieve: ptr<i8> = calloc(NFACT + 1, 1)
    if freq_arr == null || sieve == null { return 1 }
    let mut i: i64 = 0
    while i <= NFACT {
        sieve[i] = 1
        i = i + 1
    }
    sieve[0] = 0
    sieve[1] = 0
    i = 2
    while i * i <= NFACT {
        if sieve[i] != 0 {
            let mut j: i64 = i * i
            while j <= NFACT {
                sieve[j] = 0
                j = j + i
            }
        }
        i = i + 1
    }
    MAX_E = 0
    let mut p: i64 = 2
    while p <= NFACT {
        if sieve[p] != 0 {
            let mut e: i64 = 0
            let mut nn: i64 = NFACT
            while nn > 0 {
                nn = nn / p
                e = e + nn
            }
            freq_arr[e] = freq_arr[e] + 1
            if e > MAX_E { MAX_E = e }
        }
        p = p + 1
    }
    free(sieve)
    FREQ1 = freq_arr[1]
    EXP_LIST = calloc(MAX_E + 1, 8)
    EXP_CNT = calloc(MAX_E + 1, 8)
    if EXP_LIST == null || EXP_CNT == null { return 1 }
    NEXP = 0
    let mut e2: i64 = 1
    while e2 <= MAX_E {
        if freq_arr[e2] != 0 {
            EXP_LIST[NEXP] = e2
            EXP_CNT[NEXP] = freq_arr[e2]
            NEXP = NEXP + 1
        }
        e2 = e2 + 1
    }
    free(freq_arr)

    STRIDE = MAX_E + 1
    R2MAX = K / 2 + 1
    let nbase: i64 = (K + 1) * R2MAX * STRIDE
    BASE = calloc(nbase, 8)
    if BASE == null { return 1 }
    let inv: ptr<i64> = calloc(MAX_E + 1, 8)
    if inv == null { return 1 }
    if MAX_E >= 1 { inv[1] = 1 }
    i = 2
    while i <= MAX_E {
        inv[i] = (MOD - (MOD / i) * inv[MOD % i] % MOD) % MOD
        i = i + 1
    }
    let mut r1: i64 = 0
    while r1 <= K {
        let off0: i64 = base_off(r1, 0)
        BASE[off0] = 1
        if r1 > 0 {
            let mut t: i64 = 1
            while t <= MAX_E {
                BASE[off0 + t] = ((BASE[off0 + t - 1] * (r1 + t - 1)) % MOD * inv[t]) % MOD
                t = t + 1
            }
        }
        let mut r2: i64 = 1
        while r2 <= K / 2 {
            let off: i64 = base_off(r1, r2)
            let prev: i64 = base_off(r1, r2 - 1)
            let mut t2: i64 = 0
            while t2 <= MAX_E {
                BASE[off + t2] = BASE[prev + t2]
                t2 = t2 + 1
            }
            t2 = 2
            while t2 <= MAX_E {
                let x: i64 = BASE[off + t2] + BASE[off + t2 - 2]
                if x >= MOD {
                    BASE[off + t2] = x - MOD
                } else {
                    BASE[off + t2] = x
                }
                t2 = t2 + 1
            }
            r2 = r2 + 1
        }
        r1 = r1 + 1
    }
    free(inv)

    INV_FACT = calloc(K + 1, 8)
    INV_POW = calloc((K + 1) * (K + 1), 8)
    PARTS = calloc(K + 1, 8)
    DP_TMP = calloc(MAX_E + 1, 8)
    if INV_FACT == null || INV_POW == null || PARTS == null || DP_TMP == null { return 1 }

    let fact: ptr<i64> = calloc(K + 1, 8)
    if fact == null { return 1 }
    fact[0] = 1
    i = 1
    while i <= K {
        fact[i] = (fact[i - 1] * i) % MOD
        i = i + 1
    }
    INV_FACT[K] = modpow(fact[K], MOD - 2)
    i = K
    while i >= 1 {
        INV_FACT[i - 1] = (INV_FACT[i] * i) % MOD
        i = i - 1
    }
    free(fact)

    let mut s: i64 = 1
    while s <= K {
        let invs: i64 = modpow(s, MOD - 2)
        let mut acc: i64 = 1
        INV_POW[s * (K + 1) + 0] = 1
        let mut c: i64 = 1
        while c <= K {
            acc = (acc * invs) % MOD
            INV_POW[s * (K + 1) + c] = acc
            c = c + 1
        }
        s = s + 1
    }

    TOTAL = 0
    BPART = 0
    gen(K, K)
    printf("%lld\n", TOTAL)

    free(DP_TMP)
    free(PARTS)
    free(INV_POW)
    free(INV_FACT)
    free(BASE)
    free(EXP_CNT)
    free(EXP_LIST)
    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 modpow_i64_i64(int64_t base0, int64_t exp0);
int64_t base_off_i64_i64(int64_t r1, int64_t r2);
void apply_parts(void);
void gen_i64_i64(int64_t rem, int64_t maxp0);
int32_t main(void);

static const int64_t MOD = 1000000007;
static const int64_t K = 30;
static const int64_t NFACT = 10000;

/* Module statics */
static int64_t MAX_E = 0;
static int64_t NEXP = 0;
static int64_t* EXP_LIST = NULL;
static int64_t* EXP_CNT = NULL;
static int64_t FREQ1 = 0;
static int64_t* BASE = NULL;
static int64_t STRIDE = 0;
static int64_t R2MAX = 0;
static int64_t* INV_FACT = NULL;
static int64_t* INV_POW = NULL;
static int64_t* PARTS = NULL;
static int64_t BPART = 0;
static int64_t TOTAL = 0;
static int64_t* DP_TMP = NULL;



int64_t modpow_i64_i64(int64_t base0, int64_t exp0) {
    int64_t r = 1;
    int64_t b = FLOW_CHECKED_MOD((base0), (MOD));
    int64_t e = exp0;
    while (e > 0) {
        if (FLOW_CHECKED_MOD((e), (2)) == 1) {
            r = FLOW_CHECKED_MOD(((r * b)), (MOD));
        }
        b = FLOW_CHECKED_MOD(((b * b)), (MOD));
        e = FLOW_CHECKED_DIV((e), (2));
    }
    return r;
}

int64_t base_off_i64_i64(int64_t r1, int64_t r2) {
    return (((r1 * R2MAX) + r2) * STRIDE);
}

void apply_parts(void) {
    int64_t r1 = 0;
    int64_t r2 = 0;
    int64_t i = (BPART - 1);
    while (i >= 0) {
        if (PARTS[i] != 1) {
            break;
        }
        r1 = (r1 + 1);
        i = (i - 1);
    }
    while (i >= 0) {
        if (PARTS[i] != 2) {
            break;
        }
        r2 = (r2 + 1);
        i = (i - 1);
    }
    if (r1 == 0) {
        if (FREQ1 != 0) {
            return;
        }
    }
    int64_t off = base_off_i64_i64(r1, r2);
    int64_t t = 0;
    while (t <= MAX_E) {
        DP_TMP[t] = BASE[(off + t)];
        t = (t + 1);
    }
    int64_t j = 0;
    while (j <= i) {
        int64_t m = PARTS[j];
        t = m;
        while (t <= MAX_E) {
            int64_t x = (DP_TMP[t] + DP_TMP[(t - m)]);
            if (x >= MOD) {
                DP_TMP[t] = (x - MOD);
            } else {
                DP_TMP[t] = x;
            }
            t = (t + 1);
        }
        j = (j + 1);
    }
    int64_t fval = 1;
    j = 0;
    while (j < NEXP) {
        int64_t c = DP_TMP[EXP_LIST[j]];
        if (c == 0) {
            fval = 0;
            break;
        }
        int64_t fe = EXP_CNT[j];
        if (fe == 1) {
            fval = FLOW_CHECKED_MOD(((fval * c)), (MOD));
        } else {
            if (c != 1) {
                fval = FLOW_CHECKED_MOD(((fval * modpow_i64_i64(c, fe))), (MOD));
            }
        }
        j = (j + 1);
    }
    if (fval == 0) {
        return;
    }
    int64_t w = 1;
    j = 0;
    while (j < BPART) {
        int64_t s = PARTS[j];
        int64_t jj = (j + 1);
        while (jj < BPART) {
            if (PARTS[jj] != s) {
                break;
            }
            jj = (jj + 1);
        }
        int64_t cnt = (jj - j);
        w = FLOW_CHECKED_MOD(((w * INV_POW[((s * (K + 1)) + cnt)])), (MOD));
        w = FLOW_CHECKED_MOD(((w * INV_FACT[cnt])), (MOD));
        j = jj;
    }
    if (((K - BPART) & 1) != 0) {
        w = (MOD - w);
    }
    TOTAL = FLOW_CHECKED_MOD(((TOTAL + (w * fval))), (MOD));
}

void gen_i64_i64(int64_t rem, int64_t maxp0) {
    if (rem == 0) {
        apply_parts();
        return;
    }
    int64_t maxp = maxp0;
    if (maxp > rem) {
        maxp = rem;
    }
    int64_t first = maxp;
    while (first >= 1) {
        PARTS[BPART] = first;
        BPART = (BPART + 1);
        gen_i64_i64((rem - first), first);
        BPART = (BPART - 1);
        first = (first - 1);
    }
}

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