基于C语言实现整数行列式

于 2024-03-17 17:58:59 发布
阅读量547 收藏 9
点赞数 12
文章标签： 线性代数 矩阵
版权声明：本文为博主原创文章，遵循 CC 4.0 BY-SA 版权协议，转载请附上原文出处链接和本声明。
本文链接：https://blog.csdn.net/2301_79847688/article/details/136785390
版权

在本文章内，将会实现行列式的建立、销毁、打印、计算四个操作。

鉴于笔者技术有限，此行列式只针对整数int型，请读者自行扩充~_~。

1.行列式的建立与销毁

我们首先建立行列式的数据类型，由于行列式规模的不确定，采用动态分配方法。

typedef struct
{
	int n;
	int* p;
}determinant;

建立函数

void build_deter(determinant *d)//创建行列式
{
	printf("请输入行列数:>");//输入行列式
	scanf("%d", &d->n);

	d->p = (int*)malloc(sizeof (int)* d->n*d->n);//为行列式分配空间
	if (d->p == NULL)
	{
		printf("%s", strerror(errno));
	}

	printf("请按行输入值:>");
	for (int i = 0; i < d->n; i++)
	{
		for (int j = 0; j < d->n; j++)
		{
			scanf("%d",(d->p)+(i*d->n+j));
		}
	}
	printf("building over!\n");
}

销毁函数

void free_deter(determinant* d)//释放空间
{
	free(d->p);
	d->p = NULL;
}

2.行列式的打印

void dis_deter(determinant d)//打印行列式
{
	for (int i = 0; i < d.n; i++)
	{
		printf("|");
		for (int j = 0; j < d.n; j++)
		{
			printf("%-3d ", *((d.p) + (i * d.n + j)));
		}
		printf("|\n");
	}
	printf("printing over!\n");
}

其中由于动态分配原因，采用*((d.p) + (i * d.n + j))来模拟二维的定位。

3.行列式的计算

贴一张行列式的计算公式

可以看出，我们只需要解决两个问题，一个是找出p1p2p3……pn的全排列，并以此找到一组a进行相加，另一个是找到p1p2p3……pn的逆序数，从而确定每一组a前面的符号。

对于全排列，我在上一文章有讲实现原理，在此不加赘述，想要了解可以自行传送~_~，https://blog.csdn.net/2301_79847688/article/details/136762462

贴个代码：

int fac(int n)//计算阶乘
{
	int result=1;
	for (int i = n; i > 0; i--)
	{
		result *= i;
	}
	return result;
}

void change(int* a, int* b)//交换函数
{
	int c;
	c = *b;
	*b = *a;
	*a = c;
}

void bubble_sort(int* arrange,int n)//排序
{
	for (int i = 0; i < n - 1; i++)
	{
		for (int j = 0; j < n - 1 - i; j++)
		{
			if (arrange[j] > arrange[j + 1])
				change(arrange+j, arrange+j+1);
		}
	}
}

void tran_arrange(int* arrange, int j, int n)//改变排列
{
	if (j == 0)//自然顺序做首项
	{
		return;
	}
	if (n == 2)
	{
		change(arrange, arrange + 1);
	}
	else if(n == 3)//等于3做奇偶不同的对换
	{
		if (j % 2 == 0)
		{
			change(arrange, arrange + 1);
		}
		else if (j % 2 == 1)
		{
			change(arrange + 1, arrange + 2);
		}
	}
	else //大于3递归
	{
		if (j % fac(n - 1) == 0)//在后面数循环结束后，调换首数和后面一个数，产生新的组合
		{
			int tmp = j;
			if (tmp > fac(n))
				tmp = tmp % fac(n);//tmp要在fac(n)内，防止（1）越界

			bubble_sort(arrange + 1,n-1);//排序方便调换
			if (tmp == 0)//自然顺序做首排列，第0次
			{
				return;
			}
			if (tmp == fac(n - 1))
				change(arrange, arrange + n - 1);//调换算法，第一次换最后一位
			else
				change(arrange, arrange + n - (tmp / fac(n - 1)) + 1);//（1）之后“/”
		}
		tran_arrange(arrange+1,j,n-1);
	}
}

