Find invertible matrix pairs where both matrices are integers

Question

Douglas Anderson 2022-1-7

1
链接

此问题的直接链接

https://ww2.mathworks.cn/matlabcentral/answers/1624660-find-invertible-matrix-pairs-where-both-matrices-are-integers

评论： John D'Errico 2022-1-8

采纳的回答： John D'Errico

Hello!

I think the summary says it all! Looking for 3x3, but assume that's not critical.

Have found solution for matrix pair where one is integer, but other is not.

Thanks!

Doug Anderson

3 个评论
显示 1更早的评论隐藏 1更早的评论

James Tursa 2022-1-7

编辑：James Tursa 2022-1-7

The cofactors comment was just by way of understanding the problem ... not a suggested method to use. If you have a matrix M with all integer elements and you know the determinant is 1, then you could probably just use round(inv(M)) as a first cut as long as the sizes of the integers are not too big. Do you already have a method for generating the original matrix of integer elements with determinant 1, or are you asking for that as well?

What is this for? I.e., what is the actual problem you are working on?

Douglas Anderson 2022-1-7

I have a method to generate a matrix of integer elements, i.e., doubles that look like integers! I know you can't invert a true integer matrix, AND "isinteger" only identifies if it's a true integer class. Is there a method to say that this is an integer "for all practical purposes"?

This is for security software. Have the matrix I posted, just looking for others. It's been working well.

Again, thanks.

Doug

请先登录，再进行评论。

请先登录，再回答此问题。

Answer 1

John D'Errico 2022-1-7

3
链接

此回答的直接链接

https://ww2.mathworks.cn/matlabcentral/answers/1624660-find-invertible-matrix-pairs-where-both-matrices-are-integers#answer_870530

编辑：John D'Errico 2022-1-7

在 MATLAB Online 中打开

It seems trivial.

A = eye(3)
A = 3×3
     1     0     0
     0     1     0
     0     0     1
inv(A)
ans = 3×3
     1     0     0
     0     1     0
     0     0     1

Both are entirely integer. As it turns out, A is idempotent, so it is its own inverse, but that was not expressed as an issue.

But surely you are asking about other matrices, that are also fully integer, yet have an all-integer inverse? This too should be entirely doable. The trick is to consider how to create them.

N = 3;
M = 5;
L = tril(randi(2*M,[N,N])-M,-1) + eye(N)
L = 3×3
     1     0     0
     1     1     0
     1    -4     1
A = L*L'
A = 3×3
     1     1     1
     1     2    -3
     1    -3    18
format long g
inv(A)
ans = 3×3
                        27                       -21                        -5
                       -21                        17                         4
                        -5                         4                         1

So A is a fully integer randomly created matrix. It has the property of a fully integer inverse. (if we ignore the crap that sometimes arise in the least significant bits of the inverse. They are just due to floating point trash.) If you wish to know why this works, look carefully at how I created the matrix. Look carefully at the matrix L.

A nice property of this creation scheme, is it will work for any size of matrix A. Just change the value of N. As well, it will be blazingly fast, requiring no solve.