【模板】高斯消元

已于 2022-02-23 21:11:36 修改
阅读量92 收藏
点赞数
分类专栏： 数论 文章标签： 高斯消元 行列式 逆矩阵 线性方程组 异或方程组
于 2021-10-16 18:11:11 首次发布
版权声明：本文为博主原创文章，遵循 CC 4.0 BY-SA 版权协议，转载请附上原文出处链接和本声明。
本文链接：https://blog.csdn.net/fighting4years/article/details/120801915
版权

高斯消元

一:求行列式
1除法

double gauss(){//高斯消元求行列式 
	int max_row;
	for(int row = 1, col = 1; row <= n && col <= n; row++, col++){
		//找到当前列绝对值最大的行 
		max_row = row;
		for(int i = row + 1; i <= n; i++) if(fabs(a[i][col]) > fabs(a[max_row][col] + eps)) max_row = i;
		if(fabs(a[max_row][col] < eps)) {return 0;}//有一个0， 行列式为0 
		//交换绝对值最大的行/ 减小除法的精度损失 
		if(row != max_row){
			for(int j = 1; j <= n; j++) swap(a[row][j], a[max_row][j]);
		}
		//消去当前列(只消下面行的)
		for(int i = row + 1; i <= n; i++){
			double temp = a[i][col] / a[row][col];
			for(int j = col; j <= n; j++){
				a[i][j] -= temp * a[row][j];
			}
		} 
	} 
	double ans = 1;
	for(int i = 1; i <= n; i++) ans = ans * a[i][i];
	return ans;
}

2逆元
ICPC 2021济南 J

#include<bits/stdc++.h>
using namespace std;
typedef long long ll;
const int N = 100 + 10;
const ll mod = 1e9 + 7;

ll a[N][N];
int n;
ll qpow(ll base, ll pow, ll mod){
    ll ans = 1;
    while(pow){
        if(pow & 1) ans = ans * base % mod;
        base = base * base % mod;
        pow >>= 1;
    }
    return ans;
}

ll gauss_det_inv(ll a[][N], int n, ll mod){//高斯消元求行列式(逆元) 
    int flag = 1;
    for(int row = 1, col = 1; row <= n && col <= n; row++, col++){
        int max_row = row;
        //排除0的情况 
        for(int i = row; i <= n; i++){
            if(a[i][col] != 0) {max_row = i; break;}
        }
        if(row != max_row){
            for(int j = 1; j <= n; j++) swap(a[row][j], a[max_row][j]);
            flag = -flag;//交换行列式两行， 行列式值取反 
        }
        //做除法
        ll inv = qpow(a[row][col], mod - 2, mod);
        for(int i = row + 1; i <= n; i++){
            ll temp = a[i][col] * inv % mod;
            for(int j = col; j <= n; j++){
                a[i][j] = ((a[i][j] - temp * a[row][j]) % mod + mod) % mod;
            }
        } 
    }
    ll ans = 1;
    for(int i = 1; i <= n; i++){
        ans = ans * a[i][i] % mod;
    }
    return (ans * flag % mod + mod) % mod;
}

char str[10010];
ll get_det(ll mod){
    ll ret = 0;
    int len = strlen(str + 1);
    for(int i = 1; i <= len; i++){
        ret = (ret * 10 + (str[i] - '0')) % mod;
    }
    return ret;
}
void solve(){
    scanf("%d", &n);
    scanf("%s", str + 1);
    ll det = get_det(mod);
    for(int i = 1; i <= n; i++){
        for(int j = 1; j <= n; j++){
            scanf("%lld", &a[i][j]);
        }
    }
    ll ans = gauss_det_inv(a, n, mod);
    if(ans == det)printf("+\n");
    else printf("-\n");
}

int main() {
    int t = 1;
    scanf("%d", &t);
    while(t--)solve();
    
    return 0;
}

3解线性方程组
1.1: 求具体的解
（1）唯一解

void gauss(double mat[N][N], double extra[N]/*矩阵的增广部分*/, int rank/*矩阵的秩*/){//化为对角矩阵
	//数据保证唯一解
	int maxRow;
	for(int row = 1, col = 1; row <= rank && col <= rank; row++, col++){
		//找到当前列绝对值最大的行 
		maxRow = row;
		for(int i = row + 1; i <= rank; i++) if(fabs(mat[i][col]) > fabs(mat[maxRow][col])) maxRow = i;
		//交换绝对值最大的行/ 减小除法的精度损失 
		if(row != maxRow){
			for(int j = 1; j <= rank; j++) swap(mat[row][j], mat[maxRow][j]);
			swap(extra[row], extra[maxRow]);
		}
		for(int i = 1; i <= rank; i++){
			if(i == row) continue;
			double rate = mat[i][col] / mat[row][col];
			for(int j = col; j <= rank; j++){
				mat[i][j] -= rate * mat[row][j];
			}
			extra[i] -= rate * extra[row];
		}
	}
	//除掉系数
	for(int i = 1; i <= rank; i++){
		extra[i] /= mat[i][i];
	}
}

