ZOJ 1337 集合、关系求pi

#include<iostream.h>
#include<math.h>
#include<iomanip.h>
//using namespace std;
int shu[100];

 

int main()
{
    int a,b,n,i,j,temp,jishu;
    double pair;
    while((cin>>n) && n!=0)
    {
        jishu=0;
        for(i=0;i<n;i++)
            cin>>shu[i];

        for(i=0;i<n;i++)
        {
            for(j=i+1;j<n;j++)
            {
                a=shu[i];
                b=shu[j];
                do
                {
                    temp=a%b;
                    a=b;
                    b=temp;
                }while(b!=0);
                if(a==1)
                    jishu++;
            }
        }
        pair=n*(n-1)/2;
        pair=sqrt(6*pair/jishu);
        if(jishu==0)
            cout<<"No estimate for this data set."<<endl;
        else
            cout<<setiosflags(ios::fixed)<<pair<<endl;
    }
    return 0;
}