对于逆序数，对于排列中一个元素，他后面有多少比他小的数，那他就构成了多少逆序数，而排列的总逆序数，就是所有元素对应逆序数的和。

大体思想是这样，看看代码：

int count_re(int* arrange,int n)//统计逆序数
{
	int count = 0;
	for (int i = 0; i < n; i++)
	{
		for (int j = i+1; j < n; j++)
		{
			if (arrange[i] > arrange[j])
				count++;
		}
	}
	if (count % 2 == 1)
		return -1;
	else
		return 1;
}

这样我们就解决了所有算法问题，只需要用循环得到每一个数再相加就可以了，看看代码：

int count_deter(determinant d)//计算行列式
{
	int result = 0;//行列式最终结果
	int one = 1;//行列式的每个积
	int *col_arrange = (int*)malloc(sizeof(int)*d.n);//创建一个记录列排列的数组
	if (col_arrange == NULL)
	{
		printf("%s", strerror(errno));
	}

	for (int i = 0; i < d.n; i++)
	{
		*(col_arrange+i) = i;
	}//赋值

	for (int r = 0; r < fac(d.n); r++)//所有行列式的结果
	{
		one = 1;
		tran_arrange(col_arrange,r,d.n);//确定一个列的排序
		for (int i = 0; i < d.n; i++)
		{
			int a = *(d.p+(i*d.n)+col_arrange[i]);//行列式每个数
			one = a*one;
		}
		result += one* count_re(col_arrange, d.n);
	}

	free(col_arrange);
	col_arrange = NULL;
	return result;
}

4.总结

这便是所有内容，感谢观看~_~,贴个代码： 

#define _CRT_SECURE_NO_WARNINGS

#include <stdio.h>
#include <stdlib.h>
#include <errno.h>
#include <string.h>

typedef struct
{
	int n;
	int* p;
}determinant;

int fac(int n)//计算阶乘
{
	int result=1;
	for (int i = n; i > 0; i--)
	{
		result *= i;
	}
	return result;
}

void change(int* a, int* b)//交换函数
{
	int c;
	c = *b;
	*b = *a;
	*a = c;
}

void bubble_sort(int* arrange,int n)//排序
{
	for (int i = 0; i < n - 1; i++)
	{
		for (int j = 0; j < n - 1 - i; j++)
		{
			if (arrange[j] > arrange[j + 1])
				change(arrange+j, arrange+j+1);
		}
	}
}

void tran_arrange(int* arrange, int j, int n)//改变排列
{
	if (j == 0)//自然顺序做首项
	{
		return;
	}
	if (n == 2)
	{
		change(arrange, arrange + 1);
	}
	else if(n == 3)//等于3做奇偶不同的对换
	{
		if (j % 2 == 0)
		{
			change(arrange, arrange + 1);
		}
		else if (j % 2 == 1)
		{
			change(arrange + 1, arrange + 2);
		}
	}
	else //大于3递归
	{
		if (j % fac(n - 1) == 0)//在后面数循环结束后，调换首数和后面一个数，产生新的组合
		{
			int tmp = j;
			if (tmp > fac(n))
				tmp = tmp % fac(n);//tmp要在fac(n)内，防止（1）越界

			bubble_sort(arrange + 1,n-1);//排序方便调换
			if (tmp == 0)//自然顺序做首排列，第0次
			{
				return;
			}
			if (tmp == fac(n - 1))
				change(arrange, arrange + n - 1);//调换算法，第一次换最后一位
			else
				change(arrange, arrange + n - (tmp / fac(n - 1)) + 1);//（1）之后“/”
		}
		tran_arrange(arrange+1,j,n-1);
	}
}


void build_deter(determinant *d)//创建行列式
{
	printf("请输入行列数:>");//输入行列式
	scanf("%d", &d->n);

	d->p = (int*)malloc(sizeof (int)* d->n*d->n);//为行列式分配空间
	if (d->p == NULL)
	{
		printf("%s", strerror(errno));
	}

	printf("请按行输入值:>");
	for (int i = 0; i < d->n; i++)
	{
		for (int j = 0; j < d->n; j++)
		{
			scanf("%d",(d->p)+(i*d->n+j));
		}
	}
	printf("building over!\n");
}

