给定圆的半径和圆心的 x、y 坐标,写一个在圆中产生均匀随机点的函数 randPoint
。
说明:
- 输入值和输出值都将是浮点数。
- 圆的半径和圆心的 x、y 坐标将作为参数传递给类的构造函数。
- 圆周上的点也认为是在圆中。
randPoint
返回一个包含随机点的x坐标和y坐标的大小为2的数组。
示例 1:
输入:
["Solution","randPoint","randPoint","randPoint"]
[[1,0,0],[],[],[]]
输出: [null,[-0.72939,-0.65505],[-0.78502,-0.28626],[-0.83119,-0.19803]]
示例 2:
输入:
["Solution","randPoint","randPoint","randPoint"]
[[10,5,-7.5],[],[],[]]
输出: [null,[11.52438,-8.33273],[2.46992,-16.21705],[11.13430,-12.42337]]
输入语法说明:
输入是两个列表:调用成员函数名和调用的参数。Solution
的构造函数有三个参数,圆的半径、圆心的 x 坐标、圆心的 y 坐标。randPoint
没有参数。输入参数是一个列表,即使参数为空,也会输入一个 [] 空列表。
思路:这道题的直观思路是生成[0,radius]的均匀分布len,生成[0,2*PI]的均匀分布degree,那么圆内的点分布就是:
x_center+len*cos(degree),y_center+len*sin(degree) 但是在这样做是不对的,对于角度degree,如果是均匀分布是没问题的,但是对于半径,如果是均匀分布是不对的,按照如上思路得出的分布图如下左图所示(而我们期望的是右图的结果):
可以看到在半径小的地方点密集,半径大的地方点稀疏,这是因为如果是均匀分布的len,在确定的len上的点的数量是一定的,换句话说,是把[0,radius]的点看成一个一个圆环的点的叠加且对于每个圆环上的点的数量是确定的,这显然是不对的,我们希望对于半径大的圆环上的点的数量多些,保证大半径和小半径的圆环的点的距离相等。
所以这里要用到几个概念:
我们希望对于确定的半径为r的圆环的圆周的点的分布数量正比于半径,即和圆周长相等:pdf(probability density function)=2*PI*r
那么对应的cdf(culmulative distribution function)=pdf的积分=PI*(r^2)
有了pdf和cdf,我们希望得到变量r的分布,即通过pdf和cdf反推出变量r的随机样本。
这里要用到:inverse transform sampling逆变换采样,先看维基百科的解释:
简单的说就是我们有个一个变量x,以及变量x对应的pdf和cdf,我们希望反推出x对应的随机样本(换句话说就是x的随机样本是怎么随着x变化的),只需要对cdf求反函数即可,对于这道题而言cdf的反函数为:y=sqrt(r),即x的随机样本与x是根号的关系才能得到我们想要的x的pdf和cdf。
参考代码:
class Solution {
public:
Solution(double radius, double x_center, double y_center) {
this->radius = radius;
this->x_center = x_center;
this->y_center = y_center;
}
vector<double> randPoint() {
double len = sqrt((double)rand() / RAND_MAX)*radius;
double degree =2* M_PI* ((double)rand() / RAND_MAX);
return { x_center + len * cos(degree),y_center + len * sin(degree) };
}
private:
double radius;
double x_center;
double y_center;
};
(这里要注意的是c++的rand函数是均匀分布,我们不能对r取根号,只能对rand函数取根号,使得rand的均匀分布变成我们想要的根号分布)
方法二:拒绝抽样,原理就是蒙特卡洛原理,利用采样点来无限逼近真实点的分布,在这道题中随机生成在圆的外切正方形的点,如果落到圆内(包括圆上)就返回值,否则重新生成点。
参考代码:
class Solution {
public:
Solution(double radius, double x_center, double y_center) {
r = radius;
centerX = x_center;
centerY = y_center;
}
vector<double> randPoint() {
while (true) {
double x = (2 * ((double)rand() / RAND_MAX) - 1)*r;
double y= (2 * ((double)rand() / RAND_MAX) - 1)*r;
if (x*x + y * y <= r*r) return {centerX+x,centerY+y};
}
}
private:
double r, centerX, centerY;
};