给定 nn 组数据 ai,bi,mi,对于每组数求出一个 xi,使其满足 ai×xi≡bi(modmi),如果无解则输出 impossible
。
输入格式
第一行包含整数 n。
接下来 n 行,每行包含一组数据 ai,bi,mi。
输出格式
输出共 n 行,每组数据输出一个整数表示一个满足条件的 xi,如果无解则输出 impossible
。
每组数据结果占一行,结果可能不唯一,输出任意一个满足条件的结果均可。
输出答案必须在 int 范围之内。
数据范围
1≤n≤105,
1≤ai,bi,mi≤2×109
输入样例:
2
2 3 6
4 3 5
输出样例:
impossible
-3
/*题目:a*x%m = b,
则可化简为:a*x = m*y+b;
a*x+m*(-y) = b (扩展欧几里得)
用扩展欧几里得求出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;
}
int d = exgcd(b,a%b,y,x);
y-=a/b*x;
return d;
}
int main()
{
int t;
cin>>t;
while(t--)
{
int a,b,m,x,y;
cin>>a>>b>>m;
int d = exgcd(a,m,x,y);
if(b%d) cout<<"impossible"<<endl;
else cout<<(ll)b/d*x%m<<endl;
/*考虑线性同余方程 ax ≡ b (mod m),其中 a、b 和 m 都是整数。根据模运算的定义,两个整数 x 和 y 满足同余关系(即 x ≡ y (mod m))当且仅当它们除以 m 的余数相等。
求出的最大公约数d一定是b的倍数,d/b是整个a*x+m*(-y) = b扩大的倍数,x应该进行相应的扩大,%m是保证其在答案在int范围*/
}
}