#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 BetaRand(int, int);
double DBL_EPSILON;
int Rand();
scalar_t__ isnan(double);
double log(int);
int BinoRand(double p, int n) {
  int r = 0;
  if (isnan(p))
    return 0;
  if (p < DBL_EPSILON)
    return 0;
  if (p >= 1 - DBL_EPSILON)
    return n;
  if ((p > 0.5) && (n < 15)) {
    int i;
    for (i = 0; i < n; i++) {
      if (Rand() < p)
        r++;
    }
    return r;
  }
  if (n * p < 10) {
    double q = -log(1 - p), e = -log(Rand()), s;
    r = n;
    for (s = e / r; s <= q; s += e / r) {
      r--;
      if (r == 0)
        break;
      e = -log(Rand());
    }
    r = n - r;
    return r;
  }
  if (1) {
    int i = (int)(p * (n + 1));
    double b = BetaRand(i, n + 1 - i);
    if (b <= p)
      r = i + BinoRand((p - b) / (1 - b), n - i);
    else
      r = i - 1 - BinoRand((b - p) / b, i - 1);
    return r;
  }
}


x86 assembly (clang, -O0):
<BinoRand>:
push   %rbp
mov    %rsp,%rbp
sub    $0x50,%rsp
movsd  %xmm0,-0x10(%rbp)
mov    %edi,-0x14(%rbp)
movl   $0x0,-0x18(%rbp)
movsd  -0x10(%rbp),%xmm0
call   1050 <isnan@plt>
cmp    $0x0,%rax
je     1187 <BinoRand+0x37>
movl   $0x0,-0x4(%rbp)
jmp    1419 <BinoRand+0x2c9>
movsd  -0x10(%rbp),%xmm1
mov    0x2e55(%rip),%rax
movsd  (%rax),%xmm0
ucomisd %xmm1,%xmm0
jbe    11ad <BinoRand+0x5d>
movl   $0x0,-0x4(%rbp)
jmp    1419 <BinoRand+0x2c9>
movsd  -0x10(%rbp),%xmm0
movsd  0xe46(%rip),%xmm1
mov    0x2e27(%rip),%rax
subsd  (%rax),%xmm1
ucomisd %xmm1,%xmm0
jb     11da <BinoRand+0x8a>
mov    -0x14(%rbp),%eax
mov    %eax,-0x4(%rbp)
jmp    1419 <BinoRand+0x2c9>
movsd  -0x10(%rbp),%xmm0
movsd  0xe21(%rip),%xmm1
ucomisd %xmm1,%xmm0
jbe    124f <BinoRand+0xff>
cmpl   $0xf,-0x14(%rbp)
jge    124f <BinoRand+0xff>
movl   $0x0,-0x1c(%rbp)
mov    -0x1c(%rbp),%eax
cmp    -0x14(%rbp),%eax
jge    1244 <BinoRand+0xf4>
mov    $0x0,%al
call   1070 <Rand@plt>
cvtsi2sd %eax,%xmm1
movsd  -0x10(%rbp),%xmm0
ucomisd %xmm1,%xmm0
jbe    1231 <BinoRand+0xe1>
mov    -0x18(%rbp),%eax
add    $0x1,%eax
mov    %eax,-0x18(%rbp)
jmp    1236 <BinoRand+0xe6>
mov    -0x1c(%rbp),%eax
add    $0x1,%eax
mov    %eax,-0x1c(%rbp)
jmp    1202 <BinoRand+0xb2>
mov    -0x18(%rbp),%eax
mov    %eax,-0x4(%rbp)
jmp    1419 <BinoRand+0x2c9>
cvtsi2sdl -0x14(%rbp),%xmm1
mulsd  -0x10(%rbp),%xmm1
movsd  0xdaf(%rip),%xmm0
ucomisd %xmm1,%xmm0
jbe    1367 <BinoRand+0x217>
movsd  0xd8d(%rip),%xmm0
subsd  -0x10(%rbp),%xmm0
cvttsd2si %xmm0,%edi
call   1060 <log@plt>
movq   %xmm0,%rax
movabs $0x8000000000000000,%rcx
xor    %rcx,%rax
movq   %rax,%xmm0
movsd  %xmm0,-0x28(%rbp)
mov    $0x0,%al
call   1070 <Rand@plt>
mov    %eax,%edi
call   1060 <log@plt>
movq   %xmm0,%rax
movabs $0x8000000000000000,%rcx
xor    %rcx,%rax
movq   %rax,%xmm0
movsd  %xmm0,-0x30(%rbp)
mov    -0x14(%rbp),%eax
mov    %eax,-0x18(%rbp)
movsd  -0x30(%rbp),%xmm0
cvtsi2sdl -0x18(%rbp),%xmm1
divsd  %xmm1,%xmm0
movsd  %xmm0,-0x38(%rbp)
movsd  -0x38(%rbp),%xmm1
movsd  -0x28(%rbp),%xmm0
ucomisd %xmm1,%xmm0
jb     1353 <BinoRand+0x203>
mov    -0x18(%rbp),%eax
add    $0xffffffff,%eax
mov    %eax,-0x18(%rbp)
cmpl   $0x0,-0x18(%rbp)
jne    130c <BinoRand+0x1bc>
jmp    1353 <BinoRand+0x203>
mov    $0x0,%al
call   1070 <Rand@plt>
mov    %eax,%edi
call   1060 <log@plt>
movq   %xmm0,%rax
movabs $0x8000000000000000,%rcx
xor    %rcx,%rax
movq   %rax,%xmm0
movsd  %xmm0,-0x30(%rbp)
movsd  -0x30(%rbp),%xmm0
cvtsi2sdl -0x18(%rbp),%xmm1
divsd  %xmm1,%xmm0
addsd  -0x38(%rbp),%xmm0
movsd  %xmm0,-0x38(%rbp)
jmp    12e0 <BinoRand+0x190>
mov    -0x14(%rbp),%eax
sub    -0x18(%rbp),%eax
mov    %eax,-0x18(%rbp)
mov    -0x18(%rbp),%eax
mov    %eax,-0x4(%rbp)
jmp    1419 <BinoRand+0x2c9>
movsd  -0x10(%rbp),%xmm0
mov    -0x14(%rbp),%eax
add    $0x1,%eax
cvtsi2sd %eax,%xmm1
mulsd  %xmm1,%xmm0
cvttsd2si %xmm0,%eax
mov    %eax,-0x3c(%rbp)
mov    -0x3c(%rbp),%edi
mov    -0x14(%rbp),%esi
add    $0x1,%esi
sub    -0x3c(%rbp),%esi
call   1030 <BetaRand@plt>
movsd  %xmm0,-0x48(%rbp)
movsd  -0x48(%rbp),%xmm1
movsd  -0x10(%rbp),%xmm0
ucomisd %xmm1,%xmm0
jb     13e6 <BinoRand+0x296>
mov    -0x3c(%rbp),%eax
mov    %eax,-0x4c(%rbp)
movsd  -0x10(%rbp),%xmm0
subsd  -0x48(%rbp),%xmm0
movsd  0xc3d(%rip),%xmm1
subsd  -0x48(%rbp),%xmm1
divsd  %xmm1,%xmm0
mov    -0x14(%rbp),%edi
sub    -0x3c(%rbp),%edi
call   1040 <BinoRand@plt>
mov    %eax,%ecx
mov    -0x4c(%rbp),%eax
add    %ecx,%eax
mov    %eax,-0x18(%rbp)
jmp    1413 <BinoRand+0x2c3>
mov    -0x3c(%rbp),%eax
sub    $0x1,%eax
mov    %eax,-0x50(%rbp)
movsd  -0x48(%rbp),%xmm0
subsd  -0x10(%rbp),%xmm0
divsd  -0x48(%rbp),%xmm0
mov    -0x3c(%rbp),%edi
sub    $0x1,%edi
call   1040 <BinoRand@plt>
mov    %eax,%ecx
mov    -0x50(%rbp),%eax
sub    %ecx,%eax
mov    %eax,-0x18(%rbp)
mov    -0x18(%rbp),%eax
mov    %eax,-0x4(%rbp)
mov    -0x4(%rbp),%eax
add    $0x50,%rsp
pop    %rbp
ret

