快比赛了还不会555，矩阵快速幂如何实现

请给出方法和代码，主要是矩阵不会写
必采纳！谢谢！

import numpy as np
nkl = input().split()
n=int(nkl[0])
k=int(nkl[1])
A=np.ones((n,n))

for i in range(n):
    vs=input().split()
    for j in range(n):
        A[i][j]=int(vs[j])
        A[i][j]=pow(A[i][j],k)%(pow(10,9)+7)

print(A)