Potenzen von Diagonalmatrizen berechnen sich über Potenzen der einzelnen Diagonalelemente.
Potenzen anderer berechnest Du am effektivsten über binäre Darstellung, d.h. 2er-Potenzen
\( 2017 = 1+32+64+128+256+512+1024 \)
\( A^{2017} = A^{1}\cdot A^{32}\cdot A^{64}\cdot A^{128}\cdot A^{256}\cdot A^{512}\cdot A^{1024} \)
Du startest mit
\(A = A \), \( B = A = A^1 \)
\( \circledast \) Wenn Du \(B\) brauchst (hier ja), dann gilt \( A = A\cdot B \)
Dann immer \( B = B\cdot B \) und weiter mit \( \circledast \).
\(B\) wird dadurch immer wieder quadriert, wobei Du wegen Effizienz nicht \(B^2\), sondern \(B\cdot B\) rechnest. Du baust damit eine Reihe auf: \( B^1 \to B^2 \to B^4 \to B^8 \to \dots \).
Immwer wenn Du eine davon brauchst (das sagt Dir die Binärdarstellung von \(2017 \) wird sie multipliziert.
Grüße,
M.B.