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算法解析相关推荐

  1. python终结一个循环额_Python语言入门之内存管理方式和垃圾回收算法解析

    本文主要向大家介绍了Python语言入门之内存管理方式和垃圾回收算法解析,通过具体的内容向大家展示,希望对大家学习Python语言有所帮助. 在列表,元组,实例,类,字典和函数中存在循环引用问题.有 ...

  2. AI学习笔记(七)图像滤波器、OpenCV算法解析

    AI学习笔记之图像滤波器.OpenCV算法解析 图像滤波器 图像噪声 噪声的产生 信噪比 高斯噪声 椒盐噪声 其他噪声 图像滤波 滤波的目的 滤波的要求 各种滤波器 均值滤波 中值滤波 最大最小值滤波 ...

  3. 回归算法分类,常用回归算法解析

    回归算法分类,常用回归算法解析 回归是数学建模.分类和预测中最古老但功能非常强大的工具之一.回归在工程.物理学.生物学.金融.社会科学等各个领域都有应用,是数据科学家常用的基本工具. 回归通常是机器学 ...

  4. CVPR目标检测与实例分割算法解析:FCOS(2019),Mask R-CNN(2019),PolarMask(2020)

    CVPR目标检测与实例分割算法解析:FCOS(2019),Mask R-CNN(2019),PolarMask(2020) 目标检测:FCOS(CVPR 2019) 目标检测算法FCOS(FCOS: ...

  5. 10没有基于策略的qos_分布式QoS算法解析

    QoS对于服务多租户多业务的整体系统来说,不管对网络还是存储,都格外重要,没有QoS,会造成不同租户及业务之间对资源的抢占,用户A用爽了,用户B却遭了殃,频频投诉,这是系统管理员最头疼的事情.我们今天 ...

  6. 脑机接口主流算法解析课程视频汇总

    目录 讲座1--SSVEP算法解析 讲座2--ERP/P300算法解析 讲座3--运动想象算法解析 讲座4--情感脑机接口算法解析 本分享为脑机学习者Rose整理发表于公众号:脑机接口社区 .QQ交流 ...

  7. yolo类检测算法解析——yolo v3

    原文:https://www.cnblogs.com/cvtoEyes/p/8608205.html yolo类检测算法解析--yolo v3 计算机视觉的发展史可谓很长了,它的分支很多,而且理论那是 ...

  8. 【机器学习】通俗的元胞自动机算法解析和应用

    [机器学习]通俗的元胞自动机算法解析和应用 文章目录 1 元胞自动机的定义 2 元胞自动机的组成 3 元胞自动机的特征 4 Python实现元胞自动机(生命游戏) 5 总结 6 Github(华盛顿州 ...

  9. 【机器学习】树回归和聚类算法解析和应用

    [机器学习]树回归和聚类算法解析和应用 文章目录 1 树回归 2 CART ( Classification And Regression Tree) 分类回归树 3 K-means3.1 合理选择 ...

  10. 【机器学习】通俗的k-近邻算法算法解析和应用

    [机器学习]通俗的k-近邻算法算法解析和应用 文章目录 1 概述 2 KNN 场景 3 KNN 原理 4 实例:改进约会网站的配对效果 5 算法总结 6 KNN算法的优缺点 7 图像分类应用 1 概述 ...

最新文章

  1. 通俗理解并行性和并发性的区别
  2. 转:HTTP 协议 GET 、POST之不同
  3. 两高发布司法解释 依法严惩涉地下钱庄犯罪
  4. 网络时代:我只是一个用户
  5. DIV+CSS规范命名集合
  6. js 加入收藏,设为首页
  7. 计算机网络组建校园网课设报告,计算机网络课设校园网组建.doc
  8. winfrom DotNetBar sideNav控件使用问题
  9. Android自定义控件(状态提示图表,我先收藏为敬
  10. oracle改了监听 12514,ORA-12514:TNS 监听问题
  11. r语言把两个折线图图像放到一个图里_OpenCV计算机视觉学习(10)——图像变换(傅里叶变换,高通滤波,低通滤波)...
  12. iOS 相机调用后为英文的解决办法
  13. MIMO技术原理、概念、现状简介
  14. python机械臂仿真_机械臂 python
  15. 我见过的郭弃疾先生(兰亭集势CEO)
  16. ORACLE updata是提示违反唯一约束条件
  17. 北京内推 | 字节跳动AML机器学习系统团队招聘机器学习训练框架研发实习生
  18. C++ Primer 5 ,众里寻他千百度,下载的大虾别忘了后面买一本正版的支持作者
  19. xbox360游戏下载_如何将游戏文件移动到另一个Xbox
  20. C语言回溯法九宫格数独问题

热门文章

  1. matlab传热模拟计算,MATLAB生物化工计算与模拟
  2. 阶段5 3.微服务项目【学成在线】_day04 页面静态化_23-页面预览-页面预览开发
  3. 阶段5 3.微服务项目【学成在线】_day02 CMS前端开发_18-CMS前端页面查询开发-页面原型-创建页面和定义路由...
  4. 阶段1 语言基础+高级_1-3-Java语言高级_06-File类与IO流_09 序列化流_2_对象的序列化流_ObjectOutputStream...
  5. 阶段1 语言基础+高级_1-3-Java语言高级_05-异常与多线程_第2节 线程实现方式_14_匿名内部类方式实现线程的创建...
  6. java自动化学习笔记
  7. eclipse中xml文件报错异常处理
  8. cnblog之初来乍到
  9. 77-CCI,Commodity Channel Index,商品通道指标.(2015.7.1)
  10. php微型mvc框架创建步骤