挑战程序设计竞赛:青蛙的约会

题目大意

在这里插入图片描述

解题思路

为了迎合公式,我们改变变量记号: 令 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) (m1n1)t(y1x1)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=m1n1,x=t,n=y1x1,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) axn(modb)
  • x b a s e = x ∗ n / d x_{base} = x * n /d xbase=xn/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
  • 0
    点赞
  • 0
    收藏
    觉得还不错? 一键收藏
  • 0
    评论

“相关推荐”对你有帮助么?

  • 非常没帮助
  • 没帮助
  • 一般
  • 有帮助
  • 非常有帮助
提交
评论
添加红包

请填写红包祝福语或标题

红包个数最小为10个

红包金额最低5元

当前余额3.43前往充值 >
需支付:10.00
成就一亿技术人!
领取后你会自动成为博主和红包主的粉丝 规则
hope_wisdom
发出的红包
实付
使用余额支付
点击重新获取
扫码支付
钱包余额 0

抵扣说明:

1.余额是钱包充值的虚拟货币,按照1:1的比例进行支付金额的抵扣。
2.余额无法直接购买下载,可以购买VIP、付费专栏及课程。

余额充值