[Add] Kronecker Product

DIRECTORY.md

@@ -970,6 +970,7 @@
   * [Count Paths](matrix/count_paths.py)
   * [Cramers Rule 2X2](matrix/cramers_rule_2x2.py)
   * [Inverse Of Matrix](matrix/inverse_of_matrix.py)
+  * [Kronecker Product](matrix/kronecker_product.py)
   * [Largest Square Area In Matrix](matrix/largest_square_area_in_matrix.py)
   * [Matrix Based Game](matrix/matrix_based_game.py)
   * [Matrix Class](matrix/matrix_class.py)
@@ -1510,6 +1511,7 @@
   * [Autocomplete Using Trie](strings/autocomplete_using_trie.py)
   * [Barcode Validator](strings/barcode_validator.py)
   * [Bitap String Match](strings/bitap_string_match.py)
+  * [Booths Algorithm](strings/booths_algorithm.py)
   * [Boyer Moore Horspool](strings/boyer_moore_horspool.py)
   * [Boyer Moore Search](strings/boyer_moore_search.py)
   * [Camel Case To Snake Case](strings/camel_case_to_snake_case.py)
@@ -1536,6 +1538,7 @@
   * [Jaro Winkler](strings/jaro_winkler.py)
   * [Join](strings/join.py)
   * [Knuth Morris Pratt](strings/knuth_morris_pratt.py)
+  * [Largest Smallest Words](strings/largest_smallest_words.py)
   * [Levenshtein Distance](strings/levenshtein_distance.py)
   * [Lower](strings/lower.py)
   * [Manacher](strings/manacher.py)

matrix/kronecker_product.py

@@ -0,0 +1,75 @@
+# @Author  : jay2219
+# @File    : kronecker_product.py
+# @Date    : 13/10/2024
+
+"""
+Perform Kronecker product of two matrices.
+https://en.wikipedia.org/wiki/Kronecker_product
+"""
+
+
+def is_2d(matrix: list[list[int]]) -> bool:
+    """
+    >>> is_2d([])
+    True
+    >>> is_2d([1, 2])
+    False
+    >>> is_2d([[1, 2], [3, 4]])
+    True
+    """
+
+    return all(isinstance(matrix, list) and (isinstance(i, list) for i in matrix))
+
+
+def kronecker_product(
+    matrix_a: list[list[int]], matrix_b: list[list[int]]
+) -> list[list[int]]:
+    """
+    :param matrix_a: A 2-D Matrix with dimension m x n
+    :param matrix_b: Another 2-D Matrix with dimension p x q
+    :return: Result of matrix_a ⊗ matrix_b
+    :raises ValueError: If the matrices are not 2-D.
+
+    >>> kronecker_product([[1, 2]], [[5, 6], [7, 8]])
+    [[5, 6, 10, 12], [7, 8, 14, 16]]
+
+    >>> kronecker_product([[1, 2], [4, 5]], [[5, 6], [7, 8]])
+    [[5, 6, 10, 12], [7, 8, 14, 16], [20, 24, 25, 30], [28, 32, 35, 40]]
+
+    >>> kronecker_product([1, 2], [[5, 6], [7, 8]])
+    Traceback (most recent call last):
+        ...
+    ValueError: Input matrices must be 2-D.
+    """
+
+    # Check if the input matrices are valid
+    if not all((is_2d(matrix_a), is_2d(matrix_b))):
+        raise ValueError("Input matrices must be 2-D.")
+
+    if not matrix_a or not matrix_b:
+        return []
+
+    rows_matrix_a, cols_matrix_a = len(matrix_a), len(matrix_a[0])
+    rows_matrix_b, cols_matrix_b = len(matrix_b), len(matrix_b[0])
+
+    # Resultant matrix dimensions
+    result = [
+        [0] * (cols_matrix_a * cols_matrix_b)
+        for _ in range(rows_matrix_a * rows_matrix_b)
+    ]
+
+    for r_index_a in range(rows_matrix_a):
+        for c_index_a in range(cols_matrix_a):
+            for r_index_b in range(rows_matrix_b):
+                for c_index_b in range(cols_matrix_b):
+                    result[r_index_a * rows_matrix_b + r_index_b][
+                        c_index_a * cols_matrix_b + c_index_b
+                    ] = matrix_a[r_index_a][c_index_a] * matrix_b[r_index_b][c_index_b]
+
+    return result
+
+
+if __name__ == "__main__":
+    import doctest
+
+    doctest.testmod()