void dis_deter(determinant d)//打印行列式
{
	for (int i = 0; i < d.n; i++)
	{
		printf("|");
		for (int j = 0; j < d.n; j++)
		{
			printf("%-3d ", *((d.p) + (i * d.n + j)));
		}
		printf("|\n");
	}
	printf("printing over!\n");
}

int count_re(int* arrange,int n)//统计逆序数
{
	int count = 0;
	for (int i = 0; i < n; i++)
	{
		for (int j = i+1; j < n; j++)
		{
			if (arrange[i] > arrange[j])
				count++;
		}
	}
	if (count % 2 == 1)
		return -1;
	else
		return 1;
}


int count_deter(determinant d)//计算行列式
{
	int result = 0;//行列式最终结果
	int one = 1;//行列式的每个积
	int *col_arrange = (int*)malloc(sizeof(int)*d.n);//创建一个记录列排列的数组
	if (col_arrange == NULL)
	{
		printf("%s", strerror(errno));
	}

	for (int i = 0; i < d.n; i++)
	{
		*(col_arrange+i) = i;
	}//赋值

	for (int r = 0; r < fac(d.n); r++)//所有行列式的结果
	{
		one = 1;
		tran_arrange(col_arrange,r,d.n);//确定一个列的排序
		for (int i = 0; i < d.n; i++)
		{
			int a = *(d.p+(i*d.n)+col_arrange[i]);//行列式每个数
			one = a*one;
		}
		result += one* count_re(col_arrange, d.n);
	}

	free(col_arrange);
	col_arrange = NULL;
	return result;
}


void free_deter(determinant* d)//释放空间
{
	free(d->p);
	d->p = NULL;
}

