矩阵模板

//矩阵大小上限
const int SIZ=100;
int MOD=100;

//矩阵大小为n*m，初始化全部为0
struct mat
{
    int n,m;
    int ar[SIZ][SIZ];
    mat()
    {
        memset(ar,0,sizeof(ar));
        n=m=SIZ;
    };
};

//矩阵乘法
mat operator *(mat a,mat b)
{
    mat c;
    c=mat();
    c.n=a.n;
    c.m=b.m;
    for(int i=1;i<=a.n;i++)
        for(int j=1;j<=b.n;j++)
            if(a.ar[i][j]!=0)
            for(int k=1;k<=a.m;k++)
            {
                c.ar[i][k]+=(a.ar[i][j]*b.ar[j][k])%MOD;
                c.ar[i][k]%=MOD;
            }
    return c;
}

//矩阵加法
mat operator +(mat a,mat b)
{
    mat c;
    c=mat();
    c.n=a.n;
    c.m=a.m;
    for(int i=1;i<=a.n;i++)
        for(int j=1;j<a.m;j++)
            c.ar[i][j]=a.ar[i][j]+b.ar[i][j];
    return c;
}

//矩阵快速幂
mat operator ^(mat a,int k)
{
    mat c;
    c=mat();
    c.n=a.n;
    c.m=a.m;
    for(int i=1;i<=a.n;i++)
        c.ar[i][i]=1;
    while(k)
    {
        if(k&1)
            c=c*a;
        a=a*a;
        k/=2;
    }
    return c;
}