线性同余方程组模板

void exgcd(int a,int b,int &d, int &x,int &y){///扩展欧几里德算法
    if(!b){x=1;y=0;d=a;}
    else{exgcd(b,a%b,d,y,x);  y-=x*(a/b);}
}
int n,m;///m代表方程组的数量
int a[11],b[11];///x%a[i]=b[i]
///求解同余方程组
bool Chi_r(int &minnum,int &mod){///minnum表示最小正整数,mod表示模
    int m1=a[0],b1=b[0],x,y,d;
    for(int i=1;i<m;++i){
        int m2=a[i],b2=b[i];
        exgcd(m1,m2,d,x,y);
        if((b2-b1)%d!=0) return false;///不存在解
        int t=m2/d;
        x=(x*((b2-b1)/d)%t+t)%t;
        b1=m1*x+b1;
        m1=m1*(m2/d);
    }
    minnum=b1;
    mod=m1;
    return true;///存在解
}

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

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

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

请填写红包祝福语或标题

红包个数最小为10个

红包金额最低5元

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

抵扣说明:

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

余额充值