社区讨论

萌新刚学OI,不是妹子,高斯消元似乎被卡精度求助

P2455[SDOI2006] 线性方程组参与者 8已保存回复 8

讨论操作

快速查看讨论及其快照的属性,并进行相关操作。

当前回复
8 条
当前快照
1 份
快照标识符
@mi7clvv8
此快照首次捕获于
2025/11/20 19:28
4 个月前
此快照最后确认于
2025/11/20 19:28
4 个月前
查看原帖
板子A过洛谷P3389 【模板】高斯消元法,复制过来加入无穷多解的判断就WA
最后两组数据是要输出无解但是我一直得出无穷多解
手动造数据测了好像也没什么问题QAQ
CPP
#include<iostream>
#include<cmath>
#include<algorithm>
#include<queue>
#include<cstring>
#include<cstdio>
using namespace std;
typedef double dd;

int read()
{
    int f=1,x=0;
    char ss=getchar();
    while(ss<'0'||ss>'9'){if(ss=='-')f=-1;ss=getchar();}
    while(ss>='0'&&ss<='9'){x=x*10+ss-'0';ss=getchar();}
    return f*x;
}

const int maxn=410;
const dd eps=1e-8;
int n;
dd a[maxn][maxn];
dd ans[maxn];

int gauss()
{
    for(int i=1;i<=n;++i)
    {
        int r=i;
        for(int j=i+1;j<=n;++j)
        if(fabs(a[j][i])>fabs(a[r][i])) r=j;
        
        if(fabs(a[r][i])<eps) return 0;
        if(r!=i) swap(a[i],a[r]);
        
        dd f=a[i][i];
        for(int j=i;j<=n+1;j++) a[i][j]/=f;
        for(int j=i+1;j<=n;j++)
        {
            f=a[j][i];
            for(int k=i;k<=n+1;k++)
            a[j][k]-=a[i][k]*f;
        }
    }
    
    int judge=0;
    for(int i=1;i<=n;++i) 
    {
        int cnt=0;
        for(int j=1;j<=n;++j)
        if(fabs(a[i][j])<eps) cnt++;
        
        if(cnt==n&&fabs(a[i][n+1])>eps) return -1;
        if(cnt==n) judge=1;
    }
    if(judge) return 0;
    
    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]);
    }
    return 1;
}

int main()
{
    n=read();
    for(int i=1;i<=n;++i)
    for(int j=1;j<=n+1;++j)
    scanf("%lf",&a[i][j]);
    
    int tt=gauss();
    if(tt!=1) printf("%d",tt);
    else{
    	for(int i=1;i<=n;++i)
        if(fabs(ans[i])<eps) printf("x%d=0\n",i);
        else printf("x%d=%.2lf\n",i,ans[i]);
    }
    return 0;
}

回复

8 条回复,欢迎继续交流。

正在加载回复...