#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 */
scalar_t__ EOF;
int /*<<< orphan*/ free(int *);
int gcd(int, int);
int *malloc(int);
int /*<<< orphan*/ printf(char *, ...);
scalar_t__ scanf(char *, int *);
double sqrt(int);
int main() {
  int n, i, j, coprime_count;
  int *ary;
  while (scanf("%d", &n) != EOF) {
    if (n == 0)
      break;
    ary = malloc(n * sizeof(int));
    for (i = 0; i < n; i++)
      scanf("%d", ary + i);
    coprime_count = 0;
    for (i = 0; i < n - 1; i++)
      for (j = i + 1; j < n; j++)
        if (gcd(ary[i], ary[j]) == 1)
          coprime_count++;
    if (coprime_count == 0)
      printf("No estimate for this data set.\n");
    else
      printf("%f\n", sqrt(6 * (n * (n - 1) / 2) / (float)coprime_count));
    free(ary);
  }
  return 0;
}


x86 assembly (clang, -O0):
<main>:
push   %rbp
mov    %rsp,%rbp
sub    $0x20,%rsp
movl   $0x0,-0x4(%rbp)
lea    0xe8a(%rip),%rdi
lea    -0x8(%rbp),%rsi
call   1060 <scanf@plt>
mov    0x2e62(%rip),%rcx
cmp    (%rcx),%rax
je     12da <main+0x17a>
cmpl   $0x0,-0x8(%rbp)
jne    119e <main+0x3e>
jmp    12da <main+0x17a>
movslq -0x8(%rbp),%rax
shl    $0x2,%rax
mov    %eax,%edi
call   1050 <malloc@plt>
mov    %rax,-0x20(%rbp)
movl   $0x0,-0xc(%rbp)
mov    -0xc(%rbp),%eax
cmp    -0x8(%rbp),%eax
jge    11ed <main+0x8d>
mov    -0x20(%rbp),%rsi
movslq -0xc(%rbp),%rax
shl    $0x2,%rax
add    %rax,%rsi
lea    0xe26(%rip),%rdi
call   1060 <scanf@plt>
mov    -0xc(%rbp),%eax
add    $0x1,%eax
mov    %eax,-0xc(%rbp)
jmp    11b8 <main+0x58>
movl   $0x0,-0x14(%rbp)
movl   $0x0,-0xc(%rbp)
mov    -0xc(%rbp),%eax
mov    -0x8(%rbp),%ecx
sub    $0x1,%ecx
cmp    %ecx,%eax
jge    1274 <main+0x114>
mov    -0xc(%rbp),%eax
add    $0x1,%eax
mov    %eax,-0x10(%rbp)
mov    -0x10(%rbp),%eax
cmp    -0x8(%rbp),%eax
jge    1261 <main+0x101>
mov    -0x20(%rbp),%rax
movslq -0xc(%rbp),%rcx
mov    (%rax,%rcx,4),%edi
mov    -0x20(%rbp),%rax
movslq -0x10(%rbp),%rcx
mov    (%rax,%rcx,4),%esi
call   1070 <gcd@plt>
cmp    $0x1,%eax
jne    124e <main+0xee>
mov    -0x14(%rbp),%eax
add    $0x1,%eax
mov    %eax,-0x14(%rbp)
jmp    1253 <main+0xf3>
mov    -0x10(%rbp),%eax
add    $0x1,%eax
mov    %eax,-0x10(%rbp)
jmp    1215 <main+0xb5>
jmp    1266 <main+0x106>
mov    -0xc(%rbp),%eax
add    $0x1,%eax
mov    %eax,-0xc(%rbp)
jmp    11fb <main+0x9b>
cmpl   $0x0,-0x14(%rbp)
jne    1291 <main+0x131>
lea    0xd7e(%rip),%rdi
mov    $0x0,%al
call   1040 <printf@plt>
jmp    12cc <main+0x16c>
mov    -0x8(%rbp),%eax
mov    -0x8(%rbp),%ecx
sub    $0x1,%ecx
imul   %ecx,%eax
mov    $0x2,%ecx
cltd
idiv   %ecx
imul   $0x6,%eax,%eax
cvtsi2ss %eax,%xmm0
cvtsi2ssl -0x14(%rbp),%xmm1
divss  %xmm1,%xmm0
cvttss2si %xmm0,%edi
call   1080 <sqrt@plt>
lea    0xd5e(%rip),%rdi
mov    $0x1,%al
call   1040 <printf@plt>
mov    -0x20(%rbp),%rdi
call   1030 <free@plt>
jmp    116f <main+0xf>
xor    %eax,%eax
add    $0x20,%rsp
pop    %rbp
ret

