高斯消元法模板

最新推荐文章于 2024-03-14 17:09:47 发布
最新推荐文章于 2024-03-14 17:09:47 发布
阅读量1k 收藏
点赞数
分类专栏： 数学
版权声明：本文为博主原创文章，遵循 CC 4.0 BY-SA 版权协议，转载请附上原文出处链接和本声明。
本文链接：https://blog.csdn.net/xianxingwuguan1/article/details/20716105
版权

浮点数版（以hdu 3359为例）：

#pragma comment(linker, "/STACK:102400000,102400000")
#include <stdio.h>
#include <iostream>
#include <algorithm>
#include <sstream>
#include <stdlib.h>
#include <string.h>
#include <limits.h>
#include <string>
#include <time.h>
#include <math.h>
#include <queue>
#include <stack>
#include <set>
#include <map>
using namespace std;
#define INF 0x3f3f3f3f
#define eps 1e-8
#define pi acos(-1.0)
typedef long long ll;
double s[20][20],a[200][200],x[200];
int equ,var;
int Guass(){//初始化时等号右边的常数放到x数组里，计算结果保存在x数组里。
	int i,j,k,col,max_r;
	for(k=0,col=0;k<equ&&col<var;k++,col++){
		max_r=k;
		for(i=k+1;i<equ;i++)
			if(fabs(a[i][col])>fabs(a[max_r][col]))
				max_r=i;
		if(fabs(a[max_r][col])<eps)return 0;
		if(k!=max_r){
			for(j=col;j<var;j++)
				swap(a[k][j],a[max_r][j]);
			swap(x[k],x[max_r]);
		}
		x[k]/=a[k][col];
		for(j=col+1;j<var;j++)a[k][j]/=a[k][col];
		a[k][col]=1;
		for(i=0;i<equ;i++)
			if(i!=k){
				x[i]-=x[k]*a[i][k];
				for(j=col+1;j<var;j++)
					a[i][j]-=a[k][j]*a[i][col];
				a[i][col]=0;
			}
	}
	return 1;
}
int main()
{
     //freopen("data.in","r",stdin);
     //freopen("data.out","w",stdout);
	 int d,tt=0,m,n;
     while(~scanf("%d%d%d",&m,&n,&d)&&(n+m+d)){
		 if(tt++)puts("");
		 for(int i=0;i<n;i++)
			 for(int j=0;j<m;j++)
				 scanf("%lf",&s[i][j]);
		 memset(a,0,sizeof(a));
		 equ=var=n*m;
		 for(int i=0;i<n;i++)
			 for(int j=0;j<m;j++){
				 int cnt=0;
				 for(int k=0;k<n;k++)
					 for(int p=0;p<m;p++)
						 if(abs(i-k)+abs(j-p)<=d){
							 a[i*m+j][k*m+p]=1;
							 cnt++;
						 }
				 x[i*m+j]=cnt*s[i][j];
			 }
		 Guass();
		// cout<<ha<<endl;
		 for(int i=0;i<n;i++){
			 for(int j=0;j<m;j++) 
				 printf ("%8.2lf", x[i*m+j]);
			 puts("");
		 }
	 }
     return 0;
}


整数版：

