先上干货:

定理1:

如果d = gcd(a,b),则必能找到正的或负的整数k和l,使ax + by = d.

(参考exgcd:http://www.cnblogs.com/dilthey/p/6804137.html

定理2:

一元线性同余方程ax ≡ n (mod b) 有解,当且仅当gcd(a,b)|n.

也就是说,解出了ax+by=gcd(a,b),就相当于解出了ax≡n(mod b) (而且只要满足gcd(a,b)|n,就一定有解)

定理3:

若gcd(a,b) = 1,则方程ax ≡ n (mod b)在[0, b-1]上有唯一解.

由上可知ax ≡ n (mod b) 即方程 ax + by = n,由定理1可知必然存在x,y满足 ax + by = gcd(a,b) = 1,则必然存在nx,ny满足 ax + by = n

则必然ax ≡ n (mod b)有解nx,由于nx+kb都是该方程的解,则在[0, b-1]上必然出现一个解。

再证唯一性:

唯一满足mod b = 0的条件的,只有x1 = x2

定理4:

若gcd(a, b) = d,则方程ax ≡ n (mod b)在[0, b/d - 1]上有唯一解.

这样,如果我有一个方程ax ≡ n (mod b)的任意解X,那么我先X mod (b/d),使得其在范围[-b/d+1,b/d-1]范围内,

再加上b/d就使其在范围[1,2(b/d)-1]范围内,再mod (b/d)就可以使其处在[0,b/d -1]范围内了。

即在[0,b/d-1]范围内的解x=[ X mod (b/d) + (b/d) ] mod (b/d).

分割线

然后正式关于题目:

题目链接:http://poj.org/problem?id=1061

题解:

青蛙A:a = ( x + mt ) mod L

青蛙B:b = ( y + nt ) mod L

现在要使得a=b,即( x + mt ) = ( y + nt ) (mod L),即( x - y ) + ( m - n ) t = kL;

即求最小的正整数t满足 ( m - n ) t + KL = y - x (K=-k),即求一元线性同余方程( m - n ) x ≡ ( y - x ) (mod L) 的解;

那么就像上面那样,先算( m - n ) * x + L * K = gcd( m - n , L ),得到一个x的值;

然后判断( y - x ) mod gcd( m - n , L ) == 0 ?,若能整除,则有解;

然后就只要先把x *= ( y - x ) / gcd( m - n , L ) ,再根据上面求出解集中的最小正整数即为答案;

PS.当然,最后一步求[0,b/d-1]范围内的唯一解时,由于要mod (b/d),我们要保证b/d为正,在本题中,即保证L/d为正,故d要为正。

AC代码:

 #include<cstdio>
#include<algorithm>
typedef long long ll;
using namespace std;
ll x,y,m,n,L;
ll d,ans,K;
void exgcd(ll a,ll b,ll &d,ll &x,ll &y){
if(!b){d=a;x=;y=;}
else {exgcd(b,a%b,d,y,x);y-=(a/b)*x;}
}
int main()
{
scanf("%lld%lld%lld%lld%lld",&x,&y,&m,&n,&L);
if(m<n)
{
swap(m,n);
swap(x,y);
}
exgcd(m-n,L,d,ans,K);
if(d== || (x-y)%d!=) printf("Impossible\n");
else
{
ans*=(y-x)/d;
ans=(ans%(L/d)+(L/d))%(L/d);
printf("%I64d\n",ans);
}
}

PS.曾经刚开始搞ACM的时候,就因为这题是POJ第一页上为数不多的中文题而钻了很久,不过以当时的实力,想想就知道怎么磕都磕不出来的,今天总算把这道题给补上了,感觉还是不错的

POJ 1061 - 青蛙的约会 - [exgcd求解一元线性同余方程]的更多相关文章

  1. POJ 1061青蛙的约会。求解(x+mT)%L=(y+nT)%L的最小步数T。

    因为是同余,所以就是(x+mT)%L-(y+nT)%L=0.可以写成(x-y+(m-n)T)%L=0.就是这个数是L的倍数啦.那么我可以这样x-y+(m-n)T + Ls = 0.就可以了,s可正可负 ...

  2. NEFU 84 - 五指山 - [exgcd求解一元线性同余方程]

    题目链接:http://acm.nefu.edu.cn/JudgeOnline/problemShow.php?problem_id=84 Time Limit:1000ms Memory Limit ...

  3. poj 1061 青蛙的约会 拓展欧几里得模板

    // poj 1061 青蛙的约会 拓展欧几里得模板 // 注意进行exgcd时,保证a,b是正数,最后的答案如果是负数,要加上一个膜 #include <cstdio> #include ...

  4. ACM: POJ 1061 青蛙的约会 -数论专题-扩展欧几里德

    POJ 1061 青蛙的约会 Time Limit:1000MS     Memory Limit:10000KB     64bit IO Format:%lld & %llu  Descr ...

  5. POJ.1061 青蛙的约会 (拓展欧几里得)

    POJ.1061 青蛙的约会 (拓展欧几里得) 题意分析 我们设两只小青蛙每只都跳了X次,由于他们相遇,可以得出他们同余,则有: 代码总览 #include <iostream> #inc ...

  6. POJ 1061 青蛙的约会(exgcd)

    嗯... 题目链接:http://poj.org/problem?id=1061 两只青蛙相遇时肯定满足:x+k*m≡y+k*n(mod L) x+k*m-(y+k*n)=L*s k*(n-m)-s* ...

  7. poj 1061 青蛙的约会 (扩展欧几里得模板)

    青蛙的约会 Time Limit:1000MS     Memory Limit:10000KB     64bit IO Format:%I64d & %I64u Submit Status ...

  8. poj 1061青蛙的约会

    青蛙的约会 Time Limit: 1000MS   Memory Limit: 10000K Total Submissions: 90083   Accepted: 16257 Descripti ...

  9. POJ 1061 青蛙的约会 题解

    青蛙的约会 Time Limit: 1000MS   Memory Limit: 10000K Total Submissions: 139755   Accepted: 31803 Descript ...

随机推荐

  1. [AX]AX2012 Interaction class

    Ax2012 Client的form如果属性FormTemplate设置为DetailsPage或者ListPage,则必须同时设置属性InteractionClass为相应的Interaction类 ...

  2. ios开发之--调整UISearchBar的输入框的背景颜色

    遍历UISearchBar的子视图,找到输入框坐在的view,添加背景颜色即可. 代码如下: UISearchBar *searchBar = [[UISearchBar alloc] initWit ...

  3. ios开发之--把秒转换为天时分秒

    把秒转换成时分秒: - (NSString *)timeFormatted:(int)totalSeconds { ; ) % ; ; return [NSString stringWithForma ...

  4. Scala单例对象和伴生对象

    1.Scala单例对象 Scala单例对象是十分重要的,没有像在Java一样,有静态类.静态成员.静态方法,但是Scala提供了object对象,这个object对象类似于Java的静态类,它的成员. ...

  5. codeforces水题100道 第三题 Codeforces Beta Round #47 A. Domino piling (math)

    题目链接:http://www.codeforces.com/problemset/problem/50/A题意:一个NxM的举行中最多能放多少个1x2的矩形.C++代码: #include < ...

  6. PHP关于按位取反结果的推导过程

    哎呀几年过去,都快把大学学的计算机导论的知识给忘完了,现在来回顾一下按位去反的流程: <?php /** 首先来补充一下基础知识: php中有4个位运算,分别是&与 |或 ^异或 ~取反 ...

  7. Glide加载图片缓存库出现——You cannot start a load for a destroyed activity

    请记住一句话:不要再非主线程里面使用Glide加载图片,如果真的使用了,请把context参数换成getApplicationContext.

  8. QT之 Hello World

    下载……   我下载的Qt creater 版本为4.2.1,Qt版本为5.8.0 打开QT Creater 1. 新建项目 New Project -> Application -> Q ...

  9. Android设计和开发系列第二篇:Action Bar(Develop—Training)

    Adding the Action Bar GET STARTED DEPENDENCIES AND PREREQUISITES Android 2.1 or higher YOU SHOULD AL ...

  10. Esper学习之十一:EPL语法(七)

    上一篇说到了EPL如何访问关系型数据库这种数据源,实际上别的数据源,比如:webservice.分布式缓存.非关系型数据库等等,Esper提供了统一的数据访问接口.然后今天会讲解如何创建另外一种事件类 ...