poj3233Matrix Power Series

链接

也是矩阵经典题目  二分递归求解

a+a^2+a^3+..+a^(k/2)+a^(k/2+1)+...+a^k = a+a^2+..+a^k/2+a^k/2(a^1+a^2+..+a^k/2)(偶数）

a+a^2+a^3+..+a^(k/2)+a^(k/2+1)+...+a^k = a+a^2+..+a^k/2+a^k/2(a^1+a^2+..+a^k/2)+a^k。 奇数

 #include <iostream>
 #include<cstdio>
 #include<cstring>
 #include<algorithm>
 #include<stdlib.h>
 #include<vector>
 #include<cmath>
 #include<queue>
 #include<set>
 using namespace std;
 #define N 1e9
 #define LL long long
 #define INF 0xfffffff
 const double eps = 1e-;
 const double pi = acos(-1.0);
 const double inf = ~0u>>;
 struct Mat
 {
     int mat[][];
 };
 int n,mod;
 Mat operator + (Mat a,Mat b)
 {
     Mat c;
     int i,j;
     for(i =  ; i < n ;i++)
         for(j =  ;j < n ;j++)
         {
             if(a.mat[i][j]+b.mat[i][j]>mod)
             c.mat[i][j] = (a.mat[i][j]+b.mat[i][j])%mod;
             else
             c.mat[i][j] = a.mat[i][j]+b.mat[i][j];
         }
     return c;
 }
 Mat operator * (Mat a,Mat b)
 {
     Mat c;
     memset(c.mat,,sizeof(c.mat));
     int i,j,k;
     for(k = ; k < n ; k++)
     {
         for(i =  ; i < n ;i++)
         {
             if(a.mat[i][k]==) continue;
             for(j =  ;j < n ;j++)
             {
                 if(b.mat[k][j]==) continue;
                 c.mat[i][j] = (c.mat[i][j]+a.mat[i][k]*b.mat[k][j])%mod;
             }
         }
     }
     return c;
 }
 Mat operator ^(Mat a,int k)
 {
     Mat c;
     int i,j;
     for(i = ; i < n ;i++)
         for(j = ; j < n ;j++)
         c.mat[i][j] = (i==j);
     for(; k ;k >>= )
     {
         if(k&) c = c*a;
         a = a*a;
     }
     return c;
 }
 Mat solve(Mat x,int k)
 {
     if(k==) return x;
     Mat c ;
     c = x^k;
     Mat a = solve(x,k/);
     Mat b = x^(k/);
     if(k&) c = a+b*a+c;
     else c = a+b*a;
     return c;
 }
 int main()
 {
     int t;
     int i,j;
     while(scanf("%d%d%d",&n,&t,&mod)!=EOF)
     {
         Mat x;
         for(i =  ; i < n ;i++)
             for(j = ; j < n ;j++)
             scanf("%d",&x.mat[i][j]);
         x = solve(x,t);
         for(i = ; i < n ;i++)
         {
             for(j = ; j < n-; j++)
             printf("%d ",x.mat[i][j]%mod);
             printf("%d\n",x.mat[i][n-]%mod);
         }
     }
     return ;
 }