suyiheng @ 2019-07-11 20:24:51
RT,高斯消元WA90,求大佬帮忙查错QAQ
#include<iostream>
#include<cstdio>
#include<algorithm>
#include<cmath>
using namespace std;
struct data{
double f[101];
}t[101];
int n,q[101],cnt=0;
double eps=0.0000001;
int main(){
scanf("%d",&n);
for(int i=0;i<n;i++)for(int j=0;j<=n;j++)scanf("%lf",&t[i].f[j]);
for(int i=0;i<n;i++){
int p=i;
for(int j=i+1;j<n;j++)if(fabs(t[j].f[i])-fabs(t[p].f[i])>eps)p=j;
swap(t[p],t[i]);
if(fabs(t[i].f[i])<eps){
cnt++;
q[cnt]=i;
continue;
}
for(int j=i+1;j<=n;j++)t[i].f[j]/=t[i].f[i];
for(int j=0;j<n;j++)if(i!=j)for(int k=i+1;k<=n;k++)t[j].f[k]-=t[j].f[i]*t[i].f[k];
}
if(cnt==0){
for(int i=0;i<n;i++){
if(fabs(t[i].f[n])<eps)printf("x%d=0\n",i+1);
else printf("x%d=%.2lf\n",i+1,t[i].f[n]);
}
}else{
for(int i=1;i<cnt;i++){
if(fabs(t[q[i]].f[n])>eps){
printf("-1");
return 0;
}
}
printf("0");
return 0;
}
}
by Gary818 @ 2019-07-11 20:31:00
@suyiheng
第一你不是蒟蒻,第二你不是刚学,第三我猜你不是妹子,所以,不予帮助
我才不会告诉你我不会
by Smile_Cindy @ 2019-07-11 20:56:01
第一你不是蒟蒻,第二你不是刚学,第三我猜你不是妹子,所以,不予帮助
by LightningUZ @ 2019-08-11 20:21:54
捕捉Alpha%%%
by LightningUZ @ 2019-08-11 20:22:11
第一你不是蒟蒻,第二你不是刚学,第三我猜你不是妹子,所以,不予帮助
by cavalry @ 2019-08-11 20:26:17
捕捉女装大佬Alpha%%%