二维数组或更高维数组的传递问题
今天打一条简单矩阵题竟然发现了自己对二维数组作为形参传递的漏洞…一直以为可以直接用int**
丢进函数里,不过想下这也不太科学…所以就做了各种实验…
下面是传递的方法
四种方法:
-
直接传递二维数组,但是必须标明第二维的值,因为如果只是传递arr[][], 编译器无法分配这样的数组,所以要这样传int arr[][3],然后运用的时候直接吧数组名丢进去就好了,当然你喜欢int arr[3][3]也是可以的。
-
方法二.是传递指针数组方式,即int (*a)[3]。1.0版本用了方法二,传递指针数组…一比直接传递数组方便,二又能直接用a[i][j]这种方便的记号查找特定元素。
-
方法三.是传递一维指针方法。传入int arr,使用的时候,二维:(arr + i*iCols + j) 表示 arr[i][j]。但传入前要先强制转换成一维指针,1.5版本的测试这种方法出现了问题。
-
也是可以丢int **arr进去的,看2.0版本。
题目是要做一个 求行列式 的练习
先贴1.0版本
#include <iostream>
#include <cmath>
#define MAX_N 4
int det(int (*m)[MAX_N], int n);
int main() {
int n;
int matrix[MAX_N][MAX_N];
std::cin >> n;
for (int i = 0; i < n; i++) {
for (int j = 0; j < n; j++) {
std::cin >> matrix[i][j];
}
}
std::cout << det(matrix, n) << std::endl;
}
int det(int (*m)[MAX_N], int n){
if (n > 2) {
int sum = 0;
int i, j, col;
int tM[MAX_N][MAX_N];
for (col = 0; col < n; col++) {
for (i = 0; i < n - 1; i++) {
for (j = 0; j < col; j++) {
tM[i][j] = m[i + 1][j];
}
}
for (i = 0; i < n - 1; i++) {
for (j = col + 1; j < n; j++) {
tM[i][j - 1] = m[i + 1][j];
}
}
sum += pow(-1, col) * m[0][col] * det(tM, n - 1);
}
return sum;
}
else if (n == 2) {
return (m[0][0] * m[1][1] - m[0][1] * m[1][0]);
}
else if (n == 1) {
return m[0][0];
}
else
return 0;
}
这里宏定义了一下矩阵维数便于对程序进行测试, 但n受宏定义的限制,表示不爽,决定用new动态动态创建数组。
1.5.问题版本
#include <iostream>
#include <cmath>
int det(int *m, int n);
int main() {
int n;
int **matrix;
std::cin >> n;
matrix = new int*[n];
for (int i = 0; i < n; i++) {
matrix[i] = new int[n];
for (int j = 0; j < n; j++) {
std::cin >> matrix[i][j];
}
}
std::cout << det((int*)matrix, n) << std::endl;
for (int i = 0; i < n; i++) {
delete []matrix[i];
}
delete []matrix;
}
int det(int *m, int n){
if (n > 2) {
int sum = 0;
int **tM;
int i, j, col;
tM = new int*[n - 1];
for (i = 0; i < n - 1; i++) {
tM[i] = new int[n - 1];
}
for (col = 0; col < n; col++) {
for (i = 0; i < n - 1; i++) {
for (j = 0; j < col; j++) {
tM[i][j] = *(m + (1 + i)*n + j);
}
}
for (i = 0; i < n - 1; i++) {
for (j = col + 1; j < n; j++) {
tM[i][j - 1] = *(m + (1 + i)*n + j);
}
}
sum += pow(-1, col) * *(m + col) * det((int*)tM, n - 1);
}
for (int i = 0; i < n - 1; i++) {
delete []tM[i];
}
delete []tM;
return sum;
}
else if (n == 2) {
return ((*m) * (*(m + n + 1)) - (*(m + 1)) * (*(m + n)));
}
else if (n == 1) {
return *m;
}
else
return 0;
}
这种输出方法莫名奇妙地输出指针地址…并没有达到方法三预期的效果(有待改良)。
2.0版本(n不受限制)
#include <iostream>
#include <cmath>
int det(int **m, int n);
int main() {
int n;
int **matrix;
std::cin >> n;
matrix = new int*[n];
for (int i = 0; i < n; i++) {
matrix[i] = new int[n];
for (int j = 0; j < n; j++) {
std::cin >> matrix[i][j];
}
}
std::cout << det(matrix, n) << std::endl;
for (int i = 0; i < n; i++) {
delete []matrix[i];
}
delete []matrix;
}
int det(int **m, int n) {
if (n > 2) {
int sum = 0;
int **tM;
int i, j, col;
tM = new int*[n - 1];
for (i = 0; i < n - 1; i++) {
tM[i] = new int[n - 1];
}
for (col = 0; col < n; col++) {
for (i = 0; i < n - 1; i++) {
for (j = 0; j < col; j++) {
tM[i][j] = m[i + 1][j];
}
}
for (i = 0; i < n - 1; i++) {
for (j = col + 1; j < n; j++) {
tM[i][j - 1] = m[i + 1][j];
}
}
sum += pow(-1, col) * m[0][col] * det(tM, n - 1);
}
for (int i = 0; i < n - 1; i++) {
delete []tM[i];
}
delete []tM;
return sum;
}
else if (n == 2) {
return (m[0][0] * m[1][1] - m[0][1] * m[1][0]);
}
else if (n == 1) {
return m[0][0];
}
else
return 0;
}
文章部分内容参考了:
http://blog.csdn.net/liuzhanchen1987/article/details/7712640