高斯消元法(高斯·约当消元法)(整数)

整数高斯消元作用在于最后得出的解可以判断是否为整数。

首先看高斯消元法(浮点)
整数的高斯消元法,在消元时通过等式乘法,将同一项系数统一成它的最小公倍数,然后再进行消元。

核心代码:

int n,A[MAXN][MAXN*2];
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;}
void Gauss()
{
    int r,c,mxr,n=N,m=N*2+1,l,ta,tb;
    for(r=1,c=1;r<=n&&c<m;r++,c++)
    {
        mxr=r;
        for(int i=r+1;i<=n;i++)
            if(abs(A[i][c])>abs(A[mxr][c]))
                mxr=i;
        if(A[mxr][c]==0)
        {r--;continue;}
        if(mxr!=r)swap(A[mxr],A[r]);
        for(int i=1;i<=n;i++)
            if(i!=r&&A[i][c])
            {
                l=lcm(abs(A[i][c]),abs(A[r][c]));
                ta=l/A[i][c];
                tb=l/A[r][c];
                for(int j=c;j<=m;j++)
                    A[i][j]=A[i][j]*ta-A[r][j]*tb;
                if(i<c)
                    A[i][i]*=ta;
            }
    }
}

因为没找到整数高斯消元裸题,所以没有完整代码。。。

  • 1
    点赞
  • 2
    收藏
    觉得还不错? 一键收藏
  • 0
    评论

“相关推荐”对你有帮助么?

  • 非常没帮助
  • 没帮助
  • 一般
  • 有帮助
  • 非常有帮助
提交
评论
添加红包

请填写红包祝福语或标题

红包个数最小为10个

红包金额最低5元

当前余额3.43前往充值 >
需支付:10.00
成就一亿技术人!
领取后你会自动成为博主和红包主的粉丝 规则
hope_wisdom
发出的红包
实付
使用余额支付
点击重新获取
扫码支付
钱包余额 0

抵扣说明:

1.余额是钱包充值的虚拟货币,按照1:1的比例进行支付金额的抵扣。
2.余额无法直接购买下载,可以购买VIP、付费专栏及课程。

余额充值