【数论】线性模方程组

由于被AK大神误导 - - 导致之前写了一个完全错误的总结 特此修改之

 

求方程组 x=a[i](mod r[i]) 的解x

即为求解方程 x=a[i]+r[i]*y

联想普通方程组的解法 一直将两个方程合并成一个 最后只剩一个方程 求解即可

设待合并的两个方程为 x=a1+r1*y、x=a2+r2*z

合并之

a1+r1*y=a2+r2*z

->r1*y+r2*z=a2-a1

拓展欧几里得求出(y,z)

则x最小为a1+r1*y

这时两方程合并为 x=a1+r1*y(mod lcm(r1,r2))

就这样两两合并 如果某次(y,z)无解 则方程组无解

 

代码:

 1 scanf("%d%d",&r1,&a1);
 2 for (i=2;i<=n;i++){
 3     scanf("%d%d",&r2,&a2);
 4     ll gc=extgcd(x,y,r1,r2);
 5     if ((a1-a2)%gc) bo=0;
 6     ll mo=r2/gc;
 7     x=(x*(a2-a1)/gc%mo+mo)%mo;
 8     a1+=r*x;
 9     r1*=r2/gc;
10 }
View Code

 

转载于:https://www.cnblogs.com/g-word/p/3372752.html

评论
添加红包

请填写红包祝福语或标题

红包个数最小为10个

红包金额最低5元

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

抵扣说明:

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

余额充值