蒟蒻刚学OI,你猜是不是妹子,求助!

P2455 [SDOI2006] 线性方程组

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%%%


|