#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);
int probBinomialDistribution(int n, double p, double prob[]) {
  int k;
  double y;
  prob[0] = y = pow(1 - p, (double)n);
  for (k = 1; k <= n; k++)
    prob[k] = y *= ((n - k + 1) * p) / (k * (1 - p));
  return 0;
}


x86 assembly (clang, -O0):
<probBinomialDistribution>:
push   %rbp
mov    %rsp,%rbp
sub    $0x30,%rsp
mov    %edi,-0x4(%rbp)
movsd  %xmm0,-0x10(%rbp)
mov    %rsi,-0x18(%rbp)
movsd  0xed4(%rip),%xmm0
subsd  -0x10(%rbp),%xmm0
cvttsd2si %xmm0,%edi
cvtsi2sdl -0x4(%rbp),%xmm0
call   1030 <pow@plt>
movsd  %xmm0,-0x28(%rbp)
mov    -0x18(%rbp),%rax
movsd  %xmm0,(%rax)
movl   $0x1,-0x1c(%rbp)
mov    -0x1c(%rbp),%eax
cmp    -0x4(%rbp),%eax
jg     11b0 <probBinomialDistribution+0xa0>
mov    -0x4(%rbp),%eax
sub    -0x1c(%rbp),%eax
add    $0x1,%eax
cvtsi2sd %eax,%xmm0
mulsd  -0x10(%rbp),%xmm0
cvtsi2sdl -0x1c(%rbp),%xmm1
movsd  0xe82(%rip),%xmm2
subsd  -0x10(%rbp),%xmm2
mulsd  %xmm2,%xmm1
divsd  %xmm1,%xmm0
mulsd  -0x28(%rbp),%xmm0
movsd  %xmm0,-0x28(%rbp)
mov    -0x18(%rbp),%rax
movslq -0x1c(%rbp),%rcx
movsd  %xmm0,(%rax,%rcx,8)
mov    -0x1c(%rbp),%eax
add    $0x1,%eax
mov    %eax,-0x1c(%rbp)
jmp    1153 <probBinomialDistribution+0x43>
xor    %eax,%eax
add    $0x30,%rsp
pop    %rbp
ret

x86 assembly (clang, -O2):
<probBinomialDistribution>:
push   %r14
push   %rbx
sub    $0x18,%rsp
mov    %rsi,%r14
mov    %edi,%ebx
movsd  0xedc(%rip),%xmm1
movsd  %xmm0,0x8(%rsp)
subsd  %xmm0,%xmm1
movsd  %xmm1,0x10(%rsp)
cvttsd2si %xmm1,%edi
xorps  %xmm0,%xmm0
cvtsi2sd %ebx,%xmm0
call   1030 <pow@plt>
movsd  %xmm0,(%r14)
test   %ebx,%ebx
jle    120f <probBinomialDistribution+0xff>
lea    0x1(%rbx),%ecx
mov    %ebx,%edx
mov    $0x1,%eax
cmp    $0x2,%ecx
je     11dc <probBinomialDistribution+0xcc>
mov    %edx,%esi
and    $0xfffffffe,%esi
mov    $0x1,%eax
movsd  0x10(%rsp),%xmm3
movsd  0x8(%rsp),%xmm4
nopw   0x0(%rax,%rax,1)
xorps  %xmm1,%xmm1
cvtsi2sd %ebx,%xmm1
mulsd  %xmm4,%xmm1
xorps  %xmm2,%xmm2
cvtsi2sd %eax,%xmm2
mulsd  %xmm3,%xmm2
divsd  %xmm2,%xmm1
mulsd  %xmm0,%xmm1
lea    -0x1(%rbx),%edi
xorps  %xmm0,%xmm0
cvtsi2sd %edi,%xmm0
movsd  %xmm1,(%r14,%rax,8)
lea    0x1(%rax),%rdi
xorps  %xmm2,%xmm2
cvtsi2sd %edi,%xmm2
mulsd  %xmm4,%xmm0
mulsd  %xmm3,%xmm2
divsd  %xmm2,%xmm0
mulsd  %xmm1,%xmm0
movsd  %xmm0,0x8(%r14,%rax,8)
add    $0x2,%rax
add    $0xfffffffe,%ebx
cmp    %rsi,%rdi
jne    1180 <probBinomialDistribution+0x70>
test   $0x1,%dl
je     120f <probBinomialDistribution+0xff>
sub    %eax,%ecx
xorps  %xmm1,%xmm1
cvtsi2sd %ecx,%xmm1
mulsd  0x8(%rsp),%xmm1
xorps  %xmm2,%xmm2
cvtsi2sd %eax,%xmm2
movsd  0x10(%rsp),%xmm3
mulsd  %xmm2,%xmm3
divsd  %xmm3,%xmm1
mulsd  %xmm1,%xmm0
movsd  %xmm0,(%r14,%rax,8)
xor    %eax,%eax
add    $0x18,%rsp
pop    %rbx
pop    %r14
ret

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

