[NOIP 2012] 同余方程 · 拓展欧几里得

就是给你a,b两个值，求ax+by=1的一组解(x,y)，使得x是所有解中最小的正整数，然后输出x。

另外求出的x可能为负数，这时该怎么办？只要把x不断加上b直到x>0为止。这个问题的解决方法也是笔者当初一个没有理解的地方，后来无意中推算出来了：ax+by=ax+by+ab-ab=a(x+b)+b(y-a)，(据某神犇说其实是ab/gcd(a,b)，但因为a,b互质，所以gcd(a,)=1).a,b值是固定的，所以这样就构成了一个对应关系：对于确定a,b的方程ax+by=gcd(a,b),如果(x,y)是一组解，那么(x+b,y-a)也是一组解，而此时不用考虑y的值，所以就得到了上文那个方法，在程序中可以简写为: x=(x%b+b)%b。

#include<stdio.h>
#include<iostream>
#include<string.h>
#include<algorithm>
#include<time.h>
#include<stdlib.h>
using namespace std;
typedef long long ll;
int a,b;
int exgcd(int a,int b,int &x,int &y)
{
    if(b==0)
    {
        x=1;
        y=0;
        return a;
    }
    int t,p;
    t=exgcd(b,a%b,x,y);
    p=x;
    x=y;
    y=p-y*(a/b);
    return t;
}
int main()
{
    int x,y;
    cin>>a>>b;
    exgcd(a,b,x,y);
    cout<<(x%b+b)%b<<endl;//把x调整到和合适答案。
    return 0;
}