long unsigned sqrtllu(long long unsigned t) {
  long long unsigned s;
  long long unsigned b;
  for (b = 0, s = t; b++, (s >>= 1) != 0;)
    ;
  s = 1LL << (b >> 1);
  if (b & 1)
    s += s >> 1;
  do {
    b = t / s;
    s = (s + b) >> 1;
  } while (b < s);
  return s;
}


x86 assembly (clang, -O0):
<sqrtllu>:
push   %rbp
mov    %rsp,%rbp
mov    %rdi,-0x8(%rbp)
movq   $0x0,-0x18(%rbp)
mov    -0x8(%rbp),%rax
mov    %rax,-0x10(%rbp)
mov    -0x18(%rbp),%rax
add    $0x1,%rax
mov    %rax,-0x18(%rbp)
mov    -0x10(%rbp),%rax
shr    $0x1,%rax
mov    %rax,-0x10(%rbp)
cmp    $0x0,%rax
je     113f <sqrtllu+0x3f>
jmp    1118 <sqrtllu+0x18>
mov    -0x18(%rbp),%rcx
shr    $0x1,%rcx
mov    $0x1,%eax
shl    %cl,%rax
mov    %rax,-0x10(%rbp)
mov    -0x18(%rbp),%rax
and    $0x1,%rax
cmp    $0x0,%rax
je     1175 <sqrtllu+0x75>
mov    -0x10(%rbp),%rax
shr    $0x1,%rax
add    -0x10(%rbp),%rax
mov    %rax,-0x10(%rbp)
jmp    117a <sqrtllu+0x7a>
mov    -0x8(%rbp),%rax
xor    %ecx,%ecx
mov    %ecx,%edx
divq   -0x10(%rbp)
mov    %rax,-0x18(%rbp)
mov    -0x10(%rbp),%rax
add    -0x18(%rbp),%rax
shr    $0x1,%rax
mov    %rax,-0x10(%rbp)
mov    -0x18(%rbp),%rax
cmp    -0x10(%rbp),%rax
jb     117a <sqrtllu+0x7a>
mov    -0x10(%rbp),%rax
pop    %rbp
ret

x86 assembly (clang, -O2):
<sqrtllu>:
xor    %eax,%eax
mov    %rdi,%rcx
mov    %rdi,%rdx
nopl   0x0(%rax,%rax,1)
add    $0x1,%rax
shr    %rdx
cmp    $0x1,%rcx
mov    %rdx,%rcx
ja     1110 <sqrtllu+0x10>
mov    %rax,%rcx
shr    %rcx
mov    $0x1,%edx
shl    %cl,%rdx
and    $0x1,%eax
neg    %rax
mov    %rdx,%rcx
shr    %rcx
and    %rax,%rcx
add    %rdx,%rcx
jmp    1163 <sqrtllu+0x63>
cs nopw 0x0(%rax,%rax,1)
nopl   0x0(%rax)
mov    %rdi,%rax
xor    %edx,%edx
div    %rcx
add    %rax,%rcx
shr    %rcx
cmp    %rcx,%rax
jae    1177 <sqrtllu+0x77>
mov    %rdi,%rax
or     %rcx,%rax
shr    $0x20,%rax
jne    1150 <sqrtllu+0x50>
mov    %edi,%eax
xor    %edx,%edx
div    %ecx
jmp    1158 <sqrtllu+0x58>
mov    %rcx,%rax
ret

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

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