博客
关于我
强烈建议你试试无所不能的chatGPT,快点击我
[Luogu] 高斯消元法
阅读量:5042 次
发布时间:2019-06-12

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

模拟,消元

#include 
#define DB doubleconst int N = 110;const DB eps = 1e-7;DB A[N][N], Answer[N];int n;DB read() {DB a; scanf("%lf", &a); return a;}int main() { std:: cin >> n; for(int i = 1; i <= n; i ++) for(int j = 1; j <= n + 1; j ++) A[i][j] = read(); for(int i = 1; i <= n; i ++) { int r = i; for(int j = i + 1; j <= n; j ++) if(abs(A[j][i]) > abs(A[r][i])) r = i; if(abs(A[r][i]) <= eps) {std:: cout << "No Solution"; return 0;} if(i != r) std:: swap(A[i], A[r]); DB now = A[i][i]; for(int j = i; j <= n + 1; j ++) A[i][j] /= now; for(int j = i + 1; j <= n; j ++) { DB Now = A[j][i]; for(int k = i; k <= n + 1; k ++) A[j][k] -= Now * A[i][k]; } } Answer[n] = A[n][n + 1]; for(int i = n - 1; i >= 1; i --) { Answer[i] = A[i][n + 1]; for(int j = n; j > n - (n - i); j --) Answer[i] -= Answer[j] * A[i][j]; } for(int i = 1; i <= n; i ++) printf("%.2lf\n", Answer[i]); return 0;}/*31 -2 3 64 -5 6 127 -8 10 21*/

 

转载于:https://www.cnblogs.com/shandongs1/p/9074109.html

你可能感兴趣的文章
模块与包
查看>>
mysql忘记root密码
查看>>
apache服务器中设置目录不可访问
查看>>
嵌入式Linux驱动学习之路(十)字符设备驱动-my_led
查看>>
【NOIP模拟】密码
查看>>
java容器---------手工实现Linkedlist 链表
查看>>
three.js 性能优化的几种方法
查看>>
《梦断代码》读书笔记(三)
查看>>
FreeMarker解析json数据
查看>>
Java8 Lambda表达应用 -- 单线程游戏server+异步数据库操作
查看>>
次序+“选择不重复的记录”(3)——最大记录
查看>>
Codeforces 450 C. Jzzhu and Chocolate
查看>>
[Unity3D]Unity3D游戏开发MatchTarget的作用攀登效果实现
查看>>
ACdream 1115 Salmon And Cat (找规律&amp;&amp;打表)
查看>>
JSON、JSONP、Ajax的区别
查看>>
AngularJS学习篇(一)
查看>>
【转载】 IP实时传输协议RTP/RTCP详解
查看>>
关于Xshell无法连接centos6.4的问题
查看>>
Linux系统的数据写入机制--延迟写入
查看>>
css3动画——基本准则
查看>>