LLVM IR (-O2):
define i32 @probBinomialDistribution(i32 %0, double %1, double* nocapture writeonly %2) {
%4 = fsub double 1.000000e+00, %1
%5 = fptosi double %4 to i32
%6 = sitofp i32 %0 to double
%7 = tail call double @pow(i32 %5, double %6)
store double %7, double* %2
%8 = icmp slt i32 %0, 1
br i1 %8, label %56, label %9
9: ; preds = %3
%10 = add nuw i32 %0, 1
%11 = zext i32 %0 to i64
%12 = and i64 %11, 1
%13 = icmp eq i32 %10, 2
br i1 %13, label %42, label %14
14: ; preds = %9
%15 = and i64 %11, 4294967294
br label %16
16: ; preds = %16, %14
%17 = phi i64 [ 1, %14 ], [ %39, %16 ]
%18 = phi double [ %7, %14 ], [ %37, %16 ]
%19 = phi i64 [ 0, %14 ], [ %40, %16 ]
%20 = trunc i64 %17 to i32
%21 = sub i32 %10, %20
%22 = sitofp i32 %21 to double
%23 = fmul double %22, %1
%24 = sitofp i32 %20 to double
%25 = fmul double %4, %24
%26 = fdiv double %23, %25
%27 = fmul double %18, %26
%28 = getelementptr inbounds double, double* %2, i64 %17
store double %27, double* %28
%29 = add nuw nsw i64 %17, 1
%30 = trunc i64 %29 to i32
%31 = sub i32 %10, %30
%32 = sitofp i32 %31 to double
%33 = fmul double %32, %1
%34 = sitofp i32 %30 to double
%35 = fmul double %4, %34
%36 = fdiv double %33, %35
%37 = fmul double %27, %36
%38 = getelementptr inbounds double, double* %2, i64 %29
store double %37, double* %38
%39 = add nuw nsw i64 %17, 2
%40 = add i64 %19, 2
%41 = icmp eq i64 %40, %15
br i1 %41, label %42, label %16
42: ; preds = %16, %9
%43 = phi i64 [ 1, %9 ], [ %39, %16 ]
%44 = phi double [ %7, %9 ], [ %37, %16 ]
%45 = icmp eq i64 %12, 0
br i1 %45, label %56, label %46
46: ; preds = %42
%47 = trunc i64 %43 to i32
%48 = sub i32 %10, %47
%49 = sitofp i32 %48 to double
%50 = fmul double %49, %1
%51 = sitofp i32 %47 to double
%52 = fmul double %4, %51
%53 = fdiv double %50, %52
%54 = fmul double %44, %53
%55 = getelementptr inbounds double, double* %2, i64 %43
store double %54, double* %55
br label %56
56: ; preds = %46, %42, %3
ret i32 0
}