unsigned int divisor_sum(unsigned int n) {
  unsigned int total = 1, power = 2;
  unsigned int p;
  for (; (n & 1) == 0; power <<= 1, n >>= 1) {
    total += power;
  }
  for (p = 3; p * p <= n; p += 2) {
    unsigned int sum = 1;
    for (power = p; n % p == 0; power *= p, n /= p) {
      sum += power;
    }
    total *= sum;
  }
  if (n > 1) {
    total *= n + 1;
  }
  return total;
}


x86 assembly (clang, -O0):
<divisor_sum>:
push   %rbp
mov    %rsp,%rbp
mov    %edi,-0x4(%rbp)
movl   $0x1,-0x8(%rbp)
movl   $0x2,-0xc(%rbp)
mov    -0x4(%rbp),%eax
and    $0x1,%eax
cmp    $0x0,%eax
jne    1144 <divisor_sum+0x44>
mov    -0xc(%rbp),%eax
add    -0x8(%rbp),%eax
mov    %eax,-0x8(%rbp)
mov    -0xc(%rbp),%eax
shl    $0x1,%eax
mov    %eax,-0xc(%rbp)
mov    -0x4(%rbp),%eax
shr    $0x1,%eax
mov    %eax,-0x4(%rbp)
jmp    1115 <divisor_sum+0x15>
movl   $0x3,-0x10(%rbp)
mov    -0x10(%rbp),%eax
imul   -0x10(%rbp),%eax
cmp    -0x4(%rbp),%eax
ja     11b6 <divisor_sum+0xb6>
movl   $0x1,-0x14(%rbp)
mov    -0x10(%rbp),%eax
mov    %eax,-0xc(%rbp)
mov    -0x4(%rbp),%eax
xor    %edx,%edx
divl   -0x10(%rbp)
cmp    $0x0,%edx
jne    119e <divisor_sum+0x9e>
mov    -0xc(%rbp),%eax
add    -0x14(%rbp),%eax
mov    %eax,-0x14(%rbp)
mov    -0x10(%rbp),%eax
imul   -0xc(%rbp),%eax
mov    %eax,-0xc(%rbp)
mov    -0x10(%rbp),%ecx
mov    -0x4(%rbp),%eax
xor    %edx,%edx
div    %ecx
mov    %eax,-0x4(%rbp)
jmp    1168 <divisor_sum+0x68>
mov    -0x14(%rbp),%eax
imul   -0x8(%rbp),%eax
mov    %eax,-0x8(%rbp)
mov    -0x10(%rbp),%eax
add    $0x2,%eax
mov    %eax,-0x10(%rbp)
jmp    114b <divisor_sum+0x4b>
cmpl   $0x1,-0x4(%rbp)
jbe    11cd <divisor_sum+0xcd>
mov    -0x4(%rbp),%eax
add    $0x1,%eax
imul   -0x8(%rbp),%eax
mov    %eax,-0x8(%rbp)
mov    -0x8(%rbp),%eax
pop    %rbp
ret

x86 assembly (clang, -O2):
<divisor_sum>:
mov    $0x1,%r8d
test   $0x1,%dil
jne    1136 <divisor_sum+0x36>
mov    $0x1,%r8d
mov    $0x2,%eax
mov    %edi,%ecx
nopl   0x0(%rax)
add    %eax,%r8d
add    %eax,%eax
shr    %ecx
test   $0x2,%dil
mov    %ecx,%edi
je     1120 <divisor_sum+0x20>
cmp    $0x9,%ecx
jb     113d <divisor_sum+0x3d>
jmp    1150 <divisor_sum+0x50>
mov    %edi,%ecx
cmp    $0x9,%ecx
jae    1150 <divisor_sum+0x50>
lea    0x1(%rcx),%edx
cmp    $0x2,%ecx
mov    $0x1,%eax
cmovae %edx,%eax
imul   %r8d,%eax
ret
mov    $0x3,%r9d
jmp    1173 <divisor_sum+0x73>
nopl   0x0(%rax,%rax,1)
imul   %esi,%r8d
add    $0x2,%r9d
mov    %r9d,%eax
imul   %r9d,%eax
cmp    %ecx,%eax
ja     113d <divisor_sum+0x3d>
mov    %ecx,%eax
xor    %edx,%edx
div    %r9d
mov    $0x1,%esi
test   %edx,%edx
jne    1160 <divisor_sum+0x60>
mov    $0x1,%esi
mov    %r9d,%edi
nopl   0x0(%rax,%rax,1)
add    %edi,%esi
imul   %r9d,%edi
mov    %ecx,%eax
xor    %edx,%edx
div    %r9d
mov    %eax,%ecx
xor    %edx,%edx
div    %r9d
test   %edx,%edx
je     1190 <divisor_sum+0x90>
jmp    1160 <divisor_sum+0x60>

