Focus_on @ 2022-04-20 18:32:04
rt,所有讨论区的hack好像都过了,就是不知道第八个点为什么不对
#include<cstdio>
#include<cstring>
#include<algorithm>
#include<cmath>
using namespace std;
const int N=61;
const int inf=1e9;
const double eps=1e-8;
double a[N][N],b[N],ans[N];
int n,fl;
int main(){
scanf("%d",&n);
for(int i=1;i<=n;i++){
for(int j=1;j<=n;j++) scanf("%lf",&a[i][j]);
scanf("%lf",&b[i]);
}
fl=-1;
for(int i=1;i<=n;i++){
int now=i;
for(int j=i+1;j<=n;j++)
if(a[j][i]>a[i][i]) now=j;
for(int j=1;j<=n;j++) swap(a[i][j],a[now][j]);
swap(b[i],b[now]);
if(abs(a[i][i])<eps){
// fl=i;
continue;
}
for(int j=i+1;j<=n;j++){
double p=a[j][i]/a[i][i];
for(int k=i;k<=n;k++)
a[j][k]-=a[i][k]*p;
b[j]-=b[i]*p;
}
}
for(int i=n;i>=1;i--){
if(abs(a[i][i])<eps){
if(abs(b[i])>eps){
puts("-1");
return 0;
}
else{
puts("0");
return 0;
}
}
ans[i]=b[i]/a[i][i];
for(int j=1;j<i;j++)
b[j]-=a[j][i]*ans[i];
}
for(int i=1;i<=n;i++)
printf("x%d=%.2lf\n",i,ans[i]);
return 0;
}
by Focus_on @ 2022-04-20 18:36:08
lz是辣鸡whk狗,很可能要到周末才能回复。恳请各位大佬帮忙看看啊!!!!!!!
by Alex_Wei @ 2022-04-20 18:36:39
看返回信息应该是精度问题。
by Focus_on @ 2022-04-20 18:37:26
@Alex_Wei 那请问除了eps外还有什么解决方法吗/kel
by Alex_Wei @ 2022-04-20 18:39:12
你的高斯消元写得好像有问题。
by Alex_Wei @ 2022-04-20 18:41:36
@Focus_on 高斯消元当
by Alex_Wei @ 2022-04-20 18:42:36
诶我想一想,好像这里没有问题。
by Focus_on @ 2022-04-20 18:43:47
感觉上跳行好像没问题??/kel
by Alex_Wei @ 2022-04-20 18:45:59
好怪哦
by Focus_on @ 2022-04-20 18:48:47
煞笔楼主要卷whk了,对大神不胜感激~~
周末再看
by Alex_Wei @ 2022-04-20 19:03:23
@Focus_on 直接跳过是有问题的啊,这一行的信息就没了 ……