double modpow(long base, long exponent, long modulus) {
  long result = 1;
  while (exponent > 0) {
    if (exponent & 1) {
      result = (result * base) % modulus;
    }
    exponent >>= 1;
    base = (base * base) % modulus;
  }
  return (double)result;
}


x86 assembly (clang, -O0):
<modpow>:
push   %rbp
mov    %rsp,%rbp
mov    %rdi,-0x8(%rbp)
mov    %rsi,-0x10(%rbp)
mov    %rdx,-0x18(%rbp)
movq   $0x1,-0x20(%rbp)
cmpq   $0x0,-0x10(%rbp)
jle    116c <modpow+0x6c>
mov    -0x10(%rbp),%rax
and    $0x1,%rax
cmp    $0x0,%rax
je     1148 <modpow+0x48>
mov    -0x20(%rbp),%rax
imul   -0x8(%rbp),%rax
cqto
idivq  -0x18(%rbp)
mov    %rdx,-0x20(%rbp)
mov    -0x10(%rbp),%rax
sar    $0x1,%rax
mov    %rax,-0x10(%rbp)
mov    -0x8(%rbp),%rax
imul   -0x8(%rbp),%rax
cqto
idivq  -0x18(%rbp)
mov    %rdx,-0x8(%rbp)
jmp    1118 <modpow+0x18>
cvtsi2sdq -0x20(%rbp),%xmm0
pop    %rbp
ret

x86 assembly (clang, -O2):
<modpow>:
mov    $0x1,%r8d
test   %rsi,%rsi
jle    117a <modpow+0x7a>
mov    %rdx,%rcx
jmp    1145 <modpow+0x45>
mov    %r8,%rax
cqto
idiv   %rcx
mov    %rdx,%r8
imul   %rdi,%rdi
mov    %rdi,%rax
or     %rcx,%rax
shr    $0x20,%rax
je     1170 <modpow+0x70>
mov    %rdi,%rax
cqto
idiv   %rcx
mov    %rdx,%rdi
mov    %rsi,%rax
shr    %rax
cmp    $0x1,%rsi
mov    %rax,%rsi
jbe    117a <modpow+0x7a>
test   $0x1,%sil
je     111b <modpow+0x1b>
imul   %rdi,%r8
mov    %r8,%rax
or     %rcx,%rax
shr    $0x20,%rax
jne    1110 <modpow+0x10>
mov    %r8d,%eax
xor    %edx,%edx
div    %ecx
mov    %edx,%r8d
jmp    111b <modpow+0x1b>
nopw   0x0(%rax,%rax,1)
mov    %edi,%eax
xor    %edx,%edx
div    %ecx
mov    %edx,%edi
jmp    1136 <modpow+0x36>
cvtsi2sd %r8,%xmm0
ret

LLVM IR (-O0):
define double @modpow(i64 %0, i64 %1, i64 %2) {
%4 = alloca i64
%5 = alloca i64
%6 = alloca i64
%7 = alloca i64
store i64 %0, i64* %4
store i64 %1, i64* %5
store i64 %2, i64* %6
store i64 1, i64* %7
br label %8
8: ; preds = %21, %3
%9 = load i64, i64* %5
%10 = icmp sgt i64 %9, 0
br i1 %10, label %11, label %29
11: ; preds = %8
%12 = load i64, i64* %5
%13 = and i64 %12, 1
%14 = icmp ne i64 %13, 0
br i1 %14, label %15, label %21
15: ; preds = %11
%16 = load i64, i64* %7
%17 = load i64, i64* %4
%18 = mul nsw i64 %16, %17
%19 = load i64, i64* %6
%20 = srem i64 %18, %19
store i64 %20, i64* %7
br label %21
21: ; preds = %15, %11
%22 = load i64, i64* %5
%23 = ashr i64 %22, 1
store i64 %23, i64* %5
%24 = load i64, i64* %4
%25 = load i64, i64* %4
%26 = mul nsw i64 %24, %25
%27 = load i64, i64* %6
%28 = srem i64 %26, %27
store i64 %28, i64* %4
br label %8
29: ; preds = %8
%30 = load i64, i64* %7
%31 = sitofp i64 %30 to double
ret double %31
}

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