poj 1061 青蛙的约会

最新推荐文章于 2022-04-14 00:39:57 发布

mars_ch

最新推荐文章于 2022-04-14 00:39:57 发布

阅读量569

点赞数 1

分类专栏：数论文章标签： poj

本文链接：https://blog.csdn.net/mars_ch/article/details/51307358

版权

数论专栏收录该内容

4 篇文章 0 订阅

订阅专栏

每一道中文题，都不是！天！下！掉！馅饼！！

这是一道数论题，我现在在写他的时候也还是极晕的。。

题目：
poj 1061

题意：
略。

分析：
看到这道题，想到的是数学问题：
很容易可以列出个关系式：
(x+m*T)-(y+n*T)=p*l【其中T是答案，p是任意整数代表几圈】】
化简一下：
x-y=p*l+(n-m)*T;
woc一看这个式子我激动啊，这不就是那个，就是内个
ax+by=c;吗。不是内个【【扩展欧几里得吗！！】】

然后我就懵逼了。woc这数据范围，woc这玩意怎么求！！

然后让我们来吃个饭，之后慢慢分析一下。

首先：
扩展欧几里得：对于两个不全为0的整数a、b，必存在一组解x,y，使得ax+by==gcd(a,b)；//借这个式子求出两个基本解，x0,y0,但还不能保证是不是正确答案！【【需要说明的是原题中的C可不一定是最大公约数哦。没准乘过k；
但是
同理：
若是c=k*gcd(a,b);
则有 axk+byk=c;
x1=x*k; y1=y*k;

我们可以由此得到x,y的通式：
x=x0*b/d;y=y0*b/d【【可以自己证明一下】】【【类比原式】】

然后就要求扩展欧几里得的gcd：
b*x+a%b*y=d;2)
变形:b*x+a*y-(a/b)*b*y=d;
再变:a*y+b(x-a/b*y)=d; 3)
与2）式比较:
得：
当a=b;b=a%b时：x=y;y=x-a/b*y;
调用过程中的x,y就是对应的a,b的解；
当回到顶层时，a,b就是最初的a，b所以此时的x，y就是所求解；

剩下的就是求最小值，这就很好说了；

#include<cstdio>
#include<algorithm>
#include<cstring>
#define LL long long
//by mars_ch 
using namespace std;
LL x,y,m,n,l,a,b,c,d,X,Y;
LL gcd(LL w,LL z)
{
    LL p,q;
    if(z == 0)
    {
        x=1;
        y=0;
        return w;
    }
    q=gcd(z,w%z);
    p=x;
    x=y;
    y=p-(w/z)*y;
    return q;
}
int main()
{
    scanf("%lld%lld%lld%lld%lld",&X,&Y,&m,&n,&l);
    a=n-m;
    b=l;
    c=X-Y;
    d=gcd(a,b);
    if(c%d !=0)
    {
        printf("Impossible\n");
    }
    else{
        x=x*(c/d);
        y=y*(c/d);
        LL ans=x*d/b; // x/(b/d)  整形出发 
        ans=x-ans*b/d;
        if(ans<0)            //可能是负数，不加样例就是-1，所以这样例给的还是很尽职尽责的。 
        {
            ans+=b/d;
        }
        printf("%lld\n",ans);
    }
}