扩展欧几里德 noip 2012 day2 同余方程

博主回来啦！（？）

开始一年两度的集训啦，刷一些noip的题，当然我还没参加过普及组嘿。。。

所以树状数组的更新可能会延迟。。

好啦，今天更新的是noip提高组（我也不知道为什么刷提高组）2012 day2同余方程，听起来特别高大上，但是并没有。和第二题借教室一比，题目的长短也是惊人，这么看来，这道题一定是考数学的，题目描述如下：求关于x的同余方程ax≡ 1 (mod b)的最小正整数解。

拿到题目之后我有点懵。。。

等等，那个像三的东西是什么，一脸茫然，我开始方了。。同学看了看，诶，这不就是恒等于吗，这你都不会。。。我真的不会啊，欲哭无泪。其实这个所谓的恒等于和等于的区别也没有很大。反正这道题里，没区别。这道题讲的就是a*x模b=1，对，没错，听起来就这么简单。当然，所有题都可以暴力，只是看你分得的多少了，而这道题十分凶残的要用扩欧，当然，扩欧是什么？？？

严格来讲，叫扩展欧几里得算法，欧几里德算法，对于大家一定太过熟悉了，辗转相除求gcd，扩展欧几里得就是a*x+b*y=gcd(a,b)。

那么和这道题又有什么联系呢？

gcd(a,b)=gcd(b,a%b)=...=gcd(c,0)=c=ax+by

gcd(a,b)=gcd(b,a%b)=c

ax+by=bx'+a%by'

ax+by=bx+(a-a/b*b)y'

     =ay+b(x'-a/by')

x=y'

y=x'-a/b*y'

推导出的算法，当然光堆公式是不行的，首先，gcd到最后，后面的b项一定是0，最后辗转相除的结果是a项，然后a项一定是c。然后可以推导出前面的一步，然后一直推导下去，gcd(a,b)=gcd(b,a%b)，在每步设一个x'和y'，学过代数的一定是懂的。最后就可以推出最后的答案，也就是原始不带撇的x和y。好好看看上面的数学推导式！

下面是ac代码：

#include<iostream>
using namespace std;
void ss(long long a,long long b,long long&x,long long&y)
{
    if(b==0)
        x=1,y=0;
    else
        ss(b,a%b,y,x),y-=x*(a/b);
}
int main()
{   long long a,b,x,y;
    scanf("%lld%lld",&a,&b);
    ss(a,b,x,y);
    if(a*x+b*y==1) printf("%lld",(x%b+b)%b);
    else printf("%lld",(y%b+b)%b);
    return 0;
}

对没错，提高组day2的第一题代码就这么短，只要把代码弄懂就行了，如果扩欧的思想的推导没弄懂也没关系，因为也算是纯数学了，看题目描述就能看出来，所以只要记住模板，加油！！

下篇博文再见ヾ(￣▽￣)Bye~Bye~