/* ***********************************************
Author :rabbit
Created Time :2014/3/7 15:19:27
File Name :11.cpp
************************************************ */
#pragma comment(linker, "/STACK:102400000,102400000")
#include <stdio.h>
#include <iostream>
#include <algorithm>
#include <sstream>
#include <stdlib.h>
#include <string.h>
#include <limits.h>
#include <string>
#include <time.h>
#include <math.h>
#include <queue>
#include <stack>
#include <set>
#include <map>
using namespace std;
#define INF 0x3f3f3f3f
#define eps 1e-8
#define pi acos(-1.0)
typedef long long ll;
const int maxn=50;
int a[maxn][maxn];//增广矩阵。
int x[maxn];//解集。
bool free_x[maxn];//标记是否为不确定的变元。
int gcd(int a,int b){
	if(a==0)return b;
	return gcd(b%a,a);
}
int lcm(int a,int b){
	return a/gcd(a,b)*b;
}
int Guass(int equ,int var){
	int max_r;//当前这一列最大的行。
	int ta,tb,LCM,temp,free_x_num,free_index,i,j,k,col;
	for(i=0;i<=var;i++)x[i]=0,free_x[i]=true;
	for(col=0,k=0;k<equ&&col<var;k++,col++){
		max_r=k;
		for(i=k+1;i<equ;i++)
			if(abs(a[i][col])>abs(a[max_r][col]))max_r=i;
		if(max_r!=k){
			for(j=k;j<var+1;j++)swap(a[k][j],a[max_r][j]);
		}
		if(a[k][col]==0){
			k--;continue;//col列第k行之下全是0.处理当前行的下一列。
		}
		for(i=k+1;i<equ;i++)
			if(a[i][col]){
				LCM=lcm(abs(a[i][col]),abs(a[k][col]));
				ta=LCM/abs(a[i][col]);
				tb=LCM/abs(a[k][col]);
				if(a[i][col]*a[k][col]<0)tb=-tb;
				for(j=col;j<var+1;j++)
					a[i][j]=a[i][j]*ta-a[k][j]*tb;
			}
	}
	for(i=k;i<equ;i++)if(a[i][col])return -1;//无解。
	if(k<var){
		for(i=k-1;i>=0;i--){
			free_x_num=0;
			for(j=0;j<var;j++)
				if(a[i][j]&&free_x[j])free_x_num++,free_index=j;
			if(free_x_num>1)continue;//无法求解出确定变元。
			temp=a[i][var];
			for(j=0;j<var;j++)
				if(a[i][j]&&j!=free_index)temp-=a[i][j]*x[j];
			x[free_index]=temp/a[i][free_index];//求出变元。
			free_x[free_index]=0;//该变元是确定的。
		}
		return var-k;//自由变元有var-k个。
	}
	for(i=var-1;i>=0;i--){
		temp=a[i][var];
		for(j=i+1;j<var;j++)
			if(a[i][j])temp-=a[i][j]*x[j];
		if(temp%a[i][i])return -2;//有浮点数解，但是无整数解。
		x[i]=temp/a[i][i];
	}
	return 0;
}
int main()
{
    int i, j;
    int equ,var;
    while (scanf("%d %d", &equ, &var) != EOF){
        memset(a, 0, sizeof(a));
        for (i = 0; i < equ; i++)
            for (j = 0; j < var + 1; j++)
                scanf("%d", &a[i][j]);
        int free_num = Guass(equ,var);
        if (free_num == -1) printf("无解!\n");
        else if (free_num == -2) printf("有浮点数解，无整数解!\n");
        else if (free_num > 0){
            printf("无穷多解! 自由变元个数为%d\n", free_num);
            for (i = 0; i < var; i++){
                if (free_x[i]) printf("x%d 是不确定的\n", i + 1);
                else printf("x%d: %d\n", i + 1, x[i]);
            }
        }
        else{
            for (i = 0; i < var; i++)
                printf("x%d: %d\n", i + 1, x[i]);
        }
        printf("\n");
    }
    return 0;
}
/*
4 4
2 1 -5 1 8
1 -3 0 -6 9
0 2 -1 2 -5
1 4 -7 6 0

4 4
1 1 1 1 5
1 2 -1 4 -2
2 -3 -1 -5 -2
3 1 2 11 0
*/

位运算高斯消元（以poj  1222为例）：

代码：