x86 assembly (clang, -O2):
<main>:
push   %rbp
push   %r15
push   %r14
push   %r13
push   %r12
push   %rbx
push   %rax
lea    0xe7e(%rip),%rdi
lea    0x4(%rsp),%rsi
call   1070 <scanf@plt>
mov    0x2e55(%rip),%r15
cmp    (%r15),%rax
je     130e <main+0x19e>
mov    0x4(%rsp),%ebp
test   %ebp,%ebp
je     130e <main+0x19e>
lea    0xe51(%rip),%r12
lea    0xe51(%rip),%r13
cs nopw 0x0(%rax,%rax,1)
lea    0x0(,%rbp,4),%edi
call   1060 <malloc@plt>
mov    %rax,%r14
test   %ebp,%ebp
jle    12e0 <main+0x170>
mov    %r14,%rbp
xor    %ebx,%ebx
nopl   0x0(%rax)
mov    %r12,%rdi
mov    %rbp,%rsi
call   1070 <scanf@plt>
add    $0x1,%rbx
movslq 0x4(%rsp),%rax
add    $0x4,%rbp
cmp    %rax,%rbx
jl     11e0 <main+0x70>
cmp    $0x2,%eax
jl     12e0 <main+0x170>
xor    %ebx,%ebx
mov    $0x1,%r13d
xor    %ebp,%ebp
jmp    122f <main+0xbf>
cs nopw 0x0(%rax,%rax,1)
nopl   0x0(%rax)
lea    -0x1(%rax),%ecx
movslq %ecx,%rdx
add    $0x1,%r13
cmp    %rdx,%rbx
jge    1280 <main+0x110>
mov    %rbx,%r12
add    $0x1,%rbx
movslq %eax,%rcx
cmp    %rcx,%rbx
jge    1220 <main+0xb0>
mov    %r13,%r15
cs nopw 0x0(%rax,%rax,1)
nopl   0x0(%rax,%rax,1)
mov    (%r14,%r12,4),%edi
mov    (%r14,%r15,4),%esi
call   1080 <gcd@plt>
xor    %ecx,%ecx
cmp    $0x1,%eax
sete   %cl
add    %ecx,%ebp
add    $0x1,%r15
mov    0x4(%rsp),%eax
cmp    %r15d,%eax
jg     1250 <main+0xe0>
jmp    1220 <main+0xb0>
cs nopw 0x0(%rax,%rax,1)
test   %ebp,%ebp
mov    0x2d5f(%rip),%r15
lea    0xd70(%rip),%r12
lea    0xd70(%rip),%r13
je     12e0 <main+0x170>
imul   %eax,%ecx
mov    %ecx,%eax
shr    $0x1f,%eax
add    %ecx,%eax
and    $0xfffffffe,%eax
lea    (%rax,%rax,2),%eax
xorps  %xmm0,%xmm0
cvtsi2ss %eax,%xmm0
xorps  %xmm1,%xmm1
cvtsi2ss %ebp,%xmm1
divss  %xmm1,%xmm0
cvttss2si %xmm0,%edi
call   1090 <sqrt@plt>
lea    0xd38(%rip),%rdi
mov    $0x1,%al
call   1050 <printf@plt>
jmp    12e8 <main+0x178>
cs nopw 0x0(%rax,%rax,1)
xchg   %ax,%ax
mov    %r13,%rdi
call   1040 <puts@plt>
mov    %r14,%rdi
call   1030 <free@plt>
mov    %r12,%rdi
lea    0x4(%rsp),%rsi
call   1070 <scanf@plt>
cmp    (%r15),%rax
je     130e <main+0x19e>
mov    0x4(%rsp),%ebp
test   %ebp,%ebp
jne    11c0 <main+0x50>
xor    %eax,%eax
add    $0x8,%rsp
pop    %rbx
pop    %r12
pop    %r13
pop    %r14
pop    %r15
pop    %rbp
ret

