二元一次不定方程@整数解问题

news2024/11/16 10:37:03

文章目录

    • 二元一次不定方程|整数解
    • 定理1整数解存在充要条件
    • 定理2 通解
    • 特解
      • 知识回顾
      • 利用辗转相除法求
      • 使用表达式凑出通解

二元一次不定方程|整数解

二元一次不定方程的一般形式为

a x + b y = c ax + by = c ax+by=c(1)

其中 a a a b b b c c c 是整数,且 a a a b b b 都不为零。

解不定方程 (1) 就是要求出 (1) 式的所有整数解

定理1整数解存在充要条件

方程 (1) 有整数解充分必要条件是 ( a , b ) ∣ c (a,b)|c (a,b)c

证明:

条件: ( a , b ) ∣ c (a,b)|c (a,b)c,

结论:方程(1): a x + b y = c ax + by = c ax+by=c

必要性是明显的。

( a , b ) ∣ a (a,b)|a (a,b)a, ( a , b ) ∣ b (a,b)|b (a,b)b,且 x , y ∈ Z x,y\in\mathbb{Z} x,yZ,则由整除的性质有 ( a , b ) ∣ ( a x + b y ) (a,b)|(ax+by) (a,b)(ax+by),即 ( a ∣ b ) ∣ c (a|b)|c (ab)c

充分性:

由裴蜀恒等式知存在整数 u 0 u_0 u0 v 0 v_0 v0,使

a u 0 + b v 0 = ( a , b ) . au_0 + bv_0 = (a,b). au0+bv0=(a,b).(1-1),将方程(1-1)两边同时乘以 c ( a , b ) \frac{c}{(a,b)} (a,b)c,得 a u 0 c ( a , b ) au_{0}\frac{c}{(a,b)} au0(a,b)c+ b v 0 c ( a , b ) bv_{0}\frac{c}{(a,b)} bv0(a,b)c= c c c

由于 ( a , b ) ∣ c (a,b)|c (a,b)c,于是 c ( a , b ) ∈ Z \frac{c}{(a,b)}\in\mathbb{Z} (a,b)cZ, x 0 = u 0 c ( a , b ) x_0 = u_0 \frac{c}{(a,b)} x0=u0(a,b)c y 0 = v 0 c ( a , b ) y_0 = v_0 \frac{c}{(a,b)} y0=v0(a,b)c 就是 (1) 式的一组整数解

定理2 通解

若知道 (1) 式的一个特解 x 0 ( x 0 , y 0 ) \mathbf{x}_0(x_0, y_0) x0(x0,y0),则立刻就可以写出 (1) 式的全部解