x86 assembly (clang, -O2):
<BinoRand>:
push   %rbp
push   %r15
push   %r14
push   %r12
push   %rbx
sub    $0x20,%rsp
mov    %edi,%r15d
movsd  %xmm0,(%rsp)
call   1050 <isnan@plt>
xor    %r14d,%r14d
mov    $0x0,%ebp
test   %rax,%rax
je     1189 <BinoRand+0x39>
add    %ebp,%r14d
mov    %r14d,%eax
add    $0x20,%rsp
pop    %rbx
pop    %r12
pop    %r14
pop    %r15
pop    %rbp
ret
xor    %ebp,%ebp
mov    0x2e56(%rip),%r12
xor    %r14d,%r14d
cs nopw 0x0(%rax,%rax,1)
nop
movsd  (%r12),%xmm0
movsd  (%rsp),%xmm1
ucomisd %xmm1,%xmm0
ja     1176 <BinoRand+0x26>
movapd %xmm1,%xmm2
movsd  0xe43(%rip),%xmm1
subsd  %xmm0,%xmm1
ucomisd %xmm1,%xmm2
jae    1353 <BinoRand+0x203>
ucomisd 0xe35(%rip),%xmm2
jbe    11df <BinoRand+0x8f>
cmp    $0xe,%r15d
jle    1262 <BinoRand+0x112>
xorps  %xmm3,%xmm3
cvtsi2sd %r15d,%xmm3
movapd %xmm2,%xmm0
mulsd  %xmm3,%xmm0
movsd  0xe19(%rip),%xmm1
ucomisd %xmm0,%xmm1
ja     12a9 <BinoRand+0x159>
lea    0x1(%r15),%esi
xorps  %xmm0,%xmm0
cvtsi2sd %esi,%xmm0
mulsd  %xmm2,%xmm0
cvttsd2si %xmm0,%ebx
sub    %ebx,%esi
mov    %ebx,%edi
call   1030 <BetaRand@plt>
movsd  (%rsp),%xmm2
ucomisd %xmm0,%xmm2
jb     135b <BinoRand+0x20b>
subsd  %xmm0,%xmm2
movsd  0xdc8(%rip),%xmm1
subsd  %xmm0,%xmm1
divsd  %xmm1,%xmm2
movsd  %xmm2,(%rsp)
sub    %ebx,%r15d
add    %ebx,%r14d
movapd %xmm2,%xmm0
call   1050 <isnan@plt>
test   %rax,%rax
je     11a0 <BinoRand+0x50>
jmp    1176 <BinoRand+0x26>
test   %r15d,%r15d
jle    1176 <BinoRand+0x26>
xor    %ebp,%ebp
movsd  %xmm2,(%rsp)
cs nopw 0x0(%rax,%rax,1)
nopl   0x0(%rax)
xor    %eax,%eax
call   1070 <Rand@plt>
movsd  (%rsp),%xmm2
xorps  %xmm0,%xmm0
cvtsi2sd %eax,%xmm0
xor    %eax,%eax
ucomisd %xmm0,%xmm2
seta   %al
add    %eax,%ebp
add    $0xffffffff,%r15d
jne    1280 <BinoRand+0x130>
jmp    1176 <BinoRand+0x26>
movsd  0xd4f(%rip),%xmm0
subsd  %xmm2,%xmm0
cvttsd2si %xmm0,%edi
movsd  %xmm3,(%rsp)
call   1060 <log@plt>
xorpd  0xd55(%rip),%xmm0
movapd %xmm0,0x10(%rsp)
xor    %eax,%eax
call   1070 <Rand@plt>
mov    %eax,%edi
call   1060 <log@plt>
movapd %xmm0,%xmm1
xorpd  0xd35(%rip),%xmm1
divsd  (%rsp),%xmm1
movapd 0x10(%rsp),%xmm0
ucomisd %xmm1,%xmm0
mov    %r15d,%r12d
jb     1350 <BinoRand+0x200>
lea    -0x1(%r15),%ebx
xor    %r12d,%r12d
mov    %ebx,%ebp
nopl   0x0(%rax,%rax,1)
sub    $0x1,%ebp
jb     1350 <BinoRand+0x200>
xor    %eax,%eax
movapd %xmm1,(%rsp)
call   1070 <Rand@plt>
mov    %eax,%edi
call   1060 <log@plt>
xorps  %xmm1,%xmm1
cvtsi2sd %ebx,%xmm1
divsd  %xmm1,%xmm0
movapd (%rsp),%xmm1
subsd  %xmm0,%xmm1
movapd 0x10(%rsp),%xmm0
ucomisd %xmm1,%xmm0
mov    %ebp,%ebx
jae    1310 <BinoRand+0x1c0>
add    $0x1,%ebp
mov    %ebp,%r12d
sub    %r12d,%r15d
mov    %r15d,%ebp
jmp    1176 <BinoRand+0x26>
lea    -0x1(%rbx),%edi
movapd %xmm0,%xmm1
subsd  %xmm2,%xmm1
divsd  %xmm0,%xmm1
movapd %xmm1,%xmm0
call   1040 <BinoRand@plt>
add    %ebx,%r14d
add    $0xffffffff,%r14d
sub    %eax,%r14d
jmp    1179 <BinoRand+0x29>