例题
ACwing207 球形空间产生器

#include <bits/stdc++.h>
using namespace std;
const int N = 15;
const double eps = 1e-8;
int n, m;
double a[N][N];
double mat[N][N], extra[N];//矩阵的增广部分
void gauss(double mat[N][N], double extra[N]/*矩阵的增广部分*/, int rank/*矩阵的秩*/){//化为对角矩阵
	//数据保证唯一解
	int maxRow;
	for(int row = 1, col = 1; row <= rank && col <= rank; row++, col++){
		//找到当前列绝对值最大的行 
		maxRow = row;
		for(int i = row + 1; i <= rank; i++) if(fabs(mat[i][col]) > fabs(mat[maxRow][col])) maxRow = i;
		//交换绝对值最大的行/ 减小除法的精度损失 
		if(row != maxRow){
			for(int j = 1; j <= rank; j++) swap(mat[row][j], mat[maxRow][j]);
			swap(extra[row], extra[maxRow]);
		}
		for(int i = 1; i <= rank; i++){
			if(i == row) continue;
			double rate = mat[i][col] / mat[row][col];
			for(int j = col; j <= rank; j++){
				mat[i][j] -= rate * mat[row][j];
			}
			extra[i] -= rate * extra[row];
		}
	}
	//除掉系数
	for(int i = 1; i <= rank; i++){
		extra[i] /= mat[i][i];
	}
}

void solve(){
	cin >> n;
	for(int i = 1; i <= n + 1; i++){
		for(int j = 1; j <= n; j++){
			cin >> a[i][j];
		}
	}
	for(int i = 1; i <= n; i++){
		for(int j = 1; j <= n; j++){
			mat[i][j] = 2 * (a[i + 1][j] - a[i][j]);
			extra[i] += (a[i + 1][j] * a[i + 1][j] - a[i][j] * a[i][j]);
		}
	}
	gauss(mat, extra, n);
	for(int i = 1; i <= n; i++){
		cout << fixed << setprecision(3) << extra[i]<< " ";
	}
}
int main()
{
	ios::sync_with_stdio(false);
	cin.tie(0); cout.tie(0);
	int t = 1;
	while(t--) solve();

	return 0;
}

4异或方程组
1.1求解的个数

int gaussXor(int mat[N], int rank) {//化为对角矩阵
	//由于异或只有0和1，所以将一行的状态压缩为一个二进制数,故mat只有一维
	//将矩阵的增广部分压缩进mat的第0位
	//Q: 为什么不是最高位
	//A: 二进制数高位到低位权值递减， 只有当系数矩阵一行的所有数均为0时， 才考虑矩阵的增广部分是否为0

	//返回值： 解的个数
	int maxRow;
	for (int row = 1; row <= rank; row++) {
		//找到当前值最大的行 
		for (int i = row + 1; i <= rank; i++) {
			if (mat[i] > mat[row]) swap(mat[i], mat[row]);
		}
		//判断接下来的行存在系数全为0的行
		if (mat[row] == 0) { return 1 << (rank - row + 1); }
		if (mat[row] == 1) { return 0; }
		for (int bit = n; bit; bit--) {
			if(mat[row] >> bit & 1) {
				for (int i = 1; i <= rank; i++) {
					if(i == row) continue;
					if(mat[i] >> bit & 1) mat[i] ^= mat[row];
				}
				break;
			}
		}
	}
	return 1;
}

例题
ACwing208 开关问题

