线性同余方程

题目链接
题意：给定n组数据ai,bi,mi，对于每组数求出一个xi，使其满足ai∗xi≡bi(mod mi)，如果无解则输出impossible。
输入格式
第一行包含整数n。
接下来n行，每行包含一组数据ai,bi,mi。
输出格式
输出共n行，每组数据输出一个整数表示一个满足条件的xi，如果无解则输出impossible。
每组数据结果占一行，结果可能不唯一，输出任意一个满足条件的结果均可。
输出答案必须在int范围之内。
数据范围
1≤n≤1e5,
1≤ai,bi,mi≤2∗1e9
输入样例：
2
2 3 6
4 3 5
输出样例：
impossible
7
思路：这个题要求找满足a * x % m = b的x值，变换一下就是
ax = my + b; 即是 ax + my’ = b;已知a,m,b求x,y’满足公式，不就是扩展欧几里得的问题嘛，ax’ + my’’ = d = gcd(a, m), 如果d | b 那么表示有解，x = (b/d)*x’，否则无解。
代码实现：

#include<iostream>
using namespace std;
typedef long long ll;
int exgcd(int a, int b, int &x, int &y)
{
    if(!b){
        x = 1; y = 0;
        return a;
    }
    //为了递归得到b*y + (a%b)* x = gcd(a, b), 所有这里的x和y要颠倒一下
    int d = exgcd(b, a % b, y, x);
    //由于a%b = a - (a/b)*b,
    //所有整理一下就有 a*x +b*(y - (a/b)*x) = d;
    y -= a / b * x;
    return d;
}
int main()
{
    int n;
    cin >> n;
    while(n --){
        int a, b, x, y, m;
        scanf("%d%d%d", &a, &b, &m);
        int d = exgcd(a, m, x, y);
        if(b % d) puts("impossible");
        else printf("%d\n", (ll)x * (b / d) % m);
    }
    return 0;
}