高斯消元(完全主元法 and 部分主元法) C++代码

部分主元法高斯消元

/*
算法步骤:1.枚举每一列,找到绝对值最大的一行2.将该行和第一行交换3.将该行行首置为一4.将下面所有行第 i 列置为零
*/#include <iostream>
#include <cmath>using namespace std;
const int N = 109;
const double eps = 1e-6;
double a[N][N];
int n;// 未知数个数, 方程数量
void init()
{// cout << "未知数的个数,方程的数量: ";cin >> n;for (int i = 0; i < n; i++)for (int j = 0; j < n + 1; j++)cin >> a[i][j];
}
bool Gauss()
{int l = 0, r = 0; //行 列// step 1 and 2for (l = 0, r = 0; r < n; r++){int tt = l;for (int i = l; i < n; i++)if (fabs(a[tt][r]) < fabs(a[i][r]))tt = i;if (fabs(a[tt][r]) < eps)	continue;for (int i = r; i < n + 1; i++)swap(a[tt][i], a[l][i]);	// step 3for (int i = n; i >= r; i--)a[l][i] /= a[l][r];// step 4 将下面所有行第 r 列置为零for (int i = l + 1; i < n; i++){if (fabs(a[i][r]) < eps)	continue;for (int j = n; j >= r; j--)a[i][j] -= a[l][j] * a[i][r];}l++;}for (int i = n - 1; i >= 0; i--)for (int j = i + 1; j < n; j++)a[i][n] -= a[j][n] * a[i][j];if (l == n)	return 0;return 1;
}
int main()
{init();if (!Gauss())for (int i = 0; i < n; i++)cout << "x" <<  i  << " = "<< a[i][n] << '\n';elsecout << "无解!!!\n";
}/*
3
1 1 1 6
0 4 -1 5
2 -2 1 1
*/

完全主元法高斯消元

