數學基礎
數學:
一:整除
概念:
如果 \(A\) 能整除 \(B\) 則記之為 \(A|B\) 即存在一個數 \(k \in Z\) 使 \(Ak = B\)。即 \(B\) 是 \(A\) 的倍數。
性質:
1. 如果 \(a|b,b|c\) 那么 \(a|c\) 表明整除具有傳遞性。
2. 如果 \(a|b,a|c\) 那么對于任意整數對 \((x,y)\) 滿足 \(x \in Z,y \in Z\) 都有 \(a|bx+cy\)。
下面證明這兩個性質。
證明性質1:
因為 \(a|b\),所以 \(b = xa,x \in Z\);
因為 \(b|c\),所以 \(c = yb,y \in Z\)。
綜上所以 \(c = yb\),即 \(c = y\times (xa) = xya\)。
因為 \(x \in Z,y \in Z\),所以 \(xy \in Z\),所以 \(a|c\)。
證明性質2:
因為 \(a|b\) 所以 \(b = ka,k \in Z\)。
因為 \(a|c\) 所以 \(c = Ka,K \in Z\)。
因為對于任意整數對 \((x,y)\) 滿足 \(x \in Z,y \in Z\) 所以 \(bx \in Z,cy \in Z\)。
綜上;故:原式 \(bx + cy = kax + Kay\) 合并同類項得 \(bx + cy = a(kx + Ky)\)。
因為 \(k \in Z,x \in Z,K \in Z,y \in Z\),所以 \(kx \in Z,Ky \in Z\)。
所以 \(kx + Ky \in Z\) 。
至此原式得:\(a|a\times (kx + Ky)\)
所以,對于任意整數對 \((x,y)\) 滿足 \(x \in Z,y \in Z\) 都有 \(a|bx+cy\)。
證明題:
一:設 \(a|n,b|n\) 且存在整數 \(ax + by = 1\),證明 \(ab|n\)。
證明:
因為 \(a|n,b|n\) 所以 \(n = ka,n = Kb,k \in Z,K \in Z\)。
因為 \(ax + by = 1\),
所以 \(n = n\times (ax + by)\);
即 \(n = nax + nby\),
即 \(n = (Kb)ax + (ka)by\),
即 \(n = Kabx + kaby\),
即 \(n = ab(Kx + ky)\);
因為 \(K \in Z,x \in Z,k \in Z,y \in Z\),
所以 \(Kx \in Z,ky \in Z\),
所以 \(Kx + ky \in Z\),
所以 \(ab|n\)。
二:exgcd算法
形如求解 \(ax+by=c\) 知道 \(a,b,c\) 求解整數解 \(x,y\) 的值。
通常使用 exgcd 算法求解。
首先考慮無整數解的情況,根據貝祖定理,對于 \(a\in Z,b\in Z,c\in Z\) 如果有兩個整數 \(x,y\) 滿足 \(ax+by=c\) 那么 \(c=k\times\gcd(a,b),k\in Z\)。
即:如果 \(c\mod\gcd(a,b)!=0\) 那么就無解。
如果有解,那么:
令 \(ax+by=\gcd(a,b)\),\(bx+(a\mod b)y=\gcd(b,a\mod b)\)。
因為 \(\gcd(a,b)=\gcd(b,a\mod b)\)。
所以 \(ax+by=bx+(a\mod b)y\)。
因為 \(a\mod b=a-\lfloor a/b\rfloor\times b\)。
所以原式為 \(ax+by=bx+(a-\lfloor a/b\rfloor\times b)y\)。
化簡下得:\(ax+by=bx+ay-\lfloor a/b\rfloor by\)。
也就是 \(ax+by=b(x-\lfloor a/b\rfloor y)+ay\)。
所以 \(x=y,y=x-\lfloor a/b\rfloor y\)。
但是這么做只能求出 \(ax+by=\gcd(a,b)\) 的特解,并不能求出 \(ax+by=c\) 的特解,但是因為貝祖定理得 \(c=k\times\gcd(a,b),k\in Z\) 所以將兩邊同時乘 \(c/\gcd(a,b)\) 因為 \(c\mod\gcd(a,b)=0\) 所以 \(c/\gcd(a,b)\in Z\)。
又因為 \(x,y\in Z\) 所以 \(x\times(c/\gcd(a,b))\in Z\) 同理 \(y\times(c/\gcd(a,b))\in Z\)。所以不用擔心是不是整數的問題。
模板題:
https://www.luogu.com.cn/problem/CF7C
代碼:
#include<bits/stdc++.h>
using namespace std;
typedef __int128 ll;
ll a,b,x,y,c;
ll exgcd(ll a,ll b,ll &x,ll &y){
if(!b){
x=1,y=0;
return a;
}
ll d=exgcd(b,a%b,x,y);
ll t=x;
x=y;
y=t-(a/b)*y;
return d;
}
template<typename T>inline void rd(T&r){
r=0;char c=getchar(),m=1;
for(;!isdigit(c);c=getchar()){
if(c=='-')m=-1;
}
for(;isdigit(c);c=getchar()){
r=(r<<3)+(r<<1)+(c^48);
}
r*=m;
}
template<typename T>inline void wt(T r){
if(r<0){
putchar('-');wt(-r);return;
}
if(r>9) wt(r/10);
putchar(r%10+'0');
}
int main(){
rd(a);rd(b);rd(c);
c=-c;
ll g=exgcd(a,b,x,y);
if(c%g){
puts("-1");
return 0;
}
ll t=c/g;
wt(x*t);
putchar(' ');
wt(y*t);
putchar('\n');
return 0;
}

浙公網安備 33010602011771號