形如 \(ax\equiv b(\mod n)\) 的方程称为线性同余方程,从区间 \([0,n-1]\) 中求解 \(x\).

逆元求解。

假设 \(gcd(a,n)=1\),两边同时乘上 \(a^{-1}\) 即可。

\(g=gcd(a,n)\),左侧始终可以 被 \(g\) 整除,若右侧不可则无解。

若右侧可以被 \(g\) 整除,则将 \(a,b,n\) 同时除以 \(g\),得到 \(a’x\equiv b'(\mod n’)\).

此时 \(gcd(a’,n’)=1\),回到上面的情况。

所以解为 \(x\equiv x’+i*n'(\mod n),i\in[0,g-1]\),个数为 \(g\) 个或 \(0\) 个。

扩展欧几里得算法求解。

方程可以写成 \(ax+nk= b\),有解的充要条件是 \(gcd(a,n)|b\).

先算出 \(ax_0+nk_0=gcd(a,n)\),之后两边同时除以 \(gcd(a,n)\) 再乘上 \(b\) 即为一组解。

\(a=\frac{x_0}{gcd(a,n)}b,k=\frac{k_0}{gcd(a,n)}b\).

\(gcd(a,n)=1,ax+nk=b\) 的一组特解为 \(x_0,k_0\),则任意解可以写成 \(x=x_0+nt,k=k_0-at,t\) 为任意整数。

对于求一个最小整数解,\(x=(x\mod t+t)\mod t,t=\frac{n}{gcd(a,n)}\)

原文地址:http://www.cnblogs.com/safeng/p/16907874.html

1. 本站所有资源来源于用户上传和网络,如有侵权请邮件联系站长! 2. 分享目的仅供大家学习和交流,请务用于商业用途! 3. 如果你也有好源码或者教程,可以到用户中心发布,分享有积分奖励和额外收入! 4. 本站提供的源码、模板、插件等等其他资源,都不包含技术服务请大家谅解! 5. 如有链接无法下载、失效或广告,请联系管理员处理! 6. 本站资源售价只是赞助,收取费用仅维持本站的日常运营所需! 7. 如遇到加密压缩包,默认解压密码为"gltf",如遇到无法解压的请联系管理员! 8. 因为资源和程序源码均为可复制品,所以不支持任何理由的退款兑现,请斟酌后支付下载 声明:如果标题没有注明"已测试"或者"测试可用"等字样的资源源码均未经过站长测试.特别注意没有标注的源码不保证任何可用性