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

```asm
<julius_regular_involution>:
  endbr64
  push   %r15
  push   %r14
  push   %r13
  push   %r12
  push   %rbp
  push   %rbx
  sub    $0x58,%rsp
  mov    %rdi,%r12
  mov    %rsi,%rbp
  mov    %rdx,%rbx
  mov    %rcx,%r13
  mov    %fs:0x28,%rax
  mov    %rax,0x48(%rsp)
  xor    %eax,%eax
  lea    0x30(%rsp),%rax
  lea    0x40(%rsp),%rdx
L72:
  movb   $0x0,(%rax)
  add    $0x1,%rax
  cmp    %rdx,%rax
  jne    L72
  lea    0x30(%rsp),%r14
  mov    %r13,%rsi
  mov    %r14,%rdi
  call   <aes_128_encrypt@plt>
  mov    %rbp,%rcx
  shr    $0x4,%rcx
  mov    %rsp,%r15
  mov    %r12,%rdx
  mov    %r14,%rsi
  mov    %r15,%rdi
  call   <evaluate_by_delta@plt>
  lea    0x10(%rsp),%r14
  mov    $0x10,%edx
  mov    %r15,%rsi
  mov    %r14,%rdi
  call   <memcpy@plt>
  mov    %r13,%rsi
  mov    %r14,%rdi
  call   <aes_128_encrypt@plt>
  mov    %rbp,%rax
  sub    %rbx,%rax
  test   $0x10,%al
  jne    L57
  mov    %r14,%rsi
  lea    0x20(%rsp),%rdi
  mov    $0x10,%edx
  call   <memcpy@plt>
  cmp    %rbx,%rbp
  jbe    L79
  lea    0x20(%rsp),%r13
L36:
  lea    (%r12,%rbx,1),%rdi
  mov    %r13,%rsi
  call   <add_to_first@plt>
  lea    0x30(%rsp),%rsi
  mov    %r13,%rdi
  call   <mul_by_delta@plt>
  add    $0x10,%rbx
  cmp    %rbx,%rbp
  ja     L36
L79:
  mov    0x48(%rsp),%rax
  sub    %fs:0x28,%rax
  jne    L25
  add    $0x58,%rsp
  pop    %rbx
  pop    %rbp
  pop    %r12
  pop    %r13
  pop    %r14
  pop    %r15
  ret
L57:
  lea    0x10(%rsp),%rsi
  lea    0x20(%rsp),%r13
  mov    $0x10,%edx
  mov    %r13,%rdi
  call   <memcpy@plt>
  lea    (%r12,%rbx,1),%rdi
  mov    %r13,%rsi
  call   <add_to_first@plt>
  lea    0x10(%rbx),%r13
  cmp    %r13,%rbp
  jbe    L79
  lea    0x30(%rsp),%rsi
  lea    0x20(%rsp),%r14
  mov    %r14,%rdi
  call   <mul_by_delta@plt>
  add    %r12,%r13
  mov    %r14,%rsi
  mov    %r13,%rdi
  call   <add_to_first@plt>
  lea    0x10(%rsp),%rsi
  mov    %r13,%rdi
  call   <add_to_first@plt>
  add    $0x20,%rbx
  cmp    %rbx,%rbp
  jbe    L79
  mov    %r14,%r13
L9:
  lea    (%r12,%rbx,1),%rdi
  mov    %r13,%rsi
  call   <add_to_first@plt>
  lea    0x30(%rsp),%rsi
  mov    %r13,%rdi
  call   <mul_by_delta@plt>
  add    $0x10,%rbx
  cmp    %rbx,%rbp
  ja     L9
  jmp    L79
L25:
  call   <__stack_chk_fail@plt>
```

assistant: ```c
void julius_regular_involution(unsigned char *pad, unsigned long long padlen, unsigned long long reg_index, unsigned long expanded_key[44]) {
  int i;
  unsigned long long cur;
  unsigned char seed[16];
  unsigned char mu[16];
  unsigned char tmp_mask[16];
  unsigned char delta[16];
  for (i = 0; i < 16; i++) {
    delta[i] = 0;
  }
  aes_128_encrypt(delta, expanded_key);
  evaluate_by_delta(seed, delta, pad, padlen / 16);
  memcpy(mu, seed, 16);
  aes_128_encrypt(mu, expanded_key);
  if (((padlen - reg_index) / 16) % 2 == 0) {
    cur = reg_index;
    memcpy(tmp_mask, mu, 16);
    for (cur = reg_index; cur < padlen; cur = cur + 16) {
      add_to_first(&pad[cur], tmp_mask);
      mul_by_delta(tmp_mask, delta);
    }
  } else {
    cur = reg_index;
    memcpy(tmp_mask, mu, 16);
    add_to_first(&pad[cur], tmp_mask);
    cur = cur + 16;
    if (cur < padlen) {
      mul_by_delta(tmp_mask, delta);
      add_to_first(&pad[cur], tmp_mask);
      add_to_first(&pad[cur], mu);
      cur = cur + 16;
    }
    for (; cur < padlen; cur = cur + 16) {
      add_to_first(&pad[cur], tmp_mask);
      mul_by_delta(tmp_mask, delta);
    }
  }
}
```