#include <stdio.h>
#include <stdlib.h>
#include <math.h>
void swap(int &a, int &b)
{
int c = a;
a = b;
b = c;
}
int gcd(int a, int b)
{
if(a < b) swap(a, b);
if(a % b == 0)
return b;
else return gcd(b, a%b);
}
int main()
{
int n;
int i, j, count;
scanf("%d", &n);
while( n !=0)
{
count = 0;
int * data =(int *) malloc(n * sizeof(int));
for(i = 0; i < n; i++)
{
scanf("%d",&data[i]);
}
for(i = 0; i < n; i++)
for(j = i+1; j < n; j++)
{
if(gcd(data[i], data[j]) == 1) count++;
}
if(count == 0) printf("No estimate for this data set.\n");
else
printf("%.6f\n",sqrt((double)3 * n*(n-1) /count));
scanf("%d", &n);
}
return 0;
}
分享到:
相关推荐
POJ分类POJ分类POJ分类POJ分类POJ分类POJ分类POJ分类POJ分类POJ分类POJ分类POJ分类POJ分类POJ分类POJ分类POJ分类POJ分类POJ分类POJ分类
poj 解题报告poj 解题报告poj 解题报告poj 解题报告poj 解题报告poj 解题报告poj 解题报告poj 解题报告poj 解题报告poj 解题报告poj 解题报告poj 解题报告poj 解题报告poj 解题报告poj 解题报告poj 解题报告poj 解题...
POJ第1861题源码 POJ第1861题源码 POJ第1861题源码
北大POJ1159-Palindrome 解题报告+AC代码
C语言 poj npu 西工大 C语言Poj答案全完整打包,给有需要的朋友
poj 3414解题报告poj 3414解题报告poj 3414解题报告poj 3414解题报告
poj分类poj分类poj分类poj分类
poj 1012解题报告poj 1012解题报告poj 1012解题报告poj 1012解题报告
poj 2329解题报告poj 2329解题报告poj 2329解题报告poj 2329解题报告
北大POJ2002-Squares 解题报告+AC代码
POJ1503解答 POJ1503解答,正确答案(已通过POJ)
poj 1659解题报告poj 1659解题报告poj 1659解题报告poj 1659解题报告
POJ1048,加强版的约瑟夫问题 难度中等
POJ1083的代码,POJ1083的代码,POJ1083的代码
poj 百练 题目分类 poj 百练 题目分类
POJ上的一道题目,自己写的代码,因为想下载别人的, 所以就放上了。
poj 1001答案
POJ2968代码有用,欢迎下载,POJ代码
Poj中一些题目的源代码,里面共有二十多道题目,OI
北大POJ3253-POJ3253-Fence Repair【STL优先队列】 解题报告+AC代码