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

```asm
<bn_sqr_comba4>:
  endbr64
  push   %rbx
  mov    (%rsi),%rax
  mov    %rdi,%rdx
  mov    %rsi,%rcx
  movabs $0x1fffe0000,%r9
  movzwl %ax,%edi
  shr    $0x10,%rax
  movzwl %ax,%eax
  mov    %rdi,%rsi
  imul   %rax,%rsi
  imul   %rdi,%rdi
  imul   %rax,%rax
  mov    %rsi,%r8
  shl    $0x11,%rsi
  and    %r9,%rsi
  shr    $0xf,%r8
  add    %esi,%edi
  mov    %edi,%edi
  cmp    %rsi,%rdi
  mov    %rdi,(%rdx)
  mov    (%rcx),%rsi
  adc    %r8,%rax
  mov    %rax,%r9
  mov    0x8(%rcx),%rax
  movzwl %si,%r8d
  shr    $0x10,%rsi
  movzwl %si,%esi
  movzwl %ax,%r11d
  shr    $0x10,%rax
  movzwl %ax,%edi
  mov    %r11,%r10
  imul   %r8,%r10
  mov    %rdi,%rax
  imul   %rdi,%r8
  imul   %rsi,%r11
  imul   %rsi,%rax
  lea    (%r11,%r8,1),%esi
  cmp    %rsi,%r8
  lea    0x10000(%rax),%rdi
  cmova  %rdi,%rax
  mov    %rsi,%rdi
  shl    $0x10,%rsi
  mov    %esi,%esi
  shr    $0x10,%rdi
  add    %rsi,%r10
  mov    %r10d,%r11d
  cmp    %rsi,%r11
  adc    %rdi,%rax
  and    $0x80000000,%r10d
  mov    %rax,%rsi
  mov    %rax,%r8
  lea    (%r11,%r11,1),%rax
  add    %rsi,%rsi
  shr    $0x1f,%r8
  mov    %esi,%edi
  and    $0x1,%r8d
  cmp    $0x1,%r10
  mov    %eax,%esi
  sbb    $0xffffffffffffffff,%rdi
  add    %r9d,%eax
  mov    %eax,%eax
  cmp    %rax,%rsi
  jbe    L9
  lea    0x1(%rdi),%rsi
  mov    %rsi,%rdi
  and    $0xffffffff,%edi
  je     L77
  cmp    %rsi,%rdi
  jae    L9
L35:
  add    $0x1,%r8
L9:
  movabs $0x1fffe0000,%r10
  mov    %rax,0x8(%rdx)
  mov    0x8(%rcx),%rsi
  movzwl %si,%r9d
  shr    $0x10,%rsi
  movzwl %si,%esi
  mov    %r9,%rax
  imul   %rsi,%rax
  imul   %r9,%r9
  imul   %rsi,%rsi
  mov    %rax,%r11
  shl    $0x11,%rax
  and    %r10,%rax
  shr    $0xf,%r11
  lea    (%r9,%rax,1),%r10d
  cmp    %rax,%r10
  lea    (%r10,%rdi,1),%r9d
  mov    0x10(%rcx),%rax
  adc    %r11,%rsi
  cmp    %r10,%r9
  adc    $0x0,%rsi
  xor    %edi,%edi
  movzwl %ax,%ebx
  add    %esi,%r8d
  mov    %rbx,%r10
  mov    %r8d,%r8d
  cmp    %r8,%rsi
  mov    (%rcx),%rsi
  seta   %dil
  shr    $0x10,%rax
  movzwl %si,%r11d
  shr    $0x10,%rsi
  movzwl %ax,%eax
  movzwl %si,%esi
  imul   %r11,%r10
  imul   %rsi,%rbx
  imul   %rax,%r11
  imul   %rsi,%rax
  lea    (%rbx,%r11,1),%esi
  cmp    %rsi,%r11
  lea    0x10000(%rax),%rbx
  cmova  %rbx,%rax
  mov    %rsi,%rbx
  shl    $0x10,%rsi
  mov    %esi,%esi
  shr    $0x10,%rbx
  add    %rsi,%r10
  mov    %r10d,%r11d
  cmp    %rsi,%r11
  adc    %rbx,%rax
  mov    %rax,%rsi
  and    $0x80000000,%esi
  cmp    $0x1,%rsi
  lea    (%r11,%r11,1),%rsi
  sbb    $0xffffffffffffffff,%rdi
  add    %rax,%rax
  and    $0x80000000,%r10d
  cmp    $0x1,%r10
  mov    %eax,%eax
  mov    %esi,%r10d
  sbb    $0xffffffffffffffff,%rax
  add    %r9d,%esi
  mov    %esi,%esi
  cmp    %rsi,%r10
  jbe    L31
  add    $0x1,%rax
  cmp    $0x1,%eax
  adc    $0x0,%rdi
L31:
  add    %eax,%r8d
  mov    %rsi,0x10(%rdx)
  mov    (%rcx),%rsi
  mov    %r8d,%r8d
  cmp    %rax,%r8
  mov    0x18(%rcx),%rax
  movzwl %si,%r10d
  adc    $0x0,%rdi
  shr    $0x10,%rsi
  movzwl %ax,%r11d
  shr    $0x10,%rax
  movzwl %si,%esi
  movzwl %ax,%eax
  mov    %r11,%r9
  imul   %r10,%r9
  imul   %rsi,%r11
  imul   %rax,%r10
  imul   %rsi,%rax
  lea    (%r11,%r10,1),%esi
  cmp    %rsi,%r10
  lea    0x10000(%rax),%r11
  cmova  %r11,%rax
  mov    %rsi,%r11
  shl    $0x10,%rsi
  mov    %esi,%esi
  shr    $0x10,%r11
  add    %rsi,%r9
  mov    %r9d,%r10d
  cmp    %rsi,%r10
  adc    %r11,%rax
  and    $0x80000000,%r9d
  mov    %rax,%rsi
  add    %rax,%rax
  shr    $0x1f,%rsi
  mov    %eax,%eax
  and    $0x1,%esi
  cmp    $0x1,%r9
  lea    (%r10,%r10,1),%r9
  sbb    $0xffffffffffffffff,%rax
  add    %r9d,%r8d
  mov    %r9d,%r10d
  mov    %r8d,%r8d
  cmp    %r8,%r10
  jbe    L58
  add    $0x1,%rax
  cmp    $0x1,%eax
  adc    $0x0,%rsi
L58:
  lea    (%rax,%rdi,1),%r9d
  mov    0x8(%rcx),%rdi
  cmp    %rax,%r9
  mov    0x10(%rcx),%rax
  adc    $0x0,%rsi
  movzwl %di,%r11d
  shr    $0x10,%rdi
  movzwl %ax,%ebx
  shr    $0x10,%rax
  movzwl %di,%edi
  movzwl %ax,%eax
  mov    %rbx,%r10
  imul   %r11,%r10
  imul   %rdi,%rbx
  imul   %rax,%r11
  imul   %rdi,%rax
  lea    (%rbx,%r11,1),%edi
  cmp    %rdi,%r11
  lea    0x10000(%rax),%rbx
  cmova  %rbx,%rax
  mov    %rdi,%rbx
  shl    $0x10,%rdi
  mov    %edi,%edi
  shr    $0x10,%rbx
  add    %rdi,%r10
  mov    %r10d,%r11d
  cmp    %rdi,%r11
  adc    %rbx,%rax
  mov    %rax,%rdi
  and    $0x80000000,%edi
  cmp    $0x1,%rdi
  lea    (%r11,%r11,1),%rdi
  sbb    $0xffffffffffffffff,%rsi
  add    %rax,%rax
  and    $0x80000000,%r10d
  cmp    $0x1,%r10
  mov    %eax,%eax
  mov    %edi,%r10d
  sbb    $0xffffffffffffffff,%rax
  add    %r8d,%edi
  mov    %edi,%edi
  cmp    %rdi,%r10
  jbe    L32
  add    $0x1,%rax
  cmp    $0x1,%eax
  adc    $0x0,%rsi
L32:
  add    %eax,%r9d
  mov    %rdi,0x18(%rdx)
  movabs $0x1fffe0000,%r10
  mov    %r9d,%r9d
  cmp    %rax,%r9
  mov    0x10(%rcx),%rax
  adc    $0x0,%rsi
  movzwl %ax,%r8d
  shr    $0x10,%rax
  movzwl %ax,%eax
  mov    %r8,%rdi
  imul   %rax,%rdi
  imul   %r8,%r8
  imul   %rax,%rax
  mov    %rdi,%r11
  shl    $0x11,%rdi
  and    %r10,%rdi
  shr    $0xf,%r11
  lea    (%r8,%rdi,1),%r10d
  cmp    %rdi,%r10
  lea    (%r9,%r10,1),%r8d
  adc    %r11,%rax
  cmp    %r10,%r8
  adc    $0x0,%rax
  xor    %edi,%edi
  lea    (%rax,%rsi,1),%r9d
  mov    0x8(%rcx),%rsi
  cmp    %r9,%rax
  mov    0x18(%rcx),%rax
  seta   %dil
  movzwl %si,%r11d
  shr    $0x10,%rsi
  movzwl %ax,%ebx
  shr    $0x10,%rax
  movzwl %si,%esi
  movzwl %ax,%eax
  mov    %rbx,%r10
  imul   %r11,%r10
  imul   %rsi,%rbx
  imul   %rax,%r11
  imul   %rsi,%rax
  lea    (%rbx,%r11,1),%esi
  cmp    %rsi,%r11
  lea    0x10000(%rax),%rbx
  cmova  %rbx,%rax
  mov    %rsi,%rbx
  shl    $0x10,%rsi
  mov    %esi,%esi
  shr    $0x10,%rbx
  add    %rsi,%r10
  mov    %r10d,%r11d
  cmp    %rsi,%r11
  adc    %rbx,%rax
  mov    %rax,%rsi
  and    $0x80000000,%esi
  cmp    $0x1,%rsi
  lea    (%r11,%r11,1),%rsi
  sbb    $0xffffffffffffffff,%rdi
  add    %rax,%rax
  and    $0x80000000,%r10d
  cmp    $0x1,%r10
  mov    %eax,%eax
  mov    %esi,%r10d
  sbb    $0xffffffffffffffff,%rax
  add    %r8d,%esi
  mov    %esi,%esi
  cmp    %rsi,%r10
  jbe    L27
  add    $0x1,%rax
  cmp    $0x1,%eax
  adc    $0x0,%rdi
L27:
  lea    (%rax,%r9,1),%r8d
  mov    %rsi,0x20(%rdx)
  mov    0x10(%rcx),%rsi
  cmp    %rax,%r8
  mov    0x18(%rcx),%rax
  adc    $0x0,%rdi
  movzwl %si,%r10d
  shr    $0x10,%rsi
  movzwl %ax,%r11d
  shr    $0x10,%rax
  movzwl %si,%esi
  movzwl %ax,%eax
  mov    %r11,%r9
  imul   %r10,%r9
  imul   %rsi,%r11
  imul   %rax,%r10
  imul   %rsi,%rax
  lea    (%r11,%r10,1),%esi
  cmp    %rsi,%r10
  lea    0x10000(%rax),%r11
  cmova  %r11,%rax
  mov    %rsi,%r11
  shl    $0x10,%rsi
  mov    %esi,%esi
  shr    $0x10,%r11
  add    %rsi,%r9
  mov    %r9d,%r10d
  cmp    %rsi,%r10
  adc    %r11,%rax
  and    $0x80000000,%r9d
  mov    %rax,%rsi
  add    %rax,%rax
  shr    $0x1f,%rsi
  mov    %eax,%eax
  and    $0x1,%esi
  cmp    $0x1,%r9
  lea    (%r10,%r10,1),%r9
  sbb    $0xffffffffffffffff,%rax
  mov    %r9d,%r10d
  add    %r8d,%r9d
  mov    %r9d,%r9d
  cmp    %r9,%r10
  jbe    L14
  add    $0x1,%rax
  cmp    $0x1,%eax
  adc    $0x0,%rsi
L14:
  lea    (%rax,%rdi,1),%r8d
  mov    %r9,0x28(%rdx)
  pop    %rbx
  movabs $0x1fffe0000,%r10
  cmp    %rax,%r8
  mov    0x18(%rcx),%rax
  adc    $0x0,%rsi
  movzwl %ax,%edi
  shr    $0x10,%rax
  movzwl %ax,%eax
  mov    %rdi,%rcx
  imul   %rax,%rcx
  imul   %rdi,%rdi
  imul   %rax,%rax
  mov    %rcx,%r9
  shl    $0x11,%rcx
  and    %r10,%rcx
  shr    $0xf,%r9
  add    %ecx,%edi
  mov    %edi,%edi
  cmp    %rcx,%rdi
  lea    (%r8,%rdi,1),%ecx
  adc    %r9,%rax
  cmp    %rdi,%rcx
  mov    %rcx,0x30(%rdx)
  adc    %esi,%eax
  mov    %rax,0x38(%rdx)
  ret
  nopl   0x0(%rax,%rax,1)
L77:
  add    $0x1,%r8
  jmp    L35
```

