矩阵的逆

给一个矩阵,求它的逆。

将原矩阵A和一个单位矩阵E作成大矩阵(A, E),用初等行变换将大矩阵中的A变为E,则会得到(E, A-1)的形式。

时间复杂度:O(n3)

inline vector<double> operator * (vector<double> a, double b)
{
	int N = a.size();
	vector<double> res(N, 0);
	for(int i = 0; i < N; ++i)
		res[i] = a[i] * b;
	return res;
}

inline vector<double> operator - (vector<double> a, vector<double> b)
{
	int N = a.size();
	vector<double> res(N, 0);
	for(int i = 0; i < N; ++i)
		res[i] = a[i] - b[i];
	return res;
}

inline void inverse(vector<double> A[], vector<double> C[], int N)
{
	for(int i = 0; i < N; ++i)
		C[i] = vector<double>(N, 0);
	for(int i = 0; i < N; ++i)
		C[i][j] = 1;
	for(int i = 0; i < N; ++i)
	{
		for(int j = i; j < N; ++j)
			if(fabs(A[i][j]) > 0)
			{
				swap(A[i], A[j]);
				swap(C[i], C[j]);
				break;
			}
		C[i] = C[i] * (1 / A[i][i]);
		A[i] = A[i] * (1 / A[i][i]);
		for(int j = 0; j < N; ++j)
			if(j != i && fabs(A[i][j] > 0))
			{
				C[j] = C[j] - C[i] * A[j][i];
				A[j] = A[j] - A[i] * A[j][i];
			}
	}
}
评论 2
添加红包

请填写红包祝福语或标题

红包个数最小为10个

红包金额最低5元

当前余额3.43前往充值 >
需支付:10.00
成就一亿技术人!
领取后你会自动成为博主和红包主的粉丝 规则
hope_wisdom
发出的红包

打赏作者

水蛙菌

你的鼓励将是我创作的最大动力

¥1 ¥2 ¥4 ¥6 ¥10 ¥20
扫码支付:¥1
获取中
扫码支付

您的余额不足,请更换扫码支付或充值

打赏作者

实付
使用余额支付
点击重新获取
扫码支付
钱包余额 0

抵扣说明:

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

余额充值