思路是判定非唯一解时往后面列消元，和第一篇题解相似

#include <bits/stdc++.h>
using namespace std;
double a[55][55],b[55],c[55];
double exs=1e-8;
bool pd(double x,double y)
{
    if(x<y)swap(x,y);
    if((x-y)<exs)return 1;
    return 0;
}
signed main()
{
	int n;
    cin>>n;
    for(int i=1;i<=n;i++)
    {
    	for(int j=1;j<=n;j++)
    	 scanf("%lf",&a[i][j]);
    	scanf("%lf",&b[i]);
	}
	for(int i=1;i<=n;i++)
	{
	  int p=i;
	  while(pd(a[i][p],0))
	  {
	  	if(p>n)break;
		int j=i+1;
	  	while(pd(a[j][p],0)&&j<=n)j++;
		if(j>n){p++;continue;}
		for(int k=p;k<=n;k++)swap(a[i][k],a[j][k]);
	    swap(b[i],b[j]);
	  }
	  for(int j=i+1;j<=n;j++)
	  {
	  	double pp=a[j][p]/a[i][p];
	  	for(int k=p;k<=n;k++)
	  	a[j][k]-=pp*a[i][k];
	  	b[j]-=pp*b[i];
	  }
	} 
	bool dj=0;
	for(int i=1;i<=n;i++)
	{
	  bool p=0;
	  for(int j=1;j<=n;j++)
	  if(a[i][j])p=1;
	  if(!p)
	  {
	  	dj=1;
	  	if(b[i]){cout<<-1;return 0;}
	  }	
	}
	if(dj){cout<<0;return 0;}
    for(int i=n;i>=1;i--)
    {
		for(int j=n-i;j>=1;j--)
         b[i]-=a[i][n-j+1]*c[n-j+1];   	
    	c[i]=b[i]/a[i][i];
	}
	for(int i=1;i<=n;i++)
	{
		printf("x%d=",i);
	    printf("%.2lf\n",c[i]); 
	}
	return 0; 
}