#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 pow(int, double);
double BezierBlend(int k, double mu, int n) {
  int nn, kn, nkn;
  double blend = 1;
  nn = n;
  kn = k;
  nkn = n - k;
  while (nn >= 1) {
    blend *= nn;
    nn--;
    if (kn > 1) {
      blend /= (double)kn;
      kn--;
    }
    if (nkn > 1) {
      blend /= (double)nkn;
      nkn--;
    }
  }
  if (k > 0)
    blend *= pow(mu, (double)k);
  if (n - k > 0)
    blend *= pow(1 - mu, (double)(n - k));
  return (blend);
}


x86 assembly (clang, -O0):
<BezierBlend>:
push   %rbp
mov    %rsp,%rbp
sub    $0x30,%rsp
mov    %edi,-0x4(%rbp)
movsd  %xmm0,-0x10(%rbp)
mov    %esi,-0x14(%rbp)
movsd  0xed5(%rip),%xmm0
movsd  %xmm0,-0x28(%rbp)
mov    -0x14(%rbp),%eax
mov    %eax,-0x18(%rbp)
mov    -0x4(%rbp),%eax
mov    %eax,-0x1c(%rbp)
mov    -0x14(%rbp),%eax
sub    -0x4(%rbp),%eax
mov    %eax,-0x20(%rbp)
cmpl   $0x1,-0x18(%rbp)
jl     11b8 <BezierBlend+0xa8>
cvtsi2sdl -0x18(%rbp),%xmm0
mulsd  -0x28(%rbp),%xmm0
movsd  %xmm0,-0x28(%rbp)
mov    -0x18(%rbp),%eax
add    $0xffffffff,%eax
mov    %eax,-0x18(%rbp)
cmpl   $0x1,-0x1c(%rbp)
jle    118d <BezierBlend+0x7d>
cvtsi2sdl -0x1c(%rbp),%xmm1
movsd  -0x28(%rbp),%xmm0
divsd  %xmm1,%xmm0
movsd  %xmm0,-0x28(%rbp)
mov    -0x1c(%rbp),%eax
add    $0xffffffff,%eax
mov    %eax,-0x1c(%rbp)
cmpl   $0x1,-0x20(%rbp)
jle    11b3 <BezierBlend+0xa3>
cvtsi2sdl -0x20(%rbp),%xmm1
movsd  -0x28(%rbp),%xmm0
divsd  %xmm1,%xmm0
movsd  %xmm0,-0x28(%rbp)
mov    -0x20(%rbp),%eax
add    $0xffffffff,%eax
mov    %eax,-0x20(%rbp)
jmp    1145 <BezierBlend+0x35>
cmpl   $0x0,-0x4(%rbp)
jle    11db <BezierBlend+0xcb>
cvttsd2si -0x10(%rbp),%edi
cvtsi2sdl -0x4(%rbp),%xmm0
call   1030 <pow@plt>
mulsd  -0x28(%rbp),%xmm0
movsd  %xmm0,-0x28(%rbp)
mov    -0x14(%rbp),%eax
sub    -0x4(%rbp),%eax
cmp    $0x0,%eax
jle    1214 <BezierBlend+0x104>
movsd  0xe0e(%rip),%xmm0
subsd  -0x10(%rbp),%xmm0
cvttsd2si %xmm0,%edi
mov    -0x14(%rbp),%eax
sub    -0x4(%rbp),%eax
cvtsi2sd %eax,%xmm0
call   1030 <pow@plt>
mulsd  -0x28(%rbp),%xmm0
movsd  %xmm0,-0x28(%rbp)
movsd  -0x28(%rbp),%xmm0
add    $0x30,%rsp
pop    %rbp
ret

