高斯消元（浮点数）

#include<iostream>#include<cmath>#define eps (1e-9)using namespace std;const int MAXN=10005;double a[MAXN][MAXN],x[MAXN];int equ,var;int Gauss(){    int i,j,k,col,max_r;    for(k=0,col=0;k<equ&&col<var;k++,col++)    {        max_r=k;        for(i=k+1;i<equ;i++)            if(fabs(a[max_r][col])>fabs(a[max_r][col]))                max_r=i;        if(fabs(a[max_r][col])<eps) return 0;        if(k!=max_r)        {            for(j=col;j<var;j++)                swap(a[k][j],a[max_r][j]);            swap(x[k],x[max_r]);        }        x[k]/=a[k][col];        for(j=col+1;j<var;j++) a[k][j]/=a[k][col];        a[k][col]=1;        for(i=0;i<equ;i++)            if(i!=k)        {            x[i]-=x[k]*a[i][k];            for(j=col+1;j<var;j++) a[i][j]-=a[k][j]*a[i][col];            a[i][col]=0;        }    }    return 1;}