矩阵运算:
NumPy系统是Python的一种开源的数值计算扩展。这种工具可用来存储和处理大型矩阵,比Python自身的嵌套列表(nested list structure)结构要高效的多(该结构也可以用来表示矩阵(matrix))
import numpy as np
from numpy.linalg import * // 矩阵运算
行列式:
a=np.mat([[1,23,3],[1,6,4],[13,12,14]])
a_det=det(a)
求逆:
a_inv=a.I
伴随矩阵:
a_adju=a_det*a_inv
求两个数的最大公因数:
def gcd(a,b):
while a!=0:
a,b = b%a,a
return b
求a模m的逆:
def exgcd(a,m):
if gcd(a,m)!=1:
return None
u1,u2,u3 = 1,0,a
v1,v2,v3 = 0,1,m
while v3!=0:
q = u3//v3
v1,v2,v3,u1,u2,u3 = (u1-q*v1),(u2-q*v2),(u3-q*v3),v1,v2,v3
return u1%m
乘法:
C=a*b
加法:
C=a+b
转置:
a_tra=a.T