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

```asm
<sha3_keccakf>:
  endbr64
  push   %r13
  push   %r12
  push   %rbp
  push   %rbx
  sub    $0x158,%rsp
  mov    %fs:0x28,%rax
  mov    %rax,0x148(%rsp)
  xor    %eax,%eax
  movl   $0x1,0x20(%rsp)
  movl   $0x8082,0x24(%rsp)
  movl   $0x808a,0x28(%rsp)
  movl   $0x80008000,0x2c(%rsp)
  movl   $0x808b,0x30(%rsp)
  movl   $0x80000001,0x34(%rsp)
  movl   $0x80008081,0x38(%rsp)
  movl   $0x8009,0x3c(%rsp)
  movl   $0x8a,0x40(%rsp)
  movl   $0x88,0x44(%rsp)
  movl   $0x80008009,0x48(%rsp)
  movl   $0x8000000a,0x4c(%rsp)
  movl   $0x8000808b,0x50(%rsp)
  movl   $0x8b,0x54(%rsp)
  movl   $0x8089,0x58(%rsp)
  movl   $0x8003,0x5c(%rsp)
  movl   $0x8002,0x60(%rsp)
  movl   $0x80,0x64(%rsp)
  movl   $0x800a,0x68(%rsp)
  movl   $0x8000000a,0x6c(%rsp)
  movl   $0x80008081,0x70(%rsp)
  movl   $0x8080,0x74(%rsp)
  movl   $0x80000001,0x78(%rsp)
  movl   $0x80008008,0x7c(%rsp)
  movl   $0x1,0x80(%rsp)
  movl   $0x3,0x84(%rsp)
  movl   $0x6,0x88(%rsp)
  movl   $0xa,0x8c(%rsp)
  movl   $0xf,0x90(%rsp)
  movl   $0x15,0x94(%rsp)
  movl   $0x1c,0x98(%rsp)
  movl   $0x24,0x9c(%rsp)
  movl   $0x2d,0xa0(%rsp)
  movl   $0x37,0xa4(%rsp)
  movl   $0x2,0xa8(%rsp)
  movl   $0xe,0xac(%rsp)
  movl   $0x1b,0xb0(%rsp)
  movl   $0x29,0xb4(%rsp)
  movl   $0x38,0xb8(%rsp)
  movl   $0x8,0xbc(%rsp)
  movl   $0x19,0xc0(%rsp)
  movl   $0x2b,0xc4(%rsp)
  movl   $0x3e,0xc8(%rsp)
  movl   $0x12,0xcc(%rsp)
  movl   $0x27,0xd0(%rsp)
  movl   $0x3d,0xd4(%rsp)
  movl   $0x14,0xd8(%rsp)
  movl   $0x2c,0xdc(%rsp)
  movl   $0xa,0xe0(%rsp)
  movl   $0x7,0xe4(%rsp)
  movl   $0xb,0xe8(%rsp)
  movl   $0x11,0xec(%rsp)
  movl   $0x12,0xf0(%rsp)
  movl   $0x3,0xf4(%rsp)
  movl   $0x5,0xf8(%rsp)
  movl   $0x10,0xfc(%rsp)
  movl   $0x8,0x100(%rsp)
  movl   $0x15,0x104(%rsp)
  movl   $0x18,0x108(%rsp)
  movl   $0x4,0x10c(%rsp)
  movl   $0xf,0x110(%rsp)
  movl   $0x17,0x114(%rsp)
  movl   $0x13,0x118(%rsp)
  movl   $0xd,0x11c(%rsp)
  movl   $0xc,0x120(%rsp)
  movl   $0x2,0x124(%rsp)
  movl   $0x14,0x128(%rsp)
  movl   $0xe,0x12c(%rsp)
  movl   $0x16,0x130(%rsp)
  movl   $0x9,0x134(%rsp)
  movl   $0x6,0x138(%rsp)
  movl   $0x1,0x13c(%rsp)
  lea    0x20(%rsp),%r10
  lea    0x80(%rsp),%r11
  lea    0x14(%rdi),%r9
  mov    $0x40,%r8d
  jmp    L23
L104:
  mov    (%r10),%eax
  xor    %eax,(%rdi)
  add    $0x4,%r10
  cmp    %r11,%r10
  je     L1
L23:
  mov    %rdi,%rsi
  mov    %rsp,%rcx
  mov    %rdi,%rdx
L96:
  mov    (%rdx),%eax
  xor    0x14(%rdx),%eax
  xor    0x28(%rdx),%eax
  xor    0x3c(%rdx),%eax
  xor    0x50(%rdx),%eax
  mov    %eax,(%rcx)
  add    $0x4,%rdx
  add    $0x4,%rcx
  cmp    %r9,%rdx
  jne    L96
  lea    0x64(%rdi),%r13
  mov    $0x0,%ebx
L12:
  lea    0x4(%rbx),%edx
  movslq %edx,%rax
  imul   $0x66666667,%rax,%rax
  sar    $0x21,%rax
  mov    %edx,%ebp
  sar    $0x1f,%ebp
  sub    %ebp,%eax
  lea    (%rax,%rax,4),%eax
  sub    %eax,%edx
  movslq %edx,%rbp
  add    $0x1,%ebx
  movslq %ebx,%rax
  imul   $0x66666667,%rax,%rax
  sar    $0x21,%rax
  mov    %ebx,%edx
  sar    $0x1f,%edx
  sub    %edx,%eax
  lea    (%rax,%rax,4),%edx
  mov    %ebx,%eax
  sub    %edx,%eax
  cltq
  mov    (%rsp,%rax,4),%eax
  lea    (%rax,%rax,1),%edx
  mov    $0x3f,%ecx
  sar    %cl,%eax
  or     %eax,%edx
  xor    (%rsp,%rbp,4),%edx
  lea    -0x64(%r13),%rax
L70:
  xor    %edx,(%rax)
  add    $0x14,%rax
  cmp    %r13,%rax
  jne    L70
  add    $0x4,%r13
  cmp    $0x5,%ebx
  jne    L12
  mov    0x4(%rdi),%r12d
  mov    $0x0,%eax
L84:
  movslq 0xe0(%rsp,%rax,1),%rdx
  lea    (%rdi,%rdx,4),%rbp
  mov    %r12d,%edx
  mov    0x0(%rbp),%r12d
  mov    %r8d,%ecx
  sub    0x80(%rsp,%rax,1),%ecx
  mov    %edx,%ebx
  sar    %cl,%ebx
  movzbl 0x80(%rsp,%rax,1),%ecx
  shl    %cl,%edx
  or     %edx,%ebx
  mov    %ebx,0x0(%rbp)
  add    $0x4,%rax
  cmp    $0x60,%rax
  jne    L84
  mov    %r12d,(%rsp)
  mov    $0x0,%ebx
L75:
  mov    $0x0,%eax
L86:
  mov    (%rsi,%rax,1),%edx
  mov    %edx,(%rsp,%rax,1)
  add    $0x4,%rax
  cmp    $0x14,%rax
  jne    L86
  mov    %rsi,%rcx
  mov    $0x0,%edx
L24:
  mov    %edx,%ebp
  add    $0x1,%edx
  movslq %edx,%rax
  imul   $0x66666667,%rax,%rax
  sar    $0x21,%rax
  mov    %edx,%r12d
  sar    $0x1f,%r12d
  sub    %r12d,%eax
  lea    (%rax,%rax,4),%r12d
  mov    %edx,%eax
  sub    %r12d,%eax
  cltq
  mov    (%rsp,%rax,4),%r12d
  not    %r12d
  add    $0x2,%ebp
  movslq %ebp,%rax
  imul   $0x66666667,%rax,%rax
  sar    $0x21,%rax
  mov    %ebp,%r13d
  sar    $0x1f,%r13d
  sub    %r13d,%eax
  lea    (%rax,%rax,4),%eax
  sub    %eax,%ebp
  movslq %ebp,%rbp
  and    (%rsp,%rbp,4),%r12d
  xor    %r12d,(%rcx)
  add    $0x4,%rcx
  cmp    $0x5,%edx
  jne    L24
  add    $0x5,%ebx
  add    $0x14,%rsi
  cmp    $0x19,%ebx
  jne    L75
  jmp    L104
L1:
  mov    0x148(%rsp),%rax
  sub    %fs:0x28,%rax
  jne    L50
  add    $0x158,%rsp
  pop    %rbx
  pop    %rbp
  pop    %r12
  pop    %r13
  ret
L50:
  call   <__stack_chk_fail@plt>
```

