快速幂取模&快速乘取模

快速幂取模

即快速求出（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$已经达到了长整型,直接循环来完成的话肯定会超时的. 那么快速幂就应运而生了. 什么是快速幂呢? 利 ...

随机推荐

Mysql5.5源码安装步骤笔记记录
1.cmake软件的安装wget https://cmake.org/files/v3.5/cmake-3.5.0-rc3.tar.gztar xf cmake-3.5.0.tar.gzcd cmak ...
XAMARIN +VS2015 ANDROID 开发禁止屏幕自动转动 Portrait,Nosensor
网上有很多java的写法,但是放在C#中都无法使用,其实有时候还是安下心来认真的去看官网文档比在百度或者google来的快的多 this.RequestedOrientation = Android. ...
9.nodejs权威指南--Socket.IO
1. Socket.IO 1.1 服务器 var http = require('http'); var sio = require('socket.io'); var fs = require('f ...
HTML标记语法之表格元素
语法与语义: <table>和</table>定义表格的开始和结束 <thead>和</thead>定义表格头部的开始和结束 <tbody> ...
Xml Schema：C#访问在complextype中插入新元素
最近用c#写Xml Schema代码,找了很久也找不到如何在已有的complextype中插入新的element,最后我充分发挥自己的聪明才智,哈哈,终于从...中找到了灵感. XmlSchemaSe ...
App接口中xml方式封装通信接口
eBay_Relist（退换刊登费）
如果物品首次刊登结束时没有人中标,或是买家并没有付款完成交易,那么卖家以通过重新刊登的方法来再次销售.如果该物品在第二次刊登时成功售出,且满足eBay的重新刊登退费条件,那么eBay便会退还重新刊登的 ...
Android写入文件操作权限
<uses-permission android:name="android.permission.WRITE_EXTERNAL_STORAGE"></uses- ...
解决ionic在Android和iOS的一些样式上的冲突
//设置默认返回按钮的文字 $ionicConfigProvider.backButton.previousTitleText(false).text('返回'); // 设置全局 $http 超时 ...