你好呀！我来帮你解决这个Sage求解二元二次型只返回单个整数解的问题～

首先得明确：Sage里的BQF.solve_integer()方法本来就只设计返回任意一个有效整数解，不会自动枚举所有可能的解，这就是为什么你只拿到一个结果的原因。不过我们有好几种办法可以拿到所有正整数解，不管是在Sage内部扩展，还是用其他工具都能实现，下面给你具体说：

一、在Sage内部获取全部解的两种方法 #

1. 暴力枚举法（简单直接，适合中小规模的m）

因为我们要找的是满足x² + n y² = m的正整数解(x,y)，那y的取值范围其实是有限的：y最大不会超过sqrt(m/n)（因为x²必须是非负的）。我们可以写个小脚本遍历所有可能的y值，检查对应的m - n y²是否是完全平方数就行，逻辑非常直观。

比如针对你的例子，在Sage里运行这段代码：

def find_all_positive_solutions(m, n):
    solutions = []
    max_y = int(sqrt(m // n)) + 1  # 加1避免边界值遗漏
    for y in range(1, max_y + 1):
        remainder = m - n * y * y
        if remainder <= 0:
            continue
        x = integer_sqrt(remainder)
        if x * x == remainder:
            solutions.append((x, y))
    return solutions

# 测试你的第一个例子
print(find_all_positive_solutions(713, 22))  # 输出 [(25, 2), (19, 4)]
# 测试第二个例子
print(find_all_positive_solutions(713, 37))  # 输出 [(26, 1), (11, 4)]

这个方法的好处是不用纠结二次型的复杂数论性质，写起来快，对于你提到的brilliant数（两个素数乘积的数）来说，计算量完全没问题。

2. 利用二次型自同构群生成解

如果m很大，暴力枚举效率不够，我们可以用二次型的自同构群来从已知的一个解生成其他解。对于x² + n y²这种对角二次型，它的自同构变换能把一个合法解映射到另一个非等价解。不过这个方法需要一定的数论基础，步骤相对复杂；对于你用到的场景，暴力枚举的效率已经足够高了，除非m特别大，否则没必要用这个方法。

二、其他可以直接返回全部解的工具 #

如果不想在Sage里写脚本，也可以用这些工具：

• PARI/GP：它的qfbsolve()函数可以返回二次型的多个解，或者用相关枚举命令直接输出所有正整数解
• Mathematica：用Reduce函数可以直接求解并返回所有满足条件的正整数解，比如输入Reduce[x^2 + 22 y^2 == 713 && x>0 && y>0, {x,y}, Integers]，会直接给出两个解
• 普通Python脚本：不用Sage的话，用Python自带的math库也能实现刚才的暴力枚举逻辑，代码和Sage里的几乎一样，只是把integer_sqrt换成int(math.isqrt(remainder))就行

补充说明 #

你提到的那些多个解，本质上是二次型的非等价解，Sage的默认方法不会主动去枚举这些，因为它的核心目标是快速找到一个可行解，而不是穷尽所有解。

备注：内容来源于stack exchange，提问作者ReverseFlowControl