[Luogu P1516]青蛙的约会

按照题意,显然可以列出同余方程,k即为所求天数,再将其化为不定方程 ,那么对这个方程用扩展欧几里德算法即可得出k,q的一组解,但是方程有解的充要条件是(m – n) 和L不同时为零并且x – y是m – n和L的因子,扩展欧几里德算出的解才是方程的解

#include<cstdio>
#include<cmath>
#include<cstring>
#include<algorithm>
#include<string>
#include<stack>
#include<queue>
#include<vector>
#include<set>
#include<ctime>
#include<cstdlib>
#define ll long long
using namespace std;
void exgcd(ll a,ll b,ll &d,ll &x,ll &y)
{
    if(!b)
    {
        x = 1;
        y = 0;
        d = a;
        return;
    }
    exgcd(b,a % b,d,x,y);
    ll tmp = x;
    x = y;
    y = tmp - (a / b) * y;
    return;
}
int main()
{
    ll x,y,m,n,l,k,q,d;
    scanf("%lld %lld %lld %lld %lld",&x,&y,&m,&n,&l);
    if(n < m)
    {
        swap(n,m);
        swap(x,y);
    }
    ll a =  n - m;
    exgcd(a,l,d,k,q);
    if((x - y) % d != 0 || m == n) //判断方程有解的条件
    {
        printf("Impossible");
        return 0;
    }
    printf("%lld",(k * (x - y) / d % (l / d) + (l / d)) % (l/ d));//为得到一组正数解
    return 0;
}

原文地址:https://www.cnblogs.com/lijilai-oi/p/10691662.html

时间: 2024-10-18 09:07:18

[Luogu P1516]青蛙的约会的相关文章

luogu P1516 青蛙的约会(线性同余方程扩展欧几里德)

