【BZOJ】1477 青蛙的约会 扩欧

110 篇文章 0 订阅
1 篇文章 0 订阅

题目传送门

开坑开坑,从此入了数论的坑再也停不下来……话说数论什么的真的好难懂啊……

这题的解题思路就是扩展欧几里得算法。对于同余方程axc(mod b),可以转化成ax+by=c来求解。

回到这题,我们可以很快的发现题目要我们求(x+m*t)-(y-n*t)=Lp(t为步数,p为圈数)

然后我们推一下公式,可以发现(m-n)*t-Lp=y-x,然后就是扩欧的套路了:对于一个不定方程ax+by=c,可以先求出a*x0+b*y0=gcd(a,b),然后同时乘上c/gcd(a,b),得到a*x0*c/gcd(a,b)+b*y0*c/gcd(a,b)=c。这里解出来的x0*c/gcd(a,b)可能是个负数,转化成正数就行了。

然后发现这题就是一个扩欧的裸题,直接AC就行了。

附上AC代码:

#include <cstdio>
using namespace std;
 
int n,m,x,y,l,X,Y,t;
 
inline int exgcd(int a,int b,int &x,int &y){
    if (b==0) return x=1,y=0,a;
    int k=exgcd(b,a%b,x,y),t=x;
    return x=y,y=t-a/b*y,k;
}
 
int main(void){
    scanf("%d%d%d%d%d",&x,&y,&m,&n,&l);
    if ((x-y)%(t=exgcd((n-m),l,X,Y))) puts("Impossible");
    else printf("%lld\n",((long long)(x-y)/t*X%(l/t)+(l/t))%(l/t));
    return 0;
}

1、资源项目源码均已通过严格测试验证,保证能够正常运行; 2、项目问题、技术讨论,可以给博主私信或留言,博主看到后会第一时间与您进行沟通; 3、本项目比较适合计算机领域相关的毕业设计课题、课程作业等使用,尤其对于人工智能、计算机科学与技术等相关专业,更为适合; 、4下载使用后,可先查看README.md或论文文件(如有),本项目仅用作交流学习参考,请切勿用于商业用途。 5、资源来自互联网采集,如有侵权,私聊博主删除。 6、可私信博主看论文后选择购买源代码。 1、资源项目源码均已通过严格测试验证,保证能够正常运行; 2、项目问题、技术讨论,可以给博主私信或留言,博主看到后会第一时间与您进行沟通; 3、本项目比较适合计算机领域相关的毕业设计课题、课程作业等使用,尤其对于人工智能、计算机科学与技术等相关专业,更为适合;、下载 4使用后,可先查看README.md或论文文件(如有),本项目仅用作交流学习参考,请切勿用于商业用途。 5、资源来自互联网采集,如有侵权,私聊博主删除。 6、可私信博主看论文后选择购买源代码。 1、资源项目源码均已通过严格测试验证,保证能够正常运行; 2、项目问题、技术讨论,可以给博主私信或留言,博主看到后会第一时间与您进行沟通; 3、本项目比较适合计算机领域相关的毕业设计课题、课程作业等使用,尤其对于人工智能、计算机科学与技术等相关专业,更为适合;、 4下载使用后,可先查看README.md或论文文件(如有),本项目仅用作交流学习参考,请切勿用于商业用途。 5、资源来自互联网采集,如有侵权,私聊博主删除。 6、可私信博主看论文后选择购买源代码。
评论
添加红包

请填写红包祝福语或标题

红包个数最小为10个

红包金额最低5元

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

抵扣说明:

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

余额充值