#define NULL ((void *)0)
typedef unsigned long size_t; // Customize by platform.
typedef long intptr_t;
typedef unsigned long uintptr_t;
typedef long scalar_t__; // Either arithmetic or pointer type.
/* By default, we understand bool (as a convenience). */
typedef int bool;
#define false 0
#define true 1

/* Forward declarations */

/* Type definitions */

/* Variables and functions */
double log(int);
double st_comb_ln(unsigned int n, unsigned int r) {
  if (r == n || r == 0)
    return 0;
  return (n - r) * log((double)n / (n - r)) + r * log((double)n / r) -
         0.5 * log(2 * 3.1415926 * r * (n - r) / n) +
         1.0 / 12 * (1.0 / n - 1.0 / (n - r) - 1.0 / r);
}


x86 assembly (clang, -O0):
<st_comb_ln>:
push   %rbp
mov    %rsp,%rbp
sub    $0x30,%rsp
mov    %edi,-0xc(%rbp)
mov    %esi,-0x10(%rbp)
mov    -0x10(%rbp),%eax
cmp    -0xc(%rbp),%eax
je     1134 <st_comb_ln+0x24>
cmpl   $0x0,-0x10(%rbp)
jne    1141 <st_comb_ln+0x31>
xorps  %xmm0,%xmm0
movsd  %xmm0,-0x8(%rbp)
jmp    1258 <st_comb_ln+0x148>
mov    -0xc(%rbp),%eax
mov    %eax,%ecx
mov    -0x10(%rbp),%edx
sub    %edx,%ecx
cvtsi2sd %rcx,%xmm1
movsd  %xmm1,-0x20(%rbp)
cvtsi2sd %rax,%xmm0
divsd  %xmm1,%xmm0
cvttsd2si %xmm0,%edi
call   1030 <log@plt>
movsd  %xmm0,-0x28(%rbp)
mov    -0x10(%rbp),%eax
cvtsi2sd %rax,%xmm1
movsd  %xmm1,-0x30(%rbp)
mov    -0xc(%rbp),%eax
cvtsi2sd %rax,%xmm0
divsd  %xmm1,%xmm0
cvttsd2si %xmm0,%edi
call   1030 <log@plt>
movsd  -0x30(%rbp),%xmm1
movsd  -0x28(%rbp),%xmm2
movaps %xmm0,%xmm3
movsd  -0x20(%rbp),%xmm0
mulsd  %xmm3,%xmm1
mulsd  %xmm2,%xmm0
addsd  %xmm1,%xmm0
movsd  %xmm0,-0x18(%rbp)
mov    -0x10(%rbp),%eax
mov    %eax,%edx
cvtsi2sd %rax,%xmm0
movsd  0xe4d(%rip),%xmm1
mulsd  %xmm1,%xmm0
mov    -0xc(%rbp),%eax
mov    %eax,%ecx
sub    %edx,%ecx
cvtsi2sd %rcx,%xmm1
mulsd  %xmm1,%xmm0
cvtsi2sd %rax,%xmm1
divsd  %xmm1,%xmm0
cvttsd2si %xmm0,%edi
call   1030 <log@plt>
movaps %xmm0,%xmm1
movsd  -0x18(%rbp),%xmm0
movsd  0xe1f(%rip),%xmm2
mulsd  %xmm2,%xmm1
addsd  %xmm0,%xmm1
mov    -0xc(%rbp),%eax
mov    %eax,%ecx
cvtsi2sd %rax,%xmm3
movsd  0xdf5(%rip),%xmm2
movaps %xmm2,%xmm0
divsd  %xmm3,%xmm0
mov    -0x10(%rbp),%eax
mov    %eax,%edx
sub    %edx,%ecx
cvtsi2sd %rcx,%xmm3
divsd  %xmm3,%xmm2
subsd  %xmm2,%xmm0
cvtsi2sd %rax,%xmm3
movsd  0xdcd(%rip),%xmm2
divsd  %xmm3,%xmm2
subsd  %xmm2,%xmm0
movsd  0xdb5(%rip),%xmm2
mulsd  %xmm2,%xmm0
addsd  %xmm1,%xmm0
movsd  %xmm0,-0x8(%rbp)
movsd  -0x8(%rbp),%xmm0
add    $0x30,%rsp
pop    %rbp
ret

x86 assembly (clang, -O2):
<st_comb_ln>:
mov    %edi,%eax
sub    %esi,%eax
je     1211 <st_comb_ln+0x101>
test   %esi,%esi
je     1211 <st_comb_ln+0x101>
sub    $0x38,%rsp
cvtsi2sd %rax,%xmm2
movsd  %xmm2,0x8(%rsp)
movd   %esi,%xmm0
movd   %edi,%xmm1
punpcklqdq %xmm0,%xmm1
movdqa 0xebb(%rip),%xmm0
por    %xmm0,%xmm1
subpd  %xmm0,%xmm1
movapd %xmm1,0x10(%rsp)
movapd %xmm1,%xmm0
divsd  %xmm2,%xmm0
cvttsd2si %xmm0,%edi
call   1030 <log@plt>
movsd  %xmm0,(%rsp)
movapd 0x10(%rsp),%xmm0
movapd %xmm0,%xmm1
unpckhpd %xmm0,%xmm1
movapd %xmm1,0x20(%rsp)
divsd  %xmm1,%xmm0
cvttsd2si %xmm0,%edi
call   1030 <log@plt>
movapd 0x20(%rsp),%xmm2
mulsd  %xmm2,%xmm0
movsd  0x8(%rsp),%xmm1
movsd  (%rsp),%xmm3
mulsd  %xmm1,%xmm3
addsd  %xmm3,%xmm0
movsd  %xmm0,(%rsp)
mulsd  0xe6c(%rip),%xmm2
mulsd  %xmm1,%xmm2
divsd  0x10(%rsp),%xmm2
cvttsd2si %xmm2,%edi
call   1030 <log@plt>
mulsd  0xe59(%rip),%xmm0
addsd  (%rsp),%xmm0
movsd  0xe54(%rip),%xmm1
divsd  0x8(%rsp),%xmm1
movapd 0xe26(%rip),%xmm2
divpd  0x10(%rsp),%xmm2
movapd %xmm2,%xmm3
subsd  %xmm1,%xmm3
unpckhpd %xmm2,%xmm2
subsd  %xmm2,%xmm3
mulsd  0xe30(%rip),%xmm3
addsd  %xmm3,%xmm0
add    $0x38,%rsp
ret
xorps  %xmm0,%xmm0
ret

