1、设rand(s,t)返回[s,t]之间的随机小数,利用该函数在一个半径为R的圆内找随机n个点,并给出时间复杂度分析。
float GetRand( float fMin, float fMax )
{
assert( fMax > fMin );
float fRand=0.f, f1=0.f;
int n; //存(fMax-fMin)的整数部分
float f; //存(fMax-fMin)的小数部分
n = static_cast<int>(fMax-fMin);
f = (fMax-fMin) - n;
if ( n>0 )
fRand = static_cast<float>( rand() % n );
else
fRand = 0.f;
do
{
f1 = static_cast<float>( rand() ) / RAND_MAX;
f1 *= ( rand()%2==0 ? 1 : -1 );
} while ( fRand+f1<0 || fRand+f1 > (fMax-fMin) );
fRand += (fMin+f1);
return fRand;
}
void ProductNrand(float r,int n){
int count=0;
float x,y,R;
while(count<n){
x=GetRand(-r,r);
y=GetRand(-r,r);
if(sqrt(x*x+y*y)<=r){
count++;
cout<<x<<","<<y<<endl;
}
}
}