#include <bits/stdc++.h>
using namespace std;
const int N = 15;
int n, m;
int mat[N];//矩阵的增广部分
int gaussXor(int mat[N], int rank) {//化为对角矩阵
	//由于异或只有0和1，所以将一行的状态压缩为一个二进制数,故mat只有一维
	//将矩阵的增广部分压缩进mat的第0位
	//Q: 为什么不是最高位
	//A: 二进制数高位到低位权值递减， 只有当系数矩阵一行的所有数均为0时， 才考虑矩阵的增广部分是否为0

	//返回值： 解的个数
	int maxRow;
	for (int row = 1; row <= rank; row++) {
		//找到当前值最大的行 
		for (int i = row + 1; i <= rank; i++) {
			if (mat[i] > mat[row]) swap(mat[i], mat[row]);
		}
		//判断接下来的行存在系数全为0的行
		if (mat[row] == 0) { return 1 << (rank - row + 1); }
		if (mat[row] == 1) { return 0; }
		for (int bit = n; bit; bit--) {
			if(mat[row] >> bit & 1) {
				for (int i = 1; i <= rank; i++) {
					if(i == row) continue;
					if(mat[i] >> bit & 1) mat[i] ^= mat[row];
				}
				break;
			}
		}
	}
	return 1;
}

void solve(){
	cin >> n;
	for (int i = 1; i <= n; i++) {
		cin >> mat[i];
	}
	for (int i = 1; i <= n; i++) {
		int t;
		cin >> t;
		mat[i] ^= t;
		mat[i] |= 1 << i;
	}
	int x, y;
	while (cin >> x >> y && x && y) {
		mat[y] |= 1 << x;
	}
	int ans = gaussXor(mat, n);
	if(ans) cout << ans << endl;
	else cout << "Oh,it's impossible~!!" << endl;
}
int main(){
	ios::sync_with_stdio(false);
	cin.tie(0); cout.tie(0);
	int t = 1;
	cin >> t;
	while (t--) solve();

	return 0;
}

5 求一组基， 并用这组基进行一系列操作
ACwing210 异或运算

#include <bits/stdc++.h>
using namespace std;
typedef long long ll;
const int N = 510;
bool canGetZero = false;
int numOfBase;//最多能取到多少数
void gaussXor (vector<ll>& mat, int n) {
	//此题中高斯消元的主要作用是 得到一组线性无关的基
	int maxRow;
	numOfBase = n;
	for (int row = 1; row <= n; row++) {
		maxRow = 0;
		for (int i = row; i <= n; i++) {
			if (!maxRow || mat[i] > mat[maxRow]) {
				maxRow = i;
			}
		}
		if (mat[maxRow] == 0) {canGetZero = true; numOfBase = row - 1; return;}
		swap(mat[row], mat[maxRow]);
		for (int bit = 60; bit >= 0; bit--) {
			if (mat[row] >> bit & 1) {
				for (int i = 1; i <= n; i++) {
					if (i == row) continue;
					if (mat[i] >> bit & 1) mat[i] ^= mat[row];
				}
				break;
			}
		}
	}
}
void solve(int cas) {
	canGetZero = false;
	cout << "Case #" << cas << ":" << endl;
	int n;
	cin >> n;
	vector<ll> mat(n + 1);
	for (int i = 1; i <= n; i++) {
		cin >> mat[i];
	}
	gaussXor(mat, n);
	int q;
	cin >> q;
	for (int i = 1; i <= q; i++) {
		ll k;
		cin >> k;
		k--;
		if (!canGetZero) k++; 
		ll ans = 0;
		if (k >= (1ll << numOfBase)) cout << -1 << endl;
		else {
			for (int i = 0; i < numOfBase; i++) {
				if (k >= (1ll << (numOfBase - 1 - i))) {
					k -= (1ll << (numOfBase - 1 - i));
					ans ^= mat[i + 1];
				}
			}
			cout << ans << endl;
		}
	}
}
int main() {
	ios::sync_with_stdio(false);
	cin.tie(0); cout.tie(0);
	int t = 1;
	cin >> t;
	for(int i = 1; i <= t; i++) solve(i);

	return 0;
}
blank_monster
关注
  • 0
    点赞
  • 踩
  • 0
    收藏
    觉得还不错? 一键收藏
  • 0
    评论
