标签:int -- 110 printf include 高斯消
高斯消元
#include<iostream>
#include<cstdio>
#include<cmath>
using namespace std;
const int eps=1e-6;// 因为double储存有误差所以用eps;
int n;double a[110][110];
int gauss();
int main()
{
scanf("%d",&n);
for(int i=0;i<n;i++)
for(int j=0;j<n+1;j++)
scanf("%lf",&a[i][j]);
int t=gauss();
if(t==1){
printf("No Solution");
}else{
for(int i=0;i<n;i++){
printf("%0.2lf\n",a[i][n]);
}
}
///
return 0;
}
/
int gauss(){
int c,r;//C表示每一列 R表示每一行
for(c=0,r=0;c<n;c++){
int t=r;
for(int i=r;i<n;i++)
if(fabs(a[i][c])>fabs(a[t][c]))
t=i;// 每次都去找到最大的
if(fabs(a[t][c])<eps)continue;// 此处则说明该元系数为0;就可能存在吴姐或无数姐情况
/*
如果说三个方程;其中一个未知数的系数都为0;
以下两种可能:
1.舍去那个未知数,如果说对于剩下的两个未知数方程有唯一解则出现了0x=0;
就是无数姐;
2.舍去之后因为那啥吴姐,所以吴姐;
*/
for(int i=c;i<=n;i++)swap(a[t][i],a[r][i]);// 交换
for(int i=n;i>=c;i--)a[r][i]/=a[r][c];// 系数归一
for(int i=r+1;i<n;i++)
if(fabs(a[i][c])>eps)// 如果不为0
for(int j=n;j>=c;j--)
a[i][j]-=a[r][j] * a[i][c];// 那就变成0
r++;//下一行
}
if(r<n)return 1;// 吴姐或无穷解
for(int i=n-1;i>=0;i--)
for(int j=i+1;j<n;j++)
a[i][n]-=a[i][j]*a[j][n]; //结果是一个阶梯,根据梯形来处理;
return 0;
}
这个模板是徐倞舟老师教的;感谢徐倞老师;
标签:int,--,110,printf,include,高斯消 来源: https://blog.csdn.net/weixin_55270923/article/details/116749575
本站声明: 1. iCode9 技术分享网(下文简称本站)提供的所有内容,仅供技术学习、探讨和分享; 2. 关于本站的所有留言、评论、转载及引用,纯属内容发起人的个人观点,与本站观点和立场无关; 3. 关于本站的所有言论和文字,纯属内容发起人的个人观点,与本站观点和立场无关; 4. 本站文章均是网友提供,不完全保证技术分享内容的完整性、准确性、时效性、风险性和版权归属;如您发现该文章侵犯了您的权益,可联系我们第一时间进行删除; 5. 本站为非盈利性的个人网站,所有内容不会用来进行牟利,也不会利用任何形式的广告来间接获益,纯粹是为了广大技术爱好者提供技术内容和技术思想的分享性交流网站。