如何在SageMath中求两个数域的交集？

你提到的枚举二次子域判别式的方法仅适用于交集为多二次域的场景，通用性有限。这里提供一种不依赖交集类型的通用解法，核心思路是将两个数域嵌入到同一个复合域中，再直接计算子域交集：

通用步骤与代码示例 #

1. 定义目标数域
2. 构造包含两者的复合域（最小公共扩域）
3. 将原数域嵌入到复合域中，得到对应子域对象
4. 调用intersection()方法计算交集

# 示例：定义两个不同的数域
K = NumberField(x^4 - 5, 'a')  # Q(√5, i)，四次域
L = NumberField(x^6 - 5, 'b')  # Q(5^(1/6), ω)，六次域，ω为三次单位根

# 构造复合域（包含K和L的最小数域）
M = K.composite_field(L)

# 将K和L的生成元映射到复合域M中
k_gen_in_M = M(K.gen())
l_gen_in_M = M(L.gen())

# 获取K和L在M中对应的子域对象
K_sub = M.subfield(k_gen_in_M)[0]
L_sub = M.subfield(l_gen_in_M)[0]

# 计算并输出交集
intersection = K_sub.intersection(L_sub)
print("交集数域：", intersection)
print("交集次数：", intersection.degree())

方法说明 #

• 复合域的构造确保两个数域都能作为其子域存在，避免了嵌入到代数闭包时的歧义
• intersection()是Sage内置的子域计算工具，能自动处理任意次数、任意类型的数域交集，无需提前推断交集结构
• 对于简单场景（比如一个数域是另一个的子域），该方法也能直接返回正确结果（如K包含于L时，交集即为K）

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