LLVM IR (-O0):
@.str = [3 x i8] c"%d\00"
@EOF = global i64 0
@.str.1 = [32 x i8] c"No estimate for this data set.\0A\00"
@.str.2 = [4 x i8] c"%f\0A\00"
define i32 @main() {
%1 = alloca i32
%2 = alloca i32
%3 = alloca i32
%4 = alloca i32
%5 = alloca i32
%6 = alloca i32*
store i32 0, i32* %1
br label %7
7: ; preds = %90, %0
%8 = call i64 @scanf(i8* getelementptr inbounds ([3 x i8], [3 x i8]* @.str, i64 0, i64 0), i32* %2)
%9 = load i64, i64* @EOF
%10 = icmp ne i64 %8, %9
br i1 %10, label %11, label %93
11: ; preds = %7
%12 = load i32, i32* %2
%13 = icmp eq i32 %12, 0
br i1 %13, label %14, label %15
14: ; preds = %11
br label %93
15: ; preds = %11
%16 = load i32, i32* %2
%17 = sext i32 %16 to i64
%18 = mul i64 %17, 4
%19 = trunc i64 %18 to i32
%20 = call i32* @malloc(i32 %19)
store i32* %20, i32** %6
store i32 0, i32* %3
br label %21
21: ; preds = %31, %15
%22 = load i32, i32* %3
%23 = load i32, i32* %2
%24 = icmp slt i32 %22, %23
br i1 %24, label %25, label %34
25: ; preds = %21
%26 = load i32*, i32** %6
%27 = load i32, i32* %3
%28 = sext i32 %27 to i64
%29 = getelementptr inbounds i32, i32* %26, i64 %28
%30 = call i64 @scanf(i8* getelementptr inbounds ([3 x i8], [3 x i8]* @.str, i64 0, i64 0), i32* %29)
br label %31
31: ; preds = %25
%32 = load i32, i32* %3
%33 = add nsw i32 %32, 1
store i32 %33, i32* %3
br label %21
34: ; preds = %21
store i32 0, i32* %5
store i32 0, i32* %3
br label %35
35: ; preds = %68, %34
%36 = load i32, i32* %3
%37 = load i32, i32* %2
%38 = sub nsw i32 %37, 1
%39 = icmp slt i32 %36, %38
br i1 %39, label %40, label %71
40: ; preds = %35
%41 = load i32, i32* %3
%42 = add nsw i32 %41, 1
store i32 %42, i32* %4
br label %43
43: ; preds = %64, %40
%44 = load i32, i32* %4
%45 = load i32, i32* %2
%46 = icmp slt i32 %44, %45
br i1 %46, label %47, label %67
47: ; preds = %43
%48 = load i32*, i32** %6
%49 = load i32, i32* %3
%50 = sext i32 %49 to i64
%51 = getelementptr inbounds i32, i32* %48, i64 %50
%52 = load i32, i32* %51
%53 = load i32*, i32** %6
%54 = load i32, i32* %4
%55 = sext i32 %54 to i64
%56 = getelementptr inbounds i32, i32* %53, i64 %55
%57 = load i32, i32* %56
%58 = call i32 @gcd(i32 %52, i32 %57)
%59 = icmp eq i32 %58, 1
br i1 %59, label %60, label %63
60: ; preds = %47
%61 = load i32, i32* %5
%62 = add nsw i32 %61, 1
store i32 %62, i32* %5
br label %63
63: ; preds = %60, %47
br label %64
64: ; preds = %63
%65 = load i32, i32* %4
%66 = add nsw i32 %65, 1
store i32 %66, i32* %4
br label %43
67: ; preds = %43
br label %68
68: ; preds = %67
%69 = load i32, i32* %3
%70 = add nsw i32 %69, 1
store i32 %70, i32* %3
br label %35
71: ; preds = %35
%72 = load i32, i32* %5
%73 = icmp eq i32 %72, 0
br i1 %73, label %74, label %76
74: ; preds = %71
%75 = call i32 (i8*, ...) @printf(i8* getelementptr inbounds ([32 x i8], [32 x i8]* @.str.1, i64 0, i64 0))
br label %90
76: ; preds = %71
%77 = load i32, i32* %2
%78 = load i32, i32* %2
%79 = sub nsw i32 %78, 1
%80 = mul nsw i32 %77, %79
%81 = sdiv i32 %80, 2
%82 = mul nsw i32 6, %81
%83 = sitofp i32 %82 to float
%84 = load i32, i32* %5
%85 = sitofp i32 %84 to float
%86 = fdiv float %83, %85
%87 = fptosi float %86 to i32
%88 = call double @sqrt(i32 %87)
%89 = call i32 (i8*, ...) @printf(i8* getelementptr inbounds ([4 x i8], [4 x i8]* @.str.2, i64 0, i64 0), double %88)
br label %90
90: ; preds = %76, %74
%91 = load i32*, i32** %6
%92 = call i32 @free(i32* %91)
br label %7
93: ; preds = %14, %7
ret i32 0
}

