窗口显示“矩阵阶数:”,依次输入矩阵的值,输入完成后得到矩阵的行列式值。
#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;
}