user: Please show me the C implementation that matches the assembly code listed below:

```asm
<col221111_32>:
  endbr64
  push   %r15
  movslq %edx,%rax
  add    %edx,%edx
  push   %r14
  add    %rsi,%rax
  push   %r13
  push   %r12
  push   %rbp
  push   %rbx
  lea    0x100(%rdi),%rbx
  mov    %rax,-0x20(%rsp)
  lea    0x400(%rdi),%rax
  mov    %rax,-0x28(%rsp)
  movslq %edx,%rax
  mov    %rax,-0x30(%rsp)
  shl    $0x2,%rax
  mov    %rbx,-0x38(%rsp)
  lea    0x500(%rdi),%rbx
  mov    %rsi,-0x18(%rsp)
  mov    %rbx,-0x10(%rsp)
  mov    %rax,-0x8(%rsp)
L100:
  mov    -0x38(%rsp),%rax
  mov    -0x28(%rsp),%r13
  mov    -0x20(%rsp),%r15
  mov    -0x18(%rsp),%r14
  lea    -0x100(%rax),%rbx
L210:
  mov    %r15,%rdx
  mov    %r14,%rax
  xor    %r9d,%r9d
  jmp    L96
  xchg   %ax,%ax
L197:
  mov    %r10d,%r12d
  sub    %ecx,%r12d
  mov    %r12b,0x1(%rax)
  lea    (%r8,%r10,1),%r12d
  cmp    $0xff,%r12d
  ja     L187
L58:
  add    %r8d,%r10d
  mov    %r10b,0x2(%rax)
L241:
  movb   $0x0,0x3(%rax)
  mov    0x4(%rbx,%rdi,4),%edi
  lea    (%rsi,%rdi,1),%r10d
  cmp    $0xff,%r10d
  ja     L262
  lea    (%rdi,%rsi,1),%r10d
  mov    %r10b,0x4(%rax)
L275:
  mov    %edi,%r10d
  sub    %ecx,%r10d
  cmp    $0xff,%r10d
  ja     L123
  mov    %edi,%r10d
  sub    %ecx,%r10d
  mov    %r10b,0x5(%rax)
L27:
  lea    (%r8,%rdi,1),%r10d
  cmp    $0xff,%r10d
  ja     L81
  add    %r8d,%edi
  mov    %dil,0x6(%rax)
L257:
  lea    0x1(,%rbp,8),%edi
  movb   $0x0,0x7(%rax)
  lea    (%r11,%rdi,4),%r10d
  add    %r10d,%r10d
  movslq %r10d,%r10
  mov    (%rbx,%r10,4),%edi
  lea    (%rsi,%rdi,1),%r11d
  cmp    $0xff,%r11d
  ja     L212
  lea    (%rdi,%rsi,1),%r11d
  mov    %r11b,(%rdx)
L103:
  mov    %edi,%r11d
  sub    %ecx,%r11d
  cmp    $0xff,%r11d
  ja     L66
  mov    %edi,%r11d
  sub    %ecx,%r11d
  mov    %r11b,0x1(%rdx)
L276:
  lea    (%r8,%rdi,1),%r11d
  cmp    $0xff,%r11d
  ja     L115
  add    %r8d,%edi
  mov    %dil,0x2(%rdx)
L124:
  movb   $0x0,0x3(%rdx)
  mov    0x4(%rbx,%r10,4),%edi
  lea    (%rsi,%rdi,1),%r10d
  cmp    $0xff,%r10d
  ja     L219
  add    %edi,%esi
  mov    %sil,0x4(%rdx)
L137:
  mov    %edi,%esi
  sub    %ecx,%esi
  cmp    $0xff,%esi
  ja     L138
  mov    %edi,%esi
  sub    %ecx,%esi
  mov    %sil,0x5(%rdx)
L10:
  lea    (%r8,%rdi,1),%ecx
  cmp    $0xff,%ecx
  ja     L112
  lea    (%rdi,%r8,1),%ecx
  mov    %cl,0x6(%rdx)
L9:
  add    $0x1,%r9
  movb   $0x0,0x7(%rdx)
  add    $0x8,%rax
  add    $0x8,%rdx
  cmp    $0x8,%r9
  je     L193
L96:
  mov    0x100(%r13,%r9,4),%esi
  mov    0x0(%r13,%r9,4),%r8d
  mov    %r9d,%ebp
  mov    %r9d,%r11d
  sar    $0x2,%ebp
  and    $0x3,%r11d
  imul   $0x32,%r8d,%edi
  imul   $0x82,%esi,%ecx
  lea    0x80(%rdi,%rcx,1),%ecx
  mov    %ebp,%edi
  shl    $0x5,%edi
  sar    $0x8,%ecx
  add    %r11d,%edi
  add    %edi,%edi
  movslq %edi,%rdi
  mov    (%rbx,%rdi,4),%r10d
  lea    (%rsi,%r10,1),%r12d
  cmp    $0xff,%r12d
  ja     L31
  lea    (%rsi,%r10,1),%r12d
  mov    %r12b,(%rax)
L258:
  mov    %r10d,%r12d
  sub    %ecx,%r12d
  cmp    $0xff,%r12d
  jbe    L197
  not    %r12d
  sar    $0x1f,%r12d
  mov    %r12b,0x1(%rax)
  lea    (%r8,%r10,1),%r12d
  cmp    $0xff,%r12d
  jbe    L58
L187:
  not    %r12d
  sar    $0x1f,%r12d
  mov    %r12b,0x2(%rax)
  jmp    L241
  nopl   0x0(%rax)
L31:
  not    %r12d
  sar    $0x1f,%r12d
  mov    %r12b,(%rax)
  jmp    L258
  nopl   0x0(%rax)
L112:
  not    %ecx
  sar    $0x1f,%ecx
  mov    %cl,0x6(%rdx)
  jmp    L9
  nopl   (%rax)
L138:
  not    %esi
  sar    $0x1f,%esi
  mov    %sil,0x5(%rdx)
  jmp    L10
  xchg   %ax,%ax
L219:
  not    %r10d
  sar    $0x1f,%r10d
  mov    %r10b,0x4(%rdx)
  jmp    L137
L115:
  not    %r11d
  sar    $0x1f,%r11d
  mov    %r11b,0x2(%rdx)
  jmp    L124
L66:
  not    %r11d
  sar    $0x1f,%r11d
  mov    %r11b,0x1(%rdx)
  jmp    L276
L212:
  not    %r11d
  sar    $0x1f,%r11d
  mov    %r11b,(%rdx)
  jmp    L103
  nop
L81:
  not    %r10d
  sar    $0x1f,%r10d
  mov    %r10b,0x6(%rax)
  jmp    L257
L123:
  not    %r10d
  sar    $0x1f,%r10d
  mov    %r10b,0x5(%rax)
  jmp    L27
L262:
  not    %r10d
  sar    $0x1f,%r10d
  mov    %r10b,0x4(%rax)
  jmp    L275
L193:
  mov    -0x30(%rsp),%rax
  add    $0x20,%r13
  add    $0x40,%rbx
  add    %rax,%r14
  add    %rax,%r15
  cmp    -0x38(%rsp),%rbx
  jne    L210
  lea    0x200(%rbx),%rax
  mov    -0x8(%rsp),%rdx
  subq   $0xffffffffffffff80,-0x28(%rsp)
  add    %rdx,-0x18(%rsp)
  add    %rdx,-0x20(%rsp)
  mov    %rax,-0x38(%rsp)
  cmp    %rax,-0x10(%rsp)
  jne    L100
  pop    %rbx
  pop    %rbp
  pop    %r12
  pop    %r13
  pop    %r14
  pop    %r15
  ret
```

