NOIP 2012 同余方程 Ex_gcd

题目描述

求关于 x 的同余方程 ax ≡ 1 (mod b)的最小正整数解。

输入格式：

输入只有一行，包含两个正整数 a, b，用一个空格隔开。

输出格式：

输出只有一行，包含一个正整数 x0，即最小正整数解。输入数据保证一定有解。

输入样例#1：

3 10

输出样例#1：

7

题解

对于$ax≡ 1(\mod b)$可以转化为: $ax-kb= 1$,我们可以直接用$Ex_Gcd()$,来求出一组解，然后套用同解公式不断的变小就可以。
注意开$long long$
不会扩展欧几里得的戳

#include<iostream>
#include<cstdlib>
#include<cstring>
#include<cstdio>
#include<cmath>
#define ll long long 
using namespace std;
ll a,b;
ll Gcd(ll a,ll b)
{
    if(!b) return a;
    else return  Gcd(b,a%b);
}
void Ex_Gcd(ll a,ll b,ll &x,ll &y)
{
    if(!b)
    {
        x=1;y=0;
        return ;
    }
    Ex_Gcd(b,a%b,x,y);
    int temp=x;
    x=y;
    y=temp-(a/b)*y;
    return ;
}
int main()
{
    scanf("%lld%lld",&a,&b);
    ll gcd=Gcd(a,b);
    ll x,y;
    Ex_Gcd(a,b,x,y);
    ll lcm=a*b/gcd;//通解
    while(x<0) x+=lcm/a;
    while(x-(lcm/a)>0) x-=lcm/a;
    printf("%lld",x);
    return 0;
}