矩阵乘法算法时间复杂度

Question

我想出了这个矩阵乘法的算法。我在某处读到矩阵乘法的时间复杂度为 o(n^2)。 但我认为我的这个算法会给出o(n^3)。 我不知道如何计算嵌套循环的时间复杂度。所以请纠正我。

for i=1 to n
   for j=1 to n    
     c[i][j]=0
     for k=1 to n
         c[i][j] = c[i][j]+a[i][k]*b[k][j]

Answer 1

使用线性代数，存在比朴素 O(n³) 实现更好复杂度的算法。 Solvay Strassen 算法通过减少复杂度实现了 O(n^2.807)每个 2x2 子矩阵所需的乘法次数从 8 到 7。

已知最快的矩阵乘法算法是 Coppersmith-Winograd 算法，复杂度为 O(n^{2.3737）。除非矩阵很大，否则这些算法不会导致计算时间的巨大差异。在实践中，使用并行算法进行矩阵乘法更容易、更快捷。}

Answer 2

简单的算法，即按照注释中所述进行更正后得到的算法，是 O(n^3)。

确实存在可以在一定程度上减少这种情况的算法，但您不太可能找到 O(n^2) 实现。我相信最有效的实施问题仍然悬而未决。

有关详细信息，请参阅关于矩阵乘法的维基百科文章。

Answer 3

将 m×n 矩阵乘以 n×p 矩阵的标准方法的复杂度为 O(mnp)。如果所有这些对你来说都是“n”，那么它是 O(n^3)，而不是 O(n^2)。编辑：一般情况下不会是 O(n^2) 。但是对于特定类型的矩阵有更快的算法 - 如果您了解更多，您可能会做得更好。

Answer 4

在矩阵乘法中，我们使用了 3 个 for 循环，因为每个 for 循环的执行都需要时间复杂度 O(n)。因此，对于三个循环，它变为O(n^3)

Answer 5

我最近在大学作业中遇到了一个矩阵乘法问题，这就是我在 O(n^2) 中解决它的方法。

import java.util.Scanner;

public class q10 {
public static int[][] multiplyMatrices(int[][] A, int[][] B) {
    int ra = A.length; // rows in A
    int ca = A[0].length; // columns in A

    int rb = B.length; // rows in B
    int cb = B[0].length; // columns in B

    // if columns of A is not equal to rows of B, then the two matrices,
    // cannot be multiplied.
    if (ca != rb) {
        System.out.println("Incorrect order, multiplication cannot be performed");
        return A;
    } else {
        // AB is the product of A and B, and it will have rows,
        // equal to rown in A and columns equal to columns in B
        int[][] AB = new int[ra][cb];

        int k = 0; // column number of matrix B, while multiplying

        int entry; // = Aij, value in ith row and at jth index

        for (int i = 0; i < A.length; i++) {
            entry = 0;
            k = 0;

            for (int j = 0; j < A[i].length; j++) {
                // to evaluate a new Aij, clear the earlier entry
                if (j == 0) {
                    entry = 0;
                }

                int currA = A[i][j]; // number selected in matrix A
                int currB = B[j][k]; // number selected in matrix B

                entry += currA * currB; // adding to the current entry

                // if we are done with all the columns for this entry,
                // reset the loop for next one.
                if (j + 1 == ca) {
                    j = -1;
                    // put the evaluated value at its position
                    AB[i][k] = entry;

                    // increase the column number of matrix B as we are done with this one
                    k++;
                }

                // if this row is done break this loop,
                // move to next row.
                if (k == cb) {
                    j = A[i].length;
                }
            }
        }
        return AB;
    }

}

@SuppressWarnings({ "resource" })
public static void main(String[] args) {
    Scanner ip = new Scanner(System.in);

    System.out.println("Input order of first matrix (r x c):");
    int ra = ip.nextInt();
    int ca = ip.nextInt();

    System.out.println("Input order of second matrix (r x c):");
    int rb = ip.nextInt();
    int cb = ip.nextInt();

    int[][] A = new int[ra][ca];
    int[][] B = new int[rb][cb];

    System.out.println("Enter values in first matrix:");
    for (int i = 0; i < ra; i++) {
        for (int j = 0; j < ca; j++) {
            A[i][j] = ip.nextInt();
        }
    }

    System.out.println("Enter values in second matrix:");
    for (int i = 0; i < rb; i++) {
        for (int j = 0; j < cb; j++) {
            B[i][j] = ip.nextInt();
        }
    }

    int[][] AB = multiplyMatrices(A, B);

    System.out.println("The product of first and second matrix is:");
    for (int i = 0; i < AB.length; i++) {
        for (int j = 0; j < AB[i].length; j++) {
            System.out.print(AB[i][j] + " ");
        }
        System.out.println();
    }
}

}