【洛谷 p3390】模板-矩阵快速幂（数论）

题目：给定n*n的矩阵A，求A^k。

解法：利用矩阵乘法的定义和快速幂解答。注意用负数，但是数据太弱没有卡到我......（P.S.不要在 typedef long long  LL; 前使用 LL......━━(￣ー￣*|||━━）

P.S.在multi函数里，若将所有相乘的和先加起来不会爆 long long ，那就最后再模会快不少。

#include<cstdio>
#include<cstdlib>
#include<cstring>
#include<iostream>
using namespace std;
typedef long long LL;

const int N=110,D=1010;
const LL K=(LL)1e12+10,mod=(LL)1e9+7;
int n;
struct mat{LL s[N][N];}a,b,c,d,e,f;

mat multi(mat c,mat d)
{
    for (int i=1;i<=n;i++)
      for (int j=1;j<=n;j++)
      {
        e.s[i][j]=0;
        for (int k=1;k<=n;k++)
          e.s[i][j]=(e.s[i][j]+mod+(c.s[i][k]*d.s[k][j])%mod+mod)%mod;
      }
    return e;
}
mat qpow(mat a,LL k)
{
    b=a; k--;
    while (k>0)
    {
      if (k&1) b=multi(b,a);
      a=multi(a,a);
      k>>=1;
    }
    return b;
}
int main()
{
    LL k;
    scanf("%d%lld",&n,&k);
    for (int i=1;i<=n;i++)
     for (int j=1;j<=n;j++)
       scanf("%lld",&f.s[i][j]);
    f=qpow(f,k);
    for (int i=1;i<=n;i++)
    {
     for (int j=1;j<=n;j++)
       printf("%lld ",f.s[i][j]);
     printf("
");
    }
    return 0;
}