poj 1061 青蛙的约会(扩展欧几里得)

          中文题就不说题意了吧……最近打算学学数论,数学相关的真是一点都不会啊Orz,首先是一发基础的……这题以前做过,但是当时做得稀里糊涂的,连扩展欧几里得还不太明白就瞎搞搞过了,今天又做了一遍,把公式神马的都重推了一遍,感觉不错,唉,以后做题要注意了,不能没搞懂就过去了。。。。

          设两只青蛙一共跳了t次,那么就可以得到一个方程:(a+mt)-(b+nt)=kL   (k=0,1,2,3,4......),移项,可以得到:(m-n)t-Lk=y-x,t和k为未知数,剩下就是扩展欧几里得直接求方程的解了,这里求出的t不一定最优,对L取模即可。


代码:


#include<iostream>
#include<cstdio>
#include<cstring>
#include<string>
#include<algorithm>
#include<map>
#include<queue>
#include<set>
#include<stack>
#include<cmath>
#include<vector>
#define inf 0x3f3f3f3f
#define Inf 0x3FFFFFFFFFFFFFFFLL
#define eps 1e-9
#define pi acos(-1.0)
using namespace std;
typedef long long ll;
ll ext_gcd(ll a,ll b,ll &x,ll &y)
{
    if(b==0)
    {
        x=1;y=0;
        return a;
    }
    ll d=ext_gcd(b,a%b,y,x);
    y-=a/b*x;
    return d;
}
int main()
{
    //freopen("in.txt","r",stdin);
    //freopen("out.txt","w",stdout);
    ll a,b,m,n,L,x,y;
    cin>>a>>b>>m>>n>>L;
    ll A=m-n,B=-L,C=b-a;
    ll d=ext_gcd(A,B,x,y);
    if(C%d!=0)
        cout<<"Impossible"<<endl;
    else
    {
        x=x*C/d;
        if(B<0) B=-B;
        x=(x%B+B)%B;
        cout<<x<<endl;
    }
    return 0;
}

  • 0
    点赞
  • 0
    收藏
    觉得还不错? 一键收藏
  • 0
    评论

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

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

请填写红包祝福语或标题

红包个数最小为10个

红包金额最低5元

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

抵扣说明:

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

余额充值