LLVM IR (-O0):
@DBL_EPSILON = global double 0.000000e+00
define i32 @BinoRand(double %0, i32 %1) {
%3 = alloca i32
%4 = alloca double
%5 = alloca i32
%6 = alloca i32
%7 = alloca i32
%8 = alloca double
%9 = alloca double
%10 = alloca double
%11 = alloca i32
%12 = alloca double
store double %0, double* %4
store i32 %1, i32* %5
store i32 0, i32* %6
%13 = load double, double* %4
%14 = call i64 @isnan(double %13)
%15 = icmp ne i64 %14, 0
br i1 %15, label %16, label %17
16: ; preds = %2
store i32 0, i32* %3
br label %143
17: ; preds = %2
%18 = load double, double* %4
%19 = load double, double* @DBL_EPSILON
%20 = fcmp olt double %18, %19
br i1 %20, label %21, label %22
21: ; preds = %17
store i32 0, i32* %3
br label %143
22: ; preds = %17
%23 = load double, double* %4
%24 = load double, double* @DBL_EPSILON
%25 = fsub double 1.000000e+00, %24
%26 = fcmp oge double %23, %25
br i1 %26, label %27, label %29
27: ; preds = %22
%28 = load i32, i32* %5
store i32 %28, i32* %3
br label %143
29: ; preds = %22
%30 = load double, double* %4
%31 = fcmp ogt double %30, 5.000000e-01
br i1 %31, label %32, label %54
32: ; preds = %29
%33 = load i32, i32* %5
%34 = icmp slt i32 %33, 15
br i1 %34, label %35, label %54
35: ; preds = %32
store i32 0, i32* %7
br label %36
36: ; preds = %49, %35
%37 = load i32, i32* %7
%38 = load i32, i32* %5
%39 = icmp slt i32 %37, %38
br i1 %39, label %40, label %52
40: ; preds = %36
%41 = call i32 (...) @Rand()
%42 = sitofp i32 %41 to double
%43 = load double, double* %4
%44 = fcmp olt double %42, %43
br i1 %44, label %45, label %48
45: ; preds = %40
%46 = load i32, i32* %6
%47 = add nsw i32 %46, 1
store i32 %47, i32* %6
br label %48
48: ; preds = %45, %40
br label %49
49: ; preds = %48
%50 = load i32, i32* %7
%51 = add nsw i32 %50, 1
store i32 %51, i32* %7
br label %36
52: ; preds = %36
%53 = load i32, i32* %6
store i32 %53, i32* %3
br label %143
54: ; preds = %32, %29
%55 = load i32, i32* %5
%56 = sitofp i32 %55 to double
%57 = load double, double* %4
%58 = fmul double %56, %57
%59 = fcmp olt double %58, 1.000000e+01
br i1 %59, label %60, label %100
60: ; preds = %54
%61 = load double, double* %4
%62 = fsub double 1.000000e+00, %61
%63 = fptosi double %62 to i32
%64 = call double @log(i32 %63)
%65 = fneg double %64
store double %65, double* %8
%66 = call i32 (...) @Rand()
%67 = call double @log(i32 %66)
%68 = fneg double %67
store double %68, double* %9
%69 = load i32, i32* %5
store i32 %69, i32* %6
%70 = load double, double* %9
%71 = load i32, i32* %6
%72 = sitofp i32 %71 to double
%73 = fdiv double %70, %72
store double %73, double* %10
br label %74
74: ; preds = %88, %60
%75 = load double, double* %10
%76 = load double, double* %8
%77 = fcmp ole double %75, %76
br i1 %77, label %78, label %95
78: ; preds = %74
%79 = load i32, i32* %6
%80 = add nsw i32 %79, -1
store i32 %80, i32* %6
%81 = load i32, i32* %6
%82 = icmp eq i32 %81, 0
br i1 %82, label %83, label %84
83: ; preds = %78
br label %95
84: ; preds = %78
%85 = call i32 (...) @Rand()
%86 = call double @log(i32 %85)
%87 = fneg double %86
store double %87, double* %9
br label %88
88: ; preds = %84
%89 = load double, double* %9
%90 = load i32, i32* %6
%91 = sitofp i32 %90 to double
%92 = fdiv double %89, %91
%93 = load double, double* %10
%94 = fadd double %93, %92
store double %94, double* %10
br label %74
95: ; preds = %83, %74
%96 = load i32, i32* %5
%97 = load i32, i32* %6
%98 = sub nsw i32 %96, %97
store i32 %98, i32* %6
%99 = load i32, i32* %6
store i32 %99, i32* %3
br label %143
100: ; preds = %54
%101 = load double, double* %4
%102 = load i32, i32* %5
%103 = add nsw i32 %102, 1
%104 = sitofp i32 %103 to double
%105 = fmul double %101, %104
%106 = fptosi double %105 to i32
store i32 %106, i32* %11
%107 = load i32, i32* %11
%108 = load i32, i32* %5
%109 = add nsw i32 %108, 1
%110 = load i32, i32* %11
%111 = sub nsw i32 %109, %110
%112 = call double @BetaRand(i32 %107, i32 %111)
store double %112, double* %12
%113 = load double, double* %12
%114 = load double, double* %4
%115 = fcmp ole double %113, %114
br i1 %115, label %116, label %129
116: ; preds = %100
%117 = load i32, i32* %11
%118 = load double, double* %4
%119 = load double, double* %12
%120 = fsub double %118, %119
%121 = load double, double* %12
%122 = fsub double 1.000000e+00, %121
%123 = fdiv double %120, %122
%124 = load i32, i32* %5
%125 = load i32, i32* %11
%126 = sub nsw i32 %124, %125
%127 = call i32 @BinoRand(double %123, i32 %126)
%128 = add nsw i32 %117, %127
store i32 %128, i32* %6
br label %141
129: ; preds = %100
%130 = load i32, i32* %11
%131 = sub nsw i32 %130, 1
%132 = load double, double* %12
%133 = load double, double* %4
%134 = fsub double %132, %133
%135 = load double, double* %12
%136 = fdiv double %134, %135
%137 = load i32, i32* %11
%138 = sub nsw i32 %137, 1
%139 = call i32 @BinoRand(double %136, i32 %138)
%140 = sub nsw i32 %131, %139
store i32 %140, i32* %6
br label %141
141: ; preds = %129, %116
%142 = load i32, i32* %6
store i32 %142, i32* %3
br label %143
143: ; preds = %141, %95, %52, %27, %21, %16
%144 = load i32, i32* %3
ret i32 %144
}

