int sqrti(int n) {
  int p = 0;
  int q = 1;
  int r = n;
  int h = 0;
  while (q <= n)
    q = 4 * q;
  while (q != 1) {
    q = q / 4;
    h = p + q;
    p = p / 2;
    if (r >= h) {
      p = p + q;
      r = r - h;
    }
  }
  return p;
}


x86 assembly (clang, -O0):
<sqrti>:
push   %rbp
mov    %rsp,%rbp
mov    %edi,-0x4(%rbp)
movl   $0x0,-0x8(%rbp)
movl   $0x1,-0xc(%rbp)
mov    -0x4(%rbp),%eax
mov    %eax,-0x10(%rbp)
movl   $0x0,-0x14(%rbp)
mov    -0xc(%rbp),%eax
cmp    -0x4(%rbp),%eax
jg     113c <sqrti+0x3c>
mov    -0xc(%rbp),%eax
shl    $0x2,%eax
mov    %eax,-0xc(%rbp)
jmp    1122 <sqrti+0x22>
jmp    1141 <sqrti+0x41>
cmpl   $0x1,-0xc(%rbp)
je     1193 <sqrti+0x93>
mov    -0xc(%rbp),%eax
mov    $0x4,%ecx
cltd
idiv   %ecx
mov    %eax,-0xc(%rbp)
mov    -0x8(%rbp),%eax
add    -0xc(%rbp),%eax
mov    %eax,-0x14(%rbp)
mov    -0x8(%rbp),%eax
mov    $0x2,%ecx
cltd
idiv   %ecx
mov    %eax,-0x8(%rbp)
mov    -0x10(%rbp),%eax
cmp    -0x14(%rbp),%eax
jl     118e <sqrti+0x8e>
mov    -0x8(%rbp),%eax
add    -0xc(%rbp),%eax
mov    %eax,-0x8(%rbp)
mov    -0x10(%rbp),%eax
sub    -0x14(%rbp),%eax
mov    %eax,-0x10(%rbp)
jmp    1141 <sqrti+0x41>
mov    -0x8(%rbp),%eax
pop    %rbp
ret

x86 assembly (clang, -O2):
<sqrti>:
mov    $0x1,%eax
cs nopw 0x0(%rax,%rax,1)
nop
mov    %eax,%ecx
lea    0x0(,%rcx,4),%eax
cmp    %edi,%ecx
jle    1110 <sqrti+0x10>
xor    %eax,%eax
cmp    $0x1,%ecx
je     1164 <sqrti+0x64>
xor    %r8d,%r8d
xor    %eax,%eax
nopl   0x0(%rax)
lea    0x3(%rcx),%esi
test   %ecx,%ecx
cmovns %ecx,%esi
sar    $0x2,%esi
lea    (%rsi,%rax,1),%r9d
mov    %eax,%edx
shr    $0x1f,%edx
add    %eax,%edx
sar    %edx
cmp    %r9d,%edi
mov    %esi,%eax
cmovl  %r8d,%eax
cmovl  %r8d,%r9d
add    %edx,%eax
sub    %r9d,%edi
and    $0xfffffffc,%ecx
cmp    $0x4,%ecx
mov    %esi,%ecx
jne    1130 <sqrti+0x30>
ret

LLVM IR (-O0):
define i32 @sqrti(i32 %0) {
%2 = alloca i32
%3 = alloca i32
%4 = alloca i32
%5 = alloca i32
%6 = alloca i32
store i32 %0, i32* %2
store i32 0, i32* %3
store i32 1, i32* %4
%7 = load i32, i32* %2
store i32 %7, i32* %5
store i32 0, i32* %6
br label %8
8: ; preds = %12, %1
%9 = load i32, i32* %4
%10 = load i32, i32* %2
%11 = icmp sle i32 %9, %10
br i1 %11, label %12, label %15
12: ; preds = %8
%13 = load i32, i32* %4
%14 = mul nsw i32 4, %13
store i32 %14, i32* %4
br label %8
15: ; preds = %8
br label %16
16: ; preds = %37, %15
%17 = load i32, i32* %4
%18 = icmp ne i32 %17, 1
br i1 %18, label %19, label %38
19: ; preds = %16
%20 = load i32, i32* %4
%21 = sdiv i32 %20, 4
store i32 %21, i32* %4
%22 = load i32, i32* %3
%23 = load i32, i32* %4
%24 = add nsw i32 %22, %23
store i32 %24, i32* %6
%25 = load i32, i32* %3
%26 = sdiv i32 %25, 2
store i32 %26, i32* %3
%27 = load i32, i32* %5
%28 = load i32, i32* %6
%29 = icmp sge i32 %27, %28
br i1 %29, label %30, label %37
30: ; preds = %19
%31 = load i32, i32* %3
%32 = load i32, i32* %4
%33 = add nsw i32 %31, %32
store i32 %33, i32* %3
%34 = load i32, i32* %5
%35 = load i32, i32* %6
%36 = sub nsw i32 %34, %35
store i32 %36, i32* %5
br label %37
37: ; preds = %30, %19
br label %16
38: ; preds = %16
%39 = load i32, i32* %3
ret i32 %39
}

LLVM IR (-O2):
define i32 @sqrti(i32 %0) {
br label %2
2: ; preds = %2, %1
%3 = phi i32 [ 1, %1 ], [ %5, %2 ]
%4 = icmp sgt i32 %3, %0
%5 = shl nsw i32 %3, 2
br i1 %4, label %6, label %2
6: ; preds = %2
%7 = icmp eq i32 %3, 1
br i1 %7, label %22, label %8
8: ; preds = %6, %8
%9 = phi i32 [ %19, %8 ], [ %0, %6 ]
%10 = phi i32 [ %12, %8 ], [ %3, %6 ]
%11 = phi i32 [ %17, %8 ], [ 0, %6 ]
%12 = sdiv i32 %10, 4
%13 = add nsw i32 %12, %11
%14 = sdiv i32 %11, 2
%15 = icmp slt i32 %9, %13
%16 = select i1 %15, i32 0, i32 %12
%17 = add nsw i32 %16, %14
%18 = select i1 %15, i32 0, i32 %13
%19 = sub nsw i32 %9, %18
%20 = and i32 %10, -4
%21 = icmp eq i32 %20, 4
br i1 %21, label %22, label %8
22: ; preds = %8, %6
%23 = phi i32 [ 0, %6 ], [ %17, %8 ]
ret i32 %23
}