assistant: ```c
void col221111_32(int *out, unsigned char *pic, int width) {
  int i, j, k;
  unsigned char *pic0, *pic1;
  int *outy, *outc;
  int cr, cg, cb, y;
  pic0 = pic;
  pic1 = pic + width;
  outy = out;
  outc = out + 64 * 4;
  for (i = 2; i > 0; i--) {
    for (j = 4; j > 0; j--) {
      for (k = 0; k < 8; k++)
        ((cb = outc[0 + 0 * 8 + k], cr = outc[64 + 0 * 8 + k], cg = (50 * cb + 130 * cr + 128) >> 8), (y = outy[(k / 4 * 8 + 0) * 8 + (k & 3) * 2 + 0], ((pic0[(k * 2 + 0) * 4 + 0]) = (y + cr), (unsigned int)(y + cr) >= 256 ? ((pic0[(k * 2 + 0) * 4 + 0]) = (y + cr) < 0 ? 0 : 255) : 0), ((pic0[(k * 2 + 0) * 4 + 1]) = (y - cg), (unsigned int)(y - cg) >= 256 ? ((pic0[(k * 2 + 0) * 4 + 1]) = (y - cg) < 0 ? 0 : 255) : 0), ((pic0[(k * 2 + 0) * 4 + 2]) = (y + cb), (unsigned int)(y + cb) >= 256 ? ((pic0[(k * 2 + 0) * 4 + 2]) = (y + cb) < 0 ? 0 : 255) : 0), pic0[(k * 2 + 0) * 4 + 3] = 0), (y = outy[(k / 4 * 8 + 0) * 8 + (k & 3) * 2 + 1], ((pic0[(k * 2 + 1) * 4 + 0]) = (y + cr), (unsigned int)(y + cr) >= 256 ? ((pic0[(k * 2 + 1) * 4 + 0]) = (y + cr) < 0 ? 0 : 255) : 0), ((pic0[(k * 2 + 1) * 4 + 1]) = (y - cg), (unsigned int)(y - cg) >= 256 ? ((pic0[(k * 2 + 1) * 4 + 1]) = (y - cg) < 0 ? 0 : 255) : 0), ((pic0[(k * 2 + 1) * 4 + 2]) = (y + cb), (unsigned int)(y + cb) >= 256 ? ((pic0[(k * 2 + 1) * 4 + 2]) = (y + cb) < 0 ? 0 : 255) : 0), pic0[(k * 2 + 1) * 4 + 3] = 0), (y = outy[(k / 4 * 8 + 1) * 8 + (k & 3) * 2 + 0], ((pic1[(k * 2 + 0) * 4 + 0]) = (y + cr), (unsigned int)(y + cr) >= 256 ? ((pic1[(k * 2 + 0) * 4 + 0]) = (y + cr) < 0 ? 0 : 255) : 0), ((pic1[(k * 2 + 0) * 4 + 1]) = (y - cg), (unsigned int)(y - cg) >= 256 ? ((pic1[(k * 2 + 0) * 4 + 1]) = (y - cg) < 0 ? 0 : 255) : 0), ((pic1[(k * 2 + 0) * 4 + 2]) = (y + cb), (unsigned int)(y + cb) >= 256 ? ((pic1[(k * 2 + 0) * 4 + 2]) = (y + cb) < 0 ? 0 : 255) : 0), pic1[(k * 2 + 0) * 4 + 3] = 0), (y = outy[(k / 4 * 8 + 1) * 8 + (k & 3) * 2 + 1], ((pic1[(k * 2 + 1) * 4 + 0]) = (y + cr), (unsigned int)(y + cr) >= 256 ? ((pic1[(k * 2 + 1) * 4 + 0]) = (y + cr) < 0 ? 0 : 255) : 0), ((pic1[(k * 2 + 1) * 4 + 1]) = (y - cg), (unsigned int)(y - cg) >= 256 ? ((pic1[(k * 2 + 1) * 4 + 1]) = (y - cg) < 0 ? 0 : 255) : 0), ((pic1[(k * 2 + 1) * 4 + 2]) = (y + cb), (unsigned int)(y + cb) >= 256 ? ((pic1[(k * 2 + 1) * 4 + 2]) = (y + cb) < 0 ? 0 : 255) : 0), pic1[(k * 2 + 1) * 4 + 3] = 0));
      outc += 8;
      outy += 16;
      pic0 += 2 * width;
      pic1 += 2 * width;
    }
    outy += 64 * 2 - 16 * 4;
  }
}
```