智障了,智障了,水一水博客. 本来是个水题,但是for循环遍历那里写挫了... One Person Game Time Limit: 2 Seconds      Memory Limit: 65536 KB There is an interesting and simple one person game. Suppose there is a number axis under your feet. You are at point A at first and your aim is…

题意:一个人在坐标A,要前往坐标B的位置.可以往左或往右走a,b,a+b个单位,求到达B的最小步数. 分析:扩展欧几里得算法求解线性方程的套路不变.令C=fabs(A-B),c = a+b, 扩展gcd分别求 ax+by=C ; ax+cy = C : bx+cy = C的最小|x|+|y|.求min{|x|+|y|}需要一点思考. 对于线性方程ax+by=c,设d = gcd(a,b) ,若方程有解,则必须d | c,特解为 (x0,y0) = ( xx*c/d,yy*c/d) .设am =…

前言 扩展欧几里得算法是一个很好的解决同余问题的算法,非常实用. 欧几里得算法 简介 欧几里得算法,又称辗转相除法. 主要用途 求最大公因数\(gcd\). 公式 \(gcd(a,b)=gcd(b,a\%b)\) 公式证明 \(a\)可以表示成\(a=kb+a\%b\)(\(k\)为自然数). 假设\(g\)是\(a,b\)的一个公约数,则有\(g|a, g|b\). \(\because a\%b=a-kb\), \(\therefore g|(a\%b),\therefore g\)是\(b…

关于扩展欧几里得从寒假时就很迷,抄题解过了同余方程,但是原理并不理解. 今天终于把坑填上了qwq. 由于本人太菜,不会用markdown,所以这篇总结是手写的(什么).(字丑不要嫌弃嘛) ********Update9.28********** 刚刚我们求出的是一组特值,那么如何求通值? 约定:设x0,y0为一组特解,t为任意整数,设a>b(不行再交换) 那么有 x=x0+b/gcd*t y=y0-a/gcd*t ******************************* 奉上三道例题: E…

exgcd入门以及同余基础 gcd,欧几里得的智慧结晶,信息竞赛的重要算法,数论的...(编不下去了 讲exgcd之前,我们先普及一下同余的性质: 若,那么 若,,且p1,p2互质, 有了这三个式子,就不用怕在计算时溢出了. 下面我会用与分别表示a与b的最大公约数与最小公倍数. 首先会来学扩欧的同学肯定都会欧几里得算法(即辗转相除法)了吧 而通过观察发现:,先除后乘防溢出. 所以与的代码如下: inline int gcd(int a,int b) {)?a:gcd(b,a%b);} inlin…

One Person Game Time Limit: 2 Seconds      Memory Limit: 65536 KB There is an interesting and simple one person game. Suppose there is a number axis under your feet. You are at point A at first and your aim is point B. There are 6 kinds of operations…

题意:求乘法逆元最小正正数解 思路:a*x≡1(mod m),则称x 是 a 关于 m 的乘法逆元,可以通过解a*x + m*y = 1解得x.那么通过EXGcd得到特解x1,最小正解x1 = x1 % m,如果x1 <=0,x1 += m,注意m是负数时取绝对值,因为是正解,所以不能用(x1%m+m)%m. 代码: #include<iostream> #include<cstdio> #include<cstring> #include<cmath>…

贝祖定理:即如果a.b是整数,那么一定存在整数x.y使得ax+by=gcd(a,b).换句话说,如果ax+by=m有解,那么m一定是gcd(a,b)的若干倍.(可以来判断一个这样的式子有没有解)有一个直接的应用就是 如果ax+by=1有解,那么gcd(a,b)=1: int gcd(int a,int b){return b==0?a:gcd(b,a%b);} 然而这并不能告诉我们x,y解是多少. 扩欧 首先我们观察上面的式子发现一定有一个解a*1+b*0=gcd(a,b).(b%a=0) 但是…

Elementary Number Theory - Extended Euclid Algorithm Time Limit : 1 sec, Memory Limit : 65536 KB Japanese version is here Extended Euclid Algorithm Given positive integers a and b, find the integer solution (x, y) to ax+by=gcd(a,b), where gcd(a,b) is…

算法思想 我们想求得一组\(x,y\)使得 \(ax+by = \gcd(a,b)\) 根据 \(\gcd(a,b) = \gcd(b,a\bmod b)\) 如果我们现在有\(x',y'\) 使得 \(bx'+(a\bmod b)y' = \gcd(b,a\bmod b)\) 那么 \(ax+by = bx'+( a-\lfloor\frac a b\rfloor b)y'\) 移项之后 \(ax+by = ay'+b(x'-\lfloor\frac a b\rfloor y')\) 我们可以…