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)n$ n，令特解是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$

【蓝桥杯】简单数论4——丢番图方程相关推荐

蓝桥杯---简单的计算器
蓝桥杯---简单的计算器题目: 问题描述编程模拟计算器的加.减.乘.除功能,根据用户输入的运算符,对两个数进行运算.(要求switch语句) 输入格式输入只有一行,用空格隔开的运算符和两个运算数 ...
蓝桥杯简单题之十六进制转八进制（JAVA版）
首先对等待程序猿成长之路系列的小伙伴们说声抱歉,春节来临,可能要拖更一段时间,但是草稿已经在写,发布将会于春节后进行发布春节第一天闲着无聊打算找个蓝桥杯的题目练练手,诶嘿,这就找到一题,话不多说,上 ...
【蓝桥杯】计算指定日期为当年第几天
文章目录前言题目分析算法难度实战 1.创建算法 2.创建测试用例 3.测试结果总结前言蓝桥杯全国软件和信息技术专业人才大赛由工业和信息化部人才交流中心主办,每年参赛人数超过30000 ...
2019第十届蓝桥杯省赛C/C++B组题解
作为很久很久没有写过算法题的考研狗 OvO(是这样的,卑微考研),我于今年再次参加了蓝桥杯... 大家都说这次蓝桥杯简单,我也是这么觉得的(或者说是我运气好),或者说是考细节吧(虽然我也没怎么考虑细节 ...
蓝桥杯历年真题，来，看看
.标题: 世纪末的星期曾有邪教称1999年12月31日是世界末日.当然该谣言已经不攻自破. 还有人称今后的某个世纪末的12月31日,如果是星期一则会.... 有趣的是,任何一个世纪末的年份的12月3 ...
【蓝桥杯】简单数论1——GCDLCM
GCD 最大公约数Greatest Common Divisor(GCD):整数a和b的GCD是指能同时整除a和b的最大整数,记为gcd(a,b).由于-a的因子和a的因子相同,因此gcd(a, b) ...
AcWing蓝桥杯AB组辅导课08、数论
文章目录前言一.数论例题例题1:AcWing 1246. 等差数列(最大公约数,第十届蓝桥杯省赛C++B第7题) 分析题解:最大公约数例题2:AcWing 1295. X的因子链(算数基本 ...
[蓝桥杯][2013年第四届真题]买不到的数目-模拟，数论
题目描述小明开了一家糖果店.他别出心裁:把水果糖包成4颗一包和7颗一包的两种.糖果不能拆包卖. 小朋友来买糖的时候,他就用这两种包装来组合.当然有些糖果数目是无法组合出来的,比如要买 10 颗糖. ...
蓝桥杯ACM简单题－生日蜡烛 2020-12-22
蓝桥杯ACM简单题-生日蜡烛 1.问题描述某人从某年开始每年都举办一次生日派对,并且每次都要吹熄与年龄相同根数的蜡烛.现在算起来,他一共吹熄了236根蜡烛.请问,他从多少岁开始过生日派对的? 2.算 ...

【蓝桥杯】简单数论4——丢番图方程