菜鸡的初入数论（1）——欧几里得与扩展欧几里得小结

emmmm..明明是拿来自己总结的。。却总感觉想要写的像给别人看。。。突然感觉好丢脸。。好了进入正题。
欧几里得算法是指欧几里得用来求最大公因数的方法——辗转相除法。既：gcd(a,b)=gcd(b（a1),a%b(b1))。当b1为零时，则a1为a和b的最大公因数；
代码如下：int gcd(int a,int b)
{
if(b==0)
return a;
else
return gcd(b,a%b);
}这里不需要关心a和b的大小关系，因为假如a<b的话，a%b=a;下一次的时候自动就转过去了，所以不需要比大小（好吧。。我经常这么干），对于关于gcd(a,b)=gcd(b,a%b)的证明的话，如下：
设c=gcd(a,b);d=gcd(b,a%b);
因为d是b和a%b的公因数，所以有
b%d=0;（1）
(a%b)%d=0;（2）
又因为a%b=a-（a/b）*b；（/为整除，既抹去后面小数点）
得出a=(a%b)+（a/b）*b（3）
由（1）（2）（3）可得出
a%d=0；(4)
4,5可得出d为a，b公因数
又因为c为a，b的最大公因数
所以有d<=c；
b%c=0；（5）
a%c=0；（6）
由（5）（6）（3）可得出

（a%b）%c=0；（7）
5，7可得出c为a%b，b的公因数
又因为d为最大公因数
所以有c<=d
所以c=d；
欧几里得算法就讲到这里，接下来说扩展欧几里得；
扩展欧几里得，用于求解不定方程：ax+by=gcd
因为是不定方程，所以多解是肯定的（废话），顺便一提，他好像肯定有解。。。。不过我不知道怎么证明；
但他的解肯定为如下：
x=x0+(b/gcd)*t
y=y0-(a/gcd)*t
其中x0,y0为一个特定解，t为一个整数。
这块比较好解释，假设t=1（懒得打t了）
a*x+b*y=a*x0+(b/gcd)*a+y0*b-(a/gcd)*b=a*x0+b*y0=gcd （又是废话。。。不过怕以后看不懂啊。。）
至于这里为什么要使用b/gcd和a/gcd，emmmm。。。。。gcd是a和b的最大公因数，所以a/gcd和b除gcd应该是互相为素数的。。。。。所以这应该就是最小的系数了。。
接下来就是扩展欧几里得的详细解释，先上代码：void exgcd (int a1,int b1,int gcd,int &x,int &y)//a1*x+b1*y=gcd
{
if(b1==0)
{
x=1;
y=0;
}
else
{
exgcd(b1,a1%b1,gcd,x,y);
int v=x;
x=y;
y=v-(a1/b1)*y;
}//逆元时a1为需要逆元数，b1为模，x为得数；&为引用，否则无法录入。。。应该没有人和我一样不知道吧。。
}按照代码解释：
首先当b1=0的时候，a1*x+b1*y=gcd 只有一个解了应该，a1=gcd，x=1
而当b1!=0的时候：
a*x+b*y=gcd；（1
由欧几里得定理可以得知
gcd(a,b)==gcd(b,a%b);
这个证明在上面哦
ok，继续
接下来由于a*x+b*y=gcd(a,b)（这里写的全一些，容易理解），
所以有b*x0+(a%b)*y0=gcd(b,a%b)=gcd(a,b);（反正都一样了，下面就继续用gcd代替了哦）
并且由求余本质(a%b=a-b*(a/b)((这里的‘/’为整除哦)))可以得到下列式子：
b*x0+(a-b*(a/b))*y0=gcd; 
在变形，得到下列式子：
a*y0+b*(x0-(a/b)*y0)=gcd;（2
由（1，（2放在一起可以得出：

x=y0;
y=x0-(a/b)*y0;
emmmm,大功告成，接下来无限递归到b1==0的时候就可以了！
证明终了！（一直想说一次这句话呢）

好吧，还有逆元。。。等下次写逆元专题一块放吧。。。懒得动弹