我的做法是判断同构方程来判断无解和无数解

但是#8和#10一直过不了

翻了很多帖子，全部的hack都没有hack掉

在线求大佬解答。

%%%谢谢大佬storz

code:

#include<cmath>
#include<cstdio>
#include<iostream>
#include<algorithm>
#define int long long
#define ld long double
using namespace std;
const ld eps=1e-7;
const int N=111,p=998244353;
int n;
bool vis[N];
ld a[N][N],b[N];
int A[N][N],B[N];
int fpow(int x,int y){
    int s=1;
    for(;y;(y>>=1)){
        if(y&1)s=(s*x)%p;
        x=(x*x)%p;
    }
    return s;
} 
void Gause_mod(int n){
    for(int i=1;i<=n;++i){
        //因为这是在%p的意义下，所以不会出现精度问题
        //又因为第i个方程的第i项绝对不是0，所以不需要交换
        int inv=fpow(A[i][i],p-2);
        for(int j=1;j<=n;++j){
            if(i==j)continue;
            int tt=A[j][i]*inv%p;
            for(int k=i;k<=n+1;++k)A[j][k]=(A[j][k]-A[i][k]*tt%p+p)%p;
        }
    }
    for(int i=1;i<=n;++i)B[i]=A[i][n+1]*fpow(A[i][i],p-2)%p;
    return;
}
void Gause(int n){
    for(int i=1;i<=n;++i){
        int r=i;
        for(int j=i+1;j<=n;++j)if(fabs(a[r][i])<fabs(a[j][i]))r=j;
//      if(fabs(a[r][i])<eps)return;
        if(i!=r)swap(a[i],a[r]);
        ld 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[j][k]-a[i][k]*div);
        }
    }
    b[n]=a[n][n+1];
    for(int i=n-1;i>=1;--i){
        for(int j=i+1;j<=n;++j){
            a[i][n+1]-=a[i][j]*a[j][n+1];
        }
        b[i]=a[i][n+1];
    }
    return;
}
int pd(int i,int j){
    bool ok=true;
    for(int k=1;k<=n;++k){
        if(A[i][k]!=A[j][k]){
            ok=false;
            break;
        }
    }
    if(ok==false)return 0;
    if(B[i]==0||B[j]==0)for(;;);
    if((A[i][n+1]==A[j][n+1])&&(B[i]==B[j]))return 1;
    return 2;
}
signed main(){
//  printf("%lld\n",__gcd(0,0));
    scanf("%lld",&n);
    for(int i=1;i<=n;++i){
        for(int j=1;j<=n+1;++j){
            scanf("%lld",A[i]+j);
            a[i][j]=(long double)(A[i][j]);
        }
    }
    bool many=false;
    for(int i=1;i<=n;++i){
        int g=A[i][1];
        for(int j=2;j<=n;++j)g=__gcd(A[i][j],g);
        if(g!=0)for(int j=1;j<=n;++j)A[i][j]/=g;
        if(g==0&&A[i][n+1]!=0){
            puts("-1");
            return 0;
        }
        if(g==0&&A[i][n+1]==0){
            many=true;
            vis[i]=true;
        }
        else{
            B[i]=g;
            int di=__gcd(B[i],A[i][n+1]);
            B[i]/=di,A[i][n+1]/=di;
        }
    }
    for(int i=1;i<n;++i){
        if(vis[i])continue;
        for(int j=i+1;j<=n;++j){
            if(vis[j])continue;
            int sao=pd(i,j);
            if(sao==1)many=true;
            if(sao==2){
                puts("-1");
                return 0;
            }
        }
    }
    if(many){
        puts("0");
        return 0;
    }
    Gause(n);
    for(int i=1;i<=n;++i){
        if(fabs(b[i])<eps)printf("x%lld=0.00\n",i);
        else printf("x%lld=%.2Lf\n",i,b[i]);
    }
    return 0;
}
/*
3
1 2 1 1
3 6 3 3
2 4 2 3

3
-2 -1 1 1
0 1 -1 5
0 1 1 0
*/

最后再次表达深深的谢意！