int rad(int n, int *primes, int primes_n) {
  int i, r;
  for (i = 0, r = 1; n != 1 && i < primes_n; i++) {
    int p = primes[i], f;
    if (p * p > n) {
      r *= n;
      break;
    }
    for (f = 0; n % p == 0; f++) {
      n /= p;
    }
    if (f != 0) {
      r *= p;
    }
  }
  return r;
}


x86 assembly (clang, -O0):
<rad>:
push   %rbp
mov    %rsp,%rbp
mov    %edi,-0x4(%rbp)
mov    %rsi,-0x10(%rbp)
mov    %edx,-0x14(%rbp)
movl   $0x0,-0x18(%rbp)
movl   $0x1,-0x1c(%rbp)
xor    %eax,%eax
cmpl   $0x1,-0x4(%rbp)
mov    %al,-0x25(%rbp)
je     1137 <rad+0x37>
mov    -0x18(%rbp),%eax
cmp    -0x14(%rbp),%eax
setl   %al
mov    %al,-0x25(%rbp)
mov    -0x25(%rbp),%al
test   $0x1,%al
jne    1147 <rad+0x47>
jmp    11cc <rad+0xcc>
mov    -0x10(%rbp),%rax
movslq -0x18(%rbp),%rcx
mov    (%rax,%rcx,4),%eax
mov    %eax,-0x20(%rbp)
mov    -0x20(%rbp),%eax
imul   -0x20(%rbp),%eax
cmp    -0x4(%rbp),%eax
jle    1174 <rad+0x74>
mov    -0x4(%rbp),%eax
imul   -0x1c(%rbp),%eax
mov    %eax,-0x1c(%rbp)
jmp    11cc <rad+0xcc>
movl   $0x0,-0x24(%rbp)
mov    -0x4(%rbp),%eax
cltd
idivl  -0x20(%rbp)
cmp    $0x0,%edx
jne    11a5 <rad+0xa5>
mov    -0x20(%rbp),%ecx
mov    -0x4(%rbp),%eax
cltd
idiv   %ecx
mov    %eax,-0x4(%rbp)
mov    -0x24(%rbp),%eax
add    $0x1,%eax
mov    %eax,-0x24(%rbp)
jmp    117b <rad+0x7b>
cmpl   $0x0,-0x24(%rbp)
je     11b9 <rad+0xb9>
mov    -0x20(%rbp),%eax
imul   -0x1c(%rbp),%eax
mov    %eax,-0x1c(%rbp)
jmp    11be <rad+0xbe>
mov    -0x18(%rbp),%eax
add    $0x1,%eax
mov    %eax,-0x18(%rbp)
jmp    111c <rad+0x1c>
mov    -0x1c(%rbp),%eax
pop    %rbp
ret

x86 assembly (clang, -O2):
<rad>:
mov    $0x1,%r9d
cmp    $0x1,%edi
je     1160 <rad+0x60>
test   %edx,%edx
jle    1160 <rad+0x60>
mov    %edx,%r8d
mov    $0x1,%r9d
xor    %r10d,%r10d
nopl   0x0(%rax,%rax,1)
mov    (%rsi,%r10,4),%ecx
mov    %ecx,%eax
imul   %ecx,%eax
cmp    %edi,%eax
jg     1164 <rad+0x64>
mov    %edi,%eax
cltd
idiv   %ecx
test   %edx,%edx
je     1140 <rad+0x40>
mov    $0x1,%ecx
jmp    114e <rad+0x4e>
nopl   (%rax)
mov    %edi,%eax
cltd
idiv   %ecx
mov    %eax,%edi
cltd
idiv   %ecx
test   %edx,%edx
je     1140 <rad+0x40>
imul   %ecx,%r9d
cmp    $0x1,%edi
je     1160 <rad+0x60>
add    $0x1,%r10
cmp    %r8,%r10
jb     1120 <rad+0x20>
mov    %r9d,%eax
ret
imul   %edi,%r9d
mov    %r9d,%eax
ret

