double probability(int people) {
  double result = 1.0;
  for (int i = 0; i < people; i++) {
    result *= (365.0 - i) / 365.0;
  }
  return 1.0 - result;
}


x86 assembly (clang, -O0):
<probability>:
push   %rbp
mov    %rsp,%rbp
mov    %edi,-0x4(%rbp)
movsd  0xef1(%rip),%xmm0
movsd  %xmm0,-0x10(%rbp)
movl   $0x0,-0x14(%rbp)
mov    -0x14(%rbp),%eax
cmp    -0x4(%rbp),%eax
jge    115c <probability+0x5c>
cvtsi2sdl -0x14(%rbp),%xmm1
movsd  0xed4(%rip),%xmm0
subsd  %xmm1,%xmm0
movsd  0xec8(%rip),%xmm1
divsd  %xmm1,%xmm0
mulsd  -0x10(%rbp),%xmm0
movsd  %xmm0,-0x10(%rbp)
mov    -0x14(%rbp),%eax
add    $0x1,%eax
mov    %eax,-0x14(%rbp)
jmp    111b <probability+0x1b>
movsd  0xe9c(%rip),%xmm0
subsd  -0x10(%rbp),%xmm0
pop    %rbp
ret

x86 assembly (clang, -O2):
<probability>:
movsd  0xef8(%rip),%xmm0
movapd %xmm0,%xmm1
test   %edi,%edi
jle    11ba <probability+0xba>
cmp    $0x1,%edi
jne    1130 <probability+0x30>
movsd  0xedf(%rip),%xmm1
xorpd  %xmm2,%xmm2
test   $0x1,%dil
jne    11a2 <probability+0xa2>
jmp    11ba <probability+0xba>
mov    %edi,%eax
and    $0xfffffffe,%eax
movsd  0xec3(%rip),%xmm1
xor    %ecx,%ecx
xorpd  %xmm2,%xmm2
movsd  0xebd(%rip),%xmm3
movsd  0xebd(%rip),%xmm4
cs nopw 0x0(%rax,%rax,1)
nopl   (%rax)
movapd %xmm3,%xmm5
subsd  %xmm2,%xmm5
divsd  %xmm3,%xmm5
mulsd  %xmm1,%xmm5
lea    0x1(%rcx),%edx
xorps  %xmm6,%xmm6
cvtsi2sd %edx,%xmm6
movapd %xmm3,%xmm1
subsd  %xmm6,%xmm1
divsd  %xmm3,%xmm1
mulsd  %xmm5,%xmm1
addsd  %xmm4,%xmm2
add    $0x2,%ecx
cmp    %ecx,%eax
jne    1160 <probability+0x60>
xorps  %xmm2,%xmm2
cvtsi2sd %ecx,%xmm2
test   $0x1,%dil
je     11ba <probability+0xba>
movsd  0xe5e(%rip),%xmm3
movapd %xmm3,%xmm4
subsd  %xmm2,%xmm4
divsd  %xmm3,%xmm4
mulsd  %xmm4,%xmm1
subsd  %xmm1,%xmm0
ret

LLVM IR (-O0):
define double @probability(i32 %0) {
%2 = alloca i32
%3 = alloca double
%4 = alloca i32
store i32 %0, i32* %2
store double 1.000000e+00, double* %3
store i32 0, i32* %4
br label %5
5: ; preds = %16, %1
%6 = load i32, i32* %4
%7 = load i32, i32* %2
%8 = icmp slt i32 %6, %7
br i1 %8, label %9, label %19
9: ; preds = %5
%10 = load i32, i32* %4
%11 = sitofp i32 %10 to double
%12 = fsub double 3.650000e+02, %11
%13 = fdiv double %12, 3.650000e+02
%14 = load double, double* %3
%15 = fmul double %14, %13
store double %15, double* %3
br label %16
16: ; preds = %9
%17 = load i32, i32* %4
%18 = add nsw i32 %17, 1
store i32 %18, i32* %4
br label %5
19: ; preds = %5
%20 = load double, double* %3
%21 = fsub double 1.000000e+00, %20
ret double %21
}

LLVM IR (-O2):
define double @probability(i32 %0) {
%2 = icmp sgt i32 %0, 0
br i1 %2, label %3, label %19
3: ; preds = %1
%4 = and i32 %0, 1
%5 = icmp eq i32 %0, 1
br i1 %5, label %10, label %6
6: ; preds = %3
%7 = and i32 %0, -2
br label %22
8: ; preds = %22
%9 = sitofp i32 %35 to double
br label %10
10: ; preds = %8, %3
%11 = phi double [ undef, %3 ], [ %34, %8 ]
%12 = phi double [ 0.000000e+00, %3 ], [ %9, %8 ]
%13 = phi double [ 1.000000e+00, %3 ], [ %34, %8 ]
%14 = icmp eq i32 %4, 0
br i1 %14, label %19, label %15
15: ; preds = %10
%16 = fsub double 3.650000e+02, %12
%17 = fdiv double %16, 3.650000e+02
%18 = fmul double %13, %17
br label %19
19: ; preds = %15, %10, %1
%20 = phi double [ 1.000000e+00, %1 ], [ %11, %10 ], [ %18, %15 ]
%21 = fsub double 1.000000e+00, %20
ret double %21
22: ; preds = %22, %6
%23 = phi i32 [ 0, %6 ], [ %35, %22 ]
%24 = phi double [ 1.000000e+00, %6 ], [ %34, %22 ]
%25 = phi i32 [ 0, %6 ], [ %36, %22 ]
%26 = sitofp i32 %23 to double
%27 = fsub double 3.650000e+02, %26
%28 = fdiv double %27, 3.650000e+02
%29 = fmul double %24, %28
%30 = or i32 %23, 1
%31 = sitofp i32 %30 to double
%32 = fsub double 3.650000e+02, %31
%33 = fdiv double %32, 3.650000e+02
%34 = fmul double %29, %33
%35 = add nuw nsw i32 %23, 2
%36 = add i32 %25, 2
%37 = icmp eq i32 %36, %7
br i1 %37, label %8, label %22
}