SageMath中获取LLL坐标变换矩阵：有理数矩阵报错咨询

问题背景 #

对有理数矩阵调用LLL(transformation=True)以获取坐标变换矩阵时，会触发TypeError: unsupported operand parent(s) for /: '<class 'tuple'>' and 'Integer Ring'错误，但直接调用LLL()能正常返回约化基。需要在支持有理数矩阵的前提下，获取将原基转换为LLL约化基的坐标变换矩阵。

解决方案 #

方法1：整数矩阵中转处理

SageMath的LLL变换矩阵功能对整数矩阵支持完好，可先将有理数矩阵缩放为整数矩阵，处理后再转换回有理数域：

# 定义原有理数矩阵
m = matrix(QQ, 2, 2, [[5, 0], [0, 2]])
# 转换为整数矩阵（若原矩阵含分数，需先求所有元素的公分母，每行乘以公分母）
m_int = matrix(ZZ, m)
# 调用LLL获取约化基和整数变换矩阵
lll_basis, T_int = m_int.LLL(transformation=True)
# 将变换矩阵转换回有理数域
T = matrix(QQ, T_int)

# 验证正确性：原矩阵 × 变换矩阵 = LLL约化基
print("原矩阵 × 变换矩阵结果：")
print(m * T)
print("LLL直接输出的约化基：")
print(m.LLL())

方法2：手动计算变换矩阵

LLL约化基是原基的线性组合，因此变换矩阵可通过原矩阵的逆矩阵与LLL基矩阵相乘得到，无需依赖transformation参数：

# 定义含分数的有理数矩阵示例
m = matrix(QQ, 2, 2, [[1/2, 0], [0, 3/4]])
# 获取LLL约化基
lll_basis = m.LLL()
# 计算变换矩阵T：满足 lll_basis = m × T
T = m.inverse() * lll_basis

# 验证正确性
print("变换矩阵T：")
print(T)
print("验证等式 m×T == lll_basis：", m * T == lll_basis)

错误原因说明 #

SageMath的LLL算法实现最初针对整数格设计，有理数矩阵的LLL(transformation=True)参数存在兼容性bug，导致类型运算错误。上述两种方法通过绕开该参数的直接调用，可在有理数矩阵场景下正常获取变换矩阵。

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