高斯消元-整型

#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;
int a[1005][1005];//增广矩阵
int x[505];//解集
bool free_x[505];//判断是否是不确定的变元
int free_num;
int equ,var;//equ个方程var个未知数，
//增广矩阵行数为equ,分别为0-equ-1;
//列数为var+1,分别为0-var
int gcd(int a,int b)
{
    if(b==0)
        return a;
    return gcd(b,a%b);
}
int lcm(int a,int b)
{
    return a/gcd(a,b)*b;
}
int gauss()
{
    int max_r,col,k,i,j;
    int temp,free_x_num,free_index;
    for(k=0,col=0;k<equ&&col<var;k++,col++)
    {
        max_r=k;
        for(i=k+1;i<equ;i++)
        {
            if(abs(a[i][col])>abs(a[max_r][col]))
                max_r=i;
        }
        if(a[max_r][col]==0)
        {
            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(a[i][col]!=0)
            {
                int Lcm=lcm(abs(a[i][col]),abs(a[k][col]));
                int ta=Lcm/abs(a[i][col]);
                int tb=Lcm/abs(a[k][col]);
                if(a[i][col]*a[k][col]<0)
                    tb=-tb;
                for(j=col;j<=var;j++)
                {
                    a[i][j]=a[i][j]*ta-a[k][j]*tb;
                }
            }
        }
    }
    for(i=k;i<equ;i++)
        if(a[i][col]!=0)
        return -1;
    if(k<var)
    {
       //首先，自由变元有var-k个，即不确定的变元至少有var-k个
       for(i=k-1;i>=0;i--)
       {
           free_x_num=0;
           //用于判断该行中的不确定变元的个数，如果超过1个，则无法求解，它们仍然为不确定的变元
           for(j=0;j<var;j++)
           {
               if(a[i][j]!=0&&free_x[j])
               {
                   free_x_num++;
                   free_index=j;
               }
           }
           if(free_x_num>1)
            continue;
           temp=a[i][var];
           for(j=0;j<var;j++)
           {
               if(a[i][j]!=0&&j!=free_index)
               {
                   temp-=a[i][j]*x[j];
               }
           }
           x[free_index]=temp/a[i][free_index];//求出改变元
           free_x[free_index]=0;//改变元是确定的
       }
       return var-k;
    }
    for(i=var-1;i>=0;i--)
    {
        int temp=a[i][var];
        for(j=i+1;j<var;j++)
        {
            temp-=a[i][j]*x[j];
        }
        x[i]=temp/a[i][i];
        if(temp%a[i][i]!=0)
            return -2;//有浮点数解
    }
    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("%d",&a[i][j]);
        }
        scanf("%d",&a[i][n]);
    }
    ans=gauss();
    if(ans==-1)
        printf("无解\n");
    else if(ans==-2)
        printf("有浮点数解\n");
    else if(ans>0)
    {
        printf("无穷多解，自由变元的个数 %d\n",ans);
        for(i=0;i<var;i++)
        {
            if(free_x[i])
                printf("不确定的 %d\n",i+1);
            else
                printf("%d %d\n",i+1,x[i]);
        }
    }
    else
    {
        for(i=0;i<var;i++)
            printf("%d\n",x[i]);
    }
}