user: 请从给定的汇编指令推导出C源代码：

```asm
<nussinov>:
  endbr64
  push   %r15
  push   %r14
  push   %r13
  push   %r12
  push   %rbp
  push   %rbx
  sub    $0xa8,%rsp
  mov    %rdi,0x50(%rsp)
  mov    %esi,%ebp
  mov    %esi,0x30(%rsp)
  mov    %rdx,0x90(%rsp)
  mov    %rcx,0x98(%rsp)
  mov    %esi,%edi
  call   <initialize_square_table@plt>
  mov    %rax,0x28(%rsp)
  mov    %ebp,%edi
  call   <initialize_square_table@plt>
  mov    %rax,%r15
  mov    %ebp,%edi
  shl    $0x2,%edi
  call   <malloc@plt>
  mov    %rax,%r14
  mov    $0x4,%edi
  call   <malloc@plt>
  mov    %rax,0x10(%rsp)
  mov    $0x4,%edi
  call   <malloc@plt>
  test   %ebp,%ebp
  jle    L22
  movslq %ebp,%rbx
  mov    %rax,%r13
  mov    %ebp,%esi
  neg    %esi
  mov    %esi,0x64(%rsp)
  mov    %rbx,%rax
  neg    %rax
  shl    $0x2,%rax
  mov    %rax,0x88(%rsp)
  movslq %esi,%rax
  mov    %rax,%rdx
  neg    %rdx
  lea    0x0(,%rdx,4),%rsi
  mov    %rsi,0x68(%rsp)
  shl    $0x2,%rbx
  mov    %rbx,(%rsp)
  add    $0x4,%rbx
  mov    %rbx,0x78(%rsp)
  mov    %ebp,%ebx
  mov    %rbx,0x80(%rsp)
  lea    -0x4(,%rax,4),%rbx
  mov    %rbx,0x40(%rsp)
  shl    $0x2,%rax
  mov    %rax,0x48(%rsp)
  movq   $0x0,0x58(%rsp)
  movq   $0x0,0x70(%rsp)
  movq   $0x0,0x20(%rsp)
  lea    -0x1(%rbp),%eax
  mov    %eax,0x60(%rsp)
  jmp    L71
L82:
  add    0x40(%rsp),%rbx
  mov    0x30(%rsp),%edi
  sub    %edi,0x18(%rsp)
  add    0x48(%rsp),%rbp
  cmp    0xc(%rsp),%r12d
  je     L62
  mov    0x1c(%rsp),%r12d
L48:
  mov    0xc(%rsp),%esi
  sub    %r12d,%esi
  mov    0x34(%rsp),%edx
  mov    0x50(%rsp),%rdi
  call   <get_binding_score@plt>
  mov    0x38(%rsp),%rdi
  add    -0x4(%rdi,%rbp,1),%eax
  mov    %eax,(%r15,%rbp,1)
  test   %r12d,%r12d
  jle    L14
  movslq 0x18(%rsp),%rax
  add    0x20(%rsp),%rax
  lea    (%r15,%rax,4),%rcx
  lea    0x4(%r14),%rax
  lea    -0x1(%r12),%edx
  lea    0x8(%r14,%rdx,4),%rsi
L76:
  mov    (%rcx),%edx
  add    -0x4(%rbx,%rax,1),%edx
  mov    %edx,(%rax)
  add    (%rsp),%rcx
  add    $0x4,%rax
  cmp    %rsi,%rax
  jne    L76
L14:
  mov    (%r15,%rbp,1),%eax
  mov    %eax,(%r14)
  lea    0x1(%r12),%eax
  mov    %eax,0x1c(%rsp)
  mov    0x10(%rsp),%rcx
  mov    %r13,%rdx
  mov    %eax,%esi
  mov    %r14,%rdi
  call   <compare@plt>
  mov    0x0(%r13),%eax
  mov    %eax,(%r15,%rbp,1)
  mov    0x10(%rsp),%rax
  mov    (%rax),%eax
  mov    0x28(%rsp),%rdi
  mov    %eax,(%rdi,%rbp,1)
  test   %eax,%eax
  jne    L82
  cmpl   $0x0,(%r14)
  je     L82
  movl   $0xffffffff,(%rdi,%rbp,1)
  jmp    L82
L62:
  addq   $0x1,0x20(%rsp)
  mov    0x20(%rsp),%rax
  mov    0x30(%rsp),%ebx
  add    %ebx,0x64(%rsp)
  mov    0x88(%rsp),%rbx
  add    %rbx,0x68(%rsp)
  mov    (%rsp),%rbx
  add    %rbx,0x70(%rsp)
  mov    0x78(%rsp),%rbx
  add    %rbx,0x58(%rsp)
  cmp    0x80(%rsp),%rax
  je     L22
L71:
  mov    0x20(%rsp),%rbp
  mov    %ebp,0xc(%rsp)
  mov    %ebp,0x34(%rsp)
  mov    0x58(%rsp),%rax
  movl   $0x0,(%r15,%rax,1)
  mov    %ebp,%ebx
  add    $0x1,%ebx
  mov    0x60(%rsp),%esi
  mov    %ebx,%edi
  call   <min@plt>
  mov    0x30(%rsp),%r12d
  imul   %r12d,%eax
  add    %ebp,%eax
  cltq
  movl   $0x0,(%r15,%rax,4)
  mov    0x60(%rsp),%esi
  mov    %ebx,%edi
  call   <min@plt>
  imul   %r12d,%eax
  add    %ebp,%eax
  cltq
  mov    0x28(%rsp),%rdi
  movl   $0xfffffffd,(%rdi,%rax,4)
  mov    0x58(%rsp),%rax
  movl   $0xfffffffe,(%rdi,%rax,1)
  mov    %ebp,%eax
  sub    $0x1,%eax
  js     L62
  mov    0x64(%rsp),%ecx
  movslq %ecx,%rax
  lea    -0x1(%rax,%rbp,1),%rdx
  lea    (%r15,%rdx,4),%rbx
  sub    %r14,%rbx
  add    %r12d,%ecx
  mov    %ecx,0x18(%rsp)
  add    %rbp,%rax
  lea    0x0(,%rax,4),%rbp
  mov    $0x1,%r12d
  mov    0x68(%rsp),%rax
  add    0x70(%rsp),%rax
  add    %r15,%rax
  mov    %rax,0x38(%rsp)
  jmp    L48
L22:
  mov    0x98(%rsp),%rax
  mov    %r15,(%rax)
  mov    0x90(%rsp),%rax
  mov    0x28(%rsp),%rdi
  mov    %rdi,(%rax)
  add    $0xa8,%rsp
  pop    %rbx
  pop    %rbp
  pop    %r12
  pop    %r13
  pop    %r14
  pop    %r15
  ret
```

