為什么我的 C 程序應該輸出一個矩陣的 n 次方，輸出我的矩陣的 2^n 次方？

Question

我的代碼應該接受矩陣 M 並將其提升到整數 A 的冪。但是，不知何故，我的輸出始終是 M^(2^A)。 例如，如果我想找到一個矩陣的 3 次方，我會得到它的 8 次方。

#include <stdio.h> 
#include <stdlib.h> 
#include <math.h>
  void multiply(int ** p, int pwr, int dim, int ** prod) {
    int m, i, j, k;
    /*if (n<pwr){*/
    int pos = 0;
    for (m = 0; m < pwr; m++) {
      for (i = 0; i < dim; i++) {
        for (j = 0; j < dim; j++) {
          for (k = 0; k < dim; k++) {
            pos += p[i][k] * p[k][j];
          }
          prod[i][j] = pos;
          pos = 0;
        }
      }
      for (i = 0; i < dim; i++) {
        for (j = 0; j < dim; j++) {
          p[i][j] = prod[i][j];
          prod[i][j] = 0;
        }
      }
    }
    /*n=n+1;
    multiply(prod, q, pwr, dim, prod);
    }*/
  }
int main(int argc, char * argv[]) {
  FILE * fp = fopen(argv[1], "r");
  int dim, pwr, i, j;
  fscanf(fp, "%d", & dim);
  int ** matrix;
  matrix = (int ** ) malloc(dim * sizeof(int * ));
  for (i = 0; i < dim; i++) {
    matrix[i] = (int * ) malloc(dim * sizeof(int));
  }
  int ** prod;
  prod = (int ** ) malloc(dim * sizeof(int * ));
  for (i = 0; i < dim; i++) {
    prod[i] = (int * ) malloc(dim * sizeof(int));
  }
  for (i = 0; i < dim; i++) {
    for (j = 0; j < dim; j++) {
      fscanf(fp, "%d", & matrix[i][j]);
    }
  }
  fscanf(fp, "%d", & pwr);
  if (pwr == 1) {
    for (i = 0; i < dim; i++) {
      for (j = 0; j < dim; j++) {
        printf("%d ", matrix[i][j]);
      }
      printf("\n");
    }
  } else if (pwr >= 2) {
    multiply(matrix, pwr, dim, prod);
    for (i = 0; i < dim; i++) {
      for (j = 0; j < dim; j++) {
        printf("%d ", matrix[i][j]);
      }
      printf("\n");
    }
  }
  return 0;
}

Answer 1

您將矩陣乘以自身，然后將結果存儲在原始矩陣中。 然后你再做一次。

非常正常，它被供電 8 次。 您需要的是另一個臨時矩陣，您可以在其中存儲結果並保留原始矩陣以與您的結果相乘。