x86 assembly (clang, -O2):
<BezierBlend>:
push   %rbx
sub    $0x10,%rsp
mov    %esi,%ebx
sub    %edi,%ebx
test   %esi,%esi
jle    1172 <BezierBlend+0x62>
movsd  0xedb(%rip),%xmm1
mov    %ebx,%eax
mov    %edi,%ecx
jmp    113c <BezierBlend+0x2c>
nopl   0x0(%rax,%rax,1)
lea    0x1(%rsi),%edx
movapd %xmm2,%xmm1
cmp    $0x1,%edx
jbe    117a <BezierBlend+0x6a>
xorps  %xmm2,%xmm2
cvtsi2sd %esi,%xmm2
mulsd  %xmm1,%xmm2
cmp    $0x2,%ecx
jl     115a <BezierBlend+0x4a>
xorps  %xmm1,%xmm1
cvtsi2sd %ecx,%xmm1
divsd  %xmm1,%xmm2
add    $0xffffffff,%ecx
add    $0xffffffff,%esi
cmp    $0x2,%eax
jl     1130 <BezierBlend+0x20>
xorps  %xmm1,%xmm1
cvtsi2sd %eax,%xmm1
divsd  %xmm1,%xmm2
add    $0xffffffff,%eax
jmp    1130 <BezierBlend+0x20>
movsd  0xe86(%rip),%xmm2
test   %edi,%edi
movsd  %xmm0,0x8(%rsp)
jle    11a4 <BezierBlend+0x94>
cvttsd2si %xmm0,%eax
xorps  %xmm0,%xmm0
cvtsi2sd %edi,%xmm0
mov    %eax,%edi
movsd  %xmm2,(%rsp)
call   1030 <pow@plt>
movsd  (%rsp),%xmm2
mulsd  %xmm0,%xmm2
test   %ebx,%ebx
jle    11d4 <BezierBlend+0xc4>
movsd  0xe50(%rip),%xmm0
subsd  0x8(%rsp),%xmm0
cvttsd2si %xmm0,%edi
xorps  %xmm0,%xmm0
cvtsi2sd %ebx,%xmm0
movsd  %xmm2,(%rsp)
call   1030 <pow@plt>
movsd  (%rsp),%xmm2
mulsd  %xmm0,%xmm2
movapd %xmm2,%xmm0
add    $0x10,%rsp
pop    %rbx
ret

LLVM IR (-O0):
define double @BezierBlend(i32 %0, double %1, i32 %2) {
%4 = alloca i32
%5 = alloca double
%6 = alloca i32
%7 = alloca i32
%8 = alloca i32
%9 = alloca i32
%10 = alloca double
store i32 %0, i32* %4
store double %1, double* %5
store i32 %2, i32* %6
store double 1.000000e+00, double* %10
%11 = load i32, i32* %6
store i32 %11, i32* %7
%12 = load i32, i32* %4
store i32 %12, i32* %8
%13 = load i32, i32* %6
%14 = load i32, i32* %4
%15 = sub nsw i32 %13, %14
store i32 %15, i32* %9
br label %16
16: ; preds = %45, %3
%17 = load i32, i32* %7
%18 = icmp sge i32 %17, 1
br i1 %18, label %19, label %46
19: ; preds = %16
%20 = load i32, i32* %7
%21 = sitofp i32 %20 to double
%22 = load double, double* %10
%23 = fmul double %22, %21
store double %23, double* %10
%24 = load i32, i32* %7
%25 = add nsw i32 %24, -1
store i32 %25, i32* %7
%26 = load i32, i32* %8
%27 = icmp sgt i32 %26, 1
br i1 %27, label %28, label %35
28: ; preds = %19
%29 = load i32, i32* %8
%30 = sitofp i32 %29 to double
%31 = load double, double* %10
%32 = fdiv double %31, %30
store double %32, double* %10
%33 = load i32, i32* %8
%34 = add nsw i32 %33, -1
store i32 %34, i32* %8
br label %35
35: ; preds = %28, %19
%36 = load i32, i32* %9
%37 = icmp sgt i32 %36, 1
br i1 %37, label %38, label %45
38: ; preds = %35
%39 = load i32, i32* %9
%40 = sitofp i32 %39 to double
%41 = load double, double* %10
%42 = fdiv double %41, %40
store double %42, double* %10
%43 = load i32, i32* %9
%44 = add nsw i32 %43, -1
store i32 %44, i32* %9
br label %45
45: ; preds = %38, %35
br label %16
46: ; preds = %16
%47 = load i32, i32* %4
%48 = icmp sgt i32 %47, 0
br i1 %48, label %49, label %57
49: ; preds = %46
%50 = load double, double* %5
%51 = fptosi double %50 to i32
%52 = load i32, i32* %4
%53 = sitofp i32 %52 to double
%54 = call double @pow(i32 %51, double %53)
%55 = load double, double* %10
%56 = fmul double %55, %54
store double %56, double* %10
br label %57
57: ; preds = %49, %46
%58 = load i32, i32* %6
%59 = load i32, i32* %4
%60 = sub nsw i32 %58, %59
%61 = icmp sgt i32 %60, 0
br i1 %61, label %62, label %73
62: ; preds = %57
%63 = load double, double* %5
%64 = fsub double 1.000000e+00, %63
%65 = fptosi double %64 to i32
%66 = load i32, i32* %6
%67 = load i32, i32* %4
%68 = sub nsw i32 %66, %67
%69 = sitofp i32 %68 to double
%70 = call double @pow(i32 %65, double %69)
%71 = load double, double* %10
%72 = fmul double %71, %70
store double %72, double* %10
br label %73
73: ; preds = %62, %57
%74 = load double, double* %10
ret double %74
}