LLVM IR (-O2):
@.str = [3 x i8] c"%d\00"
@EOF = global i64 0
@.str.2 = [4 x i8] c"%f\0A\00"
@str = [31 x i8] c"No estimate for this data set.\00"
define i32 @main() {
%1 = alloca i32
%2 = bitcast i32* %1 to i8*
%3 = call i64 @scanf(i8* getelementptr inbounds ([3 x i8], [3 x i8]* @.str, i64 0, i64 0), i32* %1)
%4 = load i64, i64* @EOF
%5 = icmp eq i64 %3, %4
%6 = load i32, i32* %1
%7 = icmp eq i32 %6, 0
%8 = select i1 %5, i1 true, i1 %7
br i1 %8, label %77, label %9
9: ; preds = %0, %69
%10 = phi i32 [ %74, %69 ], [ %6, %0 ]
%11 = shl i32 %10, 2
%12 = call i32* @malloc(i32 %11)
%13 = icmp sgt i32 %10, 0
br i1 %13, label %16, label %57
14: ; preds = %16
%15 = icmp sgt i32 %21, 1
br i1 %15, label %31, label %57
16: ; preds = %9, %16
%17 = phi i64 [ %20, %16 ], [ 0, %9 ]
%18 = getelementptr inbounds i32, i32* %12, i64 %17
%19 = call i64 @scanf(i8* getelementptr inbounds ([3 x i8], [3 x i8]* @.str, i64 0, i64 0), i32* %18)
%20 = add nuw nsw i64 %17, 1
%21 = load i32, i32* %1
%22 = sext i32 %21 to i64
%23 = icmp slt i64 %20, %22
br i1 %23, label %16, label %14
24: ; preds = %41, %31
%25 = phi i32 [ %32, %31 ], [ %52, %41 ]
%26 = phi i32 [ %35, %31 ], [ %50, %41 ]
%27 = add nsw i32 %25, -1
%28 = sext i32 %27 to i64
%29 = icmp slt i64 %36, %28
%30 = add nuw nsw i64 %34, 1
br i1 %29, label %31, label %55
31: ; preds = %14, %24
%32 = phi i32 [ %25, %24 ], [ %21, %14 ]
%33 = phi i64 [ %36, %24 ], [ 0, %14 ]
%34 = phi i64 [ %30, %24 ], [ 1, %14 ]
%35 = phi i32 [ %26, %24 ], [ 0, %14 ]
%36 = add nuw nsw i64 %33, 1
%37 = sext i32 %32 to i64
%38 = icmp slt i64 %36, %37
br i1 %38, label %39, label %24
39: ; preds = %31
%40 = getelementptr inbounds i32, i32* %12, i64 %33
br label %41
41: ; preds = %39, %41
%42 = phi i64 [ %34, %39 ], [ %51, %41 ]
%43 = phi i32 [ %35, %39 ], [ %50, %41 ]
%44 = load i32, i32* %40
%45 = getelementptr inbounds i32, i32* %12, i64 %42
%46 = load i32, i32* %45
%47 = call i32 @gcd(i32 %44, i32 %46)
%48 = icmp eq i32 %47, 1
%49 = zext i1 %48 to i32
%50 = add nsw i32 %43, %49
%51 = add nuw nsw i64 %42, 1
%52 = load i32, i32* %1
%53 = trunc i64 %51 to i32
%54 = icmp sgt i32 %52, %53
br i1 %54, label %41, label %24
55: ; preds = %24
%56 = icmp eq i32 %26, 0
br i1 %56, label %57, label %59
57: ; preds = %9, %14, %55
%58 = call i32 @puts(i8* dereferenceable(1) getelementptr inbounds ([31 x i8], [31 x i8]* @str, i64 0, i64 0))
br label %69
59: ; preds = %55
%60 = mul nsw i32 %27, %25
%61 = sdiv i32 %60, 2
%62 = mul nsw i32 %61, 6
%63 = sitofp i32 %62 to float
%64 = sitofp i32 %26 to float
%65 = fdiv float %63, %64
%66 = fptosi float %65 to i32
%67 = call double @sqrt(i32 %66)
%68 = call i32 (i8*, ...) @printf(i8* dereferenceable(1) getelementptr inbounds ([4 x i8], [4 x i8]* @.str.2, i64 0, i64 0), double %67)
br label %69
69: ; preds = %59, %57
%70 = call i32 @free(i32* %12)
%71 = call i64 @scanf(i8* getelementptr inbounds ([3 x i8], [3 x i8]* @.str, i64 0, i64 0), i32* %1)
%72 = load i64, i64* @EOF
%73 = icmp eq i64 %71, %72
%74 = load i32, i32* %1
%75 = icmp eq i32 %74, 0
%76 = select i1 %73, i1 true, i1 %75
br i1 %76, label %77, label %9
77: ; preds = %69, %0
ret i32 0
}