矩阵快速幂

应该还蛮好理解的。。

#include <cstdio>
using namespace std;
typedef long long LL;

const LL maxn=105, prime=1e9+7;
LL n, k, mat[maxn][maxn];
LL ans[maxn][maxn], tmp[maxn][maxn], tmp2[maxn][maxn];

void plus(LL a[maxn][maxn], LL b[maxn][maxn]){
    for (LL i=0; i<n; ++i) for (LL j=0; j<n; ++j){
        tmp2[i][j]=0;
        for (LL k=0; k<n; ++k)
            tmp2[i][j]=(tmp2[i][j]+(a[i][k]*b[k][j])%prime)%prime;
    }
    for (LL i=0; i<n; ++i) for (LL j=0; j<n; ++j)
        a[i][j]=tmp2[i][j];
}

int main(){
    scanf("%lld%lld", &n, &k);
    for (LL i=0; i<n; ++i) for (LL j=0; j<n; ++j)
        scanf("%lld", &mat[i][j]), tmp[i][j]=mat[i][j];
    for (LL i=0; i<n; ++i) for (LL j=0; j<n; ++j)
        if (i==j) ans[i][j]=1;
    while (k){
        if (k&1) plus(ans, tmp);
        plus(tmp, tmp);
        k>>=1;
    }
    for (LL i=0; i<n; ++i){
        for (LL j=0; j<n; ++j)
            printf("%lld ", ans[i][j]);
        printf("
");
    }
    return 0;
}