unsigned choose(unsigned n, unsigned k) {
  unsigned ans = 1;
  k = k > n - k ? n - k : k;
  unsigned j = 1;
  for (; j <= k; j++, n--) {
    if (n % j == 0) {
      ans *= n / j;
    } else if (ans % j == 0) {
      ans = ans / j * n;
    } else {
      ans = (ans * n) / j;
    }
  }
  return ans;
}


x86 assembly (clang, -O0):
<choose>:
push   %rbp
mov    %rsp,%rbp
mov    %edi,-0x4(%rbp)
mov    %esi,-0x8(%rbp)
movl   $0x1,-0xc(%rbp)
mov    -0x8(%rbp),%eax
mov    -0x4(%rbp),%ecx
sub    -0x8(%rbp),%ecx
cmp    %ecx,%eax
jbe    1130 <choose+0x30>
mov    -0x4(%rbp),%eax
sub    -0x8(%rbp),%eax
mov    %eax,-0x14(%rbp)
jmp    1136 <choose+0x36>
mov    -0x8(%rbp),%eax
mov    %eax,-0x14(%rbp)
mov    -0x14(%rbp),%eax
mov    %eax,-0x8(%rbp)
movl   $0x1,-0x10(%rbp)
mov    -0x10(%rbp),%eax
cmp    -0x8(%rbp),%eax
ja     11c9 <choose+0xc9>
mov    -0x4(%rbp),%eax
xor    %edx,%edx
divl   -0x10(%rbp)
cmp    $0x0,%edx
jne    1174 <choose+0x74>
mov    -0x4(%rbp),%eax
xor    %edx,%edx
divl   -0x10(%rbp)
imul   -0xc(%rbp),%eax
mov    %eax,-0xc(%rbp)
jmp    11ad <choose+0xad>
mov    -0xc(%rbp),%eax
xor    %edx,%edx
divl   -0x10(%rbp)
cmp    $0x0,%edx
jne    1199 <choose+0x99>
mov    -0xc(%rbp),%eax
xor    %edx,%edx
divl   -0x10(%rbp)
imul   -0x4(%rbp),%eax
mov    %eax,-0xc(%rbp)
jmp    11a8 <choose+0xa8>
mov    -0xc(%rbp),%eax
imul   -0x4(%rbp),%eax
xor    %edx,%edx
divl   -0x10(%rbp)
mov    %eax,-0xc(%rbp)
jmp    11ad <choose+0xad>
jmp    11b2 <choose+0xb2>
mov    -0x10(%rbp),%eax
add    $0x1,%eax
mov    %eax,-0x10(%rbp)
mov    -0x4(%rbp),%eax
add    $0xffffffff,%eax
mov    %eax,-0x4(%rbp)
jmp    1143 <choose+0x43>
mov    -0xc(%rbp),%eax
pop    %rbp
ret

x86 assembly (clang, -O2):
<choose>:
mov    %edi,%r8d
sub    %esi,%r8d
cmp    %esi,%r8d
cmovae %esi,%r8d
test   %r8d,%r8d
je     117a <choose+0x7a>
neg    %r8d
mov    $0x1,%esi
mov    $0x1,%ecx
jmp    1149 <choose+0x49>
cs nopw 0x0(%rax,%rax,1)
nopl   0x0(%rax,%rax,1)
imul   %eax,%ecx
add    $0xffffffff,%edi
lea    (%r8,%rsi,1),%eax
add    $0x1,%eax
mov    %esi,%edx
add    $0x1,%edx
mov    %edx,%esi
cmp    $0x1,%eax
je     1177 <choose+0x77>
mov    %edi,%eax
xor    %edx,%edx
div    %esi
test   %edx,%edx
je     1130 <choose+0x30>
mov    %ecx,%eax
xor    %edx,%edx
div    %esi
test   %edx,%edx
je     1170 <choose+0x70>
imul   %edi,%ecx
mov    %ecx,%eax
xor    %edx,%edx
div    %esi
jmp    1173 <choose+0x73>
nopl   0x0(%rax,%rax,1)
imul   %edi,%eax
mov    %eax,%ecx
jmp    1133 <choose+0x33>
mov    %ecx,%eax
ret
mov    $0x1,%eax
ret

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

LLVM IR (-O2):
define i32 @choose(i32 %0, i32 %1) {
%3 = sub i32 %0, %1
%4 = icmp ult i32 %3, %1
%5 = select i1 %4, i32 %3, i32 %1
%6 = icmp eq i32 %5, 0
br i1 %6, label %30, label %7
7: ; preds = %2, %25
%8 = phi i32 [ %27, %25 ], [ 1, %2 ]
%9 = phi i32 [ %26, %25 ], [ 1, %2 ]
%10 = phi i32 [ %28, %25 ], [ %0, %2 ]
%11 = urem i32 %10, %8
%12 = udiv i32 %10, %8
%13 = icmp eq i32 %11, 0
br i1 %13, label %14, label %16
14: ; preds = %7
%15 = mul i32 %12, %9
br label %25
16: ; preds = %7
%17 = urem i32 %9, %8
%18 = udiv i32 %9, %8
%19 = icmp eq i32 %17, 0
br i1 %19, label %20, label %22
20: ; preds = %16
%21 = mul i32 %18, %10
br label %25
22: ; preds = %16
%23 = mul i32 %9, %10
%24 = udiv i32 %23, %8
br label %25
25: ; preds = %14, %22, %20
%26 = phi i32 [ %15, %14 ], [ %21, %20 ], [ %24, %22 ]
%27 = add nuw i32 %8, 1
%28 = add i32 %10, -1
%29 = icmp eq i32 %8, %5
br i1 %29, label %30, label %7
30: ; preds = %25, %2
%31 = phi i32 [ 1, %2 ], [ %26, %25 ]
ret i32 %31
}