int main()
{

	determinant d;
	build_deter(&d);
	dis_deter(d);
	printf("\n%d",count_deter(d));
	free_deter(&d);
	return 0;

呜呜哇哇2
关注
  • 12
    点赞
  • 踩
  • 9
    收藏
    觉得还不错? 一键收藏
  • 0
    评论
C实现行列式计算
X_0001的博客
07-19 2015
行列式是线性代数学科中最基础的部分，目前基于行列式计算，有诸多网站开发完成直接计算的在线程序，但是，码农靠寄几呀~ 直接上代码： //行列式计算 #include<stdio.h> #include<math.h> int func(int a[][10],int x); int func2(int a[][10],int y, int k,int m); int test(int a[][10]) { return a[0][0]; } int main() { ...
纯C语言行列式的计算
weixin_63129581的博客
03-23 2908
萌新写了一整天终于写出来了，其中部分算法是从B站正月点灯笼 大佬那里学来的（网址[算法教程] 全排列_哔哩哔哩_bilibili），特此感谢。 支持输入分数，负数。目前无法输入小数。 输入时，同一行的元素用一个空格隔开，别多打，没有做冗余。每行元素之间用一个回车隔开，双击回车结束输入。我设置的最多是5阶行列式，如果需要更多的话就改两个#define，第一个是阶数的上限，第二个是上限的阶乘（主要是用来储存tor的）。 大概这样，上代码： #include<stdio.h> #includ
行列式计算程序设计
肖泽云
11-22 2031
由于里面有很多公式，不便于编辑，所以整理成图片！在后面附有程序所有源代码！ 综合上面的分析，现定义一个函数CalculateDet(double[,] A, int DetN)用于计算矩阵A[n,n]的行列式，其中输入参数A表示矩阵，DetN表示矩阵的阶数。该函数的代码如下：         private double CalculateDet(double[,] A, int DetN)   
c语言计算行列式
tgkyrie的博客
10-16 5698
大一刚学了c和线性代数，尝试一下写个行列式求值程序 #include<stdio.h> #include<math.h> int main(){ int det(int b[][10],int x); //定义求行列式函数 int a[10][10]; //定义一个固定数组 int n,i=0,j=0,l; printf("输入行列式阶数\n"); scanf("%d",&n); printf("输入行列式元素空格隔开\n") for(i=0;i<n;i+
C语言求行列式的值
最新发布
amen_yanwang的博客
04-20 796
c语言来求行列式的计算值
n阶行列式计算 C语言 实现
08-17
用C语言实现的n阶行列式计算,希望大家能多提意见让程序更加完善~
C语言n阶行列式算法(递归法)
m0_62264923的博客
07-17 807
【代码】C语言n阶行列式算法(递归法)
任意阶数实数方阵的行列式的值的C语言实现详解
Inspire
12-29 1109
算法描述 对于一个N阶的方阵，则它的行列式的值为，这里我们引入余子式的概念，在线性代数的课本中它的定义如下：在N阶行列式中，把元素所在的第i行和第j列划去后，留下了的N-1阶行列式叫元素的余子式，记作。则方阵A的行列式可以表示为。对于式子里头的，它是一个N-1阶的行列式，我们可以按照上述的方法进一步的分解，得到N-2阶的行列式，然后再分解得到N-3阶行列式，一直分解到2阶的行列式，这时2阶的行列式就变的非常的好求了，即，从这里我们就可以看到这是一个递归求解的过程，因此可以非常方便的用C语言...
c语言，通过计算行最简的方式来求行列式的值
weixin_63129581的博客
10-02 848
纯c语言，化简行最简式，计算行列式的值
行列式计算的C语言小程序
06-03
行列式计算的C语言小程序，用了简单的递归算法实现的，希望大家多多指教。
C语言计算行列式 4 段程序
02-03
搜集了4个用C语言计算行列式的程序，基本算法都是用高斯消元法，不过在消元之前，先把绝对值最大的元素进行适当的位置调整。希望能对C语言编程初学者有所帮助。
C 计算行列式的值
04-06
C语言，自己编的。手动输入行列式，然后计算行列式的值。
使用C++计算行列式
11-18
使用C++实现了任意阶行列式的计算。程序简单，清晰易懂。
五阶以内行列式计算 c语言
01-02
简单的c语言程序，计算行列式，五阶以内，但可以通过修改增加阶数，新手练手作品
C语言实现行列式计算
热门推荐
Dream_wave的博客
07-28 1万+
C语言实现行列式计算 前言 行列式的计算有多种方法，其中按行（列）展开法最为直观，因此也最容易编写为程序。本文使用C语言，利用函数的递归，计算任意阶行列式的值。注意：本程序中，在计算行列式之前，需输入行列式的阶数。 C语言在数学计算领域有诸多不便之处，例如数组的长度不可更改，这一点对于行列式的计算来说是极不方便的，因为余子式的阶数比原行列式的阶数要小，而在递归的过程中涉及到数组的长度。为此，设立两个函数，第一个函数DET为计算的主要部分，第二个函数Minor专门用于计算余子式的值，而第二个函数中会调用第一个
行列式计算c语言
qq_52959938的博客
04-06 371
行列式计算c语言 #include<stdio.h> #include<math.h> int main()## 标题 { int i, j, p, q,t; int s,h, k,s,g,m; //h = a[i][j], k = b[p][q]; g = h * k s += h * k; int a[3][4]; int b[4][2]; printf("输入a\n"); //输入数组a for (
【应用C】C语言实现行列式与矩阵的运算系统（+源代码）
单眼皮根号三的博客
01-11 4219
文章目录01 - 线性代数02 - 数据结构03 - 行列式运算04 - 矩阵运算05 - 结果测试 01 - 线性代数 02 - 数据结构 03 - 行列式运算 04 - 矩阵运算 05 - 结果测试    ...
C语言实现行列式计算（展开法）
luoshiyong123的博客
06-28 5934
如下图所示行列式： 程序如下： #include<stdio.h> #include<math.h> #include<fstream> #include<iostream> using namespace std; double get_num(double ori[][5],int m)//计算m行 行列式的值 { for(int i=0;i<m;i++) { for(int j=0
c语言 计算行列式,新手作品：行列式计算C语言版
weixin_42613501的博客
05-20 164
该楼层疑似违规已被系统折叠隐藏此楼查看此楼对话 ControlHeightDecrease Shift+Up Arrow 向上调整选定的控件或对话一个对话单位对话 ControlHeightIncrease Shift+Down Arrow 向下调整选定的控件或对话一个对话单位对话 ControlMoveDown Down Arr...
评论
添加红包

请填写红包祝福语或标题

个

红包个数最小为10个

元

红包金额最低5元

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

抵扣说明：

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

余额充值