顺序消去法
//顺序高斯消元
//拿列主元随便改的
int n;
double a[111][111],ans[111];
signed main()
{
cin>>n;
rpp(i,n) rpp(j,n+1) cin>>a[i][j];
for(int i=1;i<=n;++i)//枚举当前项
{
for(int j=i+1;j<=n;++j)
{
double tmp=a[j][i]/a[i][i];
for(int k=i;k<=n+1;++k) // 第i次迭代的时候前i-1项肯定是0,直接从i开始
a[j][k]-=a[i][k]*tmp;
}
// rpp(i,n)
// {
// rpp(j,n+1)
// {
// cout<<a[i][j]<<" ";
// }
// cout<<endl;
// }
// cout<<endl;
}
for(int i=n;i>=1;--i)
{
ans[i]=a[i][n+1];
for(int j=i+1;j<=n;++j)//他下面的都剪掉
ans[i]-=a[i][j]*ans[j];
ans[i]/=a[i][i];
}
rpp(i,n) printf("%.2f\n",ans[i]);
cout<<endl;
stop;
return 0;
}
列主元消除法
//列主元消除法(加减消元+代入消元)
//和顺序消除法相比,只多了一个选择每项最大主元
//和全主元消除法相比,差别在于全主元是在右下角的矩阵内选择最大值作为主元
//列主元只在该列寻找
//判断解的情况
//对于无解的判断:某一行前 nnn 个数均为 0,最后的结果却不为 000。
//对于无数解的判断:某一行 n+1n+1n+1 个数均为 000.
int n;
double a[111][111],ans[111];
signed main()
{
cin>>n;
rpp(i,n) rpp(j,n+1) cin>>a[i][j];
for(int i=1;i<=n;++i)//枚举当前项
{
//(1)找到主元所在行,换到当前行
int mx=i;
for(int j=i+1;j<=n;++j)//找到该项最大的那一个方程
if(fabs(a[j][i])>fabs(a[mx][i])) mx=j;
if(mx!=i) swap(a[mx],a[i]);//交换这两行
if(!a[i][i])
{
cout<<"No Solution"<<endl;
return 0;
}
//(2)将该行主元系数置为1
double tmp=a[i][i];
rpp(j,n+1) a[i][j]/=tmp;
//(3)将该行以下的该项系数全部变成0
for(int j=i+1;j<=n;++j)
{
tmp=a[j][i];
for(int k=i;k<=n+1;++k) // 第i次迭代的时候前i-1项肯定是0,直接从i开始
a[j][k]-=a[i][k]*tmp;
}
}
//(4)回代
for(int i=n;i>=1;--i)
{
ans[i]=a[i][n+1];
for(int j=i+1;j<=n;++j)//他下面的都剪掉
ans[i]-=a[i][j]*ans[j];
}
rpp(i,n) printf("%.2f\n",ans[i]);
cout<<endl;
stop;
return 0;
}
全主元消除法
//全主元消除法
// 和列主元相比只改了(1),然后加一个找pos的就行
int n;
int pos[111];//i列为第pos[i]项
double a[111][111],ans[111];
signed main()
{
cin>>n;
rpp(i,n) pos[i]=i;
rpp(i,n) rpp(j,n+1) cin>>a[i][j];
for(int i=1;i<=n;++i)//枚举当前项
{
//(1)在右下角中查找主元
int r=i,s=i;
for(int _=i+1;_<=n;++_)
for(int __=i+1;__<=n;++__)
if(fabs(a[_][__])>fabs(a[r][s])) r=_,s=__;
if(r!=i) swap(a[i],a[r]);//先交换行
if(s!=i)//然后交换列,同时更新pos数组
{
swap(pos[s],pos[i]);
rpp(_,n) swap(a[_][i],a[_][s]);
}
if(!a[i][i])
{
cout<<"No Solution"<<endl;
return 0;
}
//(2)将该行主元系数置为1
double tmp=a[i][i];
rpp(j,n+1) a[i][j]/=tmp;
//(3)将该行以下的该项系数全部变成0
for(int j=i+1;j<=n;++j)
{
tmp=a[j][i];
for(int k=i;k<=n+1;++k) // 第i次迭代的时候前i-1项肯定是0,直接从i开始
a[j][k]-=a[i][k]*tmp;
}
}
//(4)回代
for(int i=n;i>=1;--i)
{
ans[i]=a[i][n+1];
for(int j=i+1;j<=n;++j)//他下面的都剪掉
ans[i]-=a[i][j]*ans[j];
}
map<int,int>mp;
rpp(i,n) mp[pos[i]]=i;//第i项在那一列
rpp(i,n) printf("%.2f\n",ans[mp[i]]);
cout<<endl;
return 0;
}
高斯–约当消元
//列主元高斯--约当消元(精度最高的(相对地)
//缺点:只能判断是否有唯一解,无法判断究竟是无解还是无数解。
// 1.选择一个尚未被选过的未知数作为主元,选择包含这个主元的方程。
// 2.将这个方程主元的系数化为1。
// 3.通过加减消元,消掉其它方程中的这个未知数。
// 4.重复以上步骤,直到把所有式子变成形如: a1+b0+c*0……=d
double a[111][111];
signed main()
{
int n;cin>>n;
rpp(i,n) rpp(j,n+1) cin>>a[i][j];
rpp(i,n)//枚举当前要消除的项(列)
{
//(1)找到一个主元
int mx=i;
for(int j=i+1;j<=n;++j)//找到该项最大的那行方程//每个方程只能有一个主元/小于i的不用比了
if(fabs(a[j][i])>fabs(a[mx][i])) mx=j;
rpp(j,n+1) swap(a[mx][j],a[i][j]);//交换这两行
if(!a[i][i])//最大值等于0则说明该列都为0,肯定无解
{
puts("No Solution");
return 0;
}
//(2)将其他行中的该项消除
rpp(j,n)
{
if(j==i) continue;
double tmp=a[j][i]/a[i][i];
for(int k=i;k<=n+1;++k)//i列即之前肯定已经消完了
a[j][k]-=tmp*a[i][k];
}
// rpp(_,n)
// {
// rpp(__,n+1) printf("%6.2f",a[_][__]);
// cout<<endl;
// }
// cout<<endl;//输出中间值
}
//(3)输出最终结果
rpp(i,n) printf("%.2f\n",a[i][n+1]/a[i][i]);
stop;
return 0;
}