快速幂-竞赛算法用法-和-例题详解---Hayden's Blog.

xiaoxiao2021-02-28  83

在ACM竞赛中经常会遇到各种的求幂问题,这种问题一般情况下对应的数字都会很大,所以快速幂就会经常用到,快速幂很好用,就算没有弄清楚快速幂的原理,只要有模板很快就可以解题。

先说几个比较重要的公式: 模运算与基本四则运算有些相似,但是除法例外。其规则如下:

(a + b) % p = (a % p + b % p) % p (1) // 这个公式等下写题目就会用到的 (a – b) % p = (a % p – b % p) % p (2) (a * b) % p = (a % p * b % p) % p (3) (a^b) % p = ((a % p)^b) % p (4) //快速幂的公式

结合律:

((a+b) % p + c) % p = (a + (b+c) % p) % p (5) ((a*b) % p * c)% p = (a * (b*c) % p) % p (6)

交换律:

(a + b) % p = (b+a) % p (7) (a * b) % p = (b * a) % p (8)

分配律:

((a +b)% p * c) % p = ((a * c) % p + (b * c) % p) % p (9)

重要定理:

若a≡b (% p),则对于任意的c,都有(a + c) ≡ (b + c) (%p);(10) 若a≡b (% p),则对于任意的c,都有(a * c) ≡ (b * c) (%p);(11) 若a≡b (% p),c≡d (% p),则 (a + c) ≡ (b + d) (%p),(a – c) ≡ (b – d) (%p), (a * c) ≡ (b * d) (%p),(a / c) ≡ (b / d) (%p); (12)

先上几个最常用的模板: 1.快速乘法 : a*b

int qmul(int a,int b){// 根据数据范围可选择long long int ans=0; while(b){ if( b&1)ans+=a;//按位与完成位数为1的判断 b>>=1;a<<=1;//位运算代替/2和*2 } return ans; }

2.快速乘法取模 (a*b)%mod

int qmul_mod(int a,int b,int mod){ int ans=0; while(b){ if((b%=mod)&1)ans+=a%=mod;//这里需要b%=mod 以及a%=mod b>>=1;a<<=1; } return ans%mod; //ans也需要对mod取模 }

3.快速幂 a^b

int qpow(int a,int b){ if(a==0)return 0;//这是个坑,校赛被坑过,很多网上的实现都没写这一点 int ans=1; while(b){ if(b&1)ans*=a;//和快速乘法的区别 b>>=1;a*=a;//区别,同上 } return ans; }

4.快速幂取模 (a^b)%mod

int qpow_mod(int a,int b,int mod){ if(a==0)return 0; int ans=1; while(b){ if(b&1)ans=(ans%mod)*(a%mod);//如果确定数据不会爆的话,可写成 ans*=a%=mod; b>>=1;a*=a%=mod;//等价于a=(a%mod)*(a%mod),且将一个模运算通过赋值代替,提高了效率 } return ans%mod;//数据不会爆的话,这里的%运算会等价于第5中不断重复的 ans%mod }

例题 POJ 1995 我是在VJ上直接拉过去的 就不放链接了 题目意思: 前面的废话省略,现在有一个游戏,每个人可以有两个数字Ai 和 Bi 现在要你把所以人的Ai^Bi求出来 然后全部加起来 再取模mod 结果输出就可以了 就是:(A1^B1+A2^B2+A3^B3+……………………An^Bn)%mod; 把上面那个式子的结果输出就可以了 数据范围不大 快速幂一遍过 但是这里用到了前面说的公式

(a + b) % p = (a % p + b % p) % p 这个对于多项也是成立的,然后利用上面的公式 将每个人的Ai 和 Bi快速幂取模 然后相加 全部加起来之后再取一次模 就拿到结果了 ok 上一下我的代码:

// POJ 1995 #include <cstdio> #include <algorithm> #include <cstring> using namespace std; typedef long long ll; ll qpow_mod(ll a,ll b,ll mod){ if(a==0) return 0; ll ans = 1; while(b){ if(b&1) ans = (ans%mod)*(a%mod); b>>=1; a*=a%=mod; } return ans%mod; } int main(int argc, char const *argv[]) { int T; int m,n; ll sum; //freopen("input.txt","r",stdin); //freopen("output.txt","w",stdout); scanf("%d",&T); while(T--){ sum = 0; scanf("%d",&m); scanf("%d",&n); for(int i=0;i<n;i++){ int a,b; scanf("%d %d",&a,&b); sum = sum + qpow_mod(a,b,m); } printf("%lld\n",sum%m); } return 0; }

很尴尬 用了141ms 虽然不是特别快 但是也不慢了至少是Ac了 好了 快速幂就讲到这里 然后我觉得 毕竟是在搞竞赛的嘛,我认为准备模板啊这些的很正常, 而且在我现在实力不是特别的雄厚的时候,准备一些模板我觉得是可取的,虽然说不能完全的理解这个模板,但是多少懂一点,然后会用,在慢慢用的过程中自己在去领悟这个算法的真谛 之前打杭电的多校联赛,清华大学的team090这只队伍四次比赛三次AK,有一个很难的一题没有A出来,很佩服他们,和大神之间有很大差距,我们可以向他们学习,至少比昨天的自己更加优秀

遭遇了不理解,枯燥,TLE,WA………这些之后还在坚持的走下去,等看到了Ac,内心的成就感,我想这就是我们学算法的一个小意义所在吧 !好了,要去吃饭了。一起加油!

转载请注明原文地址: https://www.6miu.com/read-45983.html

最新回复(0)