【蓝桥杯】简单数论4——丢番图方程
创始人
2024-05-16 20:58:27
0

1、二元线性丢番图方程

方程ax +by = c被称为二元线性丢番图方程,其中a、b、c是已知整数,x、y是变量,问是否有整数解
ax + by= c实际上是二维x-y平面上的一条直线,这条直线上如果有整数坐标点,方程就有解,如果没有整数坐标点,就无解。

 如果存在一个解,就有无穷多个解。

1.1有解的判断条件和通解的形式

定理:设a,b是整数gcd(a, b)=d。如果d不能整除c,那么方程ax + by=c没有整数解,如果d能整除c,那么存在无穷多个整数解。

解释:令a=da',b= db';有ax+by = d(a' x +b'y)=c;如果x、y、a'、b'都是整数,那么c必须是d =gcd(a, b)的倍数,才有整数解

如果(x_0,y_0)是方程的一个特解,所有的解(通解)可的形式x=x_0 +(b/d)n,y= y_0 - (a/d)n,其中n是任意整数。

 

说明: x值按b/d递增,y值按- a/d递增。设(x_0,y_0)是一个格点(格点是指x、y坐标均为整数的点),移动到直线上另一个点(x_0+\Delta x,y_0+\Delta y),有a\Delta x+b\Delta y=0。△x和Ay必须是整数,(x_0+\Delta x,y_0+\Delta y)才是另一个格点。  

\Delta x最小是多少?因为a/d与b/d互素,只有\Delta x = b/d,\Delta y =- a/d时,\Delta x\Delta y才是整数,并满足a\Delta x +b\Delta y = 0。 

定理概况为: ax + by= c有解的充分必要条件d = gcd(a, b)能整除c

例:
(1)方程18x + 3y = 7没有整数解,因为gcd(18,3) = 3,3不能整除7;

(2)方程25x + 15y = 70存在无穷个解,因为gcd(25,15)= 5且5整除70,一个特解是x_0=4,y_0 = -2,通解是x=4 + 3n,y = -2- 5n

1.2例题一:线段上的格点数量

【题目描述】在二维平面上,给定两个格点p_1=(x_1,y_1)p_2=(x_2,y_2),问线段p_1p_2上除了p_1,p_2外还有几个格点?设x_1< x_2

【思路】
首先利用p_1,p_2把线段表示为方程ax + by = c的形式,它肯定有整数解。
然后在线段范围内,根据x的通解的表达式x = x_0+ (b/d)n,当x_1<x<x_2时,求出n的取值情况有多少个,这就是线段内的格点数量。

计算步骤:

(1)、用p_1(x_1,y_1)p_2(x_2,y_2)表示线段,线段表示为:

(y_2-y_1)x + (x_1-x_2)y = y_2x_1-y_1x_2

(2)、对照ax + by = c,得:
a = y_2-y_1, b = x1_-x_2,c = y_2x_1-y_1x_2

d = gcd(a,b) = gcd(\left | y_2-y_1 \right |,\left |x1-x2 \right |)

(3)、对照通解公式x = x_0+ (b/d)nn,令特解是x,代入限制条件x_1<x<x_2,有:
x_1< x+((x_1-x_2)/d)n < x2

当-d < n< 0时满足上面的表达式,此时n有d-1种取值,即线段内有d-1个格点。

2、方程的特解与扩展欧几里得算法

求解方程ax + by = c的关键是找到一个特解
根据定理的描述,解和求GCD有关;
求特解用到了欧几里得求GCD的思路,称为扩展欧几里得算法

2.1扩展欧几里得算法

方程ax + by = gcd(a, b),根据定理,它有整数解
定理:设a, b是整数且gcd(a, b)=d。如果d不能整除c,那么方程ax + by=c没有整数解,如果d能整除c,那么存在无穷多个整数解。
扩展欧几里得算法求一个特解(x_0,y_0)的代码:

def exgcd(a,b):if b == 0:return 1, 0x,y = exgcd(b,a % b)return y, x - a // b * y    # 返回特解xo,yo
a,b = map (int,input ().split())#   试试6x+15y=3
x,y = exgcd (a,b)#计算得到特解
print(x, y)

2.2扩展欧几里得算法与方程ax+by=c的特解

用扩展欧几里得算法得到ax +by =ged(a,b)的一个特解后,再利用它求方程ax +by= c的一个特解。步骤如下:
(1)判断方程ax +by = c是否有整数解,即gcd(a,b)能整除c。记d= gcd(a,b)。
(2)用扩展欧几里得算法求ax + by = d的一个特解x_0,y_0
(3)在ax_0 + by_0= d两边同时乘以c/d,得: ax_0c/d + by_0c/d=c(目的是构造c,这样和ax + by= d就能消掉c)

(4)对照ax +by =c,得到它的一个解(x_0',y_0')是:x_0'= x_0c/d,y_0'= y_0c/d

(5)方程ax + by = c的通解x=x_0'+ (b/d)n,y =y_0' - (a/d)n

 

相关内容

热门资讯

喜欢穿一身黑的男生性格(喜欢穿... 今天百科达人给各位分享喜欢穿一身黑的男生性格的知识,其中也会对喜欢穿一身黑衣服的男人人好相处吗进行解...
发春是什么意思(思春和发春是什... 本篇文章极速百科给大家谈谈发春是什么意思,以及思春和发春是什么意思对应的知识点,希望对各位有所帮助,...
网络用语zl是什么意思(zl是... 今天给各位分享网络用语zl是什么意思的知识,其中也会对zl是啥意思是什么网络用语进行解释,如果能碰巧...
为什么酷狗音乐自己唱的歌不能下... 本篇文章极速百科小编给大家谈谈为什么酷狗音乐自己唱的歌不能下载到本地?,以及为什么酷狗下载的歌曲不是...
华为下载未安装的文件去哪找(华... 今天百科达人给各位分享华为下载未安装的文件去哪找的知识,其中也会对华为下载未安装的文件去哪找到进行解...
怎么往应用助手里添加应用(应用... 今天百科达人给各位分享怎么往应用助手里添加应用的知识,其中也会对应用助手怎么添加微信进行解释,如果能...
家里可以做假山养金鱼吗(假山能... 今天百科达人给各位分享家里可以做假山养金鱼吗的知识,其中也会对假山能放鱼缸里吗进行解释,如果能碰巧解...
四分五裂是什么生肖什么动物(四... 本篇文章极速百科小编给大家谈谈四分五裂是什么生肖什么动物,以及四分五裂打一生肖是什么对应的知识点,希...
一帆风顺二龙腾飞三阳开泰祝福语... 本篇文章极速百科给大家谈谈一帆风顺二龙腾飞三阳开泰祝福语,以及一帆风顺二龙腾飞三阳开泰祝福语结婚对应...
美团联名卡审核成功待激活(美团... 今天百科达人给各位分享美团联名卡审核成功待激活的知识,其中也会对美团联名卡审核未通过进行解释,如果能...