使用R语言遍历质数对计算欧氏距离的代码问题排查

问题描述 #

• 目的：遍历给定质数序列，按顺序将相邻两个质数配对为(x,y)（如x=2,y=3；x=5,y=7……x=89,y=97），计算每对点到原点的欧氏距离；若需计算点与点之间的距离，也可对应调整逻辑。
• 需求：通过R语言for循环实现上述逻辑，原代码输出异常，需修正逻辑错误。

原代码错误分析 #

1. range()函数误用：R中range()返回输入值的最小/最大值，而非连续序列。原代码for (j in range(length(moves)-1))仅会循环两次（取1和length(moves)-1两个值），这是输出仅两条结果的核心原因。
2. 索引配对逻辑混乱：R向量索引从1开始，原代码的奇偶判断逻辑无法正确将相邻质数配对为(x,y)。
3. 结果存储失效：append(s, euclidean(x,y))未将结果赋值回s，导致最终结果向量始终为空。
4. 手动修改循环变量：j=+1属于多余操作，for循环会自动迭代变量，手动修改会打乱循环逻辑。

修正后的代码（点到原点的距离） #

# 定义点(x,y)到原点的欧氏距离函数
euclidean_to_origin <- function(x, y) sqrt(x^2 + y^2)

# 给定质数序列
moves <- c(2, 3, 5, 7, 11, 13, 17, 19, 23, 29, 31, 37, 41, 43, 47, 53, 59, 61, 67, 71, 73, 79, 83, 89, 97)

# 初始化结果存储向量
distance_results <- c()

# 按步长2循环，每两个质数为一组配对
for (i in seq(1, length(moves)-1, by = 2)) {
  x <- moves[i]
  y <- moves[i+1]
  current_dist <- euclidean_to_origin(x, y)
  distance_results <- c(distance_results, current_dist)
  print(current_dist) # 打印单组结果
}

# 输出所有结果
cat("所有配对点到原点的欧氏距离：\n")
print(distance_results)

扩展：计算相邻配对点之间的距离 #

若需求为计算前后配对点（如(2,3)与(5,7)）之间的欧氏距离，可使用以下代码：

# 定义两点之间的欧氏距离函数
euclidean_between_points <- function(p1, p2) sqrt(sum((p1 - p2)^2))

moves <- c(2, 3, 5, 7, 11, 13, 17, 19, 23, 29, 31, 37, 41, 43, 47, 53, 59, 61, 67, 71, 73, 79, 83, 89, 97)

# 将质数序列转换为配对点矩阵
points_matrix <- matrix(moves, ncol = 2, byrow = TRUE)

distance_results <- c()
# 遍历计算相邻点之间的距离
for (i in 2:nrow(points_matrix)) {
  current_dist <- euclidean_between_points(points_matrix[i-1, ], points_matrix[i, ])
  distance_results <- c(distance_results, current_dist)
  print(current_dist)
}

cat("相邻配对点之间的欧氏距离：\n")
print(distance_results)

内容的提问来源于stack exchange，提问作者Data Science Analytics Manager