高斯消元模板

typedef int Matrix[maxn][maxn];
int gcd(int a,int b)
{
    return b==0?a:gcd(b,a%b);
}
int lcm(int a,int b)
{
    return  a/gcd(a,b)*b;
}
int gauss(Matrix A,int n,int m)
{
    int row,col;
    for(row=0,col=0; row<n&&col<m; row++,col++)
    {
        int r=row;
        for(int i=row+1; i<n; i++)
            if(abs(A[i][col])>abs(A[r][col])) r=i;
        if(r!=row) for(int j=0; j<=m; j++) swap(A[r][j],A[row][j]);
        if(A[row][col]==0)
        {
            row--;
            continue;
        }
        for(int k=row+1; k<n; k++)//消元过程
            if(A[k][col])
            {
                int LCM=lcm(abs(A[k][col]),abs(A[row][col]));
                int ra=LCM/abs(A[k][col]);
                int rb=LCM/abs(A[row][col]);
                if (A[k][col]*A[row][col]<0) rb=-rb;
                for(int j=col; j<=m; j++)
                {
                    A[k][j]=A[k][j]*ra-A[row][j]*rb;
                }
            }
    }
    for(int i=row; i<n; i++)//无解情况
        if(A[i][m]) return -1;
    if(row<m) return 0;//多解
    for(int i=m-1; i>=0; i—)//回代过程
    {
        for(int j=i+1; j<m; j++)
        {
            A[i][m]-=A[j][m]*A[i][j];
        }
        A[i][m]/=A[i][i];
    }
    return 1;
}