题意 题解 做了这道题,发现扩欧快忘了. 根据题意可以很快地列出线性同余方程. 设跳了k次 x+mkΞy+nk(mod l) (m-n)kΞ-(x-y)(mod l) 然后化一下 (m-n)k+(x-y)Ξ0(mod l) 也就是前面一坨是l的倍数 不妨设 (m-n)k+(x-y)=-tl (m-n)k+tl=-(x-y) 我们要求的就是保证t<=0(因为我们设的-t倍的l,所以t<=0),k>=0时k的最小值 发现这是一个不定方程 根据裴蜀定理(这个定理搜狗输入法上没有) 当-(x-y

[Luogu 1516] 青蛙的约会

[Luogu 1516] 青蛙的约会 两年前和初中那群队友一起当模拟做过,一起 TLE,一起骂过的题目: 一年前队友的博客上被当作例题,被我抄下来去操场都带着看的题目: 终于好好学习了 ExGCD,才发现是扩欧板子题并顺手切掉的题目- 雨中灯市欲眠,原已萧萧数年. 似有故人轻叩,又将旧题重写. 能否消得你一路而来的半生风雪? #include <cstdio> #include <cstdlib> long long a, b, c, x, y, ans; void Read(vo

【Luogu】P1516青蛙的约会(线性同余方程,扩展欧几里得)

题目链接 定理:对于方程\(ax+by=c\),等价于\(a*x=c(mod b)\),有整数解的充分必要条件是c是gcd(a,b)的整数倍. --信息学奥赛之数学一本通 避免侵权.哈哈. 两只青蛙跳到一格才行,所以说 \(x+mt=y+nt(mod l) \) \((x-y)+(m-n)t=0(mod l)\) \((m-n)t+ls=(y-x)  s属于整数集\) 令a=n-m,b=l,c=gcd(a,b),d=x-y 则有\( at+bs=d\) 扩展欧几里得求解. 设c=gcd(a,b)

P1516 青蛙的约会

题目描述 两只青蛙在网上相识了,它们聊得很开心,于是觉得很有必要见一面.它们很高兴地发现它们住在同一条纬度线上,于是它们约定各自朝西跳,直到碰面为止.可是它们出发之前忘记了一件很重要的事情,既没有问清楚对方的特征,也没有约定见面的具体位置.不过青蛙们都是很乐观的,它们觉得只要一直朝着某个方向跳下去,总能碰到对方的.但是除非这两只青蛙在同一时间跳到同一点上,不然是永远都不可能碰面的.为了帮助这两只乐观的青蛙,你被要求写一个程序来判断这两只青蛙是否能够碰面,会在什么时候碰面. 我们把这两只青蛙分别叫

P1516 青蛙的约会 洛谷

https://www.luogu.org/problem/show?pid=1516 题目描述 两只青蛙在网上相识了,它们聊得很开心,于是觉得很有必要见一面.它们很高兴地发现它们住在同一条纬度线上,于是它们约定各自朝西跳,直到碰面为止.可是它们出发之前忘记了一件很重要的事情,既没有问清楚对方的特征,也没有约定见面的具体位置.不过青蛙们都是很乐观的,它们觉得只要一直朝着某个方向跳下去,总能碰到对方的.但是除非这两只青蛙在同一时间跳到同一点上,不然是永远都不可能碰面的.为了帮助这两只乐观的青蛙,你

洛谷 P1516 青蛙的约会(中国剩余定理CRT)

题目描述 两青蛙会相遇的条件[(x-y)+k(m-n)]%l=0 令a=m-n,b=l,c=x-y,所以(c+k*a)%b=0 -> aX+bY=c 方程有解,当且仅当c%Gcd(a,b)=0 令r=Gcd(a,b) 为什么(X*(c/r)%(b/r)+b/r)%(b/r)为最小解? 目标解aX+bY==c 有解的条件是c%r=0,也就是c=K*r 那么就办成了aX+bY=K*r 那么也就是解a*(X/K)+b(Y/K)=r 然后扩展欧几里得定理就是解a*X+b*Y=r 那么解出来X之后X就要乘

题解 P1516 【青蛙的约会】

题目链接 Solution 青蛙的约会 题目大意:求解不定方程\(ax+by=c\) 分析:我们可以把原来的同余式子写成一个不定方程,这部分基本操作不讲,主要讲方程求解.看到不定方程我们就想到\(exgcd\)对吧? 但是\(exgcd\)只能适用于求解\(ax+by=g\),其中\(g=gcd(a,b)\)的情况 我们设\(exgcd\)求出的一组特解是\((x_0,y_0)\),显然对于方程\(ax+by=c\)的一组解就是\((x_oc/g,y_0c/g)\)对吧?如果\(g \nmid

POJ - 1061 青蛙的约会 (扩展欧几里得算法)

Description 两只青蛙在网上相识了,它们聊得很开心,于是觉得很有必要见一面.它们很高兴地发现它们住在同一条纬度线上,于是它们约定各自朝西跳,直到碰面为止.可是它们出发之前忘记了一件很重要的事情,既没有问清楚对方的特征,也没有约定见面的具体位置.不过青蛙们都是很乐观的,它们觉得只要一直朝着某个方向跳下去,总能碰到对方的.但是除非这两只青蛙在同一时间跳到同一点上,不然是永远都不可能碰面的.为了帮助这两只乐观的青蛙,你被要求写一个程序来判断这两只青蛙是否能够碰面,会在什么时候碰面. 我们把这

POJ1061——青蛙的约会(扩展欧几里德)

青蛙的约会 Description两只青蛙在网上相识了,它们聊得很开心,于是觉得很有必要见一面.它们很高兴地发现它们住在同一条纬度线上,于是它们约定各自朝西跳,直到碰面为止.可是它们出发之前忘记了一件很重要的事情,既没有问清楚对方的特征,也没有约定见面的具体位置.不过青蛙们都是很乐观的,它们觉得只要一直朝着某个方向跳下去,总能碰到对方的.但是除非这两只青蛙在同一时间跳到同一点上,不然是永远都不可能碰面的.为了帮助这两只乐观的青蛙,你被要求写一个程序来判断这两只青蛙是否能够碰面,会在什么时候碰面.