#include <bits/stdc++.h>
using namespace std;
#define PII pair<int, int>
const int N = 109;
const double eps = 1e-6;
double a[N][N];
int id[N];
int n;void init()
{cin >> n;for (int i = 0; i < n; i++) id[i] = i;for (int i = 0; i < n; i++)for (int j = 0; j < n + 1; j++)cin >> a[i][j];
}
PII Get_Max_Idx(int l, int r)
{PII idx;double mx = 0;for (int i = l; i < n; i++)for (int j = r; j < n; j++)if (fabs(a[i][j]) > mx)mx = fabs(a[i][j]), idx = {i, j};return idx;
}
bool Gauss()
{int l = 0, r = 0;for (l = 0, r = 0; r < n; r++){PII tt = Get_Max_Idx(l, l);//行交换for (int i = r; i < n + 1; i++)swap(a[tt.first][i], a[l][i]);tt.first = l;//列交换for (int i = 0; i < n; i++)swap(a[i][tt.second], a[i][r]);swap(id[tt.second], id[r]);// step 3for (int i = n; i >= r; i--)a[l][i] /= a[l][r];// step 4 将下面所有行第 r 列置为零for (int i = l + 1; i < n; i++){if (fabs(a[i][r]) < eps)    continue;for (int j = n; j >= r; j--)a[i][j] -= a[l][j] * a[i][r];}l++;}for (int i = n - 1; i >= 0; i--)for (int j = i + 1; j < n; j++)a[i][n] -= a[j][n] * a[i][j];if (l == n) return 0;return 1;
}
int main()
{init();if (!Gauss()){for (int i = 0; i < n; i++)cout << "x" <<  id[i]  << " = "<< a[i][n] << '\n';}else    cout << "无解!!!\n";
}

本文来自互联网用户投稿,该文观点仅代表作者本人,不代表本站立场。本站仅提供信息存储空间服务,不拥有所有权,不承担相关法律责任。如若转载,请注明出处:http://www.mzph.cn/news/160670.shtml

如若内容造成侵权/违法违规/事实不符,请联系多彩编程网进行投诉反馈email:809451989@qq.com,一经查实,立即删除!

相关文章

Linux内核的内存管理

Linux内核源码内存管理主要包括以下几个部分&#xff1a; 1. 物理内存管理&#xff1a;这部分主要负责将物理内存划分为不同的页表项&#xff0c;以便操作系统能够快速地访问和操作内存。 2. 虚拟内存管理&#xff1a;这部分主要负责将用户空间的地址映射到物理内存中&#x…

linux之进程地址空间

文章目录 1.进程地址空间回顾1.1进程地址空间划分1.2验证进程地址空间划分1.简单划分2.完整划分 2.初探进程地址空间2.1初看现象2.2Makefile的简便写法 3.进程地址空间详解3.1地址空间是什么?3.2地址空间的设计/由来3.3空间区域划分3.4如何理解地址空间?3.5解释3.2的&#x1…

警惕.locked勒索病毒,您需要知道的预防和恢复方法。

尊敬的读者&#xff1a; 随着网络技术的进步&#xff0c;勒索病毒已经成为一种极具威胁性的网络犯罪工具之一。其中&#xff0c;.locked勒索病毒是一种采用高级加密算法的恶意软件&#xff0c;目的是加密用户的文件&#xff0c;并勒索赎金以提供解密密钥。本文将介绍如何应对被…

解决No Feign Client for loadBalancing defined,修改Maven依赖

Spring微服务报错&#xff1a; java.lang.IllegalStateException:FactoryBean threw exception on object creation; nested exception is java.lang.IllegalStateException: No Feign Client for loadBalancing defined. Did you forget to include spring-cloud-starter-netf…

你不知道的库:库的种类,作用和加载方式

你不知道的库&#xff1a;库的种类&#xff0c;作用和加载方式 &#x1f4df;作者主页&#xff1a;慢热的陕西人 &#x1f334;专栏链接&#xff1a;Linux &#x1f4e3;欢迎各位大佬&#x1f44d;点赞&#x1f525;关注&#x1f693;收藏&#xff0c;&#x1f349;留言 本博客…

组件化——组件的实现原理

渲染器主要负责将虚拟 DOM 渲染为真实 DOM&#xff0c;我们只需要使用虚拟 DOM 来描述最终呈现的内容即可。但当我们编写比较复杂的页面时&#xff0c;用来描述页面结构的虚拟 DOM 的代码量会变得越来越多&#xff0c;或者说页面模板会变得越来越大。这时&#xff0c;我们就需要…

iperf3 网络测试

iperf3 测试网络的上下行带宽 下载地址 https://iperf.fr/iperf-download.php 开启服务器 开启客户端 常用命令 -c 代表客户端-s 代表服务端-u 代表 udp-r 代表数据方向是否反向 https://baijiahao.baidu.com/s?id1731514357681464971&wfrspider&forpc

C++学习 --queue

目录 1&#xff0c; 什么是queue 2&#xff0c; 创建queue 2-1&#xff0c; 标准数据类型 2-2&#xff0c; 自定义数据类型 2-3&#xff0c; 其他创建方式 3&#xff0c; 操作stack 3-1&#xff0c; 赋值 3-2&#xff0c; 插入元素(push) 3-3&#xff0c; 查询元素 3…

Python简直是万能的,这5大主要用途你一定要知道!

从2015开始国内就开始慢慢接触Python了&#xff0c;从16年开始Python就已经在国内的热度更高了&#xff0c;目前也可以算的上"全民Python"了。 众所周知小学生的教材里面已经有Python了&#xff0c;国家二级计算机证也需要学习Python了&#xff01; 因为Python简单…

编程语言发展史:布尔代数和机器语言

布尔代数是一种数学理论&#xff0c;用于描述和分析逻辑和布尔值的关系。它是由英国数学家George Boole在19世纪中期发明的&#xff0c;被认为是现代计算机科学的基础之一。布尔代数的发明使得逻辑运算可以被表示为代数运算&#xff0c;从而为计算机科学的发展奠定了基础。 在…

PTA 7-4 数列求和-加强版

7-4 数列求和-加强版 分数 20 全屏浏览题目 作者 DS课程组 单位 浙江大学 给定某数字A&#xff08;1≤A≤9&#xff09;以及非负整数N&#xff08;0≤N≤100000&#xff09;&#xff0c;求数列之和SAAAAAA⋯AA⋯A&#xff08;N个A&#xff09;。例如A1, N3时&#xff0c;S1…

Unity、UE和Godot的优劣对比

先占位。。。。。。 首先说Unity和UE这两家公司&#xff0c;是行业的两座灯塔&#xff0c;对整个游戏引擎的这个行业的发展具有这种指导性的这种作作用。这两个引擎我从2016年开始就一直在用&#xff0c;结合一下业内的共识&#xff0c;一般来说认为呢&#xff0c;Unity更擅长移…

2023全球边缘计算大会深圳站-核心PPT资料下载

一、峰会简介 边缘计算&#xff0c;是指在靠近物或数据源头的一侧&#xff0c;采用网络、计算、存储、应用核心能力为一体的开放平台&#xff0c;就近提供最近端服务。其应用程序在边缘侧发起&#xff0c;产生更快的网络服务响应&#xff0c;满足行业在实时业务、应用智能、安…

LeetCode算法题解(动态规划,背包问题)|LeetCode416. 分割等和子集

LeetCode416. 分割等和子集 题目链接&#xff1a;416. 分割等和子集 题目描述&#xff1a; 给你一个 只包含正整数 的 非空 数组 nums 。请你判断是否可以将这个数组分割成两个子集&#xff0c;使得两个子集的元素和相等。 示例 1&#xff1a; 输入&#xff1a;nums [1,5,…

Linux中的进程程序替换

Linux中的进程程序替换 1. 替换原理2. 替换函数3. 函数解释4. 命名理解程序替换的意义 1. 替换原理 替换原理 用fork创建子进程后执行的是和父进程相同的程序(但有可能执行不同的代码分支),子进程往往要调用一种exec函数以执行另一个程序。当进程调用一种exec函数时,该进程的…

[Docker]九.Docker compose讲解

docker-compose 是 docker 官方的一个开源项目&#xff0c;可以实现对 docker 容器集群的快速编排, docker-compose 通过一个 配置文件 来管理多个 Docker 容器,在配置文件中&#xff0c;所有的容器通过 services 来定义&#xff0c;然后使用 docker-compose脚本 来 启动&am…

Nuxt.js Next.js Nest.js

Nuxt.js和Next.js都是服务端渲染框架(SSR)&#xff0c;属于前端框架,Nest.js则是node框架,属于后端框架。 其中Nuxt.js是vue的ssr框架&#xff0c;Next.js是react的ssr框架。 都是比vue和react更上层的前端框架。 文章目录 1.SSR2.Nuxt2.1 Nuxt的下载2.2 Nuxt的集成2.3 Nuxt…

HuggingFace-利用BERT预训练模型实现中文情感分类(下游任务)

准备数据集 使用编码工具 首先需要加载编码工具&#xff0c;编码工具可以将抽象的文字转成数字&#xff0c;便于神经网络后续的处理&#xff0c;其代码如下&#xff1a; # 定义数据集 from transformers import BertTokenizer, BertModel, AdamW # 加载tokenizer token Ber…

cobol基本动词

cobol基本动词 基本动词用于过程部中的数据处理。每个语句总是以cobol动词开头。 input&#xff08;输入&#xff09;/output&#xff08;输出&#xff09; 输入输出动词用于从用户获取数据。并显示cobol程序的输出。 accept 用于从操作系统或者用户获取数据&#xff0c;例如日…

langchain 部署组件-LangServe

原文&#xff1a;&#x1f99c;️&#x1f3d3; LangServe | &#x1f99c;️&#x1f517; Langchain LangServe &#x1f6a9; We will be releasing a hosted version of LangServe for one-click deployments of LangChain applications. Sign up here to get on the wa…