码农知识堂 - 1000bd
  •   Python
  •   PHP
  •   JS/TS
  •   JAVA
  •   C/C++
  •   C#
  •   GO
  •   Kotlin
  •   Swift
  • 高斯消元


    1. #include
    2. #define IOS ios::sync_with_stdio(0);cin.tie(0);cout.tie(0);
    3. #define endl '\n'
    4. using namespace std;
    5. typedef pair<int, int> PII;
    6. typedef long long ll;
    7. typedef long double ld;
    8. const int N = 110;
    9. const double eps = 1e-6;
    10. int n;
    11. double a[N][N];
    12. int gauss()
    13. {
    14. int r, c;//行 列
    15. for(c = 0, r = 0; c < n; c ++)
    16. {
    17. int t = r;
    18. for(int i = r + 1; i < n; i ++)
    19. {//找到最大的行
    20. if(fabs(a[i][c]) > fabs(a[t][c]))
    21. {
    22. t = i;
    23. }
    24. }
    25. if(abs(a[t][c]) < eps)continue;
    26. for(int j = c; j <= n; j ++)
    27. {//与最大行交换
    28. swap(a[t][j], a[r][j]);
    29. }
    30. //第一个数变为1
    31. for(int j = n; j >= c; j --)a[r][j] /= a[r][c];
    32. for(int i = r + 1; i < n; i ++)
    33. {
    34. for(int j = n; j >= c; j --)
    35. {
    36. a[i][j] -= a[i][c] * a[r][j];
    37. }
    38. }
    39. r ++;
    40. }
    41. if(r < n)
    42. {
    43. for(int i = r; i < n; i ++)
    44. if(abs(a[i][n]) > eps)
    45. return 2;
    46. return 1;
    47. }
    48. for(int i = n - 1; i >= 0; i --)
    49. {
    50. for(int j = i + 1; j < n; j ++)
    51. {
    52. a[i][n] -= a[j][n] * a[i][j];
    53. }
    54. }
    55. return 0;
    56. }
    57. int main()
    58. {
    59. scanf("%d", &n);
    60. for(int i = 0; i < n; i ++)
    61. for(int j = 0; j <= n; j ++)
    62. scanf("%lf", &a[i][j]);
    63. int t = gauss();
    64. if(t == 0)
    65. {
    66. for(int i = 0; i < n; i ++)
    67. {
    68. if(fabs(a[i][n]) < eps)a[i][n] = 0;//防止出现-0.00
    69. printf("%.2lf\n", a[i][n]);
    70. }
    71. }
    72. else if(t == 1)printf("Infinite group solutions\n");
    73. else printf("No solution\n");
    74. return 0;
    75. }

  • 相关阅读:
    Spring Cloud Alibaba —— 服务注册与配置中心
    新书《深度学习与交通大数据实战》代码与数据等资源获取声明
    干货 | 测试人职场晋升“潜规则”:15 年经验资深测试经理的职场忠告
    在Postgres中分页的五种方法,从基本到异国情调
    VB开发OCX控件的属性之自定义列表项
    SpringBoot整合Redis
    大数据开发面试笔记Day1
    全国程序员高考卷曝光,你能答对几题?
    1.Linux的目录结构
    English语法_介词 - in
  • 原文地址:https://blog.csdn.net/a1695574118/article/details/133211309
  • 最新文章
  • 【JVM】编译执行与解释执行的区别是什么?JVM 使用哪种方式?
    用 Hashids 优雅解决 C 端自增 ID 暴露问题
    V8引擎 精品漫游指南--Ignition篇(上) 指令 栈帧 槽位 调用约定 内存布局 基础内容
    LLVM Pass快速入门(四):代码插桩
    milkup:桌面端 markdown AI续写和即时渲染
    基于项目工程构建SBOM(软件物料清单)的研究
    鸿蒙应用开发UI基础第二节:鸿蒙应用程序框架核心解析与实操
    .NET 中如何快速实现 List 集合去重?
    扣子Coze实战:从0到1打造抖音+小红书热点监控智能体
    浅谈数据访问层
  • 热门文章
  • 十款代码表白小特效 一个比一个浪漫 赶紧收藏起来吧!!!
    奉劝各位学弟学妹们,该打造你的技术影响力了!
    五年了,我在 CSDN 的两个一百万。
    Java俄罗斯方块,老程序员花了一个周末,连接中学年代!
    面试官都震惊,你这网络基础可以啊!
    你真的会用百度吗?我不信 — 那些不为人知的搜索引擎语法
    心情不好的时候,用 Python 画棵樱花树送给自己吧
    通宵一晚做出来的一款类似CS的第一人称射击游戏Demo!原来做游戏也不是很难,连憨憨学妹都学会了!
    13 万字 C 语言从入门到精通保姆级教程2021 年版
    10行代码集2000张美女图,Python爬虫120例,再上征途
小工具 小游戏
Copyright © 2022 侵权请联系2656653265@qq.com    京ICP备2022015340号-1

京公网安备 11010502049817号