// File: chhylp123[P]hifiasm/gfa_ut.cpp
// Function: get_best_path(ul_resolve_t*, usg_t*, unsigned long*, unsigned long*, unsigned char*, unsigned long, unsigned long, asg64_v*, asg64_v*, unsigned long*, unsigned long*, asg64_v*)
void get_best_path(ul_resolve_t *uidx, usg_t *ng, uint64_t *int_idx, uint64_t *int_a, 
uint8_t *f, uint64_t s, uint64_t e, asg64_v *b0, asg64_v *b1, uint64_t *ridx_a, 
uint64_t *ridx, asg64_v *res)
{
    uint64_t pi, v, k, res_n = res->n, *pz;
    b0->n = b1->n = 0; 
    for (k = s, pi = (uint64_t)-1; k < e; k++) {
        v = int_a[k]; 
        if((!f[v])&&(!f[v^1])) continue;
        if((pi != (uint64_t)-1) && (f[v^1])) {
            // if(((v>>1) == 2736)) {
            //     fprintf(stderr, "+[M::%s::ii[%lu, %lu)] utg%.6dl(%c), f[v^1]::%u, v^1::%lu, putg%.6dl(%c), f[pv]::%u, pv::%lu\n", 
            //     __func__, s, e, (int32_t)(int_a[k]>>1)+1, "+-"[int_a[k]&1], f[v^1], v^1, 
            //     (int32_t)(int_a[pi]>>1)+1, "+-"[int_a[pi]&1], f[int_a[pi]], int_a[pi]);
            // }
            if(is_best_path(uidx, ng, int_idx, int_a, s, e, k, v^1, ridx_a, ridx, b0, b1, 0.51)) {
                // if(((v>>1) == 2736)) {
                //     fprintf(stderr, "-[M::%s::ii[%lu, %lu)] utg%.6dl(%c), f[v^1]::%u, v^1::%lu, putg%.6dl(%c), f[pv]::%u, pv::%lu\n", 
                //     __func__, s, e, (int32_t)(int_a[k]>>1)+1, "+-"[int_a[k]&1], f[v^1], v^1, 
                //     (int32_t)(int_a[pi]>>1)+1, "+-"[int_a[pi]&1], f[int_a[pi]], int_a[pi]);
                // }
                // fprintf(stderr, "[M::%s::] utg%.6dl(%c)->utg%.6dl(%c)\n", __func__, 
                //     (int32_t)(int_a[pi]>>1)+1, "+-"[int_a[pi]&1], (int32_t)(int_a[k]>>1)+1, "+-"[int_a[k]&1]);
                pz = (res->n > res_n)? &(res->a[res->n-1]):(NULL);
                if((!pz) || (((*pz)>>32) != pi)) {//keep the shortest pi<->k 
                    kv_pushp(uint64_t, *res, &pz); 
                    *pz = pi; (*pz) <<= 32; (*pz) |= k;
                }
            } else {
                pi = (uint64_t)-1;
            }            
        }
        if(f[v]) {
            // fprintf(stderr, "-[M::%s::] utg%.6dl(%c), f[v]::%u, v::%lu\n", 
            // __func__, (int32_t)(int_a[k]>>1)+1, "+-"[int_a[k]&1], f[v], v);
            if(is_best_path(uidx, ng, int_idx, int_a, s, e, k, v, ridx_a, ridx, b0, b1, 0.51)) {
                pi = k;//keep the shortest pi<->k
            } else {
                pi = (uint64_t)-1;
            }
        }
    }
}

