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)
L142:
  mov    -0x38(%rsp),%rax
  mov    -0x28(%rsp),%r13
  mov    -0x20(%rsp),%r15
  mov    -0x18(%rsp),%r14
  lea    -0x100(%rax),%rbx
L29:
  mov    %r15,%rdx
  mov    %r14,%rax
  xor    %r9d,%r9d
  jmp    L196
  xchg   %ax,%ax
L61:
  mov    %r10d,%r12d
  sub    %ecx,%r12d
  mov    %r12b,0x1(%rax)
  lea    (%r8,%r10,1),%r12d
  cmp    $0xff,%r12d
  ja     L55
L92:
  add    %r8d,%r10d
  mov    %r10b,0x2(%rax)
L41:
  movb   $0x0,0x3(%rax)
  mov    0x4(%rbx,%rdi,4),%edi
  lea    (%rsi,%rdi,1),%r10d
  cmp    $0xff,%r10d
  ja     L106
  lea    (%rdi,%rsi,1),%r10d
  mov    %r10b,0x4(%rax)
L72:
  mov    %edi,%r10d
  sub    %ecx,%r10d
  cmp    $0xff,%r10d
  ja     L102
  mov    %edi,%r10d
  sub    %ecx,%r10d
  mov    %r10b,0x5(%rax)
L169:
  lea    (%r8,%rdi,1),%r10d
  cmp    $0xff,%r10d
  ja     L250
  add    %r8d,%edi
  mov    %dil,0x6(%rax)
L200:
  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     L266
  lea    (%rdi,%rsi,1),%r11d
  mov    %r11b,(%rdx)
L60:
  mov    %edi,%r11d
  sub    %ecx,%r11d
  cmp    $0xff,%r11d
  ja     L141
  mov    %edi,%r11d
  sub    %ecx,%r11d
  mov    %r11b,0x1(%rdx)
L100:
  lea    (%r8,%rdi,1),%r11d
  cmp    $0xff,%r11d
  ja     L197
  add    %r8d,%edi
  mov    %dil,0x2(%rdx)
L56:
  movb   $0x0,0x3(%rdx)
  mov    0x4(%rbx,%r10,4),%edi
  lea    (%rsi,%rdi,1),%r10d
  cmp    $0xff,%r10d
  ja     L254
  add    %edi,%esi
  mov    %sil,0x4(%rdx)
L248:
  mov    %edi,%esi
  sub    %ecx,%esi
  cmp    $0xff,%esi
  ja     L243
  mov    %edi,%esi
  sub    %ecx,%esi
  mov    %sil,0x5(%rdx)
L236:
  lea    (%r8,%rdi,1),%ecx
  cmp    $0xff,%ecx
  ja     L91
  lea    (%rdi,%r8,1),%ecx
L172:
  add    $0x1,%r9
  mov    %cl,0x6(%rdx)
  add    $0x8,%rax
  add    $0x8,%rdx
  movb   $0x0,-0x1(%rdx)
  cmp    $0x8,%r9
  je     L87
L196:
  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     L121
  lea    (%rsi,%r10,1),%r12d
  mov    %r12b,(%rax)
L126:
  mov    %r10d,%r12d
  sub    %ecx,%r12d
  cmp    $0xff,%r12d
  jbe    L61
  not    %r12d
  sar    $0x1f,%r12d
  mov    %r12b,0x1(%rax)
  lea    (%r8,%r10,1),%r12d
  cmp    $0xff,%r12d
  jbe    L92
L55:
  not    %r12d
  sar    $0x1f,%r12d
  mov    %r12b,0x2(%rax)
  jmp    L41
  nopl   0x0(%rax)
L121:
  not    %r12d
  sar    $0x1f,%r12d
  mov    %r12b,(%rax)
  jmp    L126
  nopl   0x0(%rax)
L91:
  not    %ecx
  sar    $0x1f,%ecx
  jmp    L172
  nopw   0x0(%rax,%rax,1)
L243:
  not    %esi
  sar    $0x1f,%esi
  mov    %sil,0x5(%rdx)
  jmp    L236
  xchg   %ax,%ax
L254:
  not    %r10d
  sar    $0x1f,%r10d
  mov    %r10b,0x4(%rdx)
  jmp    L248
L197:
  not    %r11d
  sar    $0x1f,%r11d
  mov    %r11b,0x2(%rdx)
  jmp    L56
L141:
  not    %r11d
  sar    $0x1f,%r11d
  mov    %r11b,0x1(%rdx)
  jmp    L100
L266:
  not    %r11d
  sar    $0x1f,%r11d
  mov    %r11b,(%rdx)
  jmp    L60
  nop
L250:
  not    %r10d
  sar    $0x1f,%r10d
  mov    %r10b,0x6(%rax)
  jmp    L200
L102:
  not    %r10d
  sar    $0x1f,%r10d
  mov    %r10b,0x5(%rax)
  jmp    L169
L106:
  not    %r10d
  sar    $0x1f,%r10d
  mov    %r10b,0x4(%rax)
  jmp    L72
L87:
  mov    -0x30(%rsp),%rax
  add    $0x20,%r13
  add    $0x40,%rbx
  add    %rax,%r14
  add    %rax,%r15
  cmp    -0x38(%rsp),%rbx
  jne    L29
  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    -0x10(%rsp),%rax
  jne    L142
  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;
  }
}
```