LLVM IR (-O2):
define double @BezierBlend(i32 %0, double %1, i32 %2) {
%4 = sub nsw i32 %2, %0
%5 = icmp sgt i32 %2, 0
br i1 %5, label %6, label %31
6: ; preds = %3, %27
%7 = phi double [ %29, %27 ], [ 1.000000e+00, %3 ]
%8 = phi i32 [ %28, %27 ], [ %4, %3 ]
%9 = phi i32 [ %20, %27 ], [ %0, %3 ]
%10 = phi i32 [ %13, %27 ], [ %2, %3 ]
%11 = sitofp i32 %10 to double
%12 = fmul double %7, %11
%13 = add nsw i32 %10, -1
%14 = icmp sgt i32 %9, 1
br i1 %14, label %15, label %19
15: ; preds = %6
%16 = sitofp i32 %9 to double
%17 = fdiv double %12, %16
%18 = add nsw i32 %9, -1
br label %19
19: ; preds = %15, %6
%20 = phi i32 [ %18, %15 ], [ %9, %6 ]
%21 = phi double [ %17, %15 ], [ %12, %6 ]
%22 = icmp sgt i32 %8, 1
br i1 %22, label %23, label %27
23: ; preds = %19
%24 = sitofp i32 %8 to double
%25 = fdiv double %21, %24
%26 = add nsw i32 %8, -1
br label %27
27: ; preds = %23, %19
%28 = phi i32 [ %26, %23 ], [ %8, %19 ]
%29 = phi double [ %25, %23 ], [ %21, %19 ]
%30 = icmp ugt i32 %10, 1
br i1 %30, label %6, label %31
31: ; preds = %27, %3
%32 = phi double [ 1.000000e+00, %3 ], [ %29, %27 ]
%33 = icmp sgt i32 %0, 0
br i1 %33, label %34, label %39
34: ; preds = %31
%35 = fptosi double %1 to i32
%36 = sitofp i32 %0 to double
%37 = tail call double @pow(i32 %35, double %36)
%38 = fmul double %32, %37
br label %39
39: ; preds = %34, %31
%40 = phi double [ %38, %34 ], [ %32, %31 ]
%41 = icmp sgt i32 %4, 0
br i1 %41, label %42, label %48
42: ; preds = %39
%43 = fsub double 1.000000e+00, %1
%44 = fptosi double %43 to i32
%45 = sitofp i32 %4 to double
%46 = tail call double @pow(i32 %44, double %45)
%47 = fmul double %40, %46
br label %48
48: ; preds = %42, %39
%49 = phi double [ %47, %42 ], [ %40, %39 ]
ret double %49
}