R语言：如何低改动现有双数GCD代码求1到n序列的最大公约数

最大公约数的运算天然满足结合律，即gcd(a,b,c) = gcd(gcd(a,b), c)，基于这一性质，我们几乎不需要改动原有两数GCD的核心计算逻辑，只需要做极少量的适配，就能实现对1到n这类数值序列的整体GCD计算。

方案一：少量修改函数内部逻辑，保留原有两数计算能力 #

原有核心计算代码完全不做改动，仅在函数开头增加3行左右的输入判断逻辑，既兼容原来传入两个整数的调用方式，也支持直接传入数值序列：

GCD=function(a, b=NULL){
  # 新增的序列适配逻辑：检测到输入为长度≥2的数值向量时，逐次调用两数GCD计算
  if(length(a) >= 2 && is.null(b)){
    return(Reduce(GCD, a))
  }
  # ——— 以下为原有完全未修改的代码 ———
  m=min(a, b)
  while(a%%m>0|b%%m>0){m=m-1}
  return(m)
}

修改完成后支持的调用场景：

• 兼容原有两数计算：GCD(12, 18) 返回结果6
• 支持1到n的序列计算：GCD(1:10) 直接返回整个序列的GCD结果，按照数学规则，只要序列中包含1，最终返回值即为1
• 支持任意长度数值列表：比如GCD(c(16,24,32,40)) 会逐次计算得到最终GCD为8

方案二：零修改原有代码，调用时复用原有函数 #

如果完全不想改动原有GCD函数的任何一行代码，只需要借助R内置的Reduce函数，把两数GCD的计算逐次套用到1到n的序列上即可，完全实现零改动：

# 原有代码完全保留不做任何修改
GCD=function(a, b){
m=min(a, b)
while(a%%m>0|b%%m>0){m=m-1}
return(m)}

# 计算1到n的序列GCD，比如取n=20
n <- 20
seq_gcd <- Reduce(GCD, 1:n)

小提示：原有函数采用从较小值逐次递减遍历的方式找公约数，在序列长度较短、数值不大的场景下计算性能完全满足需求；如果后续需要处理超长序列，可再单独优化核心算法，无需改动外层适配逻辑。

内容的提问来源于stack exchange，提问作者himalaya