unsigned long CSPL_Stats_combination(unsigned long N, unsigned long k) {
  unsigned long ans = 1;
  unsigned long i;
  if ((k > N) || (N == 0) || (k == 0))
    return (0);
  for (i = 0; i < ({
                __typeof__(k) _a = (k);
                __typeof__(N - k) _b = (N - k);
                _a < _b ? _a : _b;
              });
       i++)
    ans = (ans * (N - i)) / (i + 1);
  return (ans);
}


x86 assembly (clang, -O0):
<CSPL_Stats_combination>:
push   %rbp
mov    %rsp,%rbp
mov    %rdi,-0x10(%rbp)
mov    %rsi,-0x18(%rbp)
movq   $0x1,-0x20(%rbp)
mov    -0x18(%rbp),%rax
cmp    -0x10(%rbp),%rax
ja     1138 <CSPL_Stats_combination+0x38>
cmpq   $0x0,-0x10(%rbp)
je     1138 <CSPL_Stats_combination+0x38>
cmpq   $0x0,-0x18(%rbp)
jne    1145 <CSPL_Stats_combination+0x45>
movq   $0x0,-0x8(%rbp)
jmp    11dc <CSPL_Stats_combination+0xdc>
movq   $0x0,-0x28(%rbp)
mov    -0x28(%rbp),%rax
mov    %rax,-0x48(%rbp)
mov    -0x18(%rbp),%rax
mov    %rax,-0x30(%rbp)
mov    -0x10(%rbp),%rax
sub    -0x18(%rbp),%rax
mov    %rax,-0x38(%rbp)
mov    -0x30(%rbp),%rax
cmp    -0x38(%rbp),%rax
jae    1184 <CSPL_Stats_combination+0x84>
mov    -0x30(%rbp),%rax
mov    %rax,-0x50(%rbp)
jmp    118c <CSPL_Stats_combination+0x8c>
mov    -0x38(%rbp),%rax
mov    %rax,-0x50(%rbp)
mov    -0x48(%rbp),%rax
mov    -0x50(%rbp),%rcx
mov    %rcx,-0x40(%rbp)
cmp    -0x40(%rbp),%rax
jae    11d4 <CSPL_Stats_combination+0xd4>
mov    -0x20(%rbp),%rax
mov    -0x10(%rbp),%rcx
sub    -0x28(%rbp),%rcx
imul   %rcx,%rax
mov    -0x28(%rbp),%rcx
add    $0x1,%rcx
xor    %edx,%edx
div    %rcx
mov    %rax,-0x20(%rbp)
mov    -0x28(%rbp),%rax
add    $0x1,%rax
mov    %rax,-0x28(%rbp)
jmp    114d <CSPL_Stats_combination+0x4d>
mov    -0x20(%rbp),%rax
mov    %rax,-0x8(%rbp)
mov    -0x8(%rbp),%rax
pop    %rbp
ret

x86 assembly (clang, -O2):
<CSPL_Stats_combination>:
lea    -0x1(%rsi),%rax
cmp    %rdi,%rax
jae    114c <CSPL_Stats_combination+0x4c>
mov    %rdi,%r8
sub    %rsi,%r8
cmp    %rsi,%r8
cmova  %rsi,%r8
test   %r8,%r8
je     114f <CSPL_Stats_combination+0x4f>
cmp    $0x1,%r8
jne    1155 <CSPL_Stats_combination+0x55>
mov    $0x1,%eax
xor    %r10d,%r10d
test   $0x1,%r8b
je     114b <CSPL_Stats_combination+0x4b>
sub    %r10,%rdi
imul   %rdi,%rax
add    $0x1,%r10
mov    %rax,%rcx
or     %r10,%rcx
shr    $0x20,%rcx
je     11c5 <CSPL_Stats_combination+0xc5>
xor    %edx,%edx
div    %r10
ret
xor    %eax,%eax
ret
mov    $0x1,%eax
ret
mov    %r8,%r9
and    $0xfffffffffffffffe,%r9
mov    $0x1,%eax
xor    %r10d,%r10d
mov    %rdi,%rsi
jmp    117e <CSPL_Stats_combination+0x7e>
nopl   0x0(%rax)
xor    %edx,%edx
div    %r10
add    $0xfffffffffffffffe,%rsi
cmp    %r10,%r9
je     1129 <CSPL_Stats_combination+0x29>
imul   %rsi,%rax
lea    0x1(%r10),%rcx
mov    %rax,%rdx
or     %rcx,%rdx
shr    $0x20,%rdx
je     11a0 <CSPL_Stats_combination+0xa0>
xor    %edx,%edx
div    %rcx
mov    %rax,%rcx
jmp    11a6 <CSPL_Stats_combination+0xa6>
nopl   0x0(%rax)
xor    %edx,%edx
div    %ecx
mov    %eax,%ecx
lea    -0x1(%rsi),%rax
imul   %rcx,%rax
add    $0x2,%r10
mov    %rax,%rcx
or     %r10,%rcx
shr    $0x20,%rcx
jne    1170 <CSPL_Stats_combination+0x70>
xor    %edx,%edx
div    %r10d
jmp    1175 <CSPL_Stats_combination+0x75>
xor    %edx,%edx
div    %r10d
ret

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

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