/* ***********************************************
Author :xianxingwuguan
Created Time :2014/3/7 14:25:12
File Name :11.cpp
************************************************ */
#pragma comment(linker, "/STACK:102400000,102400000")
#include <stdio.h>
#include <iostream>
#include <algorithm>
#include <sstream>
#include <stdlib.h>
#include <string.h>
#include <limits.h>
#include <string>
#include <time.h>
#include <math.h>
#include <queue>
#include <stack>
#include <set>
#include <map>
using namespace std;
#define INF 0x3f3f3f3f
#define eps 1e-8
#define pi acos(-1.0)
typedef long long ll;
const int maxn=50;
int a[maxn][maxn],x[maxn],free_num,free_x[maxn];
int Guass(int equ,int var){//-1返回无解，0返回唯一解，否则返回自由变元个数，
	int i,j,col,k,free_num=0;
	memset(x,0,sizeof(x));
	for(k=0,col=0;k<equ&&col<var;k++,col++){
		int max_r=k;
		for(i=k+1;i<equ;i++)
			if(abs(a[i][col])>abs(a[max_r][col]))max_r=i;
		if(max_r!=k){
			for(j=k;j<var+1;j++)
				swap(a[k][j],a[max_r][j]);
		}
		if(a[k][col]==0){
			k--;
			free_x[free_num++]=col;//自由变元
			continue;
		}
		for(i=k+1;i<equ;i++)
			if(a[i][col]){
				for(j=col;j<var+1;j++)
					a[i][j]^=a[k][j];
			}
	}
	for(i=k;i<equ;i++)if(a[i][col])return -1;//无解
	if(k<var)return var-k;//自由变元个数
	for(int i=var-1;i>=0;i--){
		x[i]=a[i][var];
		for(int j=i+1;j<var;j++)
			x[i]^=(a[i][j]&x[j]);
	}
	return 0;
}
int main()
{
     //freopen("data.in","r",stdin);
     //freopen("data.out","w",stdout);
 	 int T,u,v,n;
	 scanf("%d",&T);
	 for(int t=1;t<=T;t++){
		 printf("PUZZLE #%d\n",t);
		 memset(a,0,sizeof(a));
		 int var=30,equ=30;
		 for(int i=0;i<5;i++)
			 for(int j=0;j<6;j++){
				 int p=i*6+j;
				 a[p][p]=1;
				 if(i>0)a[(i-1)*6+j][p]=1;
				 if(i<4)a[(i+1)*6+j][p]=1;
				 if(j>0)a[i*6+j-1][p]=1;
				 if(j<5)a[i*6+j+1][p]=1;
			 }
         for(int i = 0;i < 30;i++)scanf("%d",&a[i][30]);
         Guass(equ,var);
         for(int i = 0;i < 5;i++)
		 {
             for(int j = 0;j < 5;j++)printf("%d ",x[i*6+j]);
             printf("%d\n",x[i*6+5]);
         }
	 }
     return 0;
}