assistant: ```c
void sha3_keccakf(uint64_t st[25]) {
  const uint64_t keccakf_rndc[24] = {0x0000000000000001, 0x0000000000008082, 0x800000000000808a, 0x8000000080008000, 0x000000000000808b, 0x0000000080000001, 0x8000000080008081, 0x8000000000008009, 0x000000000000008a, 0x0000000000000088, 0x0000000080008009, 0x000000008000000a, 0x000000008000808b, 0x800000000000008b, 0x8000000000008089, 0x8000000000008003, 0x8000000000008002, 0x8000000000000080, 0x000000000000800a, 0x800000008000000a, 0x8000000080008081, 0x8000000000008080, 0x0000000080000001, 0x8000000080008008};
  const int keccakf_rotc[24] = {1, 3, 6, 10, 15, 21, 28, 36, 45, 55, 2, 14, 27, 41, 56, 8, 25, 43, 62, 18, 39, 61, 20, 44};
  const int keccakf_piln[24] = {10, 7, 11, 17, 18, 3, 5, 16, 8, 21, 24, 4, 15, 23, 19, 13, 12, 2, 20, 14, 22, 9, 6, 1};
  int i, j, r;
  uint64_t t, bc[5];
  for (r = 0; r < 24; r++) {
    for (i = 0; i < 5; i++)
      bc[i] = st[i] ^ st[i + 5] ^ st[i + 10] ^ st[i + 15] ^ st[i + 20];
    for (i = 0; i < 5; i++) {
      t = bc[(i + 4) % 5] ^ (((bc[(i + 1) % 5]) << (1)) | ((bc[(i + 1) % 5]) >> (64 - (1))));
      for (j = 0; j < 25; j += 5)
        st[j + i] ^= t;
    }
    t = st[1];
    for (i = 0; i < 24; i++) {
      j = keccakf_piln[i];
      bc[0] = st[j];
      st[j] = (((t) << (keccakf_rotc[i])) | ((t) >> (64 - (keccakf_rotc[i]))));
      t = bc[0];
    }
    for (j = 0; j < 25; j += 5) {
      for (i = 0; i < 5; i++)
        bc[i] = st[j + i];
      for (i = 0; i < 5; i++)
        st[j + i] ^= (~bc[(i + 1) % 5]) & bc[(i + 2) % 5];
    }
    st[0] ^= keccakf_rndc[r];
  }
}
```