#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 binom(int, int);
double proba_sing(int, int, int, int);
double proba_bin(int *neighb_counts, int *weights, int *nodes, int card,
                 int rev_count) {
  double res;
  double denom = 0;
  double num = 0;
  int n, C;
  for (int i = 0; i <= card; i++) {
    res = 1;
    for (int j = 0; j < rev_count; j++) {
      n = nodes[j];
      C = neighb_counts[n];
      res *= proba_sing(C - weights[n], card - i, weights[n], i);
      if (res == 0)
        break;
    }
    if (res != 0) {
      res *= binom(i, card);
      num += i * res;
      denom += res;
    }
  }
  return num / denom / card;
}


x86 assembly (clang, -O0):
<proba_bin>:
push   %rbp
mov    %rsp,%rbp
sub    $0x50,%rsp
mov    %rdi,-0x8(%rbp)
mov    %rsi,-0x10(%rbp)
mov    %rdx,-0x18(%rbp)
mov    %ecx,-0x1c(%rbp)
mov    %r8d,-0x20(%rbp)
xorps  %xmm0,%xmm0
movsd  %xmm0,-0x30(%rbp)
xorps  %xmm0,%xmm0
movsd  %xmm0,-0x38(%rbp)
movl   $0x0,-0x44(%rbp)
mov    -0x44(%rbp),%eax
cmp    -0x1c(%rbp),%eax
jg     126b <proba_bin+0x14b>
movsd  0xe9a(%rip),%xmm0
movsd  %xmm0,-0x28(%rbp)
movl   $0x0,-0x48(%rbp)
mov    -0x48(%rbp),%eax
cmp    -0x20(%rbp),%eax
jge    11fb <proba_bin+0xdb>
mov    -0x18(%rbp),%rax
movslq -0x48(%rbp),%rcx
mov    (%rax,%rcx,4),%eax
mov    %eax,-0x3c(%rbp)
mov    -0x8(%rbp),%rax
movslq -0x3c(%rbp),%rcx
mov    (%rax,%rcx,4),%eax
mov    %eax,-0x40(%rbp)
mov    -0x40(%rbp),%edi
mov    -0x10(%rbp),%rax
movslq -0x3c(%rbp),%rcx
sub    (%rax,%rcx,4),%edi
mov    -0x1c(%rbp),%esi
sub    -0x44(%rbp),%esi
mov    -0x10(%rbp),%rax
movslq -0x3c(%rbp),%rcx
mov    (%rax,%rcx,4),%edx
mov    -0x44(%rbp),%ecx
call   1040 <proba_sing@plt>
mulsd  -0x28(%rbp),%xmm0
movsd  %xmm0,-0x28(%rbp)
movsd  -0x28(%rbp),%xmm0
xorps  %xmm1,%xmm1
ucomisd %xmm1,%xmm0
jne    11e8 <proba_bin+0xc8>
jp     11e8 <proba_bin+0xc8>
jmp    11fb <proba_bin+0xdb>
jmp    11ed <proba_bin+0xcd>
mov    -0x48(%rbp),%eax
add    $0x1,%eax
mov    %eax,-0x48(%rbp)
jmp    1172 <proba_bin+0x52>
movsd  -0x28(%rbp),%xmm0
xorps  %xmm1,%xmm1
ucomisd %xmm1,%xmm0
jne    1218 <proba_bin+0xf8>
jp     1218 <proba_bin+0xf8>
jmp    1258 <proba_bin+0x138>
mov    -0x44(%rbp),%edi
mov    -0x1c(%rbp),%esi
call   1030 <binom@plt>
mulsd  -0x28(%rbp),%xmm0
movsd  %xmm0,-0x28(%rbp)
cvtsi2sdl -0x44(%rbp),%xmm0
movsd  -0x28(%rbp),%xmm2
movsd  -0x38(%rbp),%xmm1
mulsd  %xmm2,%xmm0
addsd  %xmm1,%xmm0
movsd  %xmm0,-0x38(%rbp)
movsd  -0x28(%rbp),%xmm0
addsd  -0x30(%rbp),%xmm0
movsd  %xmm0,-0x30(%rbp)
jmp    125d <proba_bin+0x13d>
mov    -0x44(%rbp),%eax
add    $0x1,%eax
mov    %eax,-0x44(%rbp)
jmp    1152 <proba_bin+0x32>
movsd  -0x38(%rbp),%xmm0
divsd  -0x30(%rbp),%xmm0
cvtsi2sdl -0x1c(%rbp),%xmm1
divsd  %xmm1,%xmm0
add    $0x50,%rsp
pop    %rbp
ret

