题目大意
解题思路
为了迎合公式,我们改变变量记号: 令 x 1 , y 1 , m 1 , n 1 , L 1 x1,y1,m1,n1,L1 x1,y1,m1,n1,L1分别表示题意中的 x , y , m , n , L x,y,m,n,L x,y,m,n,L.
- 分析题意可知
- 故得出公式: ( m 1 − n 1 ) t ≡ ( y 1 − x 1 ) m o d ( L 1 ) (m_1-n_1)t \equiv (y_1 - x_1) mod(L_1) (m1−n1)t≡(y1−x1)mod(L1)
- 故我们令: a = m 1 − n 1 , x = t , n = y 1 − x 1 , b = L 1 a = m_1-n_1, x=t,n=y_1-x_1, b=L_1 a=m1−n1,x=t,n=y1−x1,b=L1.
然后借助公式求出最小次数
代码
#include<iostream>
using namespace std;
typedef long long ll;
ll extend_ojld(ll a, ll b, ll &x, ll &y)
{
ll d;
if(b == 0)
{
d = a;
x = 1;
y = 0;
}
else
{
d = extend_ojld(b, a%b, y, x);
y -= (a/b)*x;
}
return d;
}
int main()
{
ll x1,y1,m1,n1,L1;
cin >> x1 >> y1 >> m1 >> n1 >> L1;
ll a, x, y, b, n;
b = L1;
a = ((m1 - n1) % b + b) %b; // 防止为负数
n = ((y1 - x1) % b + b) %b; // 防止为负数
ll d = extend_ojld(a, b, x, y);
if(n % d)
cout << "Impossible" << endl;
else
{
ll x_base = x * n / d;
ll s = b / d;
cout << (x_base % s + s) % s << endl;
}
}
知识点
- 位置对应工整 a x ≡ n ( m o d b ) ax \equiv n(mod b) ax≡n(modb)
- x b a s e = x ∗ n / d x_{base} = x * n /d xbase=x∗n/d
- s = b / d s = b/d s=b/d
- ( x b a s e % s + s ) % s (x_{base}\%s+s)\%s (xbase%s+s)%s