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

```asm
<sha3_keccakf>:
  endbr64
  push   %r15
  mov    %rdi,%r11
  mov    $0xcccccccd,%r9d
  push   %r14
  lea    0x14(%rdi),%r14
  push   %r13
  push   %r12
  push   %rbp
  mov    $0x40,%ebp
  push   %rbx
  lea    0x64(%rdi),%rbx
  sub    $0x178,%rsp
  mov    %fs:0x28,%rax
  mov    %rax,0x168(%rsp)
  xor    %eax,%eax
  lea    0x44(%rsp),%r15
  lea    0x20(%rsp),%r13
  movabs $0x800080000000808a,%rax
  mov    %rax,0x48(%rsp)
  lea    0xa0(%rsp),%r12
  movabs $0x800000010000808b,%rax
  mov    %rax,0x50(%rsp)
  movabs $0x800980008081,%rax
  mov    %rax,0x58(%rsp)
  movabs $0x880000008a,%rax
  mov    %rax,0x60(%rsp)
  movabs $0x8000000a80008009,%rax
  mov    %rax,0x68(%rsp)
  movabs $0x8b8000808b,%rax
  mov    %rax,0x70(%rsp)
  movabs $0x800300008089,%rax
  mov    %rax,0x78(%rsp)
  movabs $0x8000008002,%rax
  mov    %rax,0x80(%rsp)
  movabs $0x8000000a0000800a,%rax
  mov    %rax,0x88(%rsp)
  movabs $0x808080008081,%rax
  mov    %rax,0x90(%rsp)
  movabs $0x8000800880000001,%rax
  mov    %rax,0x98(%rsp)
  movabs $0xa00000006,%rax
  mov    %rax,0xa8(%rsp)
  movabs $0x150000000f,%rax
  mov    %rax,0xb0(%rsp)
  movabs $0x240000001c,%rax
  mov    %rax,0xb8(%rsp)
  movabs $0x370000002d,%rax
  mov    %rax,0xc0(%rsp)
  movabs $0xe00000002,%rax
  mov    %rax,0xc8(%rsp)
  movabs $0x290000001b,%rax
  mov    %rax,0xd0(%rsp)
  movabs $0x800000038,%rax
  mov    %rax,0xd8(%rsp)
  movabs $0x2b00000019,%rax
  mov    %rax,0xe0(%rsp)
  movabs $0x120000003e,%rax
  mov    %rax,0xe8(%rsp)
  movabs $0x3d00000027,%rax
  mov    %rax,0xf0(%rsp)
  movabs $0x2c00000014,%rax
  movl   $0x8082,0x44(%rsp)
  movl   $0x3,0xa4(%rsp)
  mov    %rax,0xf8(%rsp)
  movabs $0x110000000b,%rax
  mov    %rax,0x108(%rsp)
  movabs $0x300000012,%rax
  mov    %rax,0x110(%rsp)
  movabs $0x1000000005,%rax
  mov    %rax,0x118(%rsp)
  movabs $0x1500000008,%rax
  mov    %rax,0x120(%rsp)
  movabs $0x400000018,%rax
  mov    %rax,0x128(%rsp)
  movabs $0x170000000f,%rax
  mov    %rax,0x130(%rsp)
  movabs $0xd00000013,%rax
  mov    %rax,0x138(%rsp)
  movabs $0x20000000c,%rax
  mov    %rax,0x140(%rsp)
  movabs $0xe00000014,%rax
  mov    %rax,0x148(%rsp)
  movabs $0x900000016,%rax
  mov    %rbx,0x18(%rsp)
  lea    0x4(%rdi),%rbx
  mov    %rax,0x150(%rsp)
  movabs $0x100000006,%rax
  movl   $0x7,0x104(%rsp)
  movl   $0x1,0xc(%rsp)
  mov    %rbx,0x10(%rsp)
  mov    %rax,0x158(%rsp)
  mov    (%rdi),%eax
  nopl   0x0(%rax,%rax,1)
L72:
  mov    0x10(%rsp),%rsi
  mov    %r13,%rdi
  mov    %r11,%rbx
  mov    %eax,%ecx
L113:
  mov    0x10(%rsi),%edx
  xor    0x24(%rsi),%edx
  add    $0x4,%rdi
  xor    0x38(%rsi),%edx
  xor    0x4c(%rsi),%edx
  xor    %ecx,%edx
  mov    %edx,-0x4(%rdi)
  cmp    %rsi,%r14
  je     L128
  mov    (%rsi),%ecx
  add    $0x4,%rsi
  jmp    L113
  nopw   0x0(%rax,%rax,1)
L128:
  mov    0x18(%rsp),%rcx
  xor    %edi,%edi
  mov    %rcx,%r10
L13:
  lea    0x4(%rdi),%esi
  mov    $0x3f,%ecx
  mov    %rsi,%rdx
  imul   %r9,%rsi
  shr    $0x22,%rsi
  lea    (%rsi,%rsi,4),%esi
  sub    %esi,%edx
  lea    0x1(%rdi),%esi
  mov    %rsi,%rdi
  imul   %r9,%rsi
  movslq %edx,%rdx
  shr    $0x22,%rsi
  lea    (%rsi,%rsi,4),%r8d
  mov    %edi,%esi
  sub    %r8d,%esi
  movslq %esi,%rsi
  mov    0x20(%rsp,%rsi,4),%r8d
  lea    (%r8,%r8,1),%esi
  sar    %cl,%r8d
  mov    %r10,%rcx
  or     %r8d,%esi
  xor    0x20(%rsp,%rdx,4),%esi
  lea    -0x64(%r10),%rdx
  xor    %esi,%eax
  mov    %eax,(%rdx)
  add    $0x14,%rdx
  cmp    %rdx,%rcx
  je     L25
L86:
  mov    (%rdx),%eax
  add    $0x14,%rdx
  xor    %esi,%eax
  mov    %eax,-0x14(%rdx)
  cmp    %rdx,%rcx
  jne    L86
L25:
  mov    %rcx,%r10
  add    $0x4,%r10
  cmp    $0x5,%edi
  je     L7
  mov    -0x64(%r10),%eax
  jmp    L13
  nopl   0x0(%rax,%rax,1)
L7:
  mov    0x4(%r11),%eax
  xor    %edx,%edx
  mov    $0x1,%edi
  mov    $0xa,%ecx
  jmp    L121
  nopw   0x0(%rax,%rax,1)
L4:
  movslq 0x100(%rsp,%rdx,1),%rcx
  mov    (%r12,%rdx,1),%edi
  mov    %r10d,%eax
L121:
  lea    (%r11,%rcx,4),%r8
  mov    %ebp,%ecx
  mov    %eax,%esi
  add    $0x4,%rdx
  sub    %edi,%ecx
  mov    (%r8),%r10d
  sar    %cl,%esi
  mov    %edi,%ecx
  shl    %cl,%eax
  or     %eax,%esi
  mov    %esi,(%r8)
  cmp    $0x60,%rdx
  jne    L4
  xor    %r8d,%r8d
L32:
  movdqu (%rbx),%xmm0
  mov    %rbx,%rsi
  xor    %edx,%edx
  movaps %xmm0,0x0(%r13)
  mov    0x10(%rbx),%eax
  mov    %eax,0x10(%r13)
L104:
  lea    0x1(%rdx),%ecx
  mov    %edx,%eax
  mov    %rcx,%rdx
  imul   %r9,%rcx
  shr    $0x22,%rcx
  lea    (%rcx,%rcx,4),%edi
  mov    %edx,%ecx
  sub    %edi,%ecx
  lea    0x2(%rax),%edi
  mov    %rdi,%rax
  imul   %r9,%rdi
  movslq %ecx,%rcx
  mov    0x20(%rsp,%rcx,4),%ecx
  shr    $0x22,%rdi
  not    %ecx
  lea    (%rdi,%rdi,4),%edi
  sub    %edi,%eax
  cltq
  and    0x20(%rsp,%rax,4),%ecx
  xor    %ecx,(%rsi)
  add    $0x4,%rsi
  cmp    $0x5,%edx
  jne    L104
  add    $0x5,%r8d
  add    $0x14,%rbx
  cmp    $0x19,%r8d
  jne    L32
  mov    0xc(%rsp),%eax
  xor    (%r11),%eax
  mov    %eax,(%r11)
  cmp    %r15,%r12
  je     L42
  mov    (%r15),%ebx
  add    $0x4,%r15
  mov    %ebx,0xc(%rsp)
  jmp    L72
L42:
  mov    0x168(%rsp),%rax
  sub    %fs:0x28,%rax
  jne    L108
  add    $0x178,%rsp
  pop    %rbx
  pop    %rbp
  pop    %r12
  pop    %r13
  pop    %r14
  pop    %r15
  ret
L108:
  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];
  }
}
```