社区讨论
为什么输出错还有80分
P2455[SDOI2006] 线性方程组参与者 5已保存回复 4
讨论操作
快速查看讨论及其快照的属性,并进行相关操作。
- 当前回复
- 4 条
- 当前快照
- 1 份
- 快照标识符
- @mhj01vyt
- 此快照首次捕获于
- 2025/11/03 18:30 4 个月前
- 此快照最后确认于
- 2025/11/03 18:30 4 个月前
这是我交了80pts的代码:
CPP#include<bits/stdc++.h>
using namespace std;
const int N=100+8;
const double eps=1e-7;
int n;
double a[N][N];
double ans[N];
double dabs(double x)
{
return (x<0?-x:x);
}
int main(){
cin>>n;
for(int i=1;i<=n;i++)
for(int j=1;j<=n+1;j++)
cin>>a[i][j];
for(int i=1;i<=n;i++)
{
int id=i;
for(int j=i;j<=n;j++)
{
if(dabs(a[id][i])<dabs(a[j][i]))
id=j;
}
if(dabs(a[id][i])<eps)
{
if(dabs(a[id][n+1])<eps)
{
cout<<0;
return 0;
}
cout<<-1;
return 0;
}
if(i!=id)
swap(a[i],a[id]);
double div=a[i][i];
for(int j=i;j<=n+1;j++)
a[i][j]/=div;
for(int j=i+1;j<=n;j++)
{
div=a[j][i];
for(int k=i;k<=n+1;k++)
a[j][k]-=a[i][k]*div;
}
}
ans[n]=a[n][n+1];
for(int i=n-1;i>=1;i--)
{
ans[i]=a[i][n+1];
for(int j=i+1;j<=n;j++)
{
ans[i]-=a[i][j]*ans[j];
}
}
for(int i=1;i<=n;i++)
{
cout<<"x1="<<fixed<<setprecision(2)<<ans[i]<<'\n';
}
return 0;
}
回复
共 4 条回复,欢迎继续交流。
正在加载回复...