0
点赞
收藏
分享

微信扫一扫

【C++】求给定阶数的行列式值(定义法)

尤克乔乔 2022-03-25 阅读 116
c++

窗口显示“矩阵阶数:”,依次输入矩阵的值,输入完成后得到矩阵的行列式值。

#include<iostream>
using namespace std;
int n;
int* p = new int[n+1];//用于存放全排列结果
double sum = 0;
int* po_g = new int[n+1];//存放自然数数列
double* po = new double[n*n+1];//用于获取用户输入的信息
bool inv(const int*, int);//求逆序数的函数
void sumup(int* p, double* po)//求和函数
{
    int yes=1;
    bool in=inv(p,n);
    if(in) yes=-1;
    double pus=yes;
    for (int i = 0; i < n; i++)
    {
        pus = *(po + n * i + *(p + i))*pus;
    
    }
    sum = pus+sum;
}
void swap(int* a, int i,int p)//全排列的交换函数
{
    int temp = *(a+i);
    *(a+i) = *(a+p);
    *(a+p) = temp;
}
void sw(int* p1, int* p2)//逆序数用的交换函数
{
    int temp = *p1;
    *p1 = *p2;
    *p2 = temp;
}
bool inv(const int *p,int size)
{
    int*p_t=new int[n+1];
    for(int i=0;i<n;i++)
    {
        *(p_t+i)=*(p+i);
    }
    int count = 0;
    for (int i = 0; i < size; i++)
    {
        for (int j = 0; j < size-i-1; j++)
        {
            if (*(p_t+j) > *(p_t+j+1)) {
                sw(p_t + j, p_t + j + 1);
                count++;
            }
        }
    }
    if(p_t=NULL)    delete[] p_t;
    if(count%2==0)return false;
    if(count%2==1) return true;
}
void getArry(int a[], int n) {
    for (int i = 0; i < n; i++) {
        *(p+i)=a[i];
    }
}
void perm(int a[], int pp, int q) {//全排列算法内嵌行列式乘法
    if ( pp== q) {
        getArry(a, q + 1);
        sumup(p,po);
    }
    else {
        for (int i = pp; i <= q; i++) {
            swap(a, i, pp);
            perm(a, pp + 1, q);
            swap(a, i, pp);
        }
    }
}
int main()
{
    cout << "矩阵阶数:";
    cin >> n;
    for (int i = 0; i < n; i++)
    {
        for (int j = 0; j < n; j++)
        {
            cout << "输入第" << i + 1 << "行,第" << j + 1 << "列:";
            cin >> *(po + i * n + j);
        }
    }
    for (int i = 0; i < n; i++)
    {
        *(po_g + i) = i;
    }
    perm(po_g, 0, n - 1);
    cout << "行列式值为" << sum << endl;
    if(p == NULL)delete[] p;
    if (po_g == NULL)delete[] po_g;
    if (po == NULL)delete[] po;
    system("pause");
    return 0;
}

举报

相关推荐

0 条评论