高斯消元-浮点数

#include<iostream>
#include<stdio.h>
#include<string.h>
#include<stack>
#include<queue>
#include<map>
#include<set>
#include<vector>
#include<math.h>
#include<algorithm>
using namespace std;
typedef long long ll;
const double eps=1e-7;
const int maxn=500;
double a[maxn][maxn];
double x[maxn];//解集
bool free_x[maxn];//判断是否是不确定的变元
int equ,var;//equ个方程var个未知数，
//增广矩阵行数为equ,分别为0-equ-1;
//列数为var+1,分别为0-var
int gauss()
{
    int max_r,col,k,i,j;
    for(k=0,col=0;k<equ&&col<var;k++,col++)
    {
        max_r=k;
        for(i=k+1;i<equ;i++)
        {
            if((fabs(a[i][col])-fabs(a[max_r][col]))>eps)
                max_r=i;
        }
        if(fabs(a[k][col])<eps)
        {
            k--;
            continue;
        }
        if(max_r!=k)
        {
            for(j=col;j<=var;j++)
                swap(a[k][j],a[max_r][j]);
        }
        for(i=k+1;i<equ;i++)
        {
            if(fabs(a[i][col])>eps)
            {
                double temp=a[i][col]/a[k][col];
                a[i][col]=0.0;
                for(j=col;j<=var;j++)
                    a[i][j]-=a[k][j]*temp;
            }
        }
    }
    for(i=k;i<equ;i++)
        if(fabs(a[i][col])>eps)
        return -1;//无解
    if(k<var)//多解
    {
           int free_x_num;
           int free_index;
           double temp;
           for(i=k-1;i>=0;i--)
           {
               free_x_num=0;
               for(j=0;j<var;j++)
               {
                   if(free_x[j]&&fabs(a[i][j])>eps)
                   {
                       free_x_num++;
                       free_index=j;
                   }
               }
               if(free_x_num>1)
                continue;
               temp=a[i][var];
               for(j=0;j<var;j++)
               {
                   if(fabs(a[i][j])>eps&&j!=free_index)
                    temp-=a[i][j]*x[j];
               }
               x[free_index]=temp/a[i][free_index];
               free_x[free_index]=0;
           }
    }
    for(i=var-1;i>=0;i--)
    {
        double temp=a[i][var];
        for(j=i+1;j<var;j++)
        {
            if(fabs(a[i][j])>eps)
                temp-=a[i][j]*x[j];
        }
        x[i]=temp/a[i][i];
    }
    return 0;
}
int main()
{
    int n,m,i,j,ans;
    scanf("%d%d",&n,&m);
    equ=m;
    var=n;
    for(i=0;i<m;i++)
    {
        for(j=0;j<n;j++)
        {
            scanf("%lf",&a[i][j]);
        }
        scanf("%lf",&a[i][n]);
    }
    ans=gauss();
    if(ans==-1)
        printf("No solutions\n");
    else if(ans==1)
        printf("Many solutions\n");
    else
    {
        for(i=0;i<n;i++)
            printf("%d\n",(int)(x[i]+0.5));
    }
}