c语言二进制法快速取幂模算法,快速幂取模&快速乘取模

快速幂取模

即快速求出(a^b)mod c 的值。由于当a、b的值非常大时直接求a^b可能造成溢出，并且效率低。

思路

原理就是基于\(a*b \% c = ((a \% c)*(b \% c))\% c\),\(a^b \% c = (a \% c)^b \% c\)公式。

求解快速幂：

设指数b用二进制表示为\(b = (b_n b_{n-1}...b_2b_1b_0)_2\)，

\(b = b_0 + b_1*2^1 + b_2*2^2+...+b_{n-1}*2^{n-1} + b_n*2^n\)，

\(a^b = a^{b_0 + b_1*2^1 + b_2*2^2+...+b_{n-1}*2^{n-1} + b_n*2^n} = a^{b0}*a^{b_1*2^1}*a^{b_2*2^2} *...*a^{b_{n-1}*2^{n-1}} * a^{b_n*2^n}\)，

\(a^b \% c = a^{b0}*a^{b_1*2^1}*a^{b_2*2^2} *...*a^{b_{n-1}*2^{n-1}} * a^{b_n*2^n} \% c\)，

设\(K_n = (a^{b_n*2^n})\%c\)，求Kn的话，当bn=0时Kn=1，bn=1时\(Kn=(a^{2^n})\%c\)，因此再考虑计算\((a^{2^n})\%c\)。

\((a^{2^n})\%c = [(a^{2^{n-1}}\%c)*(a^{2^{n-1}}\%c)]\%c\)由此递推。

代码

python

def quick_powmod(a, b, c):

a = a % c

ans = 1 # 存放结果

while b != 0:

if b & 1: # 二进制与

ans = (ans * a) % c

a = (a * a) % c # 取模是防止溢出

b >>= 1 # 二进制向右移动一位

return ans

例如a=2 b=10 c=3，b的二进制表示为1010。

\(2^{10} = 2^{0+ 1*2^1+0*2^2+1*2^3}\)，b的二进制位从右往左取，为0时，累乘a，为1时，将累乘结果乘到ans里。

快速乘取模

使用二进制将乘法转换为加法。

思路

与快速幂取模类似，将一个乘数转换为二进制计算。

例如\(20*14 = 20*(1110)_2 = 20*2^0*0 + 20*2^1*1+20*2^2*1+20*2^3*1\)

代码

Python

def quick_mulmod(a, b, c):

ans = 0

a = a % c

while b != 0:

if b & 1 :

ans = (ans + a) % c

a = (2*a) % c

b >>= 1

return ans

乍一看还是很不好懂的，举个例子推导一遍就明白了。

乘方快速幂 OR 乘法快速幂

关于快速幂这个算法,已经不想多说,很早也就会了这个算法,但是原来一直靠着模板云里雾里的,最近重新学习,发现忽视了一个重要的问题,就是若取模的数大于int型,即若为__int64的时候应该怎么办,这样就 ...

Educational Codeforces Round 13——D. Iterated Linear Function(矩阵快速幂或普通快速幂水题)

D. Iterated Linear Function time limit per test 1 second memory limit per test 256 megabytes input ...

51nod 1113 矩阵快速幂( 矩阵快速幂经典模板 )

1113 矩阵快速幂 链接:传送门 思路:经典矩阵快速幂,模板题,经典矩阵快速幂模板. /******************************************************* ...

【转】C语言快速幂取模算法小结

(转自:http://www.jb51.net/article/54947.htm) 本文实例汇总了C语言实现的快速幂取模算法,是比较常见的算法.分享给大家供大家参考之用.具体如下: 首先,所谓的快速 ...

快速幂取模_C++

一.题目背景 已知底数a,指数b,取模值mo 求ans = ab % mo 二.朴素算法(已知可跳过) ans = 1,循环从 i 到 b ,每次将 ans = ans * a % mo 时间复杂度O ...

NYOJ-676小明的求助，快速幂求模，快速幂核心代码；

小明的求助 时间限制:2000 ms  |  内存限制:65535 KB 难度:2 描述 小明对数学很有兴趣,今天老师出了道作业题,让他求整数N的后M位,他瞬间感觉老师在作弄他,因为这是so easy ...

hdu4549 矩阵快速幂 + 欧拉降幂

R - M斐波那契数列 Time Limit:1000MS     Memory Limit:32768KB     64bit IO Format:%I64d & %I64u Submit  ...

矩阵快速幂——POJ3070

矩阵快速幂和普通的快速幂差不多,只不过写起来比较麻烦一点,需要重载*运算符. 模板: struct mat { int m[maxn][maxn]; }unit; mat operator * (ma ...

[学习笔记]快速幂&&快速乘

本质:二进制拆分(你说倍增我也没脾气).然后是一个配凑. 合起来就是边二进制拆分,边配凑. 快速乘(其实龟速):把乘数二进制拆分.利用乘法分配率. 用途:防止爆long long 代码: ll qk( ...

快速幂C++实现

快速幂模板题 很明显,这个题目不能用简单的\(for\)循环+\(mod\)来完成,因为指数\(p\)已经达到了长整型,直接循环来完成的话肯定会超时的. 那么快速幂就应运而生了. 什么是快速幂呢? 利 ...

随机推荐

leetcode日记 Product of Array Except Self

Given an array of n integers where n > 1, nums, return an array output such that output[i] is equ ...

javascript必知必会：面象对象编程

面象对象编程技术的核心理念:封装.继承.多态:在一些主流的高级编程语言中,比如:C#,VB.NET,JAVA,PHP等都是很容易实现的,而如果要在javascript中实现面象对象编程,可就不那么直接 ...

谈谈我的编程之路---WAMP(一)

WAMP的一些配置与使用心得(PHP) 记得第一次接触PHP的时候,我都不知道PHP为什么要大写,但是我却用它来进行工作了,有时候生活就是一场美丽的邂逅 青涩的我,在ES哥的引领下,第一次接触到了WA ...

Linux的僵尸进程产生原因及解决方法

Linux的僵尸进程产生原因及解决方法: 1. 产生原因: 在UNIX 系统中,一个进程结束了,但是他的父进程没有等待(调用wait / waitpid)他,那么他将变成一个僵尸进程.通过ps命令查看 ...

JAVA-3-水仙花

public static void main(String[] args) { // TODO 自动生成的方法存根 int i = 100; while (i < 1000) { int a, ...

让MyEclipse里的Tomcat自动reloadable

1  修改server.xml Context path="/***" docBase="XXX" reloadable="true"/&g ...

OpenGL&period;Vertex Array Object &lpar;VAO&rpar;&period;

OpenGL抛弃glEnable(),glColor(),glVertex(),glEnable()这一套流程的函数和管线以后,就需要一种新的方法来传递数据到Graphics Card来渲染几何体,我 ...

date、sleep和usleep命令

bash&shell系列文章:http://www.cnblogs.com/f-ck-need-u/p/7048359.html 1.1 date命令 date用于获取和设置操作系统的时间,还 ...

深入springboot原理——动手封装一个starter

从上一篇文章 我们已经知道springboot的起步依赖与自动配置的机制.spring-bo ...

git的安装以及入门

安装:https://blog.csdn.net/itpinpai/article/details/48105445 (1)下载文件 初始化 git init 连远程服务器 git remote ad ...