行列式模


18

什么是用于计算与系数的整数矩阵的行列式的已知有效的算法,残基的环模。数可能不是素数,但复合材料(以便计算在环执行,而不是一个场)。Zmmm

据我所知(请参阅下文),大多数算法都是对高斯消去算法的修改。问题是这些程序的计算效率。

如果碰巧有其他方法,我也对此感到好奇。

提前致谢。

更新:

让我解释一下这个问题的根源。假设,是质数。因此Z m是一个场。在这种情况下,我们可以使用小于m的数字执行所有计算,因此我们对数字的所有运算都有一些不错的上限:加法,乘法和求逆---运行高斯消除所需的所有运算。mZmm

另一方面,在不是素数的情况下,我们无法对某些数字进行求逆。因此,我们需要一些技巧来计算行列式。m

现在,我很好奇完成这项工作的已知技巧是什么,以及是否可以在书籍和论文中找到这些技巧。


3
您所说的``高效''是什么意思?问题显然在P
大卫

2
固定不变?怎么给?m
迈克尔·布朗丹

2
你什么意思?可以用一元写吗?
迈克尔·布朗丹2012年

5
我还是不明白这个问题。整数矩阵的行列式可以在多项式时间内计算,因此您可以将该值取模。无需在执行部门ž 或发现的因式分解mZmm
大卫

2
@ValeriySokolov:那是基本的线性代数。例如,请检查Christos H. Papadimitriou 的计算复杂性问题11.5.3 。
伊藤刚(Tsuyoshi Ito)

Answers:


15

如果你知道的因式分解你可以计算模每个p è 分开,然后用中国剩余相结合的结果。如果e i = 1,则计算模p e i i很容易,因为这是一个字段。对于更大的e i,可以使用Hensel提升。 m=p1e1pnenpieiei=1pieiei


谢谢!就像我一直在寻找的东西。这是行列式的惯例吗?(欢迎参考)。
Valeriy Sokolov

6
这些是来自计算机代数的标准技术。看看von zur Gathen和Gerhard撰写的Modern Computer Algebra或任何其他有关计算机代数的书。对于您的特定问题,另请参阅Pan,Yu和Stewart撰写的以下论文comet.lehman.cuny.edu/vpan/pdf/pan146.pdf
MarkusBläser2012年

17

Mahajan和Vinay有一种组合算法可以在交换环上工作:http : //cjtcs.cs.uchicago.edu/articles/1997/5/contents.html


感谢您的回复,其中包含指向非常有趣的论文的链接。
Valeriy Sokolov

我还相信,自从本文的作者解决了更一般的问题(对于任何交换环)以来,算法就更加有效。
Valeriy Sokolov

“有”是指“已知”或“存在”(但尚未找到)吗?这是一个合理的猜测,但我有点怀疑以小小的复合数为模的因果环的结构能否对您有很大帮助。如果我错了,我会觉得很有趣。
Sasho Nikolov 2012年

1
@ValeriySokolov说句公道话,因为答案确实回答了您的问题,所以您可以考虑接受它(或者如果您希望等待可能不会更合理的更好答案)
Suresh Venkat 2012年

@SashoNikolov我发现Wolfram Mathematica以某种方式对此进行了计算。在“实施说明”中,他们说:Det使用模块化方法和行归约,并使用中文余数定理构造结果。我想知道他们到底在做什么,但是快速搜索却没有给我任何帮助。至于“小合成 ”,它仅表示我想将此环中加法和乘法的复杂度视为O 1 。那就是所有像O log m 的因素都被视为O 1 mO(1)O(logm)O(1)
瓦列里·索科洛夫

11

为了解决这个问题,有一种基于史密斯范式快速确定性算法,该算法的最坏情况下的复杂度在整数模的矩阵相乘成本上限较高。对于任何矩阵A,算法都会输出其Smith正规形式,从中可以轻松计算det A mAdet(A)

ωn×nZmO(nω)Zm

AZmn×ndet(A)O(nω)Zm

当这是在1996年写的时候,没有渐近地更快的选择(本文提到了以前存在相同边界的算法,但我不知道哪些算法是概率性的)。

mm

O(log2m)O(M(logm)loglogm)M(t)tω


θω

也许,我不知道最常见的表示法。
Juan Bermejo Vega'2

我认为您是对的,我将其更改为“主流”
Juan Bermejo Vega
By using our site, you acknowledge that you have read and understand our Cookie Policy and Privacy Policy.
Licensed under cc by-sa 3.0 with attribution required.