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
A2017=A1⋅A32⋅A64⋅A128⋅A256⋅A512⋅A1024
Du startest mit
A=A, B=A=A1
⊛ Wenn Du B brauchst (hier ja), dann gilt A=A⋅B
Dann immer B=B⋅B und weiter mit ⊛.
B wird dadurch immer wieder quadriert, wobei Du wegen Effizienz nicht B2, sondern B⋅B rechnest. Du baust damit eine Reihe auf: B1→B2→B4→B8→….
Immwer wenn Du eine davon brauchst (das sagt Dir die Binärdarstellung von 2017 wird sie multipliziert.
Grüße,
M.B.