博客
关于我
强烈建议你试试无所不能的chatGPT,快点击我
AcWing - 高斯消元解线性方程组(高斯消元)
阅读量:2000 次
发布时间:2019-04-28

本文共 2034 字,大约阅读时间需要 6 分钟。

题目链接:

时/空限制:1s / 64MB

题目描述

输入一个包含n个方程n个未知数的线性方程组。

方程组中的系数为实数。

求解这个方程组。

下图为一个包含m个方程n个未知数的线性方程组示例:

输入格式

第一行包含整数n。

接下来n行,每行包含n+1个实数,表示一个方程的n个系数以及等号右侧的常数。

输出格式

如果给定线性方程组存在唯一解,则输出共n行,其中第i行输出第i个未知数的解,结果保留两位小数。

如果给定线性方程组存在无数解,则输出“Infinite group solutions”。

如果给定线性方程组无解,则输出“No solution”。

数据范围

1≤n≤100,

所有输入系数以及常数均保留两位小数,绝对值均不超过100。

输入样例

3

1.00 2.00 -1.00 -6.00
2.00 1.00 -3.00 -9.00
-1.00 -1.00 2.00 7.00

输出样例

1.00

-2.00
3.00

解题思路

题意:给你一个n个方程,n个未知数,求解这个方程组。

思路:利用高斯消元法,详见。

Accepted Code:

/*  * @Author: lzyws739307453  * @Language: C++  */#include 
using namespace std;const int MAXN = 105;const double eps = 1e-8;double a[MAXN][MAXN];int Select(double a[][MAXN], int n) { int r = 1; for (int c = 1; c <= n; c++) { int t = r; for (int i = r + 1; i <= n; i++) if (fabs(a[i][c]) > fabs(a[t][c])) t = i; if (fabs(a[t][c]) < eps)//无解或无穷多解 continue; if (t != r) for (int i = c; i <= n + 1; i++) swap(a[r][i], a[t][i]); for (int i = n + 1; i >= c; i--) a[r][i] /= a[r][c]; for (int i = r + 1; i <= n; i++) if (fabs(a[i][c]) > eps) for (int j = n + 1; j >= c; j--) a[i][j] -= a[r][j] * a[i][c]; r++; } return r;}int Gauss(double a[][MAXN], int n) { int r = Select(a, n); if (r <= n) { for (int i = r; i <= n; i++) if (fabs(a[i][n + 1]) > eps) return 0;//无解 return 2;//无穷多解 } for (int i = n; i >= 1; i--) for (int j = i + 1; j <= n; j++) a[i][n + 1] -= a[i][j] * a[j][n + 1]; return 1;//有唯一解}int main() { int n; scanf("%d", &n); for (int i = 1; i <= n; i++) for (int j = 1; j <= n + 1; j++) scanf("%lf", &a[i][j]); int Judge = Gauss(a, n); if (!Judge) printf("No solution\n"); else if (Judge > 1) printf("Infinite group solutions\n"); else { for (int i = 1; i <= n; i++) printf("%.2f\n", a[i][n + 1]); } return 0;}

转载地址:http://qubtf.baihongyu.com/

你可能感兴趣的文章
人工智能为什么这么火?看看安防江湖30年血战就知道了
查看>>
“前端智能为安防产生新的数据价值”
查看>>
高斯混合模型
查看>>
(8)CMake入门笔记--CMake语法
查看>>
3D点云图实验
查看>>
头文件中 #ifndef---#define---#endif的作用
查看>>
分析Linux内核启动过程:从start_kernel到init
查看>>
系统调用过程的理解
查看>>
Ant内置任务之whichresource
查看>>
Ant内置任务之symlink
查看>>
jface databinding:部分实现POJO对象的监测
查看>>
jface databinding:构建一个改进版的通用型数值验证器StringToNumberValidator
查看>>
深入理解python--线程、进程与协程(1)
查看>>
Java--流重点总结初稿
查看>>
Html2Servlet--Html代码转换为Servlet小程序
查看>>
HTTP认证方式
查看>>
JavaWeb:HttpServletResponse和HttpServletRequest
查看>>
图书商城:分类模块
查看>>
图书商城:订单模块
查看>>
开源全能播放器Vitamio的使用
查看>>