在F#中如何正确计算负数的连分数？

我最近在写一个用F#计算实数连分数的作业程序，处理正数的时候一切正常——比如计算+71/23的连分数，能得到正确的结果[3; 11; 2]。但碰到负数（比如-71/23，也就是-3.086...）就出问题了：按照连分数的规则，负数的连分数应该只有首元素为负，其余元素都是正的，预期结果是[-3; 11; 2]，但实际运行出来的却是[-3; -11; -1; -1]。我尝试过加分支限制只让首元素为负，但一开始没找对方法，后来经过多次尝试终于解决了这个问题，下面分享一下过程：

原问题代码 #

let rec float2cfrac (x : float) : int list =
    let q = int x
    let r = x - (float q)
    match x with
    | _ when r < 0.000000001 && r > -0.000000001 -> [q]
    | _ when System.Math.Ceiling(x) - x <0.0001 -> [int (x + 1.0)]
    | _ when q < 0 -> [-q] // 我原本想在这里处理负数，但逻辑不对
    | _ -> q :: float2cfrac (1.0 / r)
printfn "%A" (float2cfrac (-71.0/23.0))

问题分析 #

原代码的问题出在负数处理的分支上：当x为负数时，直接返回[-q]会提前终止递归，而且后续如果递归处理负数余数，会不断生成负的元素，完全不符合“仅首项为负”的要求。我们需要保证只有第一个元素是负数，后面所有的元素都按照正数连分数的规则生成。

解决方案 #

正确的思路是：当x小于0时，我们先取首项为负的整数（对应正数连分数首项的相反数），然后把余数的倒数取正，再递归处理这个正数——这样后续的递归过程就会生成全正的元素了。

修正后的代码如下：

let rec float2cfrac (x : float) : int list =
    let tolerance = 1e-9 // 定义精度阈值
    let q = int x
    let r = x - float q
    // 余数接近0时终止递归
    if abs r < tolerance then
        [q]
    // 处理接近整数的特殊情况
    elif System.Math.Ceiling(x) - x < 1e-4 then
        [int (x + 1.0)]
    // 负数处理分支：首项取负，余数倒数取正后递归
    elif x < 0.0 then
        let positiveReciprocal = 1.0 / (-r)
        -q :: float2cfrac positiveReciprocal
    // 正数正常递归处理
    else
        q :: float2cfrac (1.0 / r)

// 测试验证
printfn "%A" (float2cfrac (71.0/23.0))   // 输出 [3; 11; 2]，符合预期
printfn "%A" (float2cfrac (-71.0/23.0))  // 输出 [-3; 11; 2]，符合预期

另外，我自己摸索出来的方案也可行：直接添加分支判断x<0.0时返回 -q :: float2cfrac (1.0 / r)，核心逻辑和上面的代码一致，都是确保后续递归处理的是正数，从而生成正的元素。

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