#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 exp(double);
void pdfluria(double m, double phi, int k, double *problist) {
  double cumuphi;
  int n, j;
  for (n = 1; n <= k; n++)
    problist[n] = 0;
  problist[0] = exp(-m);
  for (n = 1; n <= k; n++) {
    cumuphi = 1.0;
    for (j = 1; j <= n; j++) {
      problist[n] =
          problist[n] + cumuphi * (1 - j * phi / (j + 1)) * problist[n - j];
      cumuphi = cumuphi * phi;
    }
    problist[n] = problist[n] * m / n;
  };
}


x86 assembly (clang, -O0):
<pdfluria>:
push   %rbp
mov    %rsp,%rbp
sub    $0x30,%rsp
movsd  %xmm0,-0x8(%rbp)
movsd  %xmm1,-0x10(%rbp)
mov    %edi,-0x14(%rbp)
mov    %rsi,-0x20(%rbp)
movl   $0x1,-0x2c(%rbp)
mov    -0x2c(%rbp),%eax
cmp    -0x14(%rbp),%eax
jg     115a <pdfluria+0x4a>
mov    -0x20(%rbp),%rax
movslq -0x2c(%rbp),%rcx
xorps  %xmm0,%xmm0
movsd  %xmm0,(%rax,%rcx,8)
mov    -0x2c(%rbp),%eax
add    $0x1,%eax
mov    %eax,-0x2c(%rbp)
jmp    1130 <pdfluria+0x20>
movsd  -0x8(%rbp),%xmm0
movq   %xmm0,%rax
movabs $0x8000000000000000,%rcx
xor    %rcx,%rax
movq   %rax,%xmm0
call   1030 <exp@plt>
mov    -0x20(%rbp),%rax
movsd  %xmm0,(%rax)
movl   $0x1,-0x2c(%rbp)
mov    -0x2c(%rbp),%eax
cmp    -0x14(%rbp),%eax
jg     126a <pdfluria+0x15a>
movsd  0xe62(%rip),%xmm0
movsd  %xmm0,-0x28(%rbp)
movl   $0x1,-0x30(%rbp)
mov    -0x30(%rbp),%eax
cmp    -0x2c(%rbp),%eax
jg     1234 <pdfluria+0x124>
mov    -0x20(%rbp),%rax
movslq -0x2c(%rbp),%rcx
movsd  (%rax,%rcx,8),%xmm1
movsd  -0x28(%rbp),%xmm0
cvtsi2sdl -0x30(%rbp),%xmm3
mulsd  -0x10(%rbp),%xmm3
mov    -0x30(%rbp),%eax
add    $0x1,%eax
cvtsi2sd %eax,%xmm2
divsd  %xmm2,%xmm3
movsd  0xe18(%rip),%xmm2
subsd  %xmm3,%xmm2
mulsd  %xmm2,%xmm0
mov    -0x20(%rbp),%rax
mov    -0x2c(%rbp),%ecx
sub    -0x30(%rbp),%ecx
movslq %ecx,%rcx
movsd  (%rax,%rcx,8),%xmm2
mulsd  %xmm2,%xmm0
addsd  %xmm1,%xmm0
mov    -0x20(%rbp),%rax
movslq -0x2c(%rbp),%rcx
movsd  %xmm0,(%rax,%rcx,8)
movsd  -0x28(%rbp),%xmm0
mulsd  -0x10(%rbp),%xmm0
movsd  %xmm0,-0x28(%rbp)
mov    -0x30(%rbp),%eax
add    $0x1,%eax
mov    %eax,-0x30(%rbp)
jmp    11aa <pdfluria+0x9a>
mov    -0x20(%rbp),%rax
movslq -0x2c(%rbp),%rcx
movsd  (%rax,%rcx,8),%xmm0
mulsd  -0x8(%rbp),%xmm0
cvtsi2sdl -0x2c(%rbp),%xmm1
divsd  %xmm1,%xmm0
mov    -0x20(%rbp),%rax
movslq -0x2c(%rbp),%rcx
movsd  %xmm0,(%rax,%rcx,8)
mov    -0x2c(%rbp),%eax
add    $0x1,%eax
mov    %eax,-0x2c(%rbp)
jmp    118a <pdfluria+0x7a>
add    $0x30,%rsp
pop    %rbp
ret

x86 assembly (clang, -O2):
<pdfluria>:
push   %r15
push   %r14
push   %rbx
sub    $0x20,%rsp
mov    %rsi,%rbx
test   %edi,%edi
jle    1210 <pdfluria+0xf0>
mov    %edi,%r14d
lea    0x8(%rbx),%rdi
mov    %r14d,%edx
shl    $0x3,%rdx
xor    %r15d,%r15d
xor    %esi,%esi
movaps %xmm0,0x10(%rsp)
movsd  %xmm1,0x8(%rsp)
call   1030 <memset@plt>
movaps 0xea2(%rip),%xmm0
xorps  0x10(%rsp),%xmm0
call   1040 <exp@plt>
movsd  0x8(%rsp),%xmm6
movapd 0x10(%rsp),%xmm5
movsd  %xmm0,(%rbx)
add    $0x1,%r14d
mov    $0x1,%eax
movsd  0xe87(%rip),%xmm0
nopl   0x0(%rax)
movsd  (%rbx,%rax,8),%xmm1
mov    $0x1,%ecx
mov    %r15,%rdx
movapd %xmm0,%xmm2
cs nopw 0x0(%rax,%rax,1)
nopl   0x0(%rax,%rax,1)
xorps  %xmm3,%xmm3
cvtsi2sd %ecx,%xmm3
mulsd  %xmm6,%xmm3
add    $0x1,%ecx
xorps  %xmm4,%xmm4
cvtsi2sd %ecx,%xmm4
divsd  %xmm4,%xmm3
movapd %xmm0,%xmm4
subsd  %xmm3,%xmm4
mulsd  %xmm2,%xmm4
mulsd  (%rbx,%rdx,8),%xmm4
addsd  %xmm4,%xmm1
movsd  %xmm1,(%rbx,%rax,8)
mulsd  %xmm6,%xmm2
add    $0xffffffffffffffff,%rdx
jb     11b0 <pdfluria+0x90>
mulsd  %xmm5,%xmm1
xorps  %xmm2,%xmm2
cvtsi2sd %eax,%xmm2
divsd  %xmm2,%xmm1
movsd  %xmm1,(%rbx,%rax,8)
add    $0x1,%rax
add    $0x1,%r15
cmp    %r14,%rax
jne    1190 <pdfluria+0x70>
jmp    1221 <pdfluria+0x101>
xorpd  0xde8(%rip),%xmm0
call   1040 <exp@plt>
movsd  %xmm0,(%rbx)
add    $0x20,%rsp
pop    %rbx
pop    %r14
pop    %r15
ret

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

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