LLVM IR (-O0):
define i32 @divisor_sum(i32 %0) {
%2 = alloca i32
%3 = alloca i32
%4 = alloca i32
%5 = alloca i32
%6 = alloca i32
store i32 %0, i32* %2
store i32 1, i32* %3
store i32 2, i32* %4
br label %7
7: ; preds = %15, %1
%8 = load i32, i32* %2
%9 = and i32 %8, 1
%10 = icmp eq i32 %9, 0
br i1 %10, label %11, label %20
11: ; preds = %7
%12 = load i32, i32* %4
%13 = load i32, i32* %3
%14 = add i32 %13, %12
store i32 %14, i32* %3
br label %15
15: ; preds = %11
%16 = load i32, i32* %4
%17 = shl i32 %16, 1
store i32 %17, i32* %4
%18 = load i32, i32* %2
%19 = lshr i32 %18, 1
store i32 %19, i32* %2
br label %7
20: ; preds = %7
store i32 3, i32* %5
br label %21
21: ; preds = %49, %20
%22 = load i32, i32* %5
%23 = load i32, i32* %5
%24 = mul i32 %22, %23
%25 = load i32, i32* %2
%26 = icmp ule i32 %24, %25
br i1 %26, label %27, label %52
27: ; preds = %21
store i32 1, i32* %6
%28 = load i32, i32* %5
store i32 %28, i32* %4
br label %29
29: ; preds = %38, %27
%30 = load i32, i32* %2
%31 = load i32, i32* %5
%32 = urem i32 %30, %31
%33 = icmp eq i32 %32, 0
br i1 %33, label %34, label %45
34: ; preds = %29
%35 = load i32, i32* %4
%36 = load i32, i32* %6
%37 = add i32 %36, %35
store i32 %37, i32* %6
br label %38
38: ; preds = %34
%39 = load i32, i32* %5
%40 = load i32, i32* %4
%41 = mul i32 %40, %39
store i32 %41, i32* %4
%42 = load i32, i32* %5
%43 = load i32, i32* %2
%44 = udiv i32 %43, %42
store i32 %44, i32* %2
br label %29
45: ; preds = %29
%46 = load i32, i32* %6
%47 = load i32, i32* %3
%48 = mul i32 %47, %46
store i32 %48, i32* %3
br label %49
49: ; preds = %45
%50 = load i32, i32* %5
%51 = add i32 %50, 2
store i32 %51, i32* %5
br label %21
52: ; preds = %21
%53 = load i32, i32* %2
%54 = icmp ugt i32 %53, 1
br i1 %54, label %55, label %60
55: ; preds = %52
%56 = load i32, i32* %2
%57 = add i32 %56, 1
%58 = load i32, i32* %3
%59 = mul i32 %58, %57
store i32 %59, i32* %3
br label %60
60: ; preds = %55, %52
%61 = load i32, i32* %3
ret i32 %61
}

LLVM IR (-O2):
define i32 @divisor_sum(i32 %0) {
%2 = and i32 %0, 1
%3 = icmp eq i32 %2, 0
br i1 %3, label %8, label %4
4: ; preds = %8, %1
%5 = phi i32 [ %0, %1 ], [ %14, %8 ]
%6 = phi i32 [ 1, %1 ], [ %12, %8 ]
%7 = icmp ult i32 %5, 9
br i1 %7, label %39, label %17
8: ; preds = %1, %8
%9 = phi i32 [ %13, %8 ], [ 2, %1 ]
%10 = phi i32 [ %12, %8 ], [ 1, %1 ]
%11 = phi i32 [ %14, %8 ], [ %0, %1 ]
%12 = add i32 %9, %10
%13 = shl i32 %9, 1
%14 = lshr i32 %11, 1
%15 = and i32 %11, 2
%16 = icmp eq i32 %15, 0
br i1 %16, label %8, label %4
17: ; preds = %4, %32
%18 = phi i32 [ %36, %32 ], [ 3, %4 ]
%19 = phi i32 [ %35, %32 ], [ %6, %4 ]
%20 = phi i32 [ %33, %32 ], [ %5, %4 ]
%21 = urem i32 %20, %18
%22 = icmp eq i32 %21, 0
br i1 %22, label %23, label %32
23: ; preds = %17, %23
%24 = phi i32 [ %27, %23 ], [ 1, %17 ]
%25 = phi i32 [ %28, %23 ], [ %18, %17 ]
%26 = phi i32 [ %29, %23 ], [ %20, %17 ]
%27 = add i32 %24, %25
%28 = mul i32 %25, %18
%29 = udiv i32 %26, %18
%30 = urem i32 %29, %18
%31 = icmp eq i32 %30, 0
br i1 %31, label %23, label %32
32: ; preds = %23, %17
%33 = phi i32 [ %20, %17 ], [ %29, %23 ]
%34 = phi i32 [ 1, %17 ], [ %27, %23 ]
%35 = mul i32 %34, %19
%36 = add i32 %18, 2
%37 = mul i32 %36, %36
%38 = icmp ugt i32 %37, %33
br i1 %38, label %39, label %17
39: ; preds = %32, %4
%40 = phi i32 [ %5, %4 ], [ %33, %32 ]
%41 = phi i32 [ %6, %4 ], [ %35, %32 ]
%42 = icmp ugt i32 %40, 1
%43 = add nuw i32 %40, 1
%44 = select i1 %42, i32 %43, i32 1
%45 = mul i32 %44, %41
ret i32 %45
}