高斯消元模板
来源:互联网 发布:淘宝宝贝描述低能买吗 编辑:程序博客网 时间:2024/06/05 21:08
#include<bits/stdc++.h>using namespace std;#define eps 1e-9const int MAXN=220;const int N=3;double a[MAXN][MAXN],x[MAXN];//方程的左边的矩阵和等式右边的值,求解之后x存的就是结果int equ,var;//方程数和未知数个数/**返回0表示无解,1表示有解*/int Gauss()//高斯消元 { int i,j,k,col,max_r; for(k=0,col=0;k<equ&&col<var;k++,col++) { max_r=k; for(i=k+1;i<equ;i++) if(fabs(a[i][col])>fabs(a[max_r][col])) max_r=i; if(fabs(a[max_r][col])<eps)return 0; if(k!=max_r) { for(j=col;j<var;j++) swap(a[k][j],a[max_r][j]); swap(x[k],x[max_r]); } x[k]/=a[k][col]; for(j=col+1;j<var;j++)a[k][j]/=a[k][col]; a[k][col]=1; for(i=0;i<equ;i++) if(i!=k) { x[i]-=x[k]*a[i][k]; for(j=col+1;j<var;j++)a[i][j]-=a[k][j]*a[i][col]; a[i][col]=0; } } return 1;}int main(){while(1){equ=var=N;for(int i=0;i<N;i++){for(int j=0;j<N;j++) {scanf("%lf",&a[i][j]);}scanf("%lf",&x[i]);}if(Gauss()) for(int i=0;i<N;i++) printf("x[%d]=%lf\n",i,x[i]);else printf("无解\n");}}
阅读全文
0 0
- 高斯消元模板
- 高斯消元模板
- 高斯消元模板
- 高斯消元模板
- 高斯消元模板
- 高斯消元模板
- 高斯消元模板
- 高斯消元模板
- 高斯消元模板
- 高斯消元模板
- 高斯消元模板
- 高斯消元模板
- 高斯消元模板
- 高斯消元 【模板】
- 高斯消元模板
- 高斯消元模板
- 高斯消元模板
- 高斯消元模板
- 【学习笔记一】- JS 类型识别 和 deepclone 深复制
- 去除List集合中的重复值(四种好用的方法)
- Android APP上线必不可少的应用市场!
- 小练习1
- 超级记忆法(4)——第二小时
- 高斯消元模板
- 随便写点
- 搭建ACE及TAO开发环境的快速方法 (Windows篇)
- Idea中项目文件管理
- JS、JSP、EL 内置对象简单整理
- 多校10_1010 hdu6180 贪心+STL:multiset
- java面试题<1>
- 【Linux学习笔记】11:其它帮助命令
- 2017 多校训练第十场 HDU 6178 Monkeys 最小点覆盖+树形DP