我正在实施Kraskov等人（2004）的Estimating Mutual Information论文中的公式8，并存在以下问题：
鉴于向量X = [X_1,...,X_N]和r = [r_1,...,r_N]，我需要计算A= [A_1,...A_N]，其中A_i是X中r_i内的点数半径为X_i。
如果r是固定数字，也就是说，如果所有点的每个点周围的半径都是固定的，我可以轻松使用rangesearch。但因为它是一个矢量（每个点的半径不同），我不知道如何快速完成。穷举搜索（或制作任何N^2大小的距离数组）并不好，因为N大约为一百万。