zoj 1337 Pi

题目要求给n个数，从中选出两个数，不考虑次序，则个数为s，这些对中两个数不存在公因子，即最大公约数为1，的个数为num，有

6/PI^2=num/s

即可求出PI。。。

#include<stdio.h>
#include<math.h>
int gcf(int x,int y)           //最大公约数
{
	if(x<y)return gcf(y,x);
	if(x%y==0)return y;
	return gcf(y,x%y);
}
int main()
{
	int n,a[55],i,j,num,s;
	while(scanf("%d",&n)!=EOF,n)
	{
		for(i=0;i<n;i++)scanf("%d",&a[i]);
		s=n%2?(n-1)/2*n:n/2*(n-1);        //c（n，2）
		num=0;
		for(i=0;i<n;i++)
			for(j=i+1;j<n;j++)
				if(gcf(a[i],a[j])==1)num++;
		if(num==0)printf("No estimate for this data set.\n");
		else
			printf("%.6lf\n",sqrt(6.0*s/num));
	}
	return 0;
}