前置条件:
ax + by = c , gcd(a, b) = d
计算:
x + y = ;
设 = a1 , = b1 ;
原式变为 a1x + b1y = ;
易知 gcd(a1, b1) = 1 ;
设 a1m +b1n = 1 ;
由扩展欧几里得算法可得 m0 为一个特解 ;
m = m0 + b1t , t 为任意整数 ;
则 a1x + b1y = 的通解为 :
x = m * ;
一个特解为 x0 = m0 * ;
得到 x = x0 + ;(这里把b1换成了)
所以 x = x0 + , t为任意整数 (因为 t 不一定是 d 的倍数, 所以不能约分)
在学CRT的时候自己琢磨的,如有错误,敬请斧正。
✳
谢谢你读到这里。
更新于 2023-01-01