社区讨论
萌新刚学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 条回复,欢迎继续交流。
正在加载回复...