二进制GCD算法解析
UPD 2018.3.30 这个好像就是更相减损术的样子emmm
UPD 2018.5.22 好像不是更相减损术而是叫Stein算法的样子emmm
蒟蒻来做个二进制GCD笔记。
为什么要写这个东西呢,因为按照ysy神犇在这次luogu夏令营的说法,常数会小很多。
我再查了一下(ysy神犇没说实现啊orz),这玩意的原理说起来大概是这样的:
因为普通的辗转相除法求gcd需要用到取模,所以常数比较慢。
我们使用另一种算法:
求gcd(a,b)。有三种情况:
1.a,b为偶数,则gcd(a,b)=2*gcd(a/2,b/2)
2.a为奇数,b为偶数,则gcd(a,b)=gcd(a,b/2)
3.a,b为奇数。假设a>=b,则gcd(a,b)=gcd((a-b)/2,b)
4.a为0,则返回b
然后以上的除法和判奇偶性可以使用位运算处理。
取较大值可以使用内联函数降低常数。
然后....emmm
下面是两种算法的执行时间。原题:BZOJ 入门P1234
可能是范围开的不够大?(一脸懵逼)
然后我自己测了一组比较大的数据。
9223372036854775800 529215046068469760
每个程序运行5次,忽略误差较大的运行,取平均值,得到的数据如下:
普通gcd算法:13.576ms
二进制gcd算法:13.044ms
总之...emm,掌握还是很简单的。
至于什么时候用...各位见仁见智吧。
代码如下:
#include<iostream> using namespace std; typedef long long ll; inline ll abs(ll x){return x<0?-x:x; } inline ll min(ll a,ll b){return a<b?a:b; } inline ll gcd(ll a,ll b){if(a==0)return b;if(b==0)return a;if(!(a&1)&&!(b&1))return 2*gcd(a>>1,b>>1);else if(!(a&1))return gcd(a>>1,b);else if(!(b&1))return gcd(a,b>>1);else return gcd(abs(a-b),min(a,b)); }int main(){ll a,b;cin>>a>>b;cout<<"gcd="<<gcd(a,b);return 0; }
转载于:https://www.cnblogs.com/acxblog/p/7275977.html
二进制GCD算法解析相关推荐
- python终结一个循环额_Python语言入门之内存管理方式和垃圾回收算法解析
本文主要向大家介绍了Python语言入门之内存管理方式和垃圾回收算法解析,通过具体的内容向大家展示,希望对大家学习Python语言有所帮助. 在列表,元组,实例,类,字典和函数中存在循环引用问题.有 ...
- AI学习笔记(七)图像滤波器、OpenCV算法解析
AI学习笔记之图像滤波器.OpenCV算法解析 图像滤波器 图像噪声 噪声的产生 信噪比 高斯噪声 椒盐噪声 其他噪声 图像滤波 滤波的目的 滤波的要求 各种滤波器 均值滤波 中值滤波 最大最小值滤波 ...
- 回归算法分类,常用回归算法解析
回归算法分类,常用回归算法解析 回归是数学建模.分类和预测中最古老但功能非常强大的工具之一.回归在工程.物理学.生物学.金融.社会科学等各个领域都有应用,是数据科学家常用的基本工具. 回归通常是机器学 ...
- CVPR目标检测与实例分割算法解析:FCOS(2019),Mask R-CNN(2019),PolarMask(2020)
CVPR目标检测与实例分割算法解析:FCOS(2019),Mask R-CNN(2019),PolarMask(2020) 目标检测:FCOS(CVPR 2019) 目标检测算法FCOS(FCOS: ...
- 10没有基于策略的qos_分布式QoS算法解析
QoS对于服务多租户多业务的整体系统来说,不管对网络还是存储,都格外重要,没有QoS,会造成不同租户及业务之间对资源的抢占,用户A用爽了,用户B却遭了殃,频频投诉,这是系统管理员最头疼的事情.我们今天 ...
- 脑机接口主流算法解析课程视频汇总
目录 讲座1--SSVEP算法解析 讲座2--ERP/P300算法解析 讲座3--运动想象算法解析 讲座4--情感脑机接口算法解析 本分享为脑机学习者Rose整理发表于公众号:脑机接口社区 .QQ交流 ...
- yolo类检测算法解析——yolo v3
原文:https://www.cnblogs.com/cvtoEyes/p/8608205.html yolo类检测算法解析--yolo v3 计算机视觉的发展史可谓很长了,它的分支很多,而且理论那是 ...
- 【机器学习】通俗的元胞自动机算法解析和应用
[机器学习]通俗的元胞自动机算法解析和应用 文章目录 1 元胞自动机的定义 2 元胞自动机的组成 3 元胞自动机的特征 4 Python实现元胞自动机(生命游戏) 5 总结 6 Github(华盛顿州 ...
- 【机器学习】树回归和聚类算法解析和应用
[机器学习]树回归和聚类算法解析和应用 文章目录 1 树回归 2 CART ( Classification And Regression Tree) 分类回归树 3 K-means3.1 合理选择 ...
- 【机器学习】通俗的k-近邻算法算法解析和应用
[机器学习]通俗的k-近邻算法算法解析和应用 文章目录 1 概述 2 KNN 场景 3 KNN 原理 4 实例:改进约会网站的配对效果 5 算法总结 6 KNN算法的优缺点 7 图像分类应用 1 概述 ...
最新文章
- 通俗理解并行性和并发性的区别
- 转:HTTP 协议 GET 、POST之不同
- 两高发布司法解释 依法严惩涉地下钱庄犯罪
- 网络时代:我只是一个用户
- DIV+CSS规范命名集合
- js 加入收藏,设为首页
- 计算机网络组建校园网课设报告,计算机网络课设校园网组建.doc
- winfrom DotNetBar sideNav控件使用问题
- Android自定义控件(状态提示图表,我先收藏为敬
- oracle改了监听 12514,ORA-12514:TNS 监听问题
- r语言把两个折线图图像放到一个图里_OpenCV计算机视觉学习(10)——图像变换(傅里叶变换,高通滤波,低通滤波)...
- iOS 相机调用后为英文的解决办法
- MIMO技术原理、概念、现状简介
- python机械臂仿真_机械臂 python
- 我见过的郭弃疾先生(兰亭集势CEO)
- ORACLE updata是提示违反唯一约束条件
- 北京内推 | 字节跳动AML机器学习系统团队招聘机器学习训练框架研发实习生
- C++ Primer 5 ,众里寻他千百度,下载的大虾别忘了后面买一本正版的支持作者
- xbox360游戏下载_如何将游戏文件移动到另一个Xbox
- C语言回溯法九宫格数独问题
热门文章
- matlab传热模拟计算,MATLAB生物化工计算与模拟
- 阶段5 3.微服务项目【学成在线】_day04 页面静态化_23-页面预览-页面预览开发
- 阶段5 3.微服务项目【学成在线】_day02 CMS前端开发_18-CMS前端页面查询开发-页面原型-创建页面和定义路由...
- 阶段1 语言基础+高级_1-3-Java语言高级_06-File类与IO流_09 序列化流_2_对象的序列化流_ObjectOutputStream...
- 阶段1 语言基础+高级_1-3-Java语言高级_05-异常与多线程_第2节 线程实现方式_14_匿名内部类方式实现线程的创建...
- java自动化学习笔记
- eclipse中xml文件报错异常处理
- cnblog之初来乍到
- 77-CCI,Commodity Channel Index,商品通道指标.(2015.7.1)
- php微型mvc框架创建步骤