Python求解行列式

COUNT=0
res = []
def perm(n,begin,end):
    global COUNT
    global res
    if begin>=end:
        # print(n,end='\t')
        # print(inverse(n))
        res.append((n.copy(),inverse(n)))
        # print("res:{}".format(res))
        # print()
        COUNT +=1
    else:
        i=end-1
        for num in range(begin,end):
            n[num],n[i]=n[i],n[num]
            perm(n,begin,end-1)
            n[num],n[i]=n[i],n[num]
 
def inverse(arr):
    count = 0
    for i,n in enumerate(arr):
        for j in arr[:i]:
            if j>n:
                count+=1
    return count


def det(matrix):
    length = len(matrix)
    perm(list(range(length)),0,length)
    Sum = 0
    for series,num in res:
        product = 1
        for i,mark in enumerate(series):
            product *= matrix[i][mark]
        #     print("{}*".format(matrix[i][mark]),end='')
        # print("={}".format(product))
        Sum += ((-1)**num)*product
        # print('Sum:{}'.format(Sum))
    return Sum

def main():
    n=[[1,1,1],[1,1,1],[6,4,8]]
    print(det(n))
    # print(res)

if __name__ == '__main__':
    main()

 

评论 1
添加红包

请填写红包祝福语或标题

红包个数最小为10个

红包金额最低5元

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

抵扣说明:

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

余额充值