c++ - 在处理多维数组时，为什么在计数器还在进行的时候循环里面的语句没有执行？

Question

我正在研究消除增广矩阵的问题，以便使用线性代数求解线性方程。

我正在使用二维数组来处理这个问题。但是我遇到了一个问题。我用来消除的函数是一个循环，但是在计数器还在进行的时候，里面的语句并没有被执行。因此整个循环无法正确完成。

这是我的代码:

#include <iostream>
#include <vector>

using namespace std;


vector <int> Row;
vector <vector <int> > M;

void swap(vector <vector <int> > &M, int i, vector <vector <int> > &N, int j);
void reduce(int rank, vector <vector <int> > &M);
int search_nonzero(int rank, vector <vector <int> > &M, int j);

int main(){
    cout << "plz enter the size of linear equation:";
    int rank;

    cin >> rank;
    int value = 0;;
    cout << "plz enter the value" << endl;
    for (int initial = 0; initial <= rank; ++initial){
        Row.push_back(value);
    }

    for (int i = 0; i <= rank-1; ++i){
        for (int j = 0; j <= rank; ++j){
            cout << "m[" << i+1 <<"][" << j+1 << "]:";
            cin >> value;
            Row[j] = value;
        }
        M.push_back(Row);
    }

    reduce(rank, M);

    cout << "the operated matrix is: " << endl;
    for (int i = 0; i <= rank-1; ++i){
        for (int j = 0; j <= rank; ++j){
            cout << M[i][j] << "   ";
        }
    cout << endl;
    }

}

---

void swap(vector <vector <int> > & M, int i, vector <vector <int> > &N, int j){
    vector <vector <int> > temp;
    temp.push_back(M[i]);
    M[i] = N[j];
    N[j] = temp[0];
}

int search_nonzero(int rank, vector <vector <int> > & M, int j){
    for (int s = rank-1; s >= 1; s--){
        if (M[s][j] != 0)
            return s;
    }
}

---

void reduce(int rank, vector <vector <int> > & M){
    for (int s = rank-1; s > 0; s--){
        for (int i = rank-s; i < rank; ++i){
            for (int j = rank-s-1; j <= rank; ++j){
                if (M[rank-s-1][rank-s-1] == 0){
                    swap(M, i-1, M, search_nonzero(rank, M, j));
                }
                else
                M[i][j] = M[i][j] - M[rank-s-1][j]*M[i][rank-s-1]/M[rank-s-1][rank-s-1];
            }
        }    
    }
}

运行是:

    plz enter the size of linear equation: 3
    plz enter the value
    m[1][1]= 4
    m[1][2]= 5
    m[1][3]= 6
    m[1][4]= 7
    m[2][1]= 5
    m[2][2]= 6
    m[2][3]= 4
    m[2][4]= 8
    m[3][1]= 1
    m[3][2]= 2
    m[3][3]= 5
    m[3][4]= 9
    the operated matrix is:
    4   5   6   7
    0   6   4   8
    0   0   5   9

我是用visual studio一步步调试的。问题出在 void 函数 reduce() 中。当j增加时，句子“M[i][j]=M[i][j]-M[rank-s-1][j]*M[i][rank-s-1]/M[rank - s - 1 ] [等级 - s - 1 ] ;"只是在j=0时执行了一次，而在j=1,j=2等时没有执行。

那么我的代码有什么问题呢？我猜这是由于引用调用问题，其中j=2时这个多维数组的值等没有改变，尽管计数器j改变了。

非常感谢!

Answer 1

错误在于，当您显然期望浮点运算时，您却在进行整数运算。将 vector 更改为 vector <vector <double> > .