模板
#include<cstdio>
#include<cmath>//pow函数,其实没啥用
using namespace std;
int n;long long k;
const int N=pow(10,9)+7;
struct node{
long long a[105][105];};
node shu,ans,mp;
//shu是输入的矩阵,ans是所求答案
node matrix(node x,node y){
for(int i=1;i<=n;i++)
for(int j=1;j<=n;j++){
mp.a[i][j]=0;
for(int p=1;p<=n;p++)
mp.a[i][j]=(mp.a[i][j]+x.a[i][p] * y.a[p][j])%N;