矩阵模版（例题3070、3744）

struct matrix
{
double A[NMax][NMax];
matrix operator*(const matrix & m)
{
matrix ret={0,0,0,0};
for(int i=0;i<NMax;i++)
for(int k=0;k<NMax;k++)
if(A[i][k]>0.0)
for(int j=0;j<NMax;j++)
ret.A[i][j]+=A[i][k]*m.A[k][j];
return ret;
}
}def={1,0,0,1};
matrix pow(matrix a,int n)
{
matrix tmp=a,ret=def;
while(n!=0)
{
if((n&1)!=0) ret=ret*tmp;
tmp=tmp*tmp;
n>>=1;
}
return ret;
}