萌新求助,朴素Gauss消元,WA on #8

P2455 [SDOI2006] 线性方程组

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 高斯消元当 a_{i, i} = 0 的时候,行不能跳啊,只有列可以加 1


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 直接跳过是有问题的啊,这一行的信息就没了 ……


| 下一页