POJ 1061 - 青蛙的约会 - [exgcd求解一元线性同余方程]
先上干货:
定理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求解一元线性同余方程]的更多相关文章
- 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可正可负 ...
- NEFU 84 - 五指山 - [exgcd求解一元线性同余方程]
题目链接:http://acm.nefu.edu.cn/JudgeOnline/problemShow.php?problem_id=84 Time Limit:1000ms Memory Limit ...
- poj 1061 青蛙的约会 拓展欧几里得模板
// poj 1061 青蛙的约会 拓展欧几里得模板 // 注意进行exgcd时,保证a,b是正数,最后的答案如果是负数,要加上一个膜 #include <cstdio> #include ...
- ACM: POJ 1061 青蛙的约会 -数论专题-扩展欧几里德
POJ 1061 青蛙的约会 Time Limit:1000MS Memory Limit:10000KB 64bit IO Format:%lld & %llu Descr ...
- POJ.1061 青蛙的约会 (拓展欧几里得)
POJ.1061 青蛙的约会 (拓展欧几里得) 题意分析 我们设两只小青蛙每只都跳了X次,由于他们相遇,可以得出他们同余,则有: 代码总览 #include <iostream> #inc ...
- 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* ...
- poj 1061 青蛙的约会 (扩展欧几里得模板)
青蛙的约会 Time Limit:1000MS Memory Limit:10000KB 64bit IO Format:%I64d & %I64u Submit Status ...
- poj 1061青蛙的约会
青蛙的约会 Time Limit: 1000MS Memory Limit: 10000K Total Submissions: 90083 Accepted: 16257 Descripti ...
- POJ 1061 青蛙的约会 题解
青蛙的约会 Time Limit: 1000MS Memory Limit: 10000K Total Submissions: 139755 Accepted: 31803 Descript ...
随机推荐
- NetBpm 目录
整理了一下网上的一些netbpm,虽然这项目现在不再更新了,还是想去学习一下,这个组件用时候很方便,比workFlow方便的多了 如果像jbpm那样一直更新就好了. 前两篇是个人的一个总结,后面一些是 ...
- 自定义tag标签-实现long类型转换成Date类型
数据库里存储的是bigint型的时间,entity实体中存放的是long类型的标签,现在想输出到jsp页面,由于使用的是jstl标签,而要显示的是可读的时间类型,找来找去有个 fmt:formatDa ...
- ios开发之--字符串局部改变颜色
改变指定位置字符的颜色,代码如下: NSString *descStr = @"楼主"; NSString *nickStr = [NSString stringWithForma ...
- D盾 v2.0.6.42 测试记录
0x01 前言 之前发了一篇博客<Bypass D盾_IIS防火墙SQL注入防御(多姿势)>,D哥第一时间联系我,对问题进行修复.这段时间与D哥聊了挺多关于D盾这款产品的话题,实在是很佩服 ...
- 利用DB实现简单去重处理
业务需要针对某文件进行判重操作,用Perl实现如下 #!/usr/bin/perl my %h; s/\s+$// and $h{$_}++ while <>; print "$ ...
- MySQL 安装与配置
Linux 安装 MySQL Windows 安装 MySQL 如何连接 MySQL 如何修改 MySQL 密码 如何重置 MySQL 密码
- Repeater的j简单使用
嘿嘿,今天没有任务,所以突然想起来我之前记得笔 记说要把repeater的使用以及获取值的详细使用总结一下,所以这就闲来无聊总结一下,虽然现在不会使用这些小知识点的,但是我感觉自己的学习还是要 不断地 ...
- 使用react进行父子组件传值
在单页面里面,父子组件传值是比较常见的,之前一直用vue开发,今天研究了一下react的父子组件传值,和vue差不多的思路,父组件向子组件传值,父通过初始state,子组件通过this.props进行 ...
- linux多行注释
1.多行注释: 1. 首先按esc进入命令行模式下,按下Ctrl + v,进入列(也叫区块)模式; 2. 在行首使用上下键选择需要注释的多行; 3. 按下键盘(大写)“I”键,进入插入模式 ...
- JSON语法介绍
官网:https://www.json.org/ JSON (JavaScript Object Notation) is a lightweight data-interchange f ...