// Assembly:
pushq %rbp
movq %rsp, %rbp
pushq %r15
pushq %r14
pushq %rbx
subq $0x98, %rsp
movq 0x38(%rbp), %rax
movq 0x30(%rbp), %rax
movq 0x28(%rbp), %rax
movq 0x20(%rbp), %rax
movq 0x18(%rbp), %rax
movq 0x10(%rbp), %rax
movq %rdi, -0x20(%rbp)
movq %rsi, -0x28(%rbp)
movq %rdx, -0x30(%rbp)
movq %rcx, -0x38(%rbp)
movq %r8, -0x40(%rbp)
movq %r9, -0x48(%rbp)
movq 0x38(%rbp), %rax
movq (%rax), %rax
movq %rax, -0x68(%rbp)
movq 0x20(%rbp), %rax
movq $0x0, (%rax)
movq 0x18(%rbp), %rax
movq $0x0, (%rax)
movq -0x48(%rbp), %rax
movq %rax, -0x60(%rbp)
movq $-0x1, -0x50(%rbp)
movq -0x60(%rbp), %rax
cmpq 0x10(%rbp), %rax
jae 0x362c1c
movq -0x38(%rbp), %rax
movq -0x60(%rbp), %rcx
movq (%rax,%rcx,8), %rax
movq %rax, -0x58(%rbp)
movq -0x40(%rbp), %rax
movq -0x58(%rbp), %rcx
cmpb $0x0, (%rax,%rcx)
jne 0x3629e4
movq -0x40(%rbp), %rax
movq -0x58(%rbp), %rcx
xorq $0x1, %rcx
cmpb $0x0, (%rax,%rcx)
jne 0x3629e4
jmp 0x362c0b
cmpq $-0x1, -0x50(%rbp)
je 0x362b87
movq -0x40(%rbp), %rax
movq -0x58(%rbp), %rcx
xorq $0x1, %rcx
cmpb $0x0, (%rax,%rcx)
je 0x362b87
movq -0x20(%rbp), %rdi
movq -0x28(%rbp), %rsi
movq -0x30(%rbp), %rdx
movq -0x38(%rbp), %rcx
movq -0x48(%rbp), %r8
movq 0x10(%rbp), %r9
movq -0x60(%rbp), %r15
movq -0x58(%rbp), %r14
xorq $0x1, %r14
movq 0x28(%rbp), %rbx
movq 0x30(%rbp), %r11
movq 0x18(%rbp), %r10
movq 0x20(%rbp), %rax
movsd 0x1f877(%rip), %xmm0    # 0x3822b8
movq %r15, (%rsp)
movq %r14, 0x8(%rsp)
movq %rbx, 0x10(%rsp)
movq %r11, 0x18(%rsp)
movq %r10, 0x20(%rsp)
movq %rax, 0x28(%rsp)
callq 0x362470
cmpq $0x0, %rax
je 0x362b7d
movq 0x38(%rbp), %rax
movq (%rax), %rax
cmpq -0x68(%rbp), %rax
jbe 0x362a9a
movq 0x38(%rbp), %rax
movq 0x10(%rax), %rax
movq 0x38(%rbp), %rcx
movq (%rcx), %rcx
subq $0x1, %rcx
shlq $0x3, %rcx
addq %rcx, %rax
movq %rax, -0x78(%rbp)
jmp 0x362aa2
xorl %eax, %eax
movq %rax, -0x78(%rbp)
jmp 0x362aa2
movq -0x78(%rbp), %rax
movq %rax, -0x70(%rbp)
cmpq $0x0, -0x70(%rbp)
je 0x362ac6
movq -0x70(%rbp), %rax
movq (%rax), %rax
shrq $0x20, %rax
cmpq -0x50(%rbp), %rax
je 0x362b7b
jmp 0x362ac8
movq 0x38(%rbp), %rax
movq (%rax), %rax
movq 0x38(%rbp), %rcx
cmpq 0x8(%rcx), %rax
jne 0x362b30
movq 0x38(%rbp), %rax
cmpq $0x0, 0x8(%rax)
je 0x362af5
movq 0x38(%rbp), %rax
movq 0x8(%rax), %rax
shlq %rax
movq %rax, -0x80(%rbp)
jmp 0x362b00
movl $0x2, %eax
movq %rax, -0x80(%rbp)
jmp 0x362b00
movq -0x80(%rbp), %rcx
movq 0x38(%rbp), %rax
movq %rcx, 0x8(%rax)
movq 0x38(%rbp), %rax
movq 0x10(%rax), %rdi
movq 0x38(%rbp), %rax
movq 0x8(%rax), %rsi
shlq $0x3, %rsi
callq 0x2360
movq %rax, %rcx
movq 0x38(%rbp), %rax
movq %rcx, 0x10(%rax)
movq 0x38(%rbp), %rax
movq 0x10(%rax), %rax
movq 0x38(%rbp), %rdx
movq (%rdx), %rcx
movq %rcx, %rsi
addq $0x1, %rsi
movq %rsi, (%rdx)
shlq $0x3, %rcx
addq %rcx, %rax
movq %rax, -0x70(%rbp)
movq -0x50(%rbp), %rcx
movq -0x70(%rbp), %rax
movq %rcx, (%rax)
movq -0x70(%rbp), %rax
movq (%rax), %rcx
shlq $0x20, %rcx
movq %rcx, (%rax)
movq -0x60(%rbp), %rcx
movq -0x70(%rbp), %rax
orq (%rax), %rcx
movq %rcx, (%rax)
jmp 0x362b85
movq $-0x1, -0x50(%rbp)
jmp 0x362b87
movq -0x40(%rbp), %rax
movq -0x58(%rbp), %rcx
cmpb $0x0, (%rax,%rcx)
je 0x362c09
movq -0x20(%rbp), %rdi
movq -0x28(%rbp), %rsi
movq -0x30(%rbp), %rdx
movq -0x38(%rbp), %rcx
movq -0x48(%rbp), %r8
movq 0x10(%rbp), %r9
movq -0x60(%rbp), %r15
movq -0x58(%rbp), %r14
movq 0x28(%rbp), %rbx
movq 0x30(%rbp), %r11
movq 0x18(%rbp), %r10
movq 0x20(%rbp), %rax
movsd 0x1f6eb(%rip), %xmm0    # 0x3822b8
movq %r15, (%rsp)
movq %r14, 0x8(%rsp)
movq %rbx, 0x10(%rsp)
movq %r11, 0x18(%rsp)
movq %r10, 0x20(%rsp)
movq %rax, 0x28(%rsp)
callq 0x362470
cmpq $0x0, %rax
je 0x362bff
movq -0x60(%rbp), %rax
movq %rax, -0x50(%rbp)
jmp 0x362c07
movq $-0x1, -0x50(%rbp)
jmp 0x362c09
jmp 0x362c0b
movq -0x60(%rbp), %rax
addq $0x1, %rax
movq %rax, -0x60(%rbp)
jmp 0x3629a1
addq $0x98, %rsp
popq %rbx
popq %r14
popq %r15
popq %rbp
retq
nopw (%rax,%rax)