c求逆矩阵的代码_二维矩阵求逆

c求逆矩阵的代码_二维矩阵求逆刷石油大刷自闭了,随便写点乱七八糟的东西放松一下。。逆矩阵=伴随矩阵除以对应行列式的值,所以分别把这两个求出来就解决了,另外为了避免产生分数,就不做进一步计算了。效果图如下。至于代码。。就是把上次写的行列式求值的博客稍微改了一下,emmm。#include<stdio.h>#include<iostream>#pragmawarning(disa…

大家好,又见面了,我是你们的朋友全栈君。如果您正在找激活码,请点击查看最新教程,关注关注公众号 “全栈程序员社区” 获取激活教程,可能之前旧版本教程已经失效.最新Idea2022.1教程亲测有效,一键激活。

Jetbrains全系列IDE使用 1年只要46元 售后保障 童叟无欺

 

c求逆矩阵的代码_二维矩阵求逆

#include<stdio.h>
#include<iostream>
#pragma warning(disable : 4996)

using namespace std;

const int maxn=105;

int juzhen[maxn][maxn],n;
int ans[maxn][maxn];
int book[maxn];
int nx(int n)
{
	int sum = 0;
	for (int i = 0; i < n; i++)
	{
		for (int j = 0; j < i; j++)
		{
			if (book[j] > book[i]) sum++;
		}
	}
	if (sum & 1) return -1;
	return 1;
}
int yuzi_ans;
void dfs(int yuzi[maxn][maxn], int sum, int pos,int n) {
	if (sum == n ) {
		int s = 1;
		for (int i = 0; i < n; i++) s *= yuzi[i][book[i]];
		yuzi_ans += nx(n) * s;
		return;
	}
	for (int i = 0; i < n; i++)
	{
		if (pos >> i & 1) continue;
		pos = pos | 1 << i;
		book[sum] = i;
		dfs(yuzi,sum + 1,pos,n);
		pos = pos ^ 1 << i;
	}
}
int f(int yuzi[maxn][maxn],int n) {
	yuzi_ans = 0;
	dfs(yuzi, 0, 0,n);
	return yuzi_ans;
}
void nijuzhen() {
	int yuzi[maxn][maxn];
	int bansui[maxn][maxn];
	for (int i = 0; i < n; i++) {
		for (int j = 0; j < n; j++) {
			int p = 0;
			for (int k = 0; k < n - 1; k++) {
				if (p == i) {
					p++;
				}
				int q = 0;
				for (int r = 0; r < n - 1; r++) {
					if (q == j) {
						q++;
					}
					yuzi[k][r] = juzhen[p][q++];
				}
				p++;
			}
			int flag = 1;
			if ((i + j) & 1) flag = -1;
			bansui[j][i] = f(yuzi,n-1)*flag;
		}
	}
	printf("伴随矩阵为:\n");
	for (int i = 0; i < n; i++) {
		for (int j = 0; j < n; j++) {
			printf("%d ", bansui[i][j]);
		}
		printf("\n");
	}
	printf("原矩阵对应的行列式的值为:\n");
	printf("%d", f(juzhen, n));
}

int main()
{
	printf("请输入矩阵阶数\n");
	scanf("%d", &n);
	for (int i = 0; i < n; i++)
		for (int j = 0; j < n; j++)
			scanf("%d", &juzhen[i][j]);
	nijuzhen();
	return 0;
}

 

版权声明:本文内容由互联网用户自发贡献,该文观点仅代表作者本人。本站仅提供信息存储空间服务,不拥有所有权,不承担相关法律责任。如发现本站有涉嫌侵权/违法违规的内容, 请发送邮件至 举报,一经查实,本站将立刻删除。

发布者:全栈程序员-用户IM,转载请注明出处:https://javaforall.cn/171665.html原文链接:https://javaforall.cn

【正版授权,激活自己账号】: Jetbrains全家桶Ide使用,1年售后保障,每天仅需1毛

【官方授权 正版激活】: 官方授权 正版激活 支持Jetbrains家族下所有IDE 使用个人JB账号...

(0)


相关推荐

发表回复

您的电子邮箱地址不会被公开。

关注全栈程序员社区公众号