x86 assembly (clang, -O2):
<proba_bin>:
push   %rbp
push   %r15
push   %r14
push   %r13
push   %r12
push   %rbx
sub    $0x28,%rsp
mov    %ecx,%r15d
mov    %rdx,0x20(%rsp)
test   %ecx,%ecx
js     1254 <proba_bin+0x134>
mov    %r8d,%r12d
mov    %rsi,%r13
mov    %rdi,%rbp
movslq %r8d,%r14
xor    %ebx,%ebx
xorpd  %xmm0,%xmm0
movsd  %xmm0,(%rsp)
xorpd  %xmm1,%xmm1
mov    %r15d,0x14(%rsp)
mov    %r8d,0x10(%rsp)
jmp    11bb <proba_bin+0x9b>
cs nopw 0x0(%rax,%rax,1)
nop
movsd  %xmm0,0x8(%rsp)
mov    %ebx,%edi
mov    %r15d,%esi
call   1030 <binom@plt>
xorps  %xmm1,%xmm1
cvtsi2sd %ebx,%xmm1
mulsd  0x8(%rsp),%xmm0
mulsd  %xmm0,%xmm1
movsd  0x18(%rsp),%xmm2
addsd  %xmm1,%xmm2
movapd %xmm2,%xmm1
movsd  (%rsp),%xmm2
addsd  %xmm0,%xmm2
movsd  %xmm2,(%rsp)
lea    0x1(%rbx),%eax
cmp    %r15d,%ebx
mov    %eax,%ebx
je     124d <proba_bin+0x12d>
movsd  %xmm1,0x18(%rsp)
movsd  0xe3f(%rip),%xmm0
test   %r12d,%r12d
jle    1170 <proba_bin+0x50>
sub    %ebx,%r15d
mov    $0x1,%r12d
movsd  0xe29(%rip),%xmm0
nop
movsd  %xmm0,0x8(%rsp)
mov    0x20(%rsp),%rax
movslq -0x4(%rax,%r12,4),%rax
mov    0x0(%rbp,%rax,4),%edi
mov    0x0(%r13,%rax,4),%edx
sub    %edx,%edi
mov    %r15d,%esi
mov    %ebx,%ecx
call   1040 <proba_sing@plt>
mulsd  0x8(%rsp),%xmm0
ucomisd 0xdfd(%rip),%xmm0
jne    1217 <proba_bin+0xf7>
jnp    1224 <proba_bin+0x104>
lea    0x1(%r12),%rax
cmp    %r14,%r12
mov    %rax,%r12
jl     11e0 <proba_bin+0xc0>
ucomisd 0xde4(%rip),%xmm0
mov    0x14(%rsp),%r15d
mov    0x10(%rsp),%r12d
jne    1170 <proba_bin+0x50>
jp     1170 <proba_bin+0x50>
movsd  0x18(%rsp),%xmm1
jmp    11ad <proba_bin+0x8d>
divsd  (%rsp),%xmm1
jmp    125c <proba_bin+0x13c>
movsd  0xda4(%rip),%xmm1
xorps  %xmm0,%xmm0
cvtsi2sd %r15d,%xmm0
divsd  %xmm0,%xmm1
movapd %xmm1,%xmm0
add    $0x28,%rsp
pop    %rbx
pop    %r12
pop    %r13
pop    %r14
pop    %r15
pop    %rbp
ret