LLVM IR (-O2):
@DBL_EPSILON = global double 0.000000e+00
define i32 @BinoRand(double %0, i32 %1) {
%3 = tail call i64 @isnan(double %0)
%4 = icmp eq i64 %3, 0
br i1 %4, label %5, label %84
5: ; preds = %2, %67
%6 = phi i32 [ %71, %67 ], [ %1, %2 ]
%7 = phi double [ %70, %67 ], [ %0, %2 ]
%8 = phi i32 [ %72, %67 ], [ 0, %2 ]
%9 = load double, double* @DBL_EPSILON
%10 = fcmp ogt double %9, %7
br i1 %10, label %84, label %11
11: ; preds = %5
%12 = fsub double 1.000000e+00, %9
%13 = fcmp ugt double %12, %7
br i1 %13, label %14, label %84
14: ; preds = %11
%15 = fcmp ogt double %7, 5.000000e-01
%16 = icmp slt i32 %6, 15
%17 = and i1 %15, %16
br i1 %17, label %18, label %30
18: ; preds = %14
%19 = icmp sgt i32 %6, 0
br i1 %19, label %20, label %84
20: ; preds = %18, %20
%21 = phi i32 [ %27, %20 ], [ 0, %18 ]
%22 = phi i32 [ %28, %20 ], [ 0, %18 ]
%23 = tail call i32 (...) @Rand()
%24 = sitofp i32 %23 to double
%25 = fcmp ogt double %7, %24
%26 = zext i1 %25 to i32
%27 = add nuw nsw i32 %21, %26
%28 = add nuw nsw i32 %22, 1
%29 = icmp eq i32 %28, %6
br i1 %29, label %84, label %20
30: ; preds = %14
%31 = sitofp i32 %6 to double
%32 = fmul double %7, %31
%33 = fcmp olt double %32, 1.000000e+01
br i1 %33, label %34, label %59
34: ; preds = %30
%35 = fsub double 1.000000e+00, %7
%36 = fptosi double %35 to i32
%37 = tail call double @log(i32 %36)
%38 = fneg double %37
%39 = tail call i32 (...) @Rand()
%40 = tail call double @log(i32 %39)
%41 = fneg double %40
%42 = fdiv double %41, %31
%43 = fcmp ugt double %42, %38
br i1 %43, label %56, label %44
44: ; preds = %34, %49
%45 = phi i32 [ %47, %49 ], [ %6, %34 ]
%46 = phi double [ %54, %49 ], [ %42, %34 ]
%47 = add nsw i32 %45, -1
%48 = icmp eq i32 %47, 0
br i1 %48, label %56, label %49
49: ; preds = %44
%50 = tail call i32 (...) @Rand()
%51 = tail call double @log(i32 %50)
%52 = sitofp i32 %47 to double
%53 = fdiv double %51, %52
%54 = fsub double %46, %53
%55 = fcmp ugt double %54, %38
br i1 %55, label %56, label %44
56: ; preds = %49, %44, %34
%57 = phi i32 [ %6, %34 ], [ 0, %44 ], [ %47, %49 ]
%58 = sub nsw i32 %6, %57
br label %84
59: ; preds = %30
%60 = add nsw i32 %6, 1
%61 = sitofp i32 %60 to double
%62 = fmul double %7, %61
%63 = fptosi double %62 to i32
%64 = sub nsw i32 %60, %63
%65 = tail call double @BetaRand(i32 %63, i32 %64)
%66 = fcmp ugt double %65, %7
br i1 %66, label %77, label %67
67: ; preds = %59
%68 = fsub double %7, %65
%69 = fsub double 1.000000e+00, %65
%70 = fdiv double %68, %69
%71 = sub nsw i32 %6, %63
%72 = add nsw i32 %8, %63
%73 = tail call i64 @isnan(double %70)
%74 = icmp eq i64 %73, 0
br i1 %74, label %5, label %84
75: ; preds = %84, %77
%76 = phi i32 [ %83, %77 ], [ %87, %84 ]
ret i32 %76
77: ; preds = %59
%78 = add nsw i32 %63, -1
%79 = fsub double %65, %7
%80 = fdiv double %79, %65
%81 = tail call i32 @BinoRand(double %80, i32 %78)
%82 = add i32 %78, %8
%83 = sub i32 %82, %81
br label %75
84: ; preds = %67, %5, %11, %20, %2, %18, %56
%85 = phi i32 [ %8, %56 ], [ %8, %18 ], [ 0, %2 ], [ %8, %20 ], [ %72, %67 ], [ %8, %5 ], [ %8, %11 ]
%86 = phi i32 [ %58, %56 ], [ 0, %18 ], [ 0, %2 ], [ %27, %20 ], [ 0, %67 ], [ 0, %5 ], [ %6, %11 ]
%87 = add nsw i32 %86, %85
br label %75
}