定义:设a,b是任意的两个整数,其中b不等于0,如果存在一个整数q是的等式
a=bq成立就称为b整除a,计做b|a
定理1:设a,b。c不等于0,若c|a,c|b,则对任意的整数s,t,有c|(sa+tb)
证明:因为c|a c|b
所以a=cq1 b=cq2所以sa+tb=scq1+scq2所以sa+tb=c(sq1+tq2)=cq
所以c|(sa+tb)
定理2:设a,b,c是三个不全为零的整数,如果a=bq+c,其中q是整数,则(a,b)=(b,c)即a,b的
最大公因数等于b,c的最大公因数
证明:设a,b,的最大公因数是d b和c的最大公因数为e
所以d|a d|b 所以由定理1得d|(a-bq)=c即d|c又因为e为b和c的最大公因数所以d|e(这个先不证明)
同理可得e|d 因为d|e且e|d所以d=e所以(a,b)=(b,c)
基本算法:对于不完全为 0 的非负整数 a,b,gcd(a,b)表示 a,b 的最大公约数,必然存在整数对 x,y ,使得 gcd(a,b)=ax+by。
证明:设 a>b。
1,显然当 b=0,gcd(a,b)=a。此时 x=1,y=0;
2,ab!=0 时
设 ax1+by1=gcd(a,b);
bx2+(a mod b)y2=gcd(b,a mod b);
根据朴素的欧几里德原理有 gcd(a,b)=gcd(b,a mod b);
则:ax1+by1=bx2+(a mod b)y2;
即:ax1+by1=bx2+(a-(a/b)*b)y2=ay2+bx2-(a/b)*by2;
根据恒等定理得:x1=y2; y1=x2-(a/b)*y2;
这样我们就得到了求解 x1,y1 的方法:x1,y1 的值基于 x2,y2.
上面的思想是以递归定义的,因为 gcd 不断的递归求解一定会有个时候 b=0,所以递归可以结束。
#include<stdio.h>
int x,y,q;
void extend_Eulid(int a,int b)
{
if(b == 0)
{
x = 1;
y = 0;
q = a;
}
else
{
extend_Eulid(b,a%b);
int temp = x;
x = y;
y = temp - a/b*y;
}
}
int main()
{
int a,b;
while(scanf("%d%d",&a,&b)!=EOF)
{
extend_Eulid(a,b);
printf("%d %d\n",x,y);
}
return 0;
}
a=bq成立就称为b整除a,计做b|a
定理1:设a,b。c不等于0,若c|a,c|b,则对任意的整数s,t,有c|(sa+tb)
证明:因为c|a c|b
所以a=cq1 b=cq2所以sa+tb=scq1+scq2所以sa+tb=c(sq1+tq2)=cq
所以c|(sa+tb)
定理2:设a,b,c是三个不全为零的整数,如果a=bq+c,其中q是整数,则(a,b)=(b,c)即a,b的
最大公因数等于b,c的最大公因数
证明:设a,b,的最大公因数是d b和c的最大公因数为e
所以d|a d|b 所以由定理1得d|(a-bq)=c即d|c又因为e为b和c的最大公因数所以d|e(这个先不证明)
同理可得e|d 因为d|e且e|d所以d=e所以(a,b)=(b,c)
基本算法:对于不完全为 0 的非负整数 a,b,gcd(a,b)表示 a,b 的最大公约数,必然存在整数对 x,y ,使得 gcd(a,b)=ax+by。
证明:设 a>b。
1,显然当 b=0,gcd(a,b)=a。此时 x=1,y=0;
2,ab!=0 时
设 ax1+by1=gcd(a,b);
bx2+(a mod b)y2=gcd(b,a mod b);
根据朴素的欧几里德原理有 gcd(a,b)=gcd(b,a mod b);
则:ax1+by1=bx2+(a mod b)y2;
即:ax1+by1=bx2+(a-(a/b)*b)y2=ay2+bx2-(a/b)*by2;
根据恒等定理得:x1=y2; y1=x2-(a/b)*y2;
这样我们就得到了求解 x1,y1 的方法:x1,y1 的值基于 x2,y2.
上面的思想是以递归定义的,因为 gcd 不断的递归求解一定会有个时候 b=0,所以递归可以结束。
#include<stdio.h>
int x,y,q;
void extend_Eulid(int a,int b)
{
if(b == 0)
{
x = 1;
y = 0;
q = a;
}
else
{
extend_Eulid(b,a%b);
int temp = x;
x = y;
y = temp - a/b*y;
}
}
int main()
{
int a,b;
while(scanf("%d%d",&a,&b)!=EOF)
{
extend_Eulid(a,b);
printf("%d %d\n",x,y);
}
return 0;
}