[板子]快速幂&矩阵快速幂

不会的来这看：https://www.cnblogs.com/CXCXCXC/p/4641812.html

简单的一说：当转换为二进制的时候有位运算这种黑科技，&相当于%2判断奇偶性。

x&1==0为偶，x&1==1为奇

&运算通常用于二进制取位操作，例如一个数 & 1 的结果就是取二进制的最末位。还可以判断奇偶x&1==0为偶，x&1==1为奇

只有在奇数情况的时候把base乘进去，每一次用base*base扩大平方，b的二进制去除一位。

int poww(int a,int b){
    int ans=1,base=a;
    while(b!=0){
        if(b&1==1)
        　　ans*=base;
        base*=base;
        b>>=1;
　 }
    return ans;
}

 接下来是矩阵快速幂。

摘自：blog.csdn.net/wust_zzwh/article/details/52058209

其中c[i][j]为A的第i行与B的第j列对应乘积的和，即:

#include<bits/stdc++.h>

#define LL long long
using namespace std;

LL n,k;

const long long pi=1e9+7;

struct ska{
    LL a[100+10][100+10];
}p,pp;

ska X(ska x,ska y){
    ska box;
    
    for(LL i=1;i<=n;i++){
        for(LL j=1;j<=n;j++){
            box.a[i][j]=0;
        }
    }
    
    for(LL i=1;i<=n;i++){
        
        for(LL j=1;j<=n;j++){
            
            for(LL k=1;k<=n;k++){
                
                box.a[i][j]=(box.a[i][j]+(x.a[i][k]*y.a[k][j])%pi)%pi;
                
            }
        }
    }
    
    return box;
}

ska quick_pow(LL kk){
    
    ska ans;
    
    for(LL i=1;i<=n;i++){
        ans.a[i][i]=1;
    }
   
    while(kk!=0){
        
        if(kk&1==1){
            
            ans=X(ans,p);
        }
        
        kk>>=1;
        
        p=X(p,p);
    
    }
    
    return ans;
}

int main(){
    scanf("%lld%lld",&n,&k);
    
    for(LL i=1;i<=n;i++){
        
        for(LL j=1;j<=n;j++){
            
            scanf("%lld",&p.a[i][j]);
            
        }
    }
    
    pp=quick_pow(k);
    
    for(LL i=1;i<=n;i++){
        
        for(LL j=1;j<=n;j++){
            
            printf("%lld ",pp.a[i][j]);
            
        }
        
        cout<<endl;
    }
    return 0;
}