Stein算法是一种计算两个数最大公约数的算法,是针对欧几里德算法在对大整数进行运算时,需要试商导致增加运算时间的缺陷而提出的改进算法. 算法思想: 由J. Stein 1961年提出的Stein算法很好的解决了欧几里德算法中的这个缺陷,Stein算法只有整数的移位和加减法,为了说明Stein算法的正确性,首先必须注意到以下结论: gcd(a,a)=a,也就是一个数和其自身的公约数仍是其自身. gcd(ka,kb)=k gcd(a,b),也就是最大公约数运算和倍乘运算可以交换.特殊地,当k=2时…

辗转相减求a,b的gcd其实可以优化的: 1.若a为偶数,b为奇数:gcd(a,b)=gcd(a/2,b) 2.若a为奇数,b为偶数:gcd(a,b)=gcd(a,b/2) 3.若a,b都是偶数:gcd(a,b)=2*gcd(a/2,b/2) 3.若a,b都是奇数:gcd(a,b)=gcd(a-b,b) (a>b) 然后就涉及到高精度乘单精度,高精度减高精度,高精度除单精度………

传送门 非常显du然liu的一道题 就是求GCD 因为数据范围... 所以要上压位高精+非递归的辗转相减 关于辗转相减: 如果 A是二的倍数,B是二的倍数   那么GCD(A,B)=2 * GCD(A,B) 如果只有A是二的倍数   那么GCD(A,B)=GCD(A/2,B) 如果只有B是二的倍数   那么GCD(A,B)=GCD(A,B/2) 十分显然的结论... 然后不停地让大的数减去小的数 最后当它们相等时就是GCD了(因为大的减小的一直减到不能减就相当于取模) int slove() {…

这道题直接写了我两个多小时-- 主要是写高精度的时候还存在着一些小毛病,调了很久 在输入这一块卡了很久. 然后注意这里用while的形式写,不然会炸 最后即使我已经是用的万进制了,但是交上去还是有两个点超时 然后就开始漫长的改进,一直过不了那两个点. 然后突然发现,貌似这道题没有涉及到乘法. 那不就可以直接开10的九次方为一位了(10的9次方是int的最大数量级) 我交上去之后就AC了,全部测试点交上去总和4.6秒 然后原来我一个数开的是2500这么大,因为题目给的10000位,之前用万进制除以…

题意简述 求两个整数a,b的最大公约数0 < a , b ≤ 10 ^ 10000. 题解思路 如果 a % 2 == 0 && b % 2 == 0 gcd(a,b) = gcd(a / 2, b / 2) * 2 如果 a % 2 == 0 && b % 2 != 0 gcd(a,b) = gcd(a / 2, b); 如果 a % 2 != 0 && b % 2 == 0 gcd(a,b) = gcd(a, b / 2); 如果 a % 2 !=…

1876: [SDOI2009]SuperGCD Time Limit: 4 Sec  Memory Limit: 64 MBSubmit: 3060  Solved: 1036[Submit][Status][Discuss] Description Sheng bill有着惊人的心算能力,甚至能用大脑计算出两个巨大的数的GCD(最大公约 数)!因此他经常和别人比赛计算GCD.有一天Sheng bill很嚣张地找到了你,并要求和你比 赛,但是输给Sheng bill岂不是很丢脸!所以你决定写一…

1876: [SDOI2009]SuperGCD Time Limit: 4 Sec  Memory Limit: 64 MBSubmit: 2384  Solved: 806[Submit][Status][Discuss] Description Sheng bill有着惊人的心算能力,甚至能用大脑计算出两个巨大的数的GCD(最大公约 数)!因此他经常和别人比赛计算GCD.有一天Sheng bill很嚣张地找到了你,并要求和你比 赛,但是输给Sheng bill岂不是很丢脸!所以你决定写一个…

更相减损,要用高精度.... --------------------------------------------------------------- #include<cstdio> #include<cstring> #include<cctype> #include<algorithm>   using namespace std;   const int maxn = 10009;   char S[maxn]; int Power[maxn]…

[BZOJ1876][SDOI2009]SuperGCD(数论,高精度) 题面 BZOJ 洛谷 题解 那些说数论只会\(gcd\)的人呢?我现在连\(gcd\)都不会,谁来教教我啊? 显然\(gcd\)除了辗转相除之外还可以辗转相减,然而辗转相减对于这题而言显然还不够优秀. 我们这样子来做. 如果当前\(a,b\)都是\(2\)的倍数,那么我们就把\(2\)直接同时除掉,直接在\(gcd\)中乘上一个\(2\).否则如果只有一个数是\(2\)的倍数,显然可以直接把这个\(2\)给除掉. 这样子可…

特别说明: 为了避免以后搬家时的麻烦, 这里的文章继续沿用csdn的风格和分类好了~ Emmmm这个题是一道高精度的模板题啊~ 既然是高精度的裸题, 那我们这些懒人当然是选择:用python啦~ 懒癌晚期 你看这不就做完了么←_← a=input() b=input() while(b): c=b b=a%b a=c print a 当然这份代码并不能在luogu上AC 应该是数据出锅了. (当然也不能算是出锅, 只是不太符合题目中说的输入格式而已...刻意卡python (当然还是有一些pyt…

题目:http://www.lydsy.com:808/JudgeOnline/problem.php?id=1876 分析: 以为辗转相减会TLE呢……但是好像没这个数据……就这么水过去了…… 辗转相减求a,b的gcd其实可以优化的: 1.若a为偶数,b为奇数:gcd(a,b)=gcd(a/2,b) 2.若a为奇数,b为偶数:gcd(a,b)=gcd(a,b/2) 3.若a,b都是偶数:gcd(a,b)=2*gcd(a/2,b/2) 3.若a,b都是奇数:gcd(a,b)=gcd(a-b,b)…

更相减损数. 上手就debug了3个小时,直接给我看哭了. 3个函数都写错了是什么感受? 乘2函数要从前往后乘,这样后面的数乘2进位以后不会干扰前面的数. 除2函数要从后往前除,这样前面的数借来的位不会除2次. 然后函数里面俩个类里面的变量名不要打混. 哭瞎 #include<cstdio> #include<algorithm> #include<cstring> using namespace std; + ; ,,,,,,,}; ]; int cnt; struc…

http://www.lydsy.com/JudgeOnline/problem.php?id=1876 高精压位GCD 对于  GCD(a, b)  a>b 若 a 为奇数,b 为偶数,GCD(a, b) = GCD(a, b / 2) 若 a 为偶数,b 为奇数,GCD(a, b) = GCD(a / 2, b) 若 a 为偶数,b 为偶数,GCD(a, b) = 2*GCD(a / 2, b / 2) 若 a 为奇数,b 为奇数,GCD(a, b) = GCD(a - b, b) vect…

题目链接 这题.高精度.恶心.难受. 那么高精度的gcd怎么做呢? 若a=b gcd(a,b)=a ①a偶b偶 gcd(a,b)=2*gcd(a/2,b/2) ②a偶b奇 gcd(a,b)=gcd(a/2,b) ③a奇b奇 gcd(a,b)=gcd(a-b,b) 嗯.这玩意就这样了. #include<cmath> #include<iostream> #include<cstdio> #include<algorithm> #include<cstr…

题目大意 求两个个高精度数的gcd 题目解析 在学习gcd的时候,书上就记载了"更相减损术"这一方法 基于这种方法,我们进行优化,使得我们能快速求出两个大数的gcd 对于 \(a,b\) 的 \(GCD(a, b)\) 有 [1]. 若 \(a\) 为奇数,\(b\) 为偶数,\(GCD(a, b) = GCD(a, b / 2)\) 表示 \(b\) 存在2这个因子而 \(a\) 不存在,则将 \(b\) 除以2,,不考虑因子2: [2]. 若 \(a\) 为偶数,\(b\) 为奇数…

http://www.lydsy.com/JudgeOnline/problem.php?id=1876 Description Sheng bill有着惊人的心算能力,甚至能用大脑计算出两个巨大的数的GCD(最大公约 数)!因此他经常和别人比 赛计算GCD.有一天Sheng bill很嚣张地找到了你,并要求和你比 赛,但是输给Sheng bill岂不是很丢脸!所以你 决定写一个程序来教训他. Input 共两行: 第一行:一个数A. 第二行:一个数B. 0 < A , B ≤ 10 ^ 100…

题目 Sheng bill有着惊人的心算能力,甚至能用大脑计算出两个巨大的数的GCD(最大公约 数)!因此他经常和别人比 赛计算GCD.有一天Sheng bill很嚣张地找到了你,并要求和你比 赛,但是输给Sheng bill岂不是很丢脸!所以你 决定写一个程序来教训他. 输入格式 共两行: 第一行:一个数A. 第二行:一个数B. 0 < A , B ≤ 10 ^ 10000. 输出格式 一行,表示A和B的最大公约数. 输入样例 12 54 输出样例 6 题解 时隔大半年,我回来A这道题啦[当初…

题目传送门:https://www.lydsy.com/JudgeOnline/problem.php?id=1876 一道简单的高精度+Stein算法(或者叫辗转相除法)求最大公约数板子题. md还要压位.. 代码: #include<cstdio> #include<cmath> #include<cstdlib> #include<cstring> #include<ctime> #include<iostream> #incl…

题面就是让你求两个超级大整数,求GCD 题解: 题目本意应该是出题人想考考高精度取膜 但是可以通过一种神奇的Stein算法来做 由J. Stein 1961年提出的Stein算法很好的解决了欧几里德算法中的这个缺陷,Stein算法只有整数的移位和加减法,为了说明Stein算法的正确性,首先必须注意到以下结论: gcd(a,a)=a,也就是一个数和其自身的公约数仍是其自身. gcd(ka,kb)=k gcd(a,b),也就是最大公约数运算和倍乘运算可以交换.特殊地,当k=2时,说明两个偶数的最大公…

点此看题面 大致题意: 给你两个长度\(\le10000\)的正整数,让你求它们的\(gcd\). Python​ 高精请绕道. 这题的正解应该是Python. 对于这种高精题,肯定是Python最方便了. 于是我就默默写了Python. 代码 n=(int)(input())#读入第一个整数 m=(int)(input())#读入第二个整数 while n!=m:#辗转相减法(不知道为什么辗转相除法RE了) if n>m:#如果n>m n-=m#将n减去m else:#否则 m-=n#将m减…

要你求两个非常大的数字的GCD. 不要想复杂,用高精度整更相减损术即可. #include <cstdio> #include <cstring> #include <algorithm> using namespace std; struct BigInt { static const int BASE = 10000, CARRY = 4, MAX_N = 10000; int A[MAX_N], Len; void Clear() { memset(A, 0, s…

Description Sheng bill有着惊人的心算能力,甚至能用大脑计算出两个巨大的数的GCD(最大公约 数)!因此他经常和别人比赛计算GCD.有一天Sheng bill很嚣张地找到了你,并要求和你比 赛,但是输给Sheng bill岂不是很丢脸!所以你决定写一个程序来教训他. Input 共两行: 第一行:一个数A. 第二行:一个数B. Output 一行,表示A和B的最大公约数. Sample 输入样例#1: 12 54 输出样例#1: 6 Hints 对于20%的数据,0 < A…

高精度+GCD 唔……高精gcd其实可以这么算: \[ GCD(a,b)= \begin{cases} a & b=0 \\ 2*GCD(\frac{a}{2},\frac{b}{2}) &  a\mod 2=0,b \mod 2=0 \\ GCD(\frac{a}{2},b) &  a\mod 2=0,b \mod 2=1 \\ GCD(a,\frac{b}{2}) &  a\mod 2=1,b \mod 2=0 \\ GCD(b,a-b) & else \end…

1226: [SDOI2009]学校食堂Dining Description 小F 的学校在城市的一个偏僻角落,所有学生都只好在学校吃饭.学校有一个食堂,虽然简陋,但食堂大厨总能做出让同学们满意的菜肴.当然,不同的人口味也不一定相同,但每个人的口味都可以用一个非负整数表示.由于人手不够,食堂每次只能为一个人做菜.做每道菜所需的时间是和前一道菜有关的,若前一道菜的对应的口味是a,这一道为b,则做这道菜所需的时间为(a or b)-(a and b),而做第一道菜是不需要计算时间的.其中,or 和a…

第一次感觉Python艹题的快感 1876: [SDOI2009]SuperGCD Time Limit: 4 Sec Memory Limit: 64 MB Submit: 2461 Solved: 831 [Submit][Status][Discuss] Description Sheng bill有着惊人的心算能力,甚至能用大脑计算出两个巨大的数的GCD(最大公约 数)!因此他经常和别人比赛计算GCD.有一天Sheng bill很嚣张地找到了你,并要求和你比 赛,但是输给Sheng bi…

不知不觉AC100题了,放眼望去好像都是水题.在这里就做一个存档吧(特别感谢各位大神尤其是云神http://hi.baidu.com/greencloud和丽洁姐http://wjmzbmr.com/的blog提供的题解) 代码什么的有时间在放. 1000: A+B Problem 好吧,这也算题,记得我第一次直接copy然后他不给我加库然后就CE了= = 1001: [BeiJing2006]狼抓兔子 一道经典的最小割对偶图转spfa,参见白书,第一次打时按边建图然后TLE得很惨,然后又重打了…

一.先是一些整除的性质: •整除:若a=bk,其中a,b,k都是整数,则b整除a,记做b|a. •也称b是a的约数(因数),a是b的倍数 •显而易见的性质: •1整除任何数,任何数都整除0 •若a|b,a|c,则a|b+c, a|b-c •若a|b,则对任意整数c,a|bc •传递性:若a|b,b|c,则a|c 例1: •例题:[CF 762A]k-th divisor •求n的第k小的约数.如果不存在输出-1 •1 ≤ n ≤ 10^15, 1 ≤ k ≤ 10^9 分析:这道题显然不能用O(…

感觉今天好颓啊,我才把昨晚那五道题a了,(但我明明一直在学啊,为啥这么慢,难道是我太笨了?) 闲话少叙,先说做法 问题 A: C Looooops 题目描述 对于C的for(i=A ; i!=B ;i +=C)循环语句,给出A,B,C和k(k表示变量是在k进制下的无符号整数),判断循环次数,不能终止输出"FOREVER". 输入 多组数据,每组一行,A,B,C,k 程序以0 0 0 0结束 输出 一行一个整数,表示循环次数,或者"FOREVER" 样例输入 3 3…

bzoj1000 A+B Problem Scanner sc=new Scanner(new BufferedInputStream(System.in)); 声明读入器 nextInt 读入整数 System.out.println 带换行输出 import java.util.*; import java.io.*; public class Main{ public static void main(String[] argc){ Scanner sc = new Scanner (ne…

Search GO 说明:输入题号直接进入相应题目,如需搜索含数字的题目,请在关键词前加单引号 Problem ID Title Source AC Submit Y 1000 A+B Problem 10983 18765 Y 1036 [ZJOI2008]树的统计Count 5293 13132 Y 1588 [HNOI2002]营业额统计 5056 13607 1001 [BeiJing2006]狼抓兔子 4526 18386 Y 2002 [Hnoi2010]Bounce 弹飞绵羊 43…