扩展欧几里得 POJ 1061

感觉这道题目的数据好水啊。。。我的代码我都觉得姿势特别奇怪。。。竟然还过了。。。

好吧，原来不是姿势奇怪，而是逆元需要用的时候是余数也需要的时候，这里的余数是不需要的，所以就AC了

就说一下碰到的问题吧

设

x = (x0 + mt) % L;

y = (y0 + ny) % l;

然后x - y = 0;得到

(n - m) * t + k * l = x0 - y0;（t表示时间，k表示跑了几圈，l是一圈的长度）

然后我们和扩展欧几里得做比较:ax+by=gcd(a, b);（后面都写成gcd）

然后我们就可以发现x0-y0是已知的，那么令n-m=a，k=b，d=x0-y0，t=x，y=l；

因为d不一定是gcd，也有可能是gcd的倍数或者毫无关系

当毫无关系的时候就是impossible

当是gcd的倍数的时候，那么我们就对欧几里得进行改进。因此我们可以得到以下的式子：

(a*x*d/c) + (b*y*d/c)=d;

因此，我们的解就是x*d/c

但是这里有一个很重要的问题就是，当x时负数的时候怎么办，那么值也就是负的了呀。

因此我们想到对原来的式子进行改进，得到如下式子：

a*(x*d/c+p*b)+b*(y*d/c-p*a)=d;

所以我们的解就是((x*d/c)%b + b)% b;

转载于:https://www.cnblogs.com/heimao5027/p/5615213.html

扩展欧几里得 POJ 1061相关推荐

POJ - 1061 青蛙的约会(扩展欧几里得)
题目链接:点击查看题目大意:两只青蛙在一个单向循环数轴上跳动,给出初始位置和每秒跳动的距离以及数轴长度,问是否可以相遇,若能相遇求出最小时间题目分析:自从第一次接触扩展欧几里得以来已经有半年时间了 ...
POJ 1061扩展欧几里得
扩展欧几里得的模板题,需要注意的是为了得到一个最小正数解我们要使ax+by=c中的a,b都是正数 #include<cstdio> #include<cstring> #inc ...
POJ 1061 青蛙的约会 (扩展欧几里得)
青蛙的约会 Time Limit: 1000MS Memory Limit: 10000K Total Submissions: 97673 Accepted: 18409 Descripti ...
POJ - 2142 The Balance(扩展欧几里得)
题目链接:点击查看题目大意:给出两种重量的砝码,我们需要利用天平称出指定重量的药物,我们需要求出这两种砝码各自的数量,使得砝码数量之和最小题目分析:我们可以先列出关系式,假设两个砝码的重量分别为a ...
POJ 2142——扩展欧几里得
题目是很裸的扩展欧几里得,但是对x,y有限制条件,要求所有x,y中abs(x)+abs(y)最小,在这个条件下要求abs(a* x)+abs(b* y)最小显然我们需要用扩展欧几里得求得一组解,问题 ...
欧几里得算法扩展欧几里得
本文包括: 1.几里德算法递归实现 2.扩展欧几里德算法递归实现 -实际上就是把别人总结的,我认为有助于自己理解的内容copy过来,再加上几句自己的理解. 欧几里得算法欧几里德算法又称辗转相除法 ...
关于扩展欧几里得求最小正整数解
关于扩展欧几里得求最小正整数解 //扩展欧几里得算法求最小正整数解这里讲的是欧几里得,不是那个洗澡测浮力的阿基米德 . 什么是欧几里得算法: 扩展欧几里得算法是用来求解一类特定的不定方程的,形如ax ...
poj-青蛙的约会（扩展欧几里得）nyoj-小柯的约会
由题意可列出方程 x + m * s - (y + n * s) = k * l; 化简得 (n-m)*s + k*l = x - y 化简得 (n-m)*s = x - y (mod l); ...
欧几里得定理与扩展欧几里得定理
欧几里的定理(辗转相除法): gcd(a,b) = gcd(b,a%b) gcd(a,b)表示a,b的最大公约数证明: 设 a > b c = a%b a = k * b + c (k为某个整 ...

扩展欧几里得 POJ 1061

扩展欧几里得 POJ 1061相关推荐

最新文章

热门文章