_rabbit
关注
专栏目录
高斯消元法原理及模板
tianyuanzhidian的博客
10-11 2519
1.线形代数中的高斯消元法 高斯消元法：简单来说，就是消减未知元，来解多个方程组的一种方法，在线性代数中经常把方程组换算成矩阵，通过行列之间的加减法，把原来的矩阵转换成一个三角矩阵，然后回代求解。 解方程组例： (1)x1+x2+x3=3 2*x1+3*x2+x3=5 3*x1+5*x2+6*x3=17 (2)第二式减去第一式的二倍，第三式减去第一式的三倍。 x1+x2+x3=3 x2-x3=-1 2*x2+3*x3=8 (3)第三式减去第2式的2倍 x1+...
高斯消元法模板--矩阵分解
M_Lter的博客
07-17 614
求解方程组的解向量：Ax=b， 令 A=LU； //求解L矩阵和U矩阵可以利用一般高斯消元法，也可以用直接法 则 LUx=b; 令 y=Ux; 则 Ly=b; 所以根据L矩阵和b向量，可求出y向量，再根据y向量和U矩阵求解x向量即可 #include&lt;cstdio&gt; #include&lt;cmath&gt; using namespace std; const int...
C语言复习6
C语言学习
03-09 1958
二维数组： 类型名 数组名【行长度】 【列长度】 #include<stdio.h> #define MAXM 6 #define MAXN 6 int main(void) { int col,row,i,j,m,n; int a[MAXM][MAXN]; printf("enter number of m ,n\n"); scanf("%d %d",&m,&n); printf("please enter number %d :",m*n); f
数组 最大值，最小值的选择
damifeng8的专栏
06-17 689
/*   * 在n行n列的数组里面按一下要求选出两个数。首先从每行选出最大的数，再从n个最大的数里面选出一个最小的数；  * 其次，从每行选出最小的数，再从选出的最小的数里面选出一个最大的数。  * */ #include #include #include #define MAXN 20 int a[MAXN][MAXN] ; void main(void)
编程求二维数组(矩阵)的子矩阵之和的最大值
m0_71272694的博客
12-07 168
编程求二维数组(矩阵)的子矩阵之和的最大值。解析：
C/C++蛇形填数
My_YueR的博客
04-21 4110
再 n x n 方阵里填入1，2，···，n x n，要求填成蛇形。例如，n = 4 时方阵为： 10 11 12 1 9 16 13 2 8 15 14 3 7 6 5 4 上面方阵中，多余的空格只是为了便于观察规律，不必严格输出。n <= 8。 分析：我们可以用二维数组来储存题目中的矩阵，我们一般在进行数组声明时都会设的比题目要求大一点，这里我们可以成"int a[50][50]"。我们数组升序是按蛇形升序，即从a[0][49]开始，到a[49][49]结束，再向左开始升序即...
3-3例题-蛇形填数c/c++语言
青木年华
05-30 815
#include&lt;cstdio&gt;#include&lt;string.h&gt;#define maxn 20int a[maxn][maxn];using namespace std;int main(){  int x,y,tot=0,n;  scanf("%d",&amp;n);  memset(a,0,sizeof(a));  tot=a[x=0][y=n-1]=1;//if ...
【高斯消元】【模板】高斯消元法
blog
10-16 134
高斯消元法模板
python高斯消元法
qq_33599874的博客
03-14 676
防止对角线元素为0的情况。
算法文档无代码用高斯消元法求解线性方程组
最新发布
04-15
高斯消元法是一种用于求解线性方程组的数值方法。该方法通过行变换将线性方程组的系数矩阵转换为阶梯形矩阵，进而简化为行最简形矩阵，从而便于求解方程组。下面详细地解析该方法的各个步骤和关键知识点。 首先，...
高斯消元法小总结和瞎总结
如不如此都如此的博客
05-24 5246
啊 搞到这个算法 ，大概已经在一周前了 ，一周前做去年西安邀请赛的题，其中有一题是要用到高斯消元法求系数，然后再用矩阵快速幂求正解 ，队里ZK命自认要学数学的我 两天搞清楚那个题。。然后我就从高斯消元法入手了，大概看了两天吧，也有一定的了解了？（其实说麻烦挺麻烦，说简单也挺简单的。。主要是我天资愚钝，查阅了好多好多相关博客，模板，题解 ，还是不是很懂。。。）时间大概过了有一周吧，凭借记忆来复述一把...
【OI学习笔记】线性代数-高斯消元法
lovely_cyqの小窝
09-02 441
线性代数中的高斯消元法学习笔记，较为详细地叙述了原理、算法步骤以及给出了代码。
dp 讲课代码汇总
daydreamer23333的博客
12-10 346
下面是本节课的所涉及取模知识和题目的代码汇总 因为头文件太长了 我就没贴上来了 需要用的话记得加上相应的头文件 1.加法取模、乘法取模、减法取模 ans=(ans+b)%mod; ans=ans*b%mod; ans=(ans-b+mod)%mod 2.最长上升子序列LIS（n^2dp做法） 题号 poj2533 int a[maxn],dp[maxn]; int main() { int n; scanf("%d",&n); for(int i=1;i<=n;i+
hdu 3359 高斯消元法
_rabbit的专栏
03-07 1485
Kind of a Blur Time Limit: 2000/1000 MS (Java/Others)    Memory Limit: 32768/32768 K (Java/Others) Total Submission(s): 1226    Accepted Submission(s): 356 Problem Description Image blurring o
POJ 3243 大步小步算法
_rabbit的专栏
04-02 1248
Clever Y Time Limit: 5000MS   Memory Limit: 65536K Total Submissions: 6195   Accepted: 1500 Description Little Y finds there is a very interesting formula in mathematics:
hdu 3723 卡特兰数
_rabbit的专栏
04-04 1248
Delta Wave Time Limit: 6000/3000 MS (Java/Others)    Memory Limit: 32768/32768 K (Java/Others) Total Submission(s): 690    Accepted Submission(s): 223 Problem Description A delta wave is a hig
HDU 4336 概率dp或者容斥原理
_rabbit的专栏
04-02 1164
Card Collector Time Limit: 2000/1000 MS (Java/Others)    Memory Limit: 32768/32768 K (Java/Others) Total Submission(s): 2087    Accepted Submission(s): 984 Special Judge Problem Description In
POJ 2084 卡特兰数+高精度模板
_rabbit的专栏
03-27 1131
Game of Connections Time Limit: 1000MS   Memory Limit: 30000K Total Submissions: 7136   Accepted: 3625 Description This is a small but ancient game. You are supposed to wri
高斯消元法(Gauss Elimination)【超详解&模板】
05-25
高斯消元法是一种线性代数中求解线性方程组的常用方法，它通过一系列的行变换将增广矩阵转化为上三角矩阵，从而求出方程组的解。 下面我们来详细介绍高斯消元法的步骤： 1. 将增广矩阵写成一个矩阵形式：$$\left[\...
评论
添加红包

请填写红包祝福语或标题

个

红包个数最小为10个

元

红包金额最低5元

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

抵扣说明：

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

余额充值