社区讨论
Wa 7 9 11 12 13 80pts
P2455[SDOI2006] 线性方程组参与者 1已保存回复 0
讨论操作
快速查看讨论及其快照的属性,并进行相关操作。
- 当前回复
- 0 条
- 当前快照
- 1 份
- 快照标识符
- @mhjgwbhr
- 此快照首次捕获于
- 2025/11/04 02:22 4 个月前
- 此快照最后确认于
- 2025/11/04 02:22 4 个月前
CPP
#include<iostream>
#include<cstdio>
#include<cmath>
using namespace std;
int n;
double a[1010][1010];
void swaper(int i,int j){
for(int k=1;k<=n+1;k++){
swap(a[i][k],a[j][k]);
}
}
void simpler(int i,double d){
for(int k=1;k<=n+1;k++){
a[i][k]=a[i][k]*1.0/d;
}
}
void minuser(int i,double d,int j){
for(int k=1;k<=n+1;k++){
a[j][k]-=a[i][k]*1.0*d;
}
}
int gause(){
for(int i=1;i<=n;i++){
//cout<<i;
int mx=i;
for(int j=1;j<i;j++){
if(fabs(a[j][j])<1e-10&&fabs(a[j][i]>=1e-10)){
if(fabs(a[j][i])>fabs(a[mx][i])){
mx=j;
}
}
}
for(int j=i+1;j<=n;j++){
if(fabs(a[j][i])>fabs(a[mx][i])){
mx=j;
}
}
//cout<<mx;
swaper(i,mx);
//cout<<'y';
if(fabs(a[i][i])<1e-10){
return 0;
}
simpler(i,a[i][i]);
for(int j=1;j<=n;j++){
if(j==i){
continue;
}
minuser(i,a[j][i]/a[i][i],j);
}
}
return 1;
}
int main(){
cin>>n;
for(int i=1;i<=n;i++){
for(int j=1;j<=n+1;j++){
cin>>a[i][j];
}
}
if(!gause()){
int ans=1;
for(int i=n;i>=1;i--){
if(fabs(a[i][i])<1e-10){
int ok=0;
for(int j=1;j<=n;j++){
if(fabs(a[i][j])>=1e-10){
ok=1;
break;
}
}
if(!ok&&fabs(a[i][n+1])>=1e-10){
cout<<-1;
return 0;
}
}
}
for(int i=n;i>=1;i--){
if(fabs(a[i][i])<1e-10){
int ok=0;
for(int j=1;j<=n;j++){
if(fabs(a[i][j])>=1e-10){
ok=1;
break;
}
}
if(fabs(a[i][n+1])<1e-10){
cout<<0;
return 0;
}
}
}
}
else{
for(int i=1;i<=n;i++){
printf("x%d=%.2lf\n",i,a[i][n+1]);
}
}
return 0;
}
回复
共 0 条回复,欢迎继续交流。
正在加载回复...