线性代数,行列式求值。

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

using Vector = vector<vector<int>>;

class D
{
protected:

	int num;
	int level;
	Vector arr;

public:
	D():num(1), level(1) {}
	int(*x)[5] = new int[120][5];
	void init()
	{
		cout << "请输入行列式长度:";
		cin >> num;
		arr.resize(num);
		system("cls");
		int temp;
		for (int i = 0; i < num; i++)
		{
			for (int j = 0; j < num; j++)
			{
				cin >> temp;
				arr[i].push_back(temp);
			}
		}
	}

	int getD()
	{
		int result = 0;
		for (int i = num; i > 1; i--)
		{
			level *= i;
		}
		int* y = new int[num]();
		vector<int> m;
		for (int i = 0; i < num; i++)
		{
			*(y + i) = i;
			m.push_back(i);
		}

		getX(m, 0, 0, level);

		for (int i = 0; i < level; i++)
		{
			result += multip(arr, x[i]);
		}
		return result;
	}
	int ab(int a[])//判断正负
	{
		int data = 0;
		for (int i = 0; i < num; i++)
		{
			for (int j = i + 1; j < num; j++)
			{
				if (a[i] > a[j])
					data++;
			}
		}
		if (data % 2 == 0)
			return 1;
		else
			return -1;
	}

	int multip(Vector vec, int a[])
	{
		int data = 1;
		for (int i = 0; i < num; i++)
		{
			data *= vec[i][a[i]];
		}
		data *= ab(a);
		return data;
	}
	void getX(std::vector<int> vec, int a, int b, int size)//求列标的所有排列组合
	{
		int size_t = 1;
		if (num - b != 0)
			size_t = size / (num - b);
		for (int i = 0; i < num - b; i++)
		{
			vector<int> temp = vec;
			for (int j = 0; j < size_t; j++)
			{
				x[a + j][b] = temp[i];
			}
			temp.erase(temp.begin() + i);
			getX(temp, a, b + 1, size_t);
			a += size_t;
		}
	}
};


int main()
{
	D d;
	d.init();
	cout << "D=" << d.getD() << endl;

	system("pause");
	return 0;
}
  • 1
    点赞
  • 0
    收藏
    觉得还不错? 一键收藏
  • 1
    评论

“相关推荐”对你有帮助么?

  • 非常没帮助
  • 没帮助
  • 一般
  • 有帮助
  • 非常有帮助
提交
评论 1
添加红包

请填写红包祝福语或标题

红包个数最小为10个

红包金额最低5元

当前余额3.43前往充值 >
需支付:10.00
成就一亿技术人!
领取后你会自动成为博主和红包主的粉丝 规则
hope_wisdom
发出的红包
实付
使用余额支付
点击重新获取
扫码支付
钱包余额 0

抵扣说明:

1.余额是钱包充值的虚拟货币,按照1:1的比例进行支付金额的抵扣。
2.余额无法直接购买下载,可以购买VIP、付费专栏及课程。

余额充值