P1082 [NOIP2012 提高组] 同余方程

---

原题链接

P1082
题目类型：$普及+/提高$
AC记录：Accepted

题目大意

求关于$x$的同余方程$ax\equiv 1(\mathrm{m}\mathrm{o}\mathrm{d}b)$的最小正整数解。

输入格式

一行，包含两个正整数$a,b$，用一个空格隔开。

输出格式

一个正整数$x_0$，即最小正整数解。输入数据保证一定有解。

$\mathbf{S}\mathbf{a}\mathbf{m}\mathbf{p}\mathbf{l}\mathbf{e}$ $\mathbf{I}\mathbf{n}\mathbf{p}\mathbf{u}\mathbf{t}$

3 10

$\mathbf{S}\mathbf{a}\mathbf{m}\mathbf{p}\mathbf{l}\mathbf{e}$ $\mathbf{O}\mathbf{u}\mathbf{t}\mathbf{p}\mathbf{u}\mathbf{t}$

7

$\mathbf{H}\mathbf{i}\mathbf{n}\mathbf{t}\&\mathbf{E}\mathbf{x}\mathbf{p}\mathbf{l}\mathbf{a}\mathbf{i}\mathbf{n}$
无

数据范围

对于$40\%$的数据，$2\le b\le 1,000$。
对于$60\%$的数据，$2\le b\le 50,000,000$。
对于$100\%$的数据，$2\le a,b\le 2,000,000,000$。

解题思路

这题可以用拓展欧几里得算法（下面简称拓欧）来解决。

如果不会欧几里得算法或拓展欧几里得算法的童鞋，可以戳我

而拓欧解的方程是这样的：$ax+by=\gcd (a,b)$，与我们的题目有很大的区别。所以我们要把题目中的同余方程变成这个样子：$ax+by=1$。其中的$y$为辅助解，没什么用。又由于题目中说了：输入数据保证一定有解，所以又可以得知：$\gcd (a,b)=1$。

---

解完方程以后，还有一个需要注意的点，就是题目要求的是最小正整数解，那么我们解出一个解后，如何找出其他的解呢？
答案是：设方程的一个解为$x$，则$x+kb$都是方程的解。
如果把$x+kb$代入原方程，可以得到$a(x+kb)\equiv 1(\mathrm{m}\mathrm{o}\mathrm{d}b)$，去括号后变成$ax+kab\equiv 1(\mathrm{m}\mathrm{o}\mathrm{d}b)$。由于$kab$一定是$b$的倍数，所以方程又变成了$ax\equiv 1(\mathrm{m}\mathrm{o}\mathrm{d}b)$，即为原方程。所以求最小正整数解就相当于把$x$变成$((x+\mathrm{m}\mathrm{o}\mathrm{d}b)+b)\mathrm{m}\mathrm{o}\mathrm{d}b$。
注意：由于解出来的$x$有可能是负数，所以上式在$\mathrm{m}\mathrm{o}\mathrm{d}b$之后还要再加上一个$b$再取模。

---

最后，祝大家早日

上代码

#include<iostream>

using namespace std;

int              x,y;

int exgcd(int a,int b)
{
    if(b==0)
    {
        x=1,y=0;
        return a;
    }
    int temp=exgcd(b,a%b);
    int tx=x;
    x=y;
    y=tx-a/b*y;
    return temp;
}

int main()
{
    int a,b;
    cin>>a>>b;
    exgcd(a,b);
    // cout<<x<<" "<<y<<endl;
    x=((x%b)+b)%b;
    cout<<x<<endl;
    return 0;
}

完美切题$\sim$