Произведение матриц

def mult_matr(a1,a2):
    r1=len(a1)
    c1=len(a1[0])
    r2=len(a2)
    c2=len(a2[0])
    if c1==r2:
        res=[]
        for i in range(r1):
            row=[]
            for j in range(c2):
                s=0
                for k in range(c1):
                    s+=a1[i][k]*a2[k][j]
                row.append(s)
            res.append(row)    
        return res
    else:
        return None
        
x=[[1,2],[3,4]]
y=[[5,6],[7,8]]
 
q=mult_matr(x,y)
 
if q is None:
    print("Умножение невозможно")
else:
    print(q)

import random
 
class Matrix:
    def __init__(self, rows, cols):
        self.rows = rows
        self.cols = cols
        self.data = [[0] * cols for _ in range(rows)]
 
    @classmethod
    def from_table(cls, table):
        rows = len(table)
        cols = len(table[0])
        matrix = cls(rows, cols)
        matrix.data = table
        return matrix
 
    @classmethod
    def random(cls, rows, cols):
        matrix = cls(rows, cols)
        for i in range(rows):
            for j in range(cols):
                matrix.data[i][j] = random.randint(1, 10)
        return matrix
 
    def __str__(self):
        output = ""
        for row in self.data:
            output += " ".join(str(element) for element in row)
            output += "\n"
        return output
 
    def __add__(self, other):
        if self.rows != other.rows or self.cols != other.cols:
            raise ValueError("Matrices must have the same dimensions")
        result = Matrix(self.rows, self.cols)
        for i in range(self.rows):
            for j in range(self.cols):
                result.data[i][j] = self.data[i][j] + other.data[i][j]
        return result
 
    def __mul__(self, other):
        if isinstance(other, Matrix):
            if self.cols != other.rows:
                raise ValueError("Number of columns in the first matrix must match the number of rows in the second matrix")
            result = Matrix(self.rows, other.cols)
            for i in range(self.rows):
                for j in range(other.cols):
                    for k in range(self.cols):
                        result.data[i][j] += self.data[i][k] * other.data[k][j]
        elif isinstance(other, (int, float)):
            result = Matrix(self.rows, self.cols)
            for i in range(self.rows):
                for j in range(self.cols):
                    result.data[i][j] = self.data[i][j] * other
        else:
            raise TypeError("Multiplication is only supported between matrices and numbers")
        return result
 
 
# Создание матрицы 2x3 и заполнение случайными числами
matrix1 = Matrix.random(2, 3)
print(matrix1)
 
# Создание матрицы 2x3 и заполнение из готовой таблицы
table = [[1, 2, 3], [4, 5, 6]]
matrix2 = Matrix.from_table(table)
print(matrix2)
 
# Сложение матриц
sum_matrix = matrix1 + matrix2
print(sum_matrix)
 
# Умножение матрицы на число
scalar = 2
mul_scalar = matrix1 * scalar
print(mul_scalar)
 
 
# Умножение матриц
matrix3 = Matrix.random(3, 2)
product_matrix = matrix2 * matrix3
print(product_matrix)

import numpy as np
n, m = 2, 3  # Размеры матриц
a = np.random.randint(0, 9, [n,m])
b = np.random.randint(0, 9, [m,n])
 
print(a)
print(b)
 
if m == n:
    print()
    print('Почленное (покоординатное) перемножение\n', a * b)
    
print('Операция умножения определена для двух матриц, таких что число \nстолбцов первой равно числу строк второй.')
print('Матричное умножение (функция dot)\n',np.dot(a, b))
# Различие с dot возникают при умножении матриц размерностью 3D и выше
print('Матричное умножение (функция matmul)\n', np.matmul(a, b))

from random import randint
from math import prod
 
def mult_matr(A,B):
    return [[sum(map(prod, zip(A[i], [row[j] for row in B]))) for j in range(len(B[0]))] for i in range(len(A))]
 
a = randint(1,3)
b = randint(1,3)
c = randint(1,3)
matr_A = [[randint(1,10) for _ in range(b)] for _ in range(a)]
matr_B = [[randint(1,10) for _ in range(c)] for _ in range(b)]
print(*matr_A, sep='\n')
print()
print(*matr_B, sep='\n')
print()
print(*mult_matr(matr_A, matr_B), sep='\n')