高斯消元模板(取余方程组)

const int maxn = 110;
const int mod  = 2;
int a[maxn][maxn],x[maxn],free_x[maxn],free_num;
int Qpow( int a , int b )
{
    int res = 1;
    while ( b )
    {
        if ( b&1 )
            res = res*a%mod;
        a = a*a%mod;
        b = b>>1;
    }
    return res;
}
int Gauss( int equ , int var )
{
    int max_r,col,k;
    free_num = 0;
    for ( k=0,col=0 ; k<equ&&col<var ; k++,col++ )
    {
        max_r = k;
        for ( int 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--;
            free_x[free_num] = col;
            continue;
        }
        if ( max_r!=k )
            for ( int j=col ; j<var+1 ; j++ )
                swap ( a[k][j] , a[max_r][j] );
        for ( int 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 ( int j=col ; j<var+1 ; j++ )
                    a[i][j] = ( ( a[i][j]*ta-a[k][j]*tb )%mod+mod )%mod;
            }
    }
    for ( int i=k ; i<equ ; i++ )
        if ( a[i][col]!=0 ) return -1;
    if ( k<var ) return var-k;
    for ( int i=var-1 ; i>=0 ; i-- )
    {
        x[i] = a[i][var];
        for ( int j=i+1 ; j<var ; j++ )
            if ( a[i][j] )
                x[i] = ( ( x[i]-a[i][j]*x[j] )%mod+mod )%mod;
        x[i] = x[i]*Qpow( a[i][i] , mod-2 )%mod;
    }
    return 0;
}

 

  • 0
    点赞
  • 0
    收藏
    觉得还不错? 一键收藏
  • 0
    评论
评论
添加红包

请填写红包祝福语或标题

红包个数最小为10个

红包金额最低5元

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

抵扣说明:

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

余额充值