概率算法求PI

代码如下：

#include<iostream>
#include<ctime>
using namespace std;
int main()
{
 int i,n,sum=0;
 double x,y;
 cout<<"输入点的数量:"<<endl;
 cin>>n;
 srand(time(NULL));
 for(i=1;i<n;i++)
 {
       x=(double)rand()/RAND_MAX;
    y=(double)rand()/RAND_MAX;
    if(x*x+y*y<=1)
     sum++;
 }
 cout<<"PI="<<4.0*sum/n<<endl;
 return 0;
}