解方程组（高斯消元模板）

题目描述：
Tom 是个品学兼优的好学生，但由于智商问题，算术学得不是很好，尤其是在解方程这个方面。虽然他解决 2x=2 这样的方程游刃有余，但是对于下面这样的方程组就束手无策了。
x+y=3
x-y=1
于是他要你来帮忙。给定一个线性多元一次方程组，请你求出所有未知数的解。
保证在 int 范围内可以处理所有问题。
输入格式：
输入文件的第一行一个数字 N（1≤N≤100），表示给定的方程组中的未知数的个数，同时也是这个方程组含有的方程个数。
第 2 到 N+1 行，每行 N+1 个数。每行的前 N 个数表示第 1 到 N 个未知数的系数。第 N+1 个数表示 N 个未知数乘以各自系数后再相加的和。
输出格式：
输出一行，有 N 个整数，表示第 1 到 N 个未知数的值(整数解)，而且数据保证有整数解。
样例输入：
2
1 1 3
1 -1 1
样例输出：
2 1
题目分析：
高斯消元模板题。
附代码：

#include
#include
#include
#include
#include
#include
#include
#include
#include
#include
#include
#include
#include
using namespace std;

const int maxn=110;
int n;
double f[maxn][maxn],ans[maxn];

void gauss()
{
    for(int i=1;i<=n;i++)
    {
        int l=i;
        for(int j=l+1;j<=n;j++)
            if(fabs(f[l][i])<fabs(f[j][i]))
                l=j;
        if(l!=i)
            for(int j=i;j<=n+1;j++)
                swap(f[l][j],f[i][j]);
        for(int j=i+1;j<=n;j++)
        {
            double tmp=f[j][i]/f[i][i];
            for(int k=i;k<=n+1;k++)
                f[j][k]=f[j][k]-f[i][k]*tmp;
        }
    }
    for(int i=n;i>=1;i--)
    {
        double tmp=f[i][n+1];
        for(int j=n;j>i;j--)
            tmp-=ans[j]*f[i][j];
        ans[i]=tmp/f[i][i];
    }
}

int main()
{
    //freopen("lx.in","r",stdin);

    scanf("%d",&n);
    for(int i=1;i<=n;i++)
        for(int j=1;j<=n+1;j++)
            scanf("%lf",&f[i][j]);
    gauss();
    for(int i=1;i<=n;i++)
        printf("%d ",int(ans[i]+0.5));

    return 0;
}