assistant: ```c
void nussinov(char *sequence, int sequence_length, int **traceback_table, int **score_table) {
  int *traceback = initialize_square_table(sequence_length);
  int *score = initialize_square_table(sequence_length);
  int i, j, k;
  int cells_to_compare;
  int *comparison_list = (int *)malloc(sequence_length * sizeof(int));
  int *arg = (int *)malloc(sizeof(int));
  int *val = (int *)malloc(sizeof(int));
  for (j = 0; j < sequence_length; j++) {
    score[j * sequence_length + j] = 0;
    score[min(j + 1, sequence_length - 1) * sequence_length + j] = 0;
    for (i = 0; i < 0; i++) {
      score[max(0, j - 1 - i) * sequence_length + j] = 0;
      traceback[max(0, j - 1 - i) * sequence_length + j] = 0;
    }
    traceback[min(j + 1, sequence_length - 1) * sequence_length + j] = -3;
    traceback[j * sequence_length + j] = -2;
    for (i = j - 1 - 0; i >= 0; i--) {
      score[i * sequence_length + j] = get_binding_score(sequence, i, j) + score[(i + 1) * sequence_length + (j - 1)];
      cells_to_compare = (j - i);
      for (k = 1; k <= cells_to_compare; k++) {
        comparison_list[k] = score[i * sequence_length + (i + k - 1)] + score[(i + k) * sequence_length + j];
      }
      comparison_list[0] = score[i * sequence_length + j];
      compare(comparison_list, cells_to_compare + 1, val, arg);
      score[i * sequence_length + j] = *val;
      traceback[i * sequence_length + j] = *arg;
      if (*arg == 0) {
        if (comparison_list[0] == 0)
          traceback[i * sequence_length + j] = 0;
        else
          traceback[i * sequence_length + j] = -1;
      }
    }
  }
  *score_table = score;
  *traceback_table = traceback;
}
```