LLVM IR (-O0):
define i32 @rad(i32 %0, i32* %1, i32 %2) {
%4 = alloca i32
%5 = alloca i32*
%6 = alloca i32
%7 = alloca i32
%8 = alloca i32
%9 = alloca i32
%10 = alloca i32
store i32 %0, i32* %4
store i32* %1, i32** %5
store i32 %2, i32* %6
store i32 0, i32* %7
store i32 1, i32* %8
br label %11
11: ; preds = %56, %3
%12 = load i32, i32* %4
%13 = icmp ne i32 %12, 1
br i1 %13, label %14, label %18
14: ; preds = %11
%15 = load i32, i32* %7
%16 = load i32, i32* %6
%17 = icmp slt i32 %15, %16
br label %18
18: ; preds = %14, %11
%19 = phi i1 [ false, %11 ], [ %17, %14 ]
br i1 %19, label %20, label %59
20: ; preds = %18
%21 = load i32*, i32** %5
%22 = load i32, i32* %7
%23 = sext i32 %22 to i64
%24 = getelementptr inbounds i32, i32* %21, i64 %23
%25 = load i32, i32* %24
store i32 %25, i32* %9
%26 = load i32, i32* %9
%27 = load i32, i32* %9
%28 = mul nsw i32 %26, %27
%29 = load i32, i32* %4
%30 = icmp sgt i32 %28, %29
br i1 %30, label %31, label %35
31: ; preds = %20
%32 = load i32, i32* %4
%33 = load i32, i32* %8
%34 = mul nsw i32 %33, %32
store i32 %34, i32* %8
br label %59
35: ; preds = %20
store i32 0, i32* %10
br label %36
36: ; preds = %45, %35
%37 = load i32, i32* %4
%38 = load i32, i32* %9
%39 = srem i32 %37, %38
%40 = icmp eq i32 %39, 0
br i1 %40, label %41, label %48
41: ; preds = %36
%42 = load i32, i32* %9
%43 = load i32, i32* %4
%44 = sdiv i32 %43, %42
store i32 %44, i32* %4
br label %45
45: ; preds = %41
%46 = load i32, i32* %10
%47 = add nsw i32 %46, 1
store i32 %47, i32* %10
br label %36
48: ; preds = %36
%49 = load i32, i32* %10
%50 = icmp ne i32 %49, 0
br i1 %50, label %51, label %55
51: ; preds = %48
%52 = load i32, i32* %9
%53 = load i32, i32* %8
%54 = mul nsw i32 %53, %52
store i32 %54, i32* %8
br label %55
55: ; preds = %51, %48
br label %56
56: ; preds = %55
%57 = load i32, i32* %7
%58 = add nsw i32 %57, 1
store i32 %58, i32* %7
br label %11
59: ; preds = %31, %18
%60 = load i32, i32* %8
ret i32 %60
}

LLVM IR (-O2):
define i32 @rad(i32 %0, i32* nocapture %1, i32 %2) {
%4 = icmp ne i32 %0, 1
%5 = icmp sgt i32 %2, 0
%6 = and i1 %4, %5
br i1 %6, label %7, label %35
7: ; preds = %3
%8 = zext i32 %2 to i64
br label %9
9: ; preds = %7, %27
%10 = phi i64 [ 0, %7 ], [ %31, %27 ]
%11 = phi i32 [ %0, %7 ], [ %29, %27 ]
%12 = phi i32 [ 1, %7 ], [ %30, %27 ]
%13 = getelementptr inbounds i32, i32* %1, i64 %10
%14 = load i32, i32* %13
%15 = mul nsw i32 %14, %14
%16 = icmp sgt i32 %15, %11
br i1 %16, label %25, label %17
17: ; preds = %9
%18 = srem i32 %11, %14
%19 = icmp eq i32 %18, 0
br i1 %19, label %20, label %27
20: ; preds = %17, %20
%21 = phi i32 [ %22, %20 ], [ %11, %17 ]
%22 = sdiv i32 %21, %14
%23 = srem i32 %22, %14
%24 = icmp eq i32 %23, 0
br i1 %24, label %20, label %27
25: ; preds = %9
%26 = mul nsw i32 %11, %12
br label %35
27: ; preds = %20, %17
%28 = phi i32 [ 1, %17 ], [ %14, %20 ]
%29 = phi i32 [ %11, %17 ], [ %22, %20 ]
%30 = mul nsw i32 %28, %12
%31 = add nuw nsw i64 %10, 1
%32 = icmp ne i32 %29, 1
%33 = icmp ult i64 %31, %8
%34 = select i1 %32, i1 %33, i1 false
br i1 %34, label %9, label %35
35: ; preds = %27, %3, %25
%36 = phi i32 [ %26, %25 ], [ 1, %3 ], [ %30, %27 ]
ret i32 %36
}