LLVM IR (-O0):
define double @st_comb_ln(i32 %0, i32 %1) {
%3 = alloca double
%4 = alloca i32
%5 = alloca i32
store i32 %0, i32* %4
store i32 %1, i32* %5
%6 = load i32, i32* %5
%7 = load i32, i32* %4
%8 = icmp eq i32 %6, %7
br i1 %8, label %12, label %9
9: ; preds = %2
%10 = load i32, i32* %5
%11 = icmp eq i32 %10, 0
br i1 %11, label %12, label %13
12: ; preds = %9, %2
store double 0.000000e+00, double* %3
br label %66
13: ; preds = %9
%14 = load i32, i32* %4
%15 = load i32, i32* %5
%16 = sub i32 %14, %15
%17 = uitofp i32 %16 to double
%18 = load i32, i32* %4
%19 = uitofp i32 %18 to double
%20 = load i32, i32* %4
%21 = load i32, i32* %5
%22 = sub i32 %20, %21
%23 = uitofp i32 %22 to double
%24 = fdiv double %19, %23
%25 = fptosi double %24 to i32
%26 = call double @log(i32 %25)
%27 = load i32, i32* %5
%28 = uitofp i32 %27 to double
%29 = load i32, i32* %4
%30 = uitofp i32 %29 to double
%31 = load i32, i32* %5
%32 = uitofp i32 %31 to double
%33 = fdiv double %30, %32
%34 = fptosi double %33 to i32
%35 = call double @log(i32 %34)
%36 = fmul double %28, %35
%37 = call double @llvm.fmuladd.f64(double %17, double %26, double %36)
%38 = load i32, i32* %5
%39 = uitofp i32 %38 to double
%40 = fmul double 0x401921FB4D12D84A, %39
%41 = load i32, i32* %4
%42 = load i32, i32* %5
%43 = sub i32 %41, %42
%44 = uitofp i32 %43 to double
%45 = fmul double %40, %44
%46 = load i32, i32* %4
%47 = uitofp i32 %46 to double
%48 = fdiv double %45, %47
%49 = fptosi double %48 to i32
%50 = call double @log(i32 %49)
%51 = call double @llvm.fmuladd.f64(double -5.000000e-01, double %50, double %37)
%52 = load i32, i32* %4
%53 = uitofp i32 %52 to double
%54 = fdiv double 1.000000e+00, %53
%55 = load i32, i32* %4
%56 = load i32, i32* %5
%57 = sub i32 %55, %56
%58 = uitofp i32 %57 to double
%59 = fdiv double 1.000000e+00, %58
%60 = fsub double %54, %59
%61 = load i32, i32* %5
%62 = uitofp i32 %61 to double
%63 = fdiv double 1.000000e+00, %62
%64 = fsub double %60, %63
%65 = call double @llvm.fmuladd.f64(double 0x3FB5555555555555, double %64, double %51)
store double %65, double* %3
br label %66
66: ; preds = %13, %12
%67 = load double, double* %3
ret double %67
}

LLVM IR (-O2):
define double @st_comb_ln(i32 %0, i32 %1) {
%3 = icmp eq i32 %0, %1
%4 = icmp eq i32 %1, 0
%5 = or i1 %3, %4
br i1 %5, label %35, label %6
6: ; preds = %2
%7 = sub i32 %0, %1
%8 = uitofp i32 %7 to double
%9 = insertelement <2 x i32> poison, i32 %0, i64 0
%10 = insertelement <2 x i32> %9, i32 %1, i64 1
%11 = uitofp <2 x i32> %10 to <2 x double>
%12 = extractelement <2 x double> %11, i64 0
%13 = fdiv double %12, %8
%14 = fptosi double %13 to i32
%15 = tail call double @log(i32 %14)
%16 = extractelement <2 x double> %11, i64 1
%17 = fdiv double %12, %16
%18 = fptosi double %17 to i32
%19 = tail call double @log(i32 %18)
%20 = fmul double %19, %16
%21 = tail call double @llvm.fmuladd.f64(double %8, double %15, double %20)
%22 = fmul double %16, 0x401921FB4D12D84A
%23 = fmul double %22, %8
%24 = fdiv double %23, %12
%25 = fptosi double %24 to i32
%26 = tail call double @log(i32 %25)
%27 = tail call double @llvm.fmuladd.f64(double %26, double -5.000000e-01, double %21)
%28 = fdiv double 1.000000e+00, %8
%29 = fdiv <2 x double> <double 1.000000e+00, double 1.000000e+00>, %11
%30 = extractelement <2 x double> %29, i64 0
%31 = fsub double %30, %28
%32 = extractelement <2 x double> %29, i64 1
%33 = fsub double %31, %32
%34 = tail call double @llvm.fmuladd.f64(double %33, double 0x3FB5555555555555, double %27)
br label %35
35: ; preds = %2, %6
%36 = phi double [ %34, %6 ], [ 0.000000e+00, %2 ]
ret double %36
}