矩阵快速幂+二分 poj3233

 #include <iostream>
 #include <cstdio>
 #include <string>
 #include <cstring>
 #include <stdlib.h>
 #include <math.h>
 #include <ctype.h>
 #include <queue>
 #include <map>
 #include <set>
 #include <algorithm>

 using namespace std;
 int mod;
 int n;
 struct matrix
 {
     int ma[][];
 }init, res;
 matrix Mult(matrix x, matrix y)
 {
     int i,j,k;
     matrix tmp;
     memset(tmp.ma,,sizeof(tmp.ma));
     for(i=;i<n;i++)
     for(j=;j<n;j++)
     for(k=;k<n;k++)
     tmp.ma[i][j]=(tmp.ma[i][j]+x.ma[i][k]*y.ma[k][j])%mod;
     return tmp;
 }
 matrix Pow(matrix x, int k)
 {
     matrix tmp;
     int i,j;
     memset(tmp.ma,,sizeof(tmp.ma));
     for(i=;i<n;i++) tmp.ma[i][i]=;
     while(k){
         if(k&) tmp=Mult(tmp,x);
         x=Mult(x,x);
         k>>=;
     }
     return tmp;
 }
 matrix add(matrix x, matrix y)
 {
     int i, j;
     matrix tmp;
     for(i=;i<n;i++)
     for(j=;j<n;j++)
     tmp.ma[i][j]=(x.ma[i][j]+y.ma[i][j])%mod;
     return tmp;
 }
 matrix sum(matrix x, int k)
 {
     matrix tmp, y;
     if(k==) return x;
     tmp=sum(x,k/);
     if(k&){
         y=Pow(x,k/+);
         tmp=add(Mult(y,tmp),tmp);
         return add(tmp,y);
     }
     else{
         y=Pow(x,k/);
         return add(Mult(y,tmp),tmp);
     }
 }
 int main()
 {
     int k,m,x,i,j;
     scanf("%d%d%d",&n,&k,&mod);
     for(i=;i<n;i++){
         for(j=;j<n;j++){
             scanf("%d",&x);
             init.ma[i][j]=x%mod;
         }
     }
     res=sum(init, k);
     for(i=;i<n;i++){
         for(j=;j<n;j++){
             printf("%d",res.ma[i][j]);
             if(j!=n-) printf(" ");
         }
         puts("");
     }
     return ;
 }