sage 中的阿贝尔群商答案

【问题标题】：Abelian group quotient in sagesage 中的阿贝尔群商
【发布时间】：2016-05-25 23:47:58
【问题描述】：

让d1 和d2 成为整数Z 上的矩阵。如何在 Sage 中计算组商ker d1 / im d2？

到目前为止，我已经能够计算内核和图像的基础，如下所示：

M24 = MatrixSpace(IntegerRing(),2,4)
d1 = M24([-1,1, 1,-1, -1,1, 1,-1])
kerd1 = d1.right_kernel().basis()

M43 = MatrixSpace(IntegerRing(),4,3)
d2 = M43([1,1,-1, 1,-1,-1, 1,-1,1, 1,1,1])
imd2 = d2.column_space().basis()

给出输出：

kerd1 = [
  (1, 0, 0, -1),
  (0, 1, 0, 1),
  (0, 0, 1, 1)
]

imd2 = [
  (1, 1, 1, 1),
  (0, 2, 0, -2),
  (0, 0, 2, 2)
]

我尝试这样计算商：

Z4.<a,b,c,d> = AbelianGroup(4, [0,0,0,0])
G = Z4.subgroup([a/d, b*d, c*d])
H = Z4.subgroup([a*b*c*d, b^2/d^2, c^2*d^2])
G.quotient(H)

但我有一个NotImplementedError。

【问题讨论】：

标签： sage

【解决方案1】：

我找到了两种方法：

d1 = matrix(ZZ,4,2, [-1,1, 1,-1, -1,1, 1,-1]).transpose()
d2 = matrix(ZZ,4,3, [1,1,-1, 1,-1,-1, 1,-1,1, 1,1,1])
(d1.right_kernel() / (d2.column_space())).invariants()
# OUTPUT: (2, 2)

ChainComplex([d2, d1]).homology()[1]
# OUTPUT: C2 x C2

【讨论】：