矩阵快速幂模板

struct mat{
  ll s[N][N];
  mat(){
    memset(s,0,sizeof s);
  }

  mat operator *(const mat c){
    mat ans ;
    memset(ans.s,0, sizeof ans.s);
    for(int i=0;i<N;i++)
      for(int j=0;j<N;j++)
        for(int t=0;t<N;t++)
          ans.s[i][j] =(ans.s[i][j]+ s[i][t]*c.s[t][j])%mod;
    return ans;
  }

  mat operator % (int mod){
    for(int i=0;i<N;i++)
      for(int j=0;j<N;j++)
        s[i][j] %= mod;
    return *this;
  }

  mat operator + (const mat c){
    for(int i=0;i<N;i++)
      for(int j=0;j<N;j++)
          s[i][j] += c.s[i][j];
    return *this;
  }

  mat operator - (const mat c){
    for(int i=0;i<N;i++)
      for(int j=0;j<N;j++)
          s[i][j] -= c.s[i][j];
    return *this;
  }

  void print(){
    for(int i=0;i<N;i++){
      for(int j=0;j<N;j++)
        cout<<s[i][j]<<' ';
      cout<<endl;
    }
  }

  mat operator ^ (ll k){
    mat res,a;
    memcpy(a.s,s,sizeof a.s);
    for(int i=0;i<N;i++) res.s[i][i]=1;
    while(k){
      if(k&1) res = res*a;
      a =a*a;
      k >>= 1;
    }
    return res; 
  }
};

上面是我的基本矩阵快速幂模板，其实矩阵快速幂难的不是你怎么写，难的是你矩阵怎么构造。

矩阵的构造，就是找递推关系。

要把需要用到的递推关系包含操作矩阵上去。找到合适的初始向量和合适的操作矩阵，你基本就可以完成题目了。

然后还要注意一下相关的优化，就是矩阵模板的一些多余的操作重载应该注释掉，实现一个较简单的结构体；优化取模，优化操作边界

附上hust自己给自己挂的一套题，自己其实也没做完，才做了5道。

矩阵套题