裴蜀定理及其证明
一、裴蜀定理
对于 x , y x,y x,y 的二元一次不定方程 a x + b y = c ax+by=c ax+by=c ,其有解的充要条件为 gcd ( a , b ) ∣ c \gcd(a,b)\mid c gcd(a,b)∣c 。
1.充分性证明
充分性 若 gcd ( a , b ) ∣ c \gcd(a,b)\mid c gcd(a,b)∣c ,则 a x + b y = c ax+by=c ax+by=c 有解。
证明
设 k k k 为 a , b a,b a,b 线性组合的最小非负解。令 q = ⌊ a k ⌋ q=\left\lfloor\dfrac{a}{k}\right\rfloor q=⌊ka⌋ ,则有
r = a − q k = a − q ( a x + b y ) = a ( 1 − q x ) + b ( − q y ) \begin{aligned} r&= a-qk = a-q(ax+by) \\[6pt]&=a(1-qx)+b(-qy) \end{aligned} r=a−qk=a−q(ax+by)=a(1−qx)+b(−qy)
显然 r r r 也为 a , b a,b a,b 线性组合的解,且 0 ≤ r ≤ k 0\le r \le k 0≤r≤k 。由于 k k k 为最小非负解,故 r = 0 r = 0 r=0 ,那么 k ∣ a k \mid a k∣a 。同理可得 k ∣ b k\mid b k∣b 。
令 d = gcd ( a , b ) d=\gcd(a,b) d=gcd(a,b) ,则 k ∣ d k\mid d k∣d 且 d ≥ k d \ge k d≥k 。
因为 d ∣ a , d ∣ b d\mid a,d\mid b d∣a,d∣b 且 s s s 为 a , b a,b a,b 线性组合的解,所以 d ∣ k d \mid k d∣k 。
因为 s > 0 s > 0 s>0 ,所以 d = k d = k d=k ,则 a x + b y = c ax+by=c ax+by=c 的最小非负解为 gcd ( a , b ) \gcd(a,b) gcd(a,b) 。
显然 ∀ c = k gcd ( a , b ) , k ∈ Z + \forall c=k\gcd(a,b),k\in \mathbb{Z}^+ ∀c=kgcd(a,b),k∈Z+ 是原方程的解。
2.必要性证明
必要性 若 a x + b y = c ax+by=c ax+by=c 有解,则 gcd ( a , b ) ∣ c \gcd(a,b)\mid c gcd(a,b)∣c 。
证明 令 d = gcd ( a , b ) d=\gcd(a,b) d=gcd(a,b) ,则 d ∣ a , d ∣ b d\mid a,d\mid b d∣a,d∣b 。
因为 a x + b y = c ax+by=c ax+by=c 有解,所以 d ∣ a x , d ∣ b y d\mid ax,d\mid by d∣ax,d∣by ,故 d ∣ a x + b y = c d\mid ax+by=c d∣ax+by=c 。
3.推广
对于不定方程 x 1 y 1 + x 2 y 2 + . . . + x n y n = k ( y i ∈ Z ) x_1y_1+x_2y_2+...+x_ny_n=k~(y_i \in \mathbb{Z}) x1y1+x2y2+...+xnyn=k (yi∈Z) ,其有解的充要条件为 gcd { x i } ∣ k \gcd\{x_i\}\mid k gcd{xi}∣k 。
二、裴蜀定理模板题
题目链接:P4549 【模板】裴蜀定理
要注意负数要转成正数再算 gcd \gcd gcd 。
代码:
#include <bits/stdc++.h>
using namespace std;
#define int long long
#define R register
int n,ans,a[25];
int gcd(R int a,R int b) {return b==0?a:gcd(b,a%b);}
signed main()
{
scanf("%lld",&n);
for(R int i=1; i<=n; i++)
{
scanf("%lld",&a[i]);
a[i]=a[i]>0?a[i]:-a[i];
ans=gcd((i==1?a[1]:ans),a[i]);
}
printf("%lld\n",ans);
return 0;
}