矩阵乘法&&矩阵快速幂

int add(long long &x,long long y)
{
    x+=y;
    x%=MOD;
}
Matrix mutil(Matrix x,Matrix y)
{
    Matrix z;
    memset(z.a,0,sizeof(z.a));
    for(int i=0;i<X;i++)
    {
        for(int k=0;k<X;k++)
        {
            if(x.a[i][k]==0) continue;
            for(int j=0;j<X;j++)
            {
                add(z.a[i][j],x.a[i][k]*y.a[k][j]);
            }
        }
    }
    return z;
}
Matrix Pow(Matrix x,int n)
{
    Matrix y;
  memset(y.a,0,sizeof(y.a)); for(int i=0;i<X;i++) y.a[i][i]=1; while(n) { if(n&1) y=mutil(x,y); x=mutil(x,x); n=n/2; } return y; }

  

posted @ 2017-03-22 20:07  Heilce  阅读(177)  评论(0编辑  收藏  举报