c++使用vector求矩阵的A(行列式)

最新推荐文章于 2022-04-09 17:19:09 发布

焦志坤

最新推荐文章于 2022-04-09 17:19:09 发布

阅读量577

点赞数

分类专栏： C++学习之路

本文链接：https://blog.csdn.net/cute_boy_/article/details/119573160

版权

C++学习之路专栏收录该内容

17 篇文章 5 订阅

订阅专栏

c++使用vector求矩阵的A(行列式)

直接上代码

#include <vector>
#include <iostream>
using namespace std;


bool Iseven(int num)
{
    //用位与运算来判断奇偶（最快的判断奇偶的方法）
    return ((num & 1) == 0);
}

bool PowerIsPosition(vector<double>& vec)
{
    //count即为逆序数，初始化为0
    int count = 0;
    for (int i = 0; i < (int)vec.size(); i++)
    {
        for (int j = i + 1; j < (int)vec.size(); j++)
        {
            if (vec[i] > vec[j])
            {
                count += 1;
            }
        }
    }
    return (Iseven(count));
}

void swap(vector<double>& vec, int i, int j)
{
    int temp = vec[i];
    vec[i] = vec[j];
    vec[j] = temp;
}

void Perm(vector<double>& vec, vector<vector<double> >& vec_seq, int current_index)
{
    if (current_index == ((int)vec.size() - 1))
    {
        vec_seq.push_back(vec);
    }
    else
    {
        for (int i = current_index; i < (int)vec.size(); i++)
        {
            swap(vec, i, current_index);
            Perm(vec, vec_seq, current_index + 1);
            swap(vec, i, current_index);
        }
    }
}

vector<double> inivec(int n)
{
    vector<double> vec;
    for (int i = 0; i < n; i++)
        vec.push_back(i);
    return vec;
}

double calculate(vector<vector<double>> array)
{
    //计算行列式
    int n = array.size();
    vector<vector<double> > vec_que;
    vector<double> vec = inivec(n);
    vector<double> vec_elem;
    Perm(vec, vec_que, 0);
    //最终结果，初始化为0
    double result = 0;
    //依次为vec_que中取出行列式
    for (int i = 0; i < (int)vec_que.size(); i++)
    {
        vec_elem = vec_que[i];
        //mi即为前面(-1)的n次幂，最后结果为-1或者1
        double mi = PowerIsPosition(vec_elem) ? 1 : (-1);
        double temp = mi;
        //row号初始化为0之后依次加1
        int row = 0;
        //col号依次从vec_elem中取出
        for (int j = 0; j < (int)vec_elem.size(); j++)
        {
            int col = vec_elem[j];
            temp *= array[row++][col];
        }
        result += temp;
    }
    return result;
}


int main() {
    vector<vector<double>> nums1 = {
        {1, 2, 3},
        {4, 15, 6},
        {7, 8, 8}
    };
    double w = calculate(nums1);
    cout << w << endl;
}