求PI

很先进的方法。。。从0~32767个随即数中取出N个数，每2个为一对，共C2n对
找出其中除1以外没其他公约数的对子，m对

6/pi^2=m/C2n

还是暴力搜索。。。用时比较大

#include  < iostream >
#include  < vector >
#include  < string >
#include  < math.h >
#include  < iomanip >
using   namespace  std;


 int  main()
 {
    int n,ran;
    vector<int> num;
    while (1)
    {
        cin>>n;
        if (n==0)
        {break;    }
        num.clear();
        int sum=n*(n-1)/2;
        while (n>0)
        {
            cin>>ran;
            num.push_back(ran);
            n--;
        }
        bool match;
        int count=0;
        for (size_t i=0;i<num.size()-1;i++)
        {
            for (size_t j=i+1;j<num.size();j++)
            {
                match=true;
                for (int k=2;k<=(num[i]>num[j]?num[j]:num[i]);k++)
                {
                    if (num[i]%k==0&&num[j]%k==0)
                    {
                        match=false;
                        break;
                    }
                }
                if (match==true)
                {count++;
                }
            }
        }
        if (count==0)
        {
            cout<<"No estimate for this data set."<<endl;
        }
        else
        {
           double pi=6*sum/(double)count;
           cout<<setprecision(7)<<sqrt(pi)<<endl;
        }
    }
}

PS：貌似也可以用求2数的公约数来做，当不为1时，就不要了
int gys(int a, int b)
{ int r; //用于ab互换
if(a<b){r=a; a=b; b=r;} //如果a比b小就互换，使a比b大
while (b!=0)
{
r=a%b;
a=b;
b=r;
}
/*这个while就是这个程序的核心，它不断的使a和b做除然后求除余数然后最大的数a就没有用了，因为他们的公因数与现在的r和b的公因数相同。而这时，b的值会比r大，为了保证b为小的那个数，所以需要将b的值给a将r的值给b。然后继续做除，直到b=0时表示。已经被整除掉了。这时的a就为他们的最大公约数。*/
return a; //将最大公约数返回。
}