矩阵模版(例题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;
}


评论
添加红包

请填写红包祝福语或标题

红包个数最小为10个

红包金额最低5元

当前余额3.43前往充值 >
需支付:10.00
成就一亿技术人!
领取后你会自动成为博主和红包主的粉丝 规则
hope_wisdom
发出的红包
实付
使用余额支付
点击重新获取
扫码支付
钱包余额 0

抵扣说明:

1.余额是钱包充值的虚拟货币,按照1:1的比例进行支付金额的抵扣。
2.余额无法直接购买下载,可以购买VIP、付费专栏及课程。

余额充值