Power of Matrix 等比数列求和 矩阵版！

#include<iostream>
#include<cstdio>
#include<cmath>
#include<cstring>
#include<algorithm>
#include<queue>
#include<vector>
#include<cmath>
#include<map>
#include<stack>
#include<set>
#include<string>
using namespace std;
typedef long long  LL;
typedef unsigned long long ULL;
#define MAXN 49
#define MOD 10000007
#define INF 1000000009
const double eps = 1e-;
//矩阵快速幂
int n, k;
struct Mat
{
    int a[MAXN][MAXN];
    Mat()
    {
        memset(a, , sizeof(a));
    }
    Mat operator *(const Mat& rhs)
    {
        Mat ret;
        for (int i = ; i < n; i++)
        {
            for (int j = ; j < n; j++)
            {
                if (a[i][j])
                {
                    for (int t = ; t < n; t++)
                    {
                        ret.a[i][t] = (ret.a[i][t] + a[i][j] * rhs.a[j][t])%;
                    }
                }
            }
        }
        return ret;
    }
    Mat operator +(const Mat& rhs)
    {
        Mat ret;
        for (int i = ; i < n; i++)
        {
            for (int j = ; j < n; j++)
            {
                ret.a[i][j] += (a[i][j] + rhs.a[i][j])%;
            }
        }
        return ret;
    }
};
Mat e;
Mat fpow(const Mat& m, int b)
{
    Mat ans, tmp = m;
    for (int i = ; i < n; i++)
        ans.a[i][i] = ;
    while (b != )
    {
        if (b & )
            ans = tmp*ans;
        tmp = tmp * tmp;
        b >>= ;
    }
    return ans;
}
Mat sum(const Mat& m, int k)
{
    if (k == ) return m;
    else if (k %  == )
    {
        return (e + fpow(m, k / )) * sum(m, k / );
    }
    else if (k %  == )
    {
        Mat tmp = fpow(m, k /  + );
        return    (e + tmp)*sum(m, k / ) + tmp;
    }
}

int main()
{
    while (scanf("%d%d", &n,&k), n)
    {
        for (int i = ; i < n; i++)
            e.a[i][i] = ;
        Mat M;
        for (int i = ; i < n; i++)
        {
            for (int j = ; j < n; j++)
            {
                scanf("%d", &M.a[i][j]);
                M.a[i][j] %= ;
            }
        }
        M = sum(M, k);
        for (int i = ; i < n; i++)
        {
            printf("%d", M.a[i][]);
            for (int j = ; j < n; j++)
            {
                printf("% d", M.a[i][j]);
            }
            printf("\n");
        }
        printf("\n");
    }
}