double binomial(int x, int y, int x2) {
  if (y > x)
    return 0;
  double answer = 1.0;
  for (int n = x - y + 1; n <= x; n++) {
    double d = n - x + x2;
    answer *= n / d;
  }
  return answer;
}


x86 assembly (clang, -O0):
<binomial>:
push   %rbp
mov    %rsp,%rbp
mov    %edi,-0xc(%rbp)
mov    %esi,-0x10(%rbp)
mov    %edx,-0x14(%rbp)
mov    -0x10(%rbp),%eax
cmp    -0xc(%rbp),%eax
jle    1126 <binomial+0x26>
xorps  %xmm0,%xmm0
movsd  %xmm0,-0x8(%rbp)
jmp    1189 <binomial+0x89>
movsd  0xed2(%rip),%xmm0
movsd  %xmm0,-0x20(%rbp)
mov    -0xc(%rbp),%eax
sub    -0x10(%rbp),%eax
add    $0x1,%eax
mov    %eax,-0x24(%rbp)
mov    -0x24(%rbp),%eax
cmp    -0xc(%rbp),%eax
jg     117f <binomial+0x7f>
mov    -0x24(%rbp),%eax
sub    -0xc(%rbp),%eax
add    -0x14(%rbp),%eax
cvtsi2sd %eax,%xmm0
movsd  %xmm0,-0x30(%rbp)
cvtsi2sdl -0x24(%rbp),%xmm0
divsd  -0x30(%rbp),%xmm0
mulsd  -0x20(%rbp),%xmm0
movsd  %xmm0,-0x20(%rbp)
mov    -0x24(%rbp),%eax
add    $0x1,%eax
mov    %eax,-0x24(%rbp)
jmp    113f <binomial+0x3f>
movsd  -0x20(%rbp),%xmm0
movsd  %xmm0,-0x8(%rbp)
movsd  -0x8(%rbp),%xmm0
pop    %rbp
ret

x86 assembly (clang, -O2):
<binomial>:
xorpd  %xmm0,%xmm0
mov    %edi,%eax
sub    %esi,%eax
jl     1140 <binomial+0x40>
test   %esi,%esi
jle    1141 <binomial+0x41>
sub    %esi,%edx
add    $0x1,%edx
movsd  0xee5(%rip),%xmm0
nopl   0x0(%rax,%rax,1)
add    $0x1,%eax
xorps  %xmm1,%xmm1
cvtsi2sd %edx,%xmm1
xorps  %xmm2,%xmm2
cvtsi2sd %eax,%xmm2
divsd  %xmm1,%xmm2
mulsd  %xmm2,%xmm0
add    $0x1,%edx
cmp    %edi,%eax
jl     1120 <binomial+0x20>
ret
movsd  0xeb7(%rip),%xmm0
ret

LLVM IR (-O0):
define double @binomial(i32 %0, i32 %1, i32 %2) {
%4 = alloca double
%5 = alloca i32
%6 = alloca i32
%7 = alloca i32
%8 = alloca double
%9 = alloca i32
%10 = alloca double
store i32 %0, i32* %5
store i32 %1, i32* %6
store i32 %2, i32* %7
%11 = load i32, i32* %6
%12 = load i32, i32* %5
%13 = icmp sgt i32 %11, %12
br i1 %13, label %14, label %15
14: ; preds = %3
store double 0.000000e+00, double* %4
br label %42
15: ; preds = %3
store double 1.000000e+00, double* %8
%16 = load i32, i32* %5
%17 = load i32, i32* %6
%18 = sub nsw i32 %16, %17
%19 = add nsw i32 %18, 1
store i32 %19, i32* %9
br label %20
20: ; preds = %37, %15
%21 = load i32, i32* %9
%22 = load i32, i32* %5
%23 = icmp sle i32 %21, %22
br i1 %23, label %24, label %40
24: ; preds = %20
%25 = load i32, i32* %9
%26 = load i32, i32* %5
%27 = sub nsw i32 %25, %26
%28 = load i32, i32* %7
%29 = add nsw i32 %27, %28
%30 = sitofp i32 %29 to double
store double %30, double* %10
%31 = load i32, i32* %9
%32 = sitofp i32 %31 to double
%33 = load double, double* %10
%34 = fdiv double %32, %33
%35 = load double, double* %8
%36 = fmul double %35, %34
store double %36, double* %8
br label %37
37: ; preds = %24
%38 = load i32, i32* %9
%39 = add nsw i32 %38, 1
store i32 %39, i32* %9
br label %20
40: ; preds = %20
%41 = load double, double* %8
store double %41, double* %4
br label %42
42: ; preds = %40, %14
%43 = load double, double* %4
ret double %43
}

LLVM IR (-O2):
define double @binomial(i32 %0, i32 %1, i32 %2) {
%4 = icmp slt i32 %0, %1
br i1 %4, label %20, label %5
5: ; preds = %3
%6 = icmp sgt i32 %1, 0
br i1 %6, label %7, label %20
7: ; preds = %5
%8 = sub nsw i32 %0, %1
%9 = sub i32 %2, %0
br label %10
10: ; preds = %7, %10
%11 = phi i32 [ %8, %7 ], [ %13, %10 ]
%12 = phi double [ 1.000000e+00, %7 ], [ %18, %10 ]
%13 = add nsw i32 %11, 1
%14 = add i32 %9, %13
%15 = sitofp i32 %14 to double
%16 = sitofp i32 %13 to double
%17 = fdiv double %16, %15
%18 = fmul double %12, %17
%19 = icmp slt i32 %13, %0
br i1 %19, label %10, label %20
20: ; preds = %10, %5, %3
%21 = phi double [ 0.000000e+00, %3 ], [ 1.000000e+00, %5 ], [ %18, %10 ]
ret double %21
}