assistant: ```c
void bn_sqr_comba4(unsigned long *r, unsigned long *a) {
  unsigned long bl, bh;
  unsigned long t1, t2;
  unsigned long c1, c2, c3;
  c1 = 0;
  c2 = 0;
  c3 = 0;
  {
    unsigned long l, h, m;
    h = ((a)[0]);
    l = ((h) & (0xffff));
    h = (((h) >> 16) & (0xffff));
    m = (l) * (h);
    l *= l;
    h *= h;
    h += (m & (0xffff8000L)) >> (16 - 1);
    m = (m & (0xffff)) << (16 + 1);
    l = (l + m) & (0xffffffffL);
    if (l < m)
      h++;
    (t1) = l;
    (t2) = h;
  };
  c1 = (c1 + t1) & (0xffffffffL);
  if ((c1) < t1)
    t2++;
  c2 = (c2 + t2) & (0xffffffffL);
  if ((c2) < t2)
    c3++;
  ;
  r[0] = c1;
  c1 = 0;
  t1 = (((a)[1]) & (0xffff));
  t2 = ((((a)[1]) >> 16) & (0xffff));
  bl = (((a)[0]) & (0xffff));
  bh = ((((a)[0]) >> 16) & (0xffff));
  {
    unsigned long m, m1, lt, ht;
    lt = t1;
    ht = t2;
    m = (bh) * (lt);
    lt = (bl) * (lt);
    m1 = (bl) * (ht);
    ht = (bh) * (ht);
    m = (m + m1) & (0xffffffffL);
    if (m < m1)
      ht += ((unsigned long)((1L) & (0xffff)) << 16);
    ht += (((m) >> 16) & (0xffff));
    m1 = ((unsigned long)((m) & (0xffff)) << 16);
    lt = (lt + m1) & (0xffffffffL);
    if (lt < m1)
      ht++;
    (t1) = lt;
    (t2) = ht;
  };
  if (t2 & (0x80000000L))
    c1++;
  t2 = (t2 + t2) & (0xffffffffL);
  if (t1 & (0x80000000L))
    t2++;
  t1 = (t1 + t1) & (0xffffffffL);
  c2 = (c2 + t1) & (0xffffffffL);
  if ((c2 < t1) && (((++t2) & (0xffffffffL)) == 0))
    c1++;
  c3 = (c3 + t2) & (0xffffffffL);
  if ((c3) < t2)
    c1++;
  ;
  r[1] = c2;
  c2 = 0;
  {
    unsigned long l, h, m;
    h = ((a)[1]);
    l = ((h) & (0xffff));
    h = (((h) >> 16) & (0xffff));
    m = (l) * (h);
    l *= l;
    h *= h;
    h += (m & (0xffff8000L)) >> (16 - 1);
    m = (m & (0xffff)) << (16 + 1);
    l = (l + m) & (0xffffffffL);
    if (l < m)
      h++;
    (t1) = l;
    (t2) = h;
  };
  c3 = (c3 + t1) & (0xffffffffL);
  if ((c3) < t1)
    t2++;
  c1 = (c1 + t2) & (0xffffffffL);
  if ((c1) < t2)
    c2++;
  ;
  t1 = (((a)[2]) & (0xffff));
  t2 = ((((a)[2]) >> 16) & (0xffff));
  bl = (((a)[0]) & (0xffff));
  bh = ((((a)[0]) >> 16) & (0xffff));
  {
    unsigned long m, m1, lt, ht;
    lt = t1;
    ht = t2;
    m = (bh) * (lt);
    lt = (bl) * (lt);
    m1 = (bl) * (ht);
    ht = (bh) * (ht);
    m = (m + m1) & (0xffffffffL);
    if (m < m1)
      ht += ((unsigned long)((1L) & (0xffff)) << 16);
    ht += (((m) >> 16) & (0xffff));
    m1 = ((unsigned long)((m) & (0xffff)) << 16);
    lt = (lt + m1) & (0xffffffffL);
    if (lt < m1)
      ht++;
    (t1) = lt;
    (t2) = ht;
  };
  if (t2 & (0x80000000L))
    c2++;
  t2 = (t2 + t2) & (0xffffffffL);
  if (t1 & (0x80000000L))
    t2++;
  t1 = (t1 + t1) & (0xffffffffL);
  c3 = (c3 + t1) & (0xffffffffL);
  if ((c3 < t1) && (((++t2) & (0xffffffffL)) == 0))
    c2++;
  c1 = (c1 + t2) & (0xffffffffL);
  if ((c1) < t2)
    c2++;
  ;
  r[2] = c3;
  c3 = 0;
  t1 = (((a)[3]) & (0xffff));
  t2 = ((((a)[3]) >> 16) & (0xffff));
  bl = (((a)[0]) & (0xffff));
  bh = ((((a)[0]) >> 16) & (0xffff));
  {
    unsigned long m, m1, lt, ht;
    lt = t1;
    ht = t2;
    m = (bh) * (lt);
    lt = (bl) * (lt);
    m1 = (bl) * (ht);
    ht = (bh) * (ht);
    m = (m + m1) & (0xffffffffL);
    if (m < m1)
      ht += ((unsigned long)((1L) & (0xffff)) << 16);
    ht += (((m) >> 16) & (0xffff));
    m1 = ((unsigned long)((m) & (0xffff)) << 16);
    lt = (lt + m1) & (0xffffffffL);
    if (lt < m1)
      ht++;
    (t1) = lt;
    (t2) = ht;
  };
  if (t2 & (0x80000000L))
    c3++;
  t2 = (t2 + t2) & (0xffffffffL);
  if (t1 & (0x80000000L))
    t2++;
  t1 = (t1 + t1) & (0xffffffffL);
  c1 = (c1 + t1) & (0xffffffffL);
  if ((c1 < t1) && (((++t2) & (0xffffffffL)) == 0))
    c3++;
  c2 = (c2 + t2) & (0xffffffffL);
  if ((c2) < t2)
    c3++;
  ;
  t1 = (((a)[2]) & (0xffff));
  t2 = ((((a)[2]) >> 16) & (0xffff));
  bl = (((a)[1]) & (0xffff));
  bh = ((((a)[1]) >> 16) & (0xffff));
  {
    unsigned long m, m1, lt, ht;
    lt = t1;
    ht = t2;
    m = (bh) * (lt);
    lt = (bl) * (lt);
    m1 = (bl) * (ht);
    ht = (bh) * (ht);
    m = (m + m1) & (0xffffffffL);
    if (m < m1)
      ht += ((unsigned long)((1L) & (0xffff)) << 16);
    ht += (((m) >> 16) & (0xffff));
    m1 = ((unsigned long)((m) & (0xffff)) << 16);
    lt = (lt + m1) & (0xffffffffL);
    if (lt < m1)
      ht++;
    (t1) = lt;
    (t2) = ht;
  };
  if (t2 & (0x80000000L))
    c3++;
  t2 = (t2 + t2) & (0xffffffffL);
  if (t1 & (0x80000000L))
    t2++;
  t1 = (t1 + t1) & (0xffffffffL);
  c1 = (c1 + t1) & (0xffffffffL);
  if ((c1 < t1) && (((++t2) & (0xffffffffL)) == 0))
    c3++;
  c2 = (c2 + t2) & (0xffffffffL);
  if ((c2) < t2)
    c3++;
  ;
  r[3] = c1;
  c1 = 0;
  {
    unsigned long l, h, m;
    h = ((a)[2]);
    l = ((h) & (0xffff));
    h = (((h) >> 16) & (0xffff));
    m = (l) * (h);
    l *= l;
    h *= h;
    h += (m & (0xffff8000L)) >> (16 - 1);
    m = (m & (0xffff)) << (16 + 1);
    l = (l + m) & (0xffffffffL);
    if (l < m)
      h++;
    (t1) = l;
    (t2) = h;
  };
  c2 = (c2 + t1) & (0xffffffffL);
  if ((c2) < t1)
    t2++;
  c3 = (c3 + t2) & (0xffffffffL);
  if ((c3) < t2)
    c1++;
  ;
  t1 = (((a)[3]) & (0xffff));
  t2 = ((((a)[3]) >> 16) & (0xffff));
  bl = (((a)[1]) & (0xffff));
  bh = ((((a)[1]) >> 16) & (0xffff));
  {
    unsigned long m, m1, lt, ht;
    lt = t1;
    ht = t2;
    m = (bh) * (lt);
    lt = (bl) * (lt);
    m1 = (bl) * (ht);
    ht = (bh) * (ht);
    m = (m + m1) & (0xffffffffL);
    if (m < m1)
      ht += ((unsigned long)((1L) & (0xffff)) << 16);
    ht += (((m) >> 16) & (0xffff));
    m1 = ((unsigned long)((m) & (0xffff)) << 16);
    lt = (lt + m1) & (0xffffffffL);
    if (lt < m1)
      ht++;
    (t1) = lt;
    (t2) = ht;
  };
  if (t2 & (0x80000000L))
    c1++;
  t2 = (t2 + t2) & (0xffffffffL);
  if (t1 & (0x80000000L))
    t2++;
  t1 = (t1 + t1) & (0xffffffffL);
  c2 = (c2 + t1) & (0xffffffffL);
  if ((c2 < t1) && (((++t2) & (0xffffffffL)) == 0))
    c1++;
  c3 = (c3 + t2) & (0xffffffffL);
  if ((c3) < t2)
    c1++;
  ;
  r[4] = c2;
  c2 = 0;
  t1 = (((a)[3]) & (0xffff));
  t2 = ((((a)[3]) >> 16) & (0xffff));
  bl = (((a)[2]) & (0xffff));
  bh = ((((a)[2]) >> 16) & (0xffff));
  {
    unsigned long m, m1, lt, ht;
    lt = t1;
    ht = t2;
    m = (bh) * (lt);
    lt = (bl) * (lt);
    m1 = (bl) * (ht);
    ht = (bh) * (ht);
    m = (m + m1) & (0xffffffffL);
    if (m < m1)
      ht += ((unsigned long)((1L) & (0xffff)) << 16);
    ht += (((m) >> 16) & (0xffff));
    m1 = ((unsigned long)((m) & (0xffff)) << 16);
    lt = (lt + m1) & (0xffffffffL);
    if (lt < m1)
      ht++;
    (t1) = lt;
    (t2) = ht;
  };
  if (t2 & (0x80000000L))
    c2++;
  t2 = (t2 + t2) & (0xffffffffL);
  if (t1 & (0x80000000L))
    t2++;
  t1 = (t1 + t1) & (0xffffffffL);
  c3 = (c3 + t1) & (0xffffffffL);
  if ((c3 < t1) && (((++t2) & (0xffffffffL)) == 0))
    c2++;
  c1 = (c1 + t2) & (0xffffffffL);
  if ((c1) < t2)
    c2++;
  ;
  r[5] = c3;
  c3 = 0;
  {
    unsigned long l, h, m;
    h = ((a)[3]);
    l = ((h) & (0xffff));
    h = (((h) >> 16) & (0xffff));
    m = (l) * (h);
    l *= l;
    h *= h;
    h += (m & (0xffff8000L)) >> (16 - 1);
    m = (m & (0xffff)) << (16 + 1);
    l = (l + m) & (0xffffffffL);
    if (l < m)
      h++;
    (t1) = l;
    (t2) = h;
  };
  c1 = (c1 + t1) & (0xffffffffL);
  if ((c1) < t1)
    t2++;
  c2 = (c2 + t2) & (0xffffffffL);
  if ((c2) < t2)
    c3++;
  ;
  r[6] = c1;
  r[7] = c2;
}
```