user: 请显示与下面列出的汇编代码相匹配的C语言实现:

```asm
<MD5_String2binary>:
  endbr64
  push   %r14
  push   %r13
  push   %r12
  push   %rbp
  push   %rbx
  sub    $0x70,%rsp
  mov    %rdi,%rbx
  mov    %rsi,%r14
  mov    %fs:0x28,%rax
  mov    %rax,0x68(%rsp)
  xor    %eax,%eax
  movl   $0x67452301,0x10(%rsp)
  movl   $0xefcdab89,0x14(%rsp)
  movl   $0x98badcfe,0x18(%rsp)
  movl   $0x10325476,0x1c(%rsp)
  call   <strlen@plt>
  mov    %rax,%r13
  cmp    $0x3f,%eax
  jle    L57
  lea    -0x40(%rax),%ebp
  shr    $0x6,%ebp
  add    $0x1,%ebp
  mov    %ebp,%ebp
  shl    $0x6,%rbp
  add    %rbx,%rbp
  lea    0x10(%rsp),%r12
L36:
  lea    0x18(%rsp),%rcx
  lea    0x14(%rsp),%rdx
  lea    0x1c(%rsp),%r8
  mov    %r12,%rsi
  mov    %rbx,%rdi
  call   <MD5_Round_Calculate@plt>
  add    $0x40,%rbx
  cmp    %rbp,%rbx
  jne    L36
L0:
  mov    %r13d,%ebx
  and    $0x3f,%ebx
  lea    0x20(%rsp),%r12
  mov    %ebx,%edx
  mov    %rbp,%rsi
  mov    %r12,%rdi
  call   <strncpy@plt>
  mov    $0x40,%edx
  sub    %ebx,%edx
  mov    %ebx,%ebp
  lea    (%r12,%rbp,1),%rdi
  mov    $0x0,%esi
  call   <memset@plt>
  orb    $0x80,0x20(%rsp,%rbp,1)
  cmp    $0x37,%ebx
  ja     L52
L75:
  lea    0x0(,%r13,8),%eax
  mov    %eax,0x8(%rsp)
  lea    0x8(%rsp),%rsi
  lea    0x58(%rsp),%rdi
  mov    $0x4,%edx
  call   <memcpy@plt>
  cmp    $0x1fffffff,%r13d
  jbe    L22
  lea    0x8(,%r13,8),%eax
  mov    %eax,0xc(%rsp)
  lea    0xc(%rsp),%rsi
  lea    0x5c(%rsp),%rdi
  mov    $0x4,%edx
  call   <memcpy@plt>
L28:
  lea    0x10(%rsp),%rbx
  lea    0x18(%rsp),%rcx
  lea    0x14(%rsp),%rdx
  lea    0x20(%rsp),%rdi
  lea    0x1c(%rsp),%r8
  mov    %rbx,%rsi
  call   <MD5_Round_Calculate@plt>
  mov    $0x10,%edx
  mov    %rbx,%rsi
  mov    %r14,%rdi
  call   <memcpy@plt>
  mov    0x68(%rsp),%rax
  sub    %fs:0x28,%rax
  jne    L5
  add    $0x70,%rsp
  pop    %rbx
  pop    %rbp
  pop    %r12
  pop    %r13
  pop    %r14
  ret
L57:
  mov    %rbx,%rbp
  jmp    L0
L52:
  lea    0x10(%rsp),%rsi
  lea    0x18(%rsp),%rcx
  lea    0x14(%rsp),%rdx
  lea    0x1c(%rsp),%r8
  mov    %r12,%rdi
  call   <MD5_Round_Calculate@plt>
  mov    $0x38,%edx
  mov    $0x0,%esi
  mov    %r12,%rdi
  call   <memset@plt>
  jmp    L75
L22:
  lea    0x5c(%rsp),%rdi
  mov    $0x4,%edx
  mov    $0x0,%esi
  call   <memset@plt>
  jmp    L28
L5:
  call   <__stack_chk_fail@plt>
```

assistant: ```c
void MD5_String2binary(const char *string, unsigned char *output) {
  unsigned char padding_message[64];
  unsigned char *pstring;
  unsigned int string_byte_len, string_bit_len, copy_len, msg_digest[4];
  unsigned int *A = &msg_digest[0], *B = &msg_digest[1], *C = &msg_digest[2], *D = &msg_digest[3];
  int i;
  *A = 0x67452301;
  *B = 0xefcdab89;
  *C = 0x98badcfe;
  *D = 0x10325476;
  string_byte_len = (unsigned int)strlen(string);
  pstring = (unsigned char *)string;
  for (i = string_byte_len; 64 <= i; i -= 64, pstring += 64)
    MD5_Round_Calculate(pstring, A, B, C, D);
  copy_len = string_byte_len % 64;
  strncpy((char *)padding_message, (char *)pstring, copy_len);
  memset(padding_message + copy_len, 0, 64 - copy_len);
  padding_message[copy_len] |= 0x80;
  if (56 <= copy_len) {
    MD5_Round_Calculate(padding_message, A, B, C, D);
    memset(padding_message, 0, 56);
  }
  string_bit_len = string_byte_len * 8;
  memcpy(&padding_message[56], &string_bit_len, 4);
  if (0xffffffff / 8 < string_byte_len) {
    unsigned int high = (string_byte_len - 0xffffffff / 8) * 8;
    memcpy(&padding_message[60], &high, 4);
  } else
    memset(&padding_message[60], 0, 4);
  MD5_Round_Calculate(padding_message, A, B, C, D);
  memcpy(output, msg_digest, 16);
}
```