LLVM IR (-O0):
define double @proba_bin(i32* %0, i32* %1, i32* %2, i32 %3, i32 %4) {
%6 = alloca i32*
%7 = alloca i32*
%8 = alloca i32*
%9 = alloca i32
%10 = alloca i32
%11 = alloca double
%12 = alloca double
%13 = alloca double
%14 = alloca i32
%15 = alloca i32
%16 = alloca i32
%17 = alloca i32
store i32* %0, i32** %6
store i32* %1, i32** %7
store i32* %2, i32** %8
store i32 %3, i32* %9
store i32 %4, i32* %10
store double 0.000000e+00, double* %12
store double 0.000000e+00, double* %13
store i32 0, i32* %16
br label %18
18: ; preds = %82, %5
%19 = load i32, i32* %16
%20 = load i32, i32* %9
%21 = icmp sle i32 %19, %20
br i1 %21, label %22, label %85
22: ; preds = %18
store double 1.000000e+00, double* %11
store i32 0, i32* %17
br label %23
23: ; preds = %61, %22
%24 = load i32, i32* %17
%25 = load i32, i32* %10
%26 = icmp slt i32 %24, %25
br i1 %26, label %27, label %64
27: ; preds = %23
%28 = load i32*, i32** %8
%29 = load i32, i32* %17
%30 = sext i32 %29 to i64
%31 = getelementptr inbounds i32, i32* %28, i64 %30
%32 = load i32, i32* %31
store i32 %32, i32* %14
%33 = load i32*, i32** %6
%34 = load i32, i32* %14
%35 = sext i32 %34 to i64
%36 = getelementptr inbounds i32, i32* %33, i64 %35
%37 = load i32, i32* %36
store i32 %37, i32* %15
%38 = load i32, i32* %15
%39 = load i32*, i32** %7
%40 = load i32, i32* %14
%41 = sext i32 %40 to i64
%42 = getelementptr inbounds i32, i32* %39, i64 %41
%43 = load i32, i32* %42
%44 = sub nsw i32 %38, %43
%45 = load i32, i32* %9
%46 = load i32, i32* %16
%47 = sub nsw i32 %45, %46
%48 = load i32*, i32** %7
%49 = load i32, i32* %14
%50 = sext i32 %49 to i64
%51 = getelementptr inbounds i32, i32* %48, i64 %50
%52 = load i32, i32* %51
%53 = load i32, i32* %16
%54 = call double @proba_sing(i32 %44, i32 %47, i32 %52, i32 %53)
%55 = load double, double* %11
%56 = fmul double %55, %54
store double %56, double* %11
%57 = load double, double* %11
%58 = fcmp oeq double %57, 0.000000e+00
br i1 %58, label %59, label %60
59: ; preds = %27
br label %64
60: ; preds = %27
br label %61
61: ; preds = %60
%62 = load i32, i32* %17
%63 = add nsw i32 %62, 1
store i32 %63, i32* %17
br label %23
64: ; preds = %59, %23
%65 = load double, double* %11
%66 = fcmp une double %65, 0.000000e+00
br i1 %66, label %67, label %81
67: ; preds = %64
%68 = load i32, i32* %16
%69 = load i32, i32* %9
%70 = call double @binom(i32 %68, i32 %69)
%71 = load double, double* %11
%72 = fmul double %71, %70
store double %72, double* %11
%73 = load i32, i32* %16
%74 = sitofp i32 %73 to double
%75 = load double, double* %11
%76 = load double, double* %13
%77 = call double @llvm.fmuladd.f64(double %74, double %75, double %76)
store double %77, double* %13
%78 = load double, double* %11
%79 = load double, double* %12
%80 = fadd double %79, %78
store double %80, double* %12
br label %81
81: ; preds = %67, %64
br label %82
82: ; preds = %81
%83 = load i32, i32* %16
%84 = add nsw i32 %83, 1
store i32 %84, i32* %16
br label %18
85: ; preds = %18
%86 = load double, double* %13
%87 = load double, double* %12
%88 = fdiv double %86, %87
%89 = load i32, i32* %9
%90 = sitofp i32 %89 to double
%91 = fdiv double %88, %90
ret double %91
}

LLVM IR (-O2):
define double @proba_bin(i32* nocapture %0, i32* nocapture %1, i32* nocapture %2, i32 %3, i32 %4) {
%6 = icmp slt i32 %3, 0
br i1 %6, label %18, label %7
7: ; preds = %5
%8 = icmp sgt i32 %4, 0
%9 = sext i32 %4 to i64
br label %10
10: ; preds = %7, %48
%11 = phi i32 [ 0, %7 ], [ %51, %48 ]
%12 = phi double [ 0.000000e+00, %7 ], [ %50, %48 ]
%13 = phi double [ 0.000000e+00, %7 ], [ %49, %48 ]
br i1 %8, label %14, label %41
14: ; preds = %10
%15 = sub nsw i32 %3, %11
br label %22
16: ; preds = %48
%17 = fdiv double %49, %50
br label %18
18: ; preds = %16, %5
%19 = phi double [ %17, %16 ], [ 0x7FF8000000000000, %5 ]
%20 = sitofp i32 %3 to double
%21 = fdiv double %19, %20
ret double %21
22: ; preds = %22, %14
%23 = phi i64 [ %36, %22 ], [ 0, %14 ]
%24 = phi double [ %34, %22 ], [ 1.000000e+00, %14 ]
%25 = getelementptr inbounds i32, i32* %2, i64 %23
%26 = load i32, i32* %25
%27 = sext i32 %26 to i64
%28 = getelementptr inbounds i32, i32* %0, i64 %27
%29 = load i32, i32* %28
%30 = getelementptr inbounds i32, i32* %1, i64 %27
%31 = load i32, i32* %30
%32 = sub nsw i32 %29, %31
%33 = tail call double @proba_sing(i32 %32, i32 %15, i32 %31, i32 %11)
%34 = fmul double %24, %33
%35 = fcmp une double %34, 0.000000e+00
%36 = add nuw nsw i64 %23, 1
%37 = icmp slt i64 %36, %9
%38 = select i1 %35, i1 %37, i1 false
br i1 %38, label %22, label %39
39: ; preds = %22
%40 = fcmp une double %34, 0.000000e+00
br i1 %40, label %41, label %48
41: ; preds = %10, %39
%42 = phi double [ %34, %39 ], [ 1.000000e+00, %10 ]
%43 = tail call double @binom(i32 %11, i32 %3)
%44 = fmul double %42, %43
%45 = sitofp i32 %11 to double
%46 = tail call double @llvm.fmuladd.f64(double %45, double %44, double %13)
%47 = fadd double %12, %44
br label %48
48: ; preds = %39, %41
%49 = phi double [ %46, %41 ], [ %13, %39 ]
%50 = phi double [ %47, %41 ], [ %12, %39 ]
%51 = add nuw i32 %11, 1
%52 = icmp eq i32 %11, %3
br i1 %52, label %16, label %10
}