专栏目录
OI模板 高斯消元/线性基
im_zyINF的博客
04-12 279
高斯消元
ACM模板高斯消元
apple_52197802的博客
09-08 71
ACM模板高斯消元
高斯消元模板
qq_35975367的博客
08-04 132
高斯消元通用 #include<stdio.h> #include<algorithm> #include<iostream> #include<string.h> #include<math.h> using namespace std; const int MAXN=50; int a[MAXN][MAXN];//增广矩阵 int x[MAXN];//解集 bool free_x[MAXN];//标记是否是不确定的变元 /* vo
南阳理工acm常用模板
05-15
### 南阳理工ACM常用模板之高斯消元详解 #### 一、基础知识与预备概念 ##### 1.1 XOR运算基础 XOR运算，即异或运算，在计算机科学和数学领域中有着广泛的应用。它是一种二进制运算，对于任何两个二进制数，如果...
建议收藏算法基础课模板大全
08-12
基础算法 —— 代码模板链接 常用代码模板1——基础算法 排序 ...高斯消元 组合计数 容斥原理 简单博弈论 动态规划 背包问题 线性DP 区间DP 计数类DP 数位统计DP 状态压缩DP 树形DP 记忆化搜索 贪心
kuangbin acm模板超级好用
09-05
2.8 高斯消元（浮点数） . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 33 2.9 FFT . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 34 2.10 高斯消元...
高斯消法元解方程组模板
09-09
2. **上三角矩阵与阶梯形矩阵**：高斯消元的目标是通过行变换将增广矩阵转换为上三角矩阵或阶梯形矩阵，这样可以直接读出方程组的解。 3. **主元选择**：为了保证算法的稳定性，通常会在每一步消元过程中选择当前列...
ACM/ICPC模板
05-13
ACM/ICPC模板 内容大概有这些 其他 --高精度模板 ...--高斯消元(用bigInteger和分数类 ,Java实现) //thanks to love8909 --矩阵运算 --欧拉函数 --初等数论(模线性方程组，gcd相关等) --素数判断
mpu6050陀螺仪使用方法开发教程文档.docx
07-30
mpu6050陀螺仪使用方法
【SCI顶级】龙格库塔算法RUN-CNN-LSTM-Multihead-Attention温度预测【含源码 5747期】.zip
07-30
CSDN海神之光上传的全部代码均可运行，亲测可用，直接替换数据即可，适合小白； 1、代码压缩包内容 主函数：Main .m； 数据； 调用函数：其他m文件；无需运行 运行结果效果图； 2、代码运行版本 Matlab 2019b；若运行有误，根据提示修改；若不会，可私信博主； 3、运行操作步骤 步骤一：将所有文件放到Matlab的当前文件夹中； 步骤二：双击打开除Main.m的其他m文件； 步骤三：点击运行，等程序运行完得到结果； 4、仿真咨询 如需其他服务，可私信博主或扫描博主博客文章底部QQ名片； 4.1 CSDN博客或资源的完整代码提供 4.2 期刊或参考文献复现 4.3 Matlab程序定制 4.4 科研合作 智能优化算法-CNN-LSTM-Mutilhead-Attention回归预测系列程序定制或科研合作方向： 4.4.1 遗传算法GA/蚁群算法ACO-CNN-LSTM-Mutilhead-Attention回归预测 4.4.2 粒子群算法PSO/蛙跳算法SFLA-CNN-LSTM-Mutilhead-Attention回归预测 4.4.3 灰狼算法GWO/狼群算法WPA-CNN-LSTM-Mutilhead-Attention回归预测 4.4.4 鲸鱼算法WOA/麻雀算法SSA-CNN-LSTM-Mutilhead-Attention回归预测 4.4.5 萤火虫算法FA/差分算法DE-CNN-LSTM-Mutilhead-Attention回归预测
大模型的数学之魂：统计学习理论精要
07-30
大型机器学习模型，通常简称为“大模型”，是指具有数十亿甚至数千亿参数的深度神经网络模型。这些模型因为其庞大的参数量和复杂的计算结构，能够处理非常复杂的任务和数据，展现出强大的表达能力和预测性能。大模型在自然语言处理、计算机视觉、语音识别和推荐系统等多个领域都有广泛的应用。 大模型的设计目的是通过训练海量数据来学习复杂的模式和特征，从而具备强大的泛化能力，能够对未见过的数据做出准确的预测。与参数较少、层数较浅的小模型相比，大模型虽然需要更多的计算资源和时间来训练和推理，但它们展现出了一种被称为“涌现能力”的特性，即从原始训练数据中自动学习并发现新的、更高层次的特征和模式。 大模型的发展历程可以分为几个阶段，从早期的卷积神经网络模型发展到以Transformer为代表的全新神经网络模型，再到当前以GPT为代表的预训练大模型阶段。Transformer架构因其自注意力机制而成为大模型的核心技术，它使得模型能够更好地理解和处理文本数据中的长距离依赖关系。 大模型的分类主要包括语言大模型（NLP）、视觉大模型（CV）和多模态大模型。语言大模型专注于处理文本数据和理解自然语言，而视觉大模型用
唯美淡雅四页.zip
最新发布
07-30
【简历模板】工作总结、商业计划书、述职报告、读书分享、家长会、主题班会、端午节、期末、夏至、中国风、卡通、小清新、岗位竞聘、公司介绍、读书分享、安全教育、文明礼仪、儿童故事、绘本、防溺水、夏季安全、科技风、商务、炫酷、企业培训、自我介绍、产品介绍、师德师风、班主任培训、神话故事、巴黎奥运会、世界献血者日、防范非法集资、3D快闪、毛玻璃、人工智能等等各种样式的ppt素材风格。 设计模板、图片素材、PPT模板、视频素材、办公文档、小报模板、表格模板、音效配乐、字体库。 广告设计：海报，易拉宝，展板，宣传单，宣传栏，画册，邀请函，优惠券，贺卡，文化墙，标语，制度，名片，舞台背景，广告牌，证书，明信片，菜单，折页，封面，节目单，门头，美陈，拱门，展架等。 电商设计：主图，直通车，详情页，PC端首页，移动端首页，钻展，优惠券，促销标签，店招，店铺公告等。 图片素材：PNG素材，背景素材，矢量素材，插画，元素，艺术字，UI设计等。 视频素材：AE模板，会声会影，PR模板，视频背景，实拍短片，音效配乐。 办公文档：工作汇报，毕业答辩，企业介绍，总结计划，教学课件，求职简历等PPT/WORD模板。
仿真验证阵列可以提高接收信号信噪比
07-30
这是一份模拟了阵列输入信号及噪声，并验证了相对于阵列接收到的信号，阵列输出信号可以将信噪比提高M倍，其中M为阵列的阵元个数的代码。 代码中可以随意修改阵元个数、阵元间距、波束指向角度、信号频率等。 代码中关键部分均含有文字注释，完全不必担心看不懂。 无论是从仿真波形，还是计算的信噪比结果均能看出阵元数为M的阵列将信号的信噪比提高了M倍。
[毕业设计]Java驱动的C语言编程练习与考试平台(源代码+论文).zip
07-30
[毕业设计]Java驱动的C语言编程练习与考试平台(源代码+论文)
基于Qt的多线程局域网聊天系统（含客户端+服务端）.zip
07-30
毕业设计是高等教育阶段学生完成学业的一个重要环节，通常在学士或硕士学业即将结束时进行。这是学生将在整个学业中所学知识和技能应用到实际问题上的机会，旨在检验学生是否能够独立思考、解决问题，并展示其专业能力的一项综合性任务。 毕业设计的主要特点包括： 独立性： 毕业设计要求学生具备独立思考和解决问题的能力。学生需要选择一个合适的课题，研究相关文献，进行实地调查或实验，并提出独立见解。 实践性： 毕业设计是将理论知识应用到实际问题中的一次实践。通过完成毕业设计，学生能够将所学的专业知识转化为实际的解决方案，加深对专业领域的理解。 综合性： 毕业设计往往要求学生运用多个学科的知识，综合各种技能。这有助于培养学生的综合素养，提高他们的综合能力。 导师指导： 学生在毕业设计过程中通常由一名指导老师或导师团队提供指导和支持。导师负责引导学生确定研究方向、制定计划、提供建议，并在整个过程中监督进展。 学术规范： 毕业设计要求学生按照学术规范完成研究，包括文献综述、研究设计、数据采集与分析、结论和讨论等环节。学生需要撰写一篇完整的毕业论文，并进行答辩。
iperf3可测量丢包率
07-30
测吞吐量和丢包率
遗传算法(GA)优化长短期记忆网络的数据分类预测，GA-LSTM分类预测，多输入单输出 多特征输入单输出的二分类及多分类模型
07-30
遗传算法(GA)优化长短期记忆网络的数据分类预测，GA-LSTM分类预测，多输入单输出。 多特征输入单输出的二分类及多分类模型。程序内注释详细，直接替换数据就可以用。 程序语言为matlab，程序可出分类效果图，迭代优化图，混淆矩阵图
XOR方程组解法与ACM竞赛模板
"南阳理工acm常用模板，包含高斯消元方法的多种应用和解题策略" 在ACM竞赛编程中，高斯消元是一种常见的线性代数方法，用于解决一系列线性方程组。它在算法竞赛中尤其有用，尤其是在处理约束条件和优化问题时。本...
评论
添加红包

请填写红包祝福语或标题

个

红包个数最小为10个

元

红包金额最低5元

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

抵扣说明：

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

余额充值