矩阵快速幂取模

#include <iostream>
using namespace std;
const long long int maxn = 101, mod = 1000000007;
struct 
{
    long long int s[maxn][maxn];
};
long long a, b;
matrix res, ans;
matrix mul(matrix A, matrix B, long long int n)
{
    matrix C;
    for (long long int i = 1; i <= n; i++)
        for (long long int j = 1; j <= n; j++)
            C.s[i][j] = 0;

    for (long long int i = 1; i <= n; i++)
        for (long long int j = 1; j <= n; j++)
            for (long long int k = 1; k <= n; k++)
                C.s[i][j] += A.s[i][k] * B.s[k][j]%mod;

    return C;
}
void quick(long long int y, long long int n)//快速幂
{

    res = ans;
    while (y)
    {
        if (y & 1)
            ans = mul(ans, res, n);
        res = mul(res, res, n);
        for (int i = 1; i <= n; i++)
        {
            for (int j = 1; j <= n; j++)
            {
                ans.s[i][j] %= mod;
                res.s[i][j] %= mod;
            }
        }
        y = y >> 1;
    }
}

int main()
{
    long long int n, y;
    scanf("%lld%lld", &n, &y);
    --y;
    matrix ans0;
    for (int i = 1; i <= n; i++)
        for (int j = 1; j <= n; j++)
            cin >> ans.s[i][j];
    quick(y, n);
    //ans0 = mul(ans, ans0, n);
    /* for(int i=1 ;i<3 ;i++)
        {
            for(int j=1 ;j<3 ;j++)
                printf("%lld ",ans0.s[i][j]);
            puts("");
        }*/
    for (int i = 1; i <= n; i++)
    {
        for (int j = 1; j <= n; j++)
            cout<< ans.s[i][j]%mod<< ' ';
        puts("");
    }
    return 0;
}