首页 > 其他分享 >线性同余方程

线性同余方程

时间:2022-11-20 08:55:05浏览次数:41  
标签:方程 nk gcd equiv 线性 ax 同余 mod

形如 \(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)}\)

标签:方程,nk,gcd,equiv,线性,ax,同余,mod
From: https://www.cnblogs.com/safeng/p/16907874.html

相关文章

  • 用SPSS估计HLM多层(层次)线性模型模型|附代码数据
    原文链接:http://tecdat.cn/?p=3230作为第一步,从一个不包含协变量的空模型开始 ( 点击文末“阅读原文”获取完整代码数据******** )。每所学校的截距,β0J,然后设置为平......
  • C语言:一元二次方程求解
    题目已知一元二次方程......
  • 【线性表】之顺序表(C语言)
    【线性表】之顺序表​​线性表​​​​顺序表​​​​结构定义​​​​初始化​​​​销毁​​​​打印​​​​扩展空间​​​​尾插​​​​头插​​​​尾删​​​​头删......
  • 【线性表】之栈(C语言)
    栈​​回顾​​​​栈​​​​结构定义​​​​初始化​​​​销毁​​​​入栈​​​​出栈​​​​返回栈顶元素​​​​返回栈中元素个数​​​​判断栈是否为空​​​​......
  • 【线性表】之队列(C语言)
    队列​​队列的概念​​​​结构定义​​​​初始化​​​​销毁​​​​队尾入​​​​队头出​​​​队头出​​​​队头数据​​​​队尾数据​​​​是否为空​​​​返......
  • 《大话数据结构》线性表代码总结
    //线性表存储的结构代码#include<stdio.h>#include<stdlib.h>#include<time.h>#defineMAXSIZE1000//静态链表部分的#defineMAX_SIZE20//最大长度#defineOK1#defineER......
  • 时间序列数据挖掘之分段线性表示(PLR)
    前言本篇博客用于记录个人在时间序列数据挖掘中进行的timeseriesrepresentation的实践。主要采用PLR(piecewiselinearrepresentation)的方式进行时间序列的降......
  • bzoj #4161. Shlw loves matrixI 常系数线性递推
    常系数线性递推题目链接做法:容易得到一个\(k^3logn\)的做法,考虑更快。需要利用矩阵特征方程f(x)(恰好前些天刚学对\(A^n\)用\(f(A)\)可以对其取模。之后就可以快速求......
  • Python解答多元线性方程组
    前几天一个在乌克兰读什么生态环境工程的小朋友在朋友圈发了一张求解多元线性方程组的图片,这老师真变态,这么大的数非得让他们动手拿笔来求解,我就给他一下标准答案,以至他......
  • 使用matplotlib模拟线性回归
    首先需要两个模块:1.numpy2.matplotlib.pylab安装命令pipinstallnumpypipinstallmatplotlib线性回归的主要作用就是用一条线性的函数或者表达式来拟合在离散空间的随机......