青蛙的约会(exgcd扩展欧几里得,思维)
创始人
2025-06-01 03:07:02
0

地址:https://www.luogu.com.cn/problem/P1516

青蛙的约会

题目描述

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

我们把这两只青蛙分别叫做青蛙 A 和青蛙 B,并且规定纬度线上东经 000 度处为原点,由东往西为正方向,单位长度 111 米,这样我们就得到了一条首尾相接的数轴。设青蛙 A 的出发点坐标是 xxx,青蛙 B 的出发点坐标是 yyy。青蛙 A 一次能跳 mmm 米,青蛙 B 一次能跳 nnn 米,两只青蛙跳一次所花费的时间相同。纬度线总长 LLL 米。现在要你求出它们跳了几次以后才会碰面。

输入格式

输入只包括一行五个整数 x,y,m,n,Lx,y,m,n,Lx,y,m,n,L。

输出格式

输出碰面所需要的次数,如果永远不可能碰面则输出一行一个字符串 Impossible

样例 #1

样例输入 #1

1 2 3 4 5

样例输出 #1

4

提示

对于 100%100\%100% 的数据,1≤x≠y≤2×1091 \le x \ne y \le 2 \times 10^{9}1≤x=y≤2×109,1≤m,n≤2×1091 \le m, n \le 2 \times 10^{9}1≤m,n≤2×109,1≤L≤2.1×1091 \le L \le 2.1 \times 10^{9}1≤L≤2.1×109。

思路:
根据题意首先可以推出(mt+x)%L=(nt+y)%L ,然后可以转化成(m-n)t+kL=(y-x).(这里的t和k是未知数),所以可以转化成ax+by=c,然后就可以进行扩展欧几里得定理,求出x,(这里的a=(m-n),b=L,c=y-x.),
利用扩展欧几里得到的数是gcd(a,b),如果这个数不能够被c整除,那么这两只青蛙无论怎么跳都不能够相遇,欧几里德的公式是ax + by = gcd(a, b),但是本题需要求的是ax + by = c, 所以就有了 (ax + by) * c / gcd = gcd * c / gcd; 所以就要x * c / gcd,而且x的通解是x +k * b / gcd,所以就要取模b / gcd才能得到最小x,为了避免x是一个负数,我们应该对(x%(b/gcd)+(b/gcd)%(b/gcd),这样得到我们想要的答案.
AC代码:

#include
#define lmw ios::sync_with_stdio(0);cin.tie(0);cout.tie(0);
using namespace std;
#define int long long
int exgcd(int a,int b,int &x,int &y){if(!b){x=1,y=0;return a;}int gcd=exgcd(b,a%b,x,y);int z=x;x=y;y=z-a/b*y;return gcd;
}
signed main(){lmw;int l,r,m,n,L;cin>>l>>r>>m>>n>>L;int a=m-n;int b=L;int c=r-l;if(a<0){a=-a;c=-c;}int x,y;int op=exgcd(a,b,x,y);if(c%op!=0){cout<<"Impossible\n";}else {x=x*c/op;b/=op;x=(x%b+b)%b;cout<

本题完结。

相关内容

热门资讯

linux入门---制作进度条 了解缓冲区 我们首先来看看下面的操作: 我们首先创建了一个文件并在这个文件里面添加了...
C++ 机房预约系统(六):学... 8、 学生模块 8.1 学生子菜单、登录和注销 实现步骤: 在Student.cpp的...
JAVA多线程知识整理 Java多线程基础 线程的创建和启动 继承Thread类来创建并启动 自定义Thread类的子类&#...
【洛谷 P1090】[NOIP... [NOIP2004 提高组] 合并果子 / [USACO06NOV] Fence Repair G ...
国民技术LPUART介绍 低功耗通用异步接收器(LPUART) 简介 低功耗通用异步收发器...
城乡供水一体化平台-助力乡村振... 城乡供水一体化管理系统建设方案 城乡供水一体化管理系统是运用云计算、大数据等信息化手段࿰...
程序的循环结构和random库...   第三个参数就是步长     引入文件时记得指明字符格式,否则读入不了 ...
中国版ChatGPT在哪些方面... 目录 一、中国巨大的市场需求 二、中国企业加速创新 三、中国的人工智能发展 四、企业愿景的推进 五、...
报名开启 | 共赴一场 Flu... 2023 年 1 月 25 日,Flutter Forward 大会在肯尼亚首都内罗毕...
汇编00-MASM 和 Vis... Qt源码解析 索引 汇编逆向--- MASM 和 Visual Studio入门 前提知识ÿ...
【简陋Web应用3】实现人脸比... 文章目录🍉 前情提要🌷 效果演示🥝 实现过程1. u...
前缀和与对数器与二分法 1. 前缀和 假设有一个数组,我们想大量频繁的去访问L到R这个区间的和,...
windows安装JDK步骤 一、 下载JDK安装包 下载地址:https://www.oracle.com/jav...
分治法实现合并排序(归并排序)... 🎊【数据结构与算法】专题正在持续更新中,各种数据结构的创建原理与运用✨...
在linux上安装配置node... 目录前言1,关于nodejs2,配置环境变量3,总结 前言...
Linux学习之端口、网络协议... 端口:设备与外界通讯交流的出口 网络协议:   网络协议是指计算机通信网...
Linux内核进程管理并发同步... 并发同步并发 是指在某一时间段内能够处理多个任务的能力,而 并行 是指同一时间能够处理...
opencv学习-HOG LO... 目录1. HOG(Histogram of Oriented Gradients,方向梯度直方图)1...
EEG微状态的功能意义 导读大脑的瞬时全局功能状态反映在其电场结构上。聚类分析方法一致地提取了四种头表面脑电场结构ÿ...
【Unity 手写PBR】Bu... 写在前面 前期积累: GAMES101作业7提高-实现微表面模型你需要了解的知识 【技...