已知 ( x 0 , y 0 ) (x_0, y_0) (x0,y0)是(1)式的一个解,即 a x 0 + b y 0 = c ax_{0}+by_{0}=c ax0+by0=c,则(1)式的全部解为以下式(2)
{ x = x 0 + b ( a , b ) t , y = y 0 − a ( a , b ) t , t = 0 , ± 1 , ± 2 , ⋯ \begin{cases} x = x_0 + \frac{b}{(a,b)}t, \\ y = y_0 - \frac{a}{(a,b)}t, \end{cases} \quad t=0,\pm 1, \pm 2, \cdots {x=x0+(a,b)bt,y=y0(a,b)at,t=0,±1,±2,

证明:

易知(2)式给出的所有 ( x , y ) (x,y) (x,y)都满足方程(1) a x + b y = c ax + by = c ax+by=c

反之,设 ( x , y ) (x,y) (x,y)是(1)式的任意一个解。证明其能够被表示为式(2)的形式

a x + b y = c ax+by=c ax+by=c a x 0 + b y 0 = c , ax_0+by_0=c, ax0+by0=c,(2-0)可得 a ( x − x 0 ) + b ( y − y 0 ) = 0 , a(x-x_0)+b(y-y_0)=0, a(xx0)+b(yy0)=0,(2-1),对其两边乘以 1 ( a , b ) \frac{1}{(a,b)} (a,b)1并移项,

进而有 a ( a , b ) ( x − x 0 ) = − b ( a , b ) ( y − y 0 ) . \frac{a}{(a,b)}(x-x_0)=-\frac{b}{(a,b)}(y-y_0). (a,b)a(xx0)=(a,b)b(yy0).(3)

可见 b ( a , b ) ∣ ( a ( a , b ) ( x − x 0 ) ) \frac{b}{(a,b)}|(\frac{a}{(a,b)}(x-x_{0})) (a,b)b((a,b)a(xx0))

( a ( a , b ) , b ( a , b ) ) = 1 \left(\frac{a}{(a,b)}, \frac{b}{(a,b)}\right)=1 ((a,b)a,(a,b)b)=1,知 b ( a , b ) ∣ ( x 0 − x ) \frac{b}{(a,b)}|(x_0-x) (a,b)b(x0x)

(注意,整除的性质:若 m ∣ a 1 a 2 m|a_{1}a_{2} ma1a2, ( m , a 1 ) = 1 (m,a_{1})=1 (m,a1)=1,则 m ∣ a 2 m|a_{2} ma2;在这里 m = b ( a , b ) m=\frac{b}{(a,b)} m=(a,b)b, a 1 a_{1} a1= a ( a , b ) \frac{a}{(a,b)} (a,b)a, a 2 = ( x − x 0 ) a_{2}=(x-x_{0}) a2=(xx0))

因此存在整数 t t t,使得 x 0 − x x_{0}-x x0x= b ( a , b ) t \frac{b}{(a,b)}{t} (a,b)bt,即 x = x 0 + b ( a , b ) t . x=x_0+\frac{b}{(a,b)}t. x=x0+(a,b)bt.(3-1)

代入(3)式(或带入(2-1)更好算),即得

y = y 0 − a ( a , b ) t y=y_0-\frac{a}{(a,b)}t y=y0(a,b)at(3-2)

特解

知识回顾

对整数 a > 0 a>0 a>0, b > 0 b>0 b>0反复运用带余除法,可得下列等式
a = b q 1 + r 1 , 0 < r 1 < b , b = r 1 q 2 + r 2 , 0 < r 2 < r 1 , r 1 = r 2 q 3 + r 3 , 0 < r 3 < r 2 , r 2 = r 3 q 4 + r 4 , 0 < r 4 < r 3 , . . . r n − 2 = r n − 1 q n + r n , 0 < r n < r n − 1 , r n − 1 = r n q n + 1 + r n + 1 ( r n + 1 = 0 ) . (1) \begin{aligned} \\&a=bq_{1}+r_{1},&&0<r_{1}<b, \\&b=r_{1}q_{2}+r_{2},&&0<r_{2}<r_{1}, \\&r_{1}=r_{2}q_{3}+r_{3},&&0<r_{3}<r_{2}, \\&r_{2}=r_{3}q_{4}+r_{4},&&0<r_{4}<r_{3}, \\&... \\&r_{n-2}=r_{n-1}q_{n}+r_{n},&&0<r_{n}<r_{n-1}, \\&r_{n-1}=r_{n}q_{n+1}+r_{n+1}&&(r_{n+1}=0).\end{aligned}\tag{1} a=bq1+r1,b=r1q2+r2,r1=r2q3+r3,r2=r3q4+r4,...rn2=rn1qn+rn,rn1=rnqn+1+rn+10<r1<b,0<r2<r1,0<r3<r2,0<r4<r3,0<rn<rn1,(rn+1=0).(1)

a a a b b b是任意两个正整数,并进行了式(1)中的辗转相除法,

则有
Q k a − P k b = ( − 1 ) k − 1 r k , 1 ⩽ k ⩽ n . Q_{k}a-P_{k}b=(-1)^{k-1}r_{k},\quad 1\leqslant k\leqslant n. QkaPkb=(1)k1rk,1kn. (2)

其中
{ P 0 = 1 , P 1 = q 1 , P k = q k P k − 1 + P k − 2 , Q 0 = 0 , Q 1 = 1 , Q k = q k Q k − 1 + Q k − 2 , 2 ⩽ k ⩽ n . (3) \left\{\begin{matrix} P_{0}=1,P_{1}=q_{1},P_{k}=q_kP_{k-1}+P_{k-2}, \\Q_{0}=0,Q_{1}=1,Q_{k}=q_kQ_{k-1}+Q_{k-2}, \end{matrix}\right. 2\leqslant k\leqslant n. \tag{3} {P0=1,P1=q1,Pk=qkPk1+Pk2,Q0=0,Q1=1,Qk=qkQk1+Qk2,2kn.(3)

利用辗转相除法求

不妨设(1)式中的 a , b a, b a,b互素,不然去解方程 a ( a , b ) x + b ( a , b ) y = c ( a , b ) \frac{a}{(a,b)}x+\frac{b}{(a,b)}y=\frac{c}{(a,b)} (a,b)ax+(a,b)by=(a,b)c即可。

首先画出如下的框图,然后将由辗转相除法求出的 q 1 , q 2 , ⋯   , q n q_1, q_2, \cdots, q_n q1,q2,,qn依次填在框图的第二行里。

接着在第三行填上 P 0 = 1 , P 1 = q 1 P_0 = 1, P_1 = q_1 P0=1,P1=q1,在第四行填上 Q 0 = 0 , Q 1 = 1 Q_0 = 0, Q_1 = 1 Q0=0,Q1=1

在这里插入图片描述

由递推公式依次求出 P 2 , Q 2 , P 3 , Q 3 , ⋯   , P n , Q n P_2, Q_2, P_3, Q_3, \cdots, P_n, Q_n P2,Q2,P3,Q3,,Pn,Qn

由 辗转相除法基本定理 (取 k = n k=n k=n) 知

Q n a − P n b = ( − 1 ) n − 1 , Q_na - P_nb = (-1)^{n-1}, QnaPnb=(1)n1

a [ ( − 1 ) n − 1 Q n c ] + b [ ( − 1 ) n P n c ] = c . a[(-1)^{n-1}Q_nc] + b[(-1)^nP_nc] = c. a[(1)n1Qnc]+b[(1)nPnc]=c.

因此将以上求得的 P n , Q n P_n, Q_n Pn,Qn 代入下式

{ x 0 = ( − 1 ) n − 1 Q n c , y 0 = ( − 1 ) n P n c , \begin{cases} x_0 = (-1)^{n-1}Q_nc, \\ y_0 = (-1)^nP_nc, \end{cases} {x0=(1)n1Qnc,y0=(1)nPnc, (4)

即得(1)式的一个特解。

例:求不定方程 7 x + 4 y = 189 7x + 4y = 189 7x+4y=189 的一个特解。

( 7 , 4 ) = 1 (7,4)=1 (7,4)=1,可以用辗转相除法求一个特解

在这里插入图片描述

由(4)式可得方程的一个特解为

{ x 0 = − 100 , y 0 = 200. \begin{cases} x_0 = -100, \\ y_0 = 200. \end{cases} {x0=100,y0=200.

由定理2 即得方程的全部解为

{ x = − 100 + 4 t , t = 0 , ± 1 , ± 2 , … y = 200 − 7 t . \begin{cases} x = -100 + 4t, & t = 0, \pm 1, \pm 2, \ldots \\ y = 200 - 7t. \end{cases} {x=100+4t,y=2007t.t=0,±1,±2,

使用表达式凑出通解

也可以用如下方法求这个不定方程的一个特解:

7 x + 4 y = 100 7x + 4y = 100 7x+4y=100,推出 x = 100 − 4 y 7 = 14 + 2 − 4 y 7 x = \frac{100-4y}{7} = 14 + \frac{2-4y}{7} x=71004y=14+724y

2 − 4 y 7 ∈ Z \frac{2-4y}{7}\in\mathbb{Z} 724yZ, y = 2 − 7 k 4 , k ∈ Z y=\frac{2-7k}{4},k\in\mathbb{Z} y=427k,kZ,显然 k k k是偶数才能使 y y y为整数,可取 k = 2 k=2 k=2,此时 y y y可取特值 y 0 = 4 y_0 = 4 y0=4,则 x 0 = 12 x_0 = 12 x0=12

带入通解公式
{ x = x 0 + b ( a , b ) t , y = y 0 − a ( a , b ) t , t = 0 , ± 1 , ± 2 , ⋯ \begin{cases} x = x_0 + \frac{b}{(a,b)}t, \\ y = y_0 - \frac{a}{(a,b)}t, \end{cases} \quad t=0,\pm 1, \pm 2, \cdots {x=x0+(a,b)bt,y=y0(a,b)at,t=0,±1,±2,
不定方程的全部解为
{ x = 12 + 4 t , t = 0 , ± 1 , ± 2 , … y = 4 − 7 t , \begin{cases} x = 12 + 4t, & t = 0, \pm 1, \pm 2, \ldots \\ y = 4 - 7t, \end{cases} {x=12+4t,y=47t,t=0,±1,±2,

虽然用以上两种方法得到的解的表达式形式上不同,但本质上是相同的。

本文来自互联网用户投稿,该文观点仅代表作者本人,不代表本站立场。本站仅提供信息存储空间服务,不拥有所有权,不承担相关法律责任。如若转载,请注明出处:http://www.coloradmin.cn/o/2241444.html

如若内容造成侵权/违法违规/事实不符,请联系多彩编程网进行投诉反馈,一经查实,立即删除!

相关文章

深入理解Flutter生命周期函数之StatefulWidget(一)

目录 前言 1.为什么需要生命周期函数 2.开发过程中常用的生命周期函数 1.initState() 2.didChangeDependencies() 3.build() 4.didUpdateWidget() 5.setState() 6.deactivate() 7.dispose() 3.Flutter生命周期总结 1.调用顺序 2.函数调用时机以及主要作用 4.生…

llama factory lora 微调 qwen2.5 7B Instruct模型

项目背景 甲方提供一台三卡4080显卡 需要进行qwen2.5 7b Instruct模型进行微调。以下为整体设计。 要使用 LLaMA-Factory 对 Qwen2.5 7B Instruct模型 进行 LoRA&#xff08;Low-Rank Adapters&#xff09;微调&#xff0c;流程与之前提到的 Qwen2 7B Instruct 模型类似。LoRA …

机器学习day2-特征工程

四.特征工程 1.概念 一般使用pandas来进行数据清洗和数据处理、使用sklearn来进行特征工程 将任意数据&#xff08;文本或图像等&#xff09;转换为数字特征&#xff0c;对特征进行相关的处理 步骤&#xff1a;1.特征提取&#xff1b;2.无量纲化&#xff08;预处理&#xf…

Llama架构及代码详解

Llama的框架图如图&#xff1a; 源码中含有大量分布式训练相关的代码&#xff0c;读起来比较晦涩难懂&#xff0c;所以我们对llama自顶向下进行了解析及复现&#xff0c;我们对其划分成三层&#xff0c;分别是顶层、中层、和底层&#xff0c;如下&#xff1a; Llama的整体组成…

stm32在linux环境下的开发与调试

环境安装 注&#xff1a;文末提供一键脚本 下载安装stm32cubeclt 下载地址为&#xff1a;https://www.st.com/en/development-tools/stm32cubeclt.html 选择 linux版本下载安装 安装好后默认在家目录st下 > $ ls ~/st/stm32cubeclt_1.16.0 …

第T7周:Tensorflow实现咖啡豆识别

&#x1f368; 本文为&#x1f517;365天深度学习训练营 中的学习记录博客&#x1f356; 原作者&#xff1a;K同学啊 目标 具体实现 &#xff08;一&#xff09;环境 语言环境&#xff1a;Python 3.10 编 译 器: PyCharm 框 架: &#xff08;二&#xff09;具体步骤 1. 使…

亲测有效:Maven3.8.1使用Tomcat8插件启动项目

我本地maven的settings.xml文件中的配置&#xff1a; <mirror><id>aliyunmaven</id><mirrorOf>central</mirrorOf><name>阿里云公共仓库</name><url>https://maven.aliyun.com/repository/public</url> </mirror>…

LLM - 使用 LLaMA-Factory 微调大模型 Qwen2-VL SFT(LoRA) 图像数据集 教程 (2)

欢迎关注我的CSDN&#xff1a;https://spike.blog.csdn.net/ 本文地址&#xff1a;https://spike.blog.csdn.net/article/details/143725947 免责声明&#xff1a;本文来源于个人知识与公开资料&#xff0c;仅用于学术交流&#xff0c;欢迎讨论&#xff0c;不支持转载。 LLaMA-…

神经网络与Transformer详解

一、模型就是一个数学公式 模型可以描述为:给定一组输入数据,经过一系列数学公式计算后,输出n个概率,分别代表该用户对话属于某分类的概率。 图中 a, b 就是模型的参数,a决定斜率,b决定截距。 二、神经网络的公式结构 举例:MNIST包含了70,000张手写数字的图像,其中…

鲸鱼机器人和乐高机器人的比较

鲸鱼机器人和乐高机器人各有其独特的优势和特点&#xff0c;家长在选择时可以根据孩子的年龄、兴趣、经济能力等因素进行综合考虑&#xff0c;选择最适合孩子的教育机器人产品。 优势 鲸鱼机器人 1&#xff09;价格亲民&#xff1a;鲸鱼机器人的产品价格相对乐高更为亲民&…

Flink Source 详解

Flink Source 详解 原文 flip-27 FLIP-27 介绍了新版本Source 接口定义及架构 相比于SourceFunction&#xff0c;新版本的Source更具灵活性&#xff0c;原因是将“splits数据获取”与真“正数据获取”逻辑进行了分离 重要部件 Source 作为工厂类&#xff0c;会创建以下两…

路漫漫其修远兮,吾将上下而求索---第一次使用github的过程记录和个人感受

文章目录 1.仓库位置2.新建仓库3.配置仓库4.克隆和上传5.推荐文章和我的感受 1.仓库位置 这个仓库的位置就是在我们的这个个人主页的右上角&#xff1b;如果是第一次注册账号的话&#xff0c;这个主页里面肯定是不存在仓库的&#xff0c;需要我们自己手动的进行创建&#xff1…

npm list -g --depth=0(用来列出全局安装的所有 npm 软件包而不显示它们的依赖项)

您提供的命令 npm list -g --depth0 是在 Node Package Manager (npm) 的上下文中使用的&#xff0c;用来列出全局安装的所有 npm 软件包而不显示它们的依赖项。 这是它的运作方式&#xff1a; npm list -g --depth0-g: 指定列表应包括全局安装的软件包。--depth0: 限制树形结…

tdengine学习笔记

官方文档&#xff1a;用 Docker 快速体验 TDengine | TDengine 文档 | 涛思数据 整体架构 TDENGINE是分布式&#xff0c;高可靠&#xff0c;支持水平扩展的架构设计 TDengine分布式架构的逻辑结构图如下 一个完整的 TDengine 系统是运行在一到多个物理节点上的&#xff0c;包含…

K8S单节点部署及集群部署

1.Minikube搭建单节点K8S 前置条件&#xff1a;安装docker&#xff0c;注意版本兼容问题 # 配置docker源 wget https://mirrors.aliyun.com/docker-ce/linux/centos/docker-ce.repo -O /etc/yum.repos.d/docker-ce.repo# 安装docker环境依赖 yum install -y yum-utils device-m…

以往运维岗本人面试真题分享

以下是本人面试运维岗的一些面试经历&#xff0c;在此做个记录分享 目录 TCP/IP三次握手 IPtables IPtables四表五链都是什么&#xff1f; nat端口如何做&#xff1f; 开放本机的80端口该如何做&#xff1f; 如何在单用户模式下引导Centos&#xff1f; nginx轮询模式都有…

STM32 串口输出调试信息

软硬件信息 CubeMX version 6.12.1Keil uVision V5.41.0.0 注意 串口有多种&#xff1a; TTL232485 串口的相关知识&#xff1a; 01-【HAL库】STM32实现串口打印&#xff08;printf方式) &#xff0c; 内含 TTL 和 232 区别。 我把 232 串口连进 STM32 串口助手收到的信息…

Python 三种方式实现自动化任务

在这篇文章中&#xff0c;我们将介绍一些用Python实现机器人过程自动化的包。机器人流程自动化&#xff08;Robotic process automation&#xff0c;简称RPA&#xff09;是指将鼠标点击和键盘按压自动化的过程&#xff0c;即模拟人类用户的操作。RPA用于各种应用程序&#xff0…

Android ART知多少?

Android 虚拟机 ART&#xff08;Android Runtime&#xff09;是 Android 平台上的应用程序运行时环境&#xff0c;用于执行应用程序的字节码。ART 自 Android 5.0&#xff08;Lollipop&#xff09;开始取代了 Dalvik&#xff0c;成为 Android 的默认运行时环境。本文将从以下几…

Vulnhub靶场 Billu_b0x 练习

目录 0x00 准备0x01 主机信息收集0x02 站点信息收集0x03 漏洞查找与利用1. 文件包含2. SQL注入3. 文件上传4. 反弹shell5. 提权&#xff08;思路1&#xff1a;ssh&#xff09;6. 提权&#xff08;思路2&#xff1a;内核&#xff09;7. 补充 0x04 总结 0x00 准备 下载链接&#…