long int binomialCoef(long int n, long int r) {
  int i = 0;
  int pascal = n;
  if (r > n) {
    return 0;
  }
  if ((r * 2) > n) {
    r = n - r;
  }
  if (r == 0) {
    return 1;
  }
  for (i = 2; i <= r; i++) {
    pascal *= (n - i + 1);
    pascal /= i;
  }
  return pascal;
}


x86 assembly (clang, -O0):
<binomialCoef>:
push   %rbp
mov    %rsp,%rbp
mov    %rdi,-0x10(%rbp)
mov    %rsi,-0x18(%rbp)
movl   $0x0,-0x1c(%rbp)
mov    -0x10(%rbp),%rax
mov    %eax,-0x20(%rbp)
mov    -0x18(%rbp),%rax
cmp    -0x10(%rbp),%rax
jle    1135 <binomialCoef+0x35>
movq   $0x0,-0x8(%rbp)
jmp    11bc <binomialCoef+0xbc>
mov    -0x18(%rbp),%rax
shl    $0x1,%rax
cmp    -0x10(%rbp),%rax
jle    1153 <binomialCoef+0x53>
mov    -0x10(%rbp),%rax
sub    -0x18(%rbp),%rax
mov    %rax,-0x18(%rbp)
cmpq   $0x0,-0x18(%rbp)
jne    116b <binomialCoef+0x6b>
movq   $0x1,-0x8(%rbp)
jmp    11bc <binomialCoef+0xbc>
movl   $0x2,-0x1c(%rbp)
movslq -0x1c(%rbp),%rax
cmp    -0x18(%rbp),%rax
jg     11b4 <binomialCoef+0xb4>
mov    -0x10(%rbp),%rcx
movslq -0x1c(%rbp),%rax
sub    %rax,%rcx
add    $0x1,%rcx
movslq -0x20(%rbp),%rax
imul   %rcx,%rax
mov    %eax,-0x20(%rbp)
mov    -0x1c(%rbp),%ecx
mov    -0x20(%rbp),%eax
cltd
idiv   %ecx
mov    %eax,-0x20(%rbp)
mov    -0x1c(%rbp),%eax
add    $0x1,%eax
mov    %eax,-0x1c(%rbp)
jmp    1172 <binomialCoef+0x72>
movslq -0x20(%rbp),%rax
mov    %rax,-0x8(%rbp)
mov    -0x8(%rbp),%rax
pop    %rbp
ret

x86 assembly (clang, -O2):
<binomialCoef>:
mov    %rdi,%rcx
sub    %rsi,%rcx
jge    110b <binomialCoef+0xb>
xor    %eax,%eax
ret
mov    %rdi,%rax
lea    (%rsi,%rsi,1),%rdx
cmp    %rdi,%rdx
cmovle %rsi,%rcx
test   %rcx,%rcx
je     1145 <binomialCoef+0x45>
cmp    $0x2,%rcx
jl     1142 <binomialCoef+0x42>
lea    -0x1(%rax),%esi
mov    $0x1,%edi
nopl   0x0(%rax)
add    $0x1,%rdi
imul   %esi,%eax
cltd
idiv   %edi
add    $0xffffffff,%esi
cmp    %rdi,%rcx
jg     1130 <binomialCoef+0x30>
cltq
ret
mov    $0x1,%eax
ret

LLVM IR (-O0):
define i64 @binomialCoef(i64 %0, i64 %1) {
%3 = alloca i64
%4 = alloca i64
%5 = alloca i64
%6 = alloca i32
%7 = alloca i32
store i64 %0, i64* %4
store i64 %1, i64* %5
store i32 0, i32* %6
%8 = load i64, i64* %4
%9 = trunc i64 %8 to i32
store i32 %9, i32* %7
%10 = load i64, i64* %5
%11 = load i64, i64* %4
%12 = icmp sgt i64 %10, %11
br i1 %12, label %13, label %14
13: ; preds = %2
store i64 0, i64* %3
br label %52
14: ; preds = %2
%15 = load i64, i64* %5
%16 = mul nsw i64 %15, 2
%17 = load i64, i64* %4
%18 = icmp sgt i64 %16, %17
br i1 %18, label %19, label %23
19: ; preds = %14
%20 = load i64, i64* %4
%21 = load i64, i64* %5
%22 = sub nsw i64 %20, %21
store i64 %22, i64* %5
br label %23
23: ; preds = %19, %14
%24 = load i64, i64* %5
%25 = icmp eq i64 %24, 0
br i1 %25, label %26, label %27
26: ; preds = %23
store i64 1, i64* %3
br label %52
27: ; preds = %23
store i32 2, i32* %6
br label %28
28: ; preds = %46, %27
%29 = load i32, i32* %6
%30 = sext i32 %29 to i64
%31 = load i64, i64* %5
%32 = icmp sle i64 %30, %31
br i1 %32, label %33, label %49
33: ; preds = %28
%34 = load i64, i64* %4
%35 = load i32, i32* %6
%36 = sext i32 %35 to i64
%37 = sub nsw i64 %34, %36
%38 = add nsw i64 %37, 1
%39 = load i32, i32* %7
%40 = sext i32 %39 to i64
%41 = mul nsw i64 %40, %38
%42 = trunc i64 %41 to i32
store i32 %42, i32* %7
%43 = load i32, i32* %6
%44 = load i32, i32* %7
%45 = sdiv i32 %44, %43
store i32 %45, i32* %7
br label %46
46: ; preds = %33
%47 = load i32, i32* %6
%48 = add nsw i32 %47, 1
store i32 %48, i32* %6
br label %28
49: ; preds = %28
%50 = load i32, i32* %7
%51 = sext i32 %50 to i64
store i64 %51, i64* %3
br label %52
52: ; preds = %49, %26, %13
%53 = load i64, i64* %3
ret i64 %53
}

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