快速幂算法


模运算:
公式:
  (a+b)mod n=(a mod n+b mod n)mod n;   
  (a-b)mod n=(a mod n-b mod n)mod n;
  (a*b)mod n=((a mod n)*(b mod n))mod n;//要保证n是整数//要知道a mod n和b mod n都是比n小的

利用这些共识可以有效地防止溢出。但是(a mo n)*(b mod n)可能会超出int ,所以中间结果有必要用long long 保存。
 
求a*b mod n
 int mul_mod(int a,int b,int n)
{
  a %=n;
  b%=n;
  return (int)((long long)a*b%n);
}
但是如果n本身超过int又在long long 范围内上述方法就不适用了。

大整数取模:

输入正整数n,m,求n mod m的值。n<=10^100,m<=10^9.

需要将大整数写成自左向右的形式,例如:1234=(((1*10)+2)*10+3)*10+4;然后对每一步进行取模;

代码实现:
#include<cstdio>
#include<cstring>
using namespace std;
int main()
{
    char n[1005];//n的值超过了整型数的范围,所以用字符串数组存储
    long long m;
    scanf("%s%lld",n,&m);
    int len =strlen(n);
    int ans=0;
    for(int i=0;i<len;i++)
    {
        ans=(int)(((long long)ans*10+n[i]-'0')%m);//对于数字来说 n='n'-'0';
    }
    printf("%d\n",ans);
    return 0;
}

快速幂算法
求 a^b mod n的值;
显然可以写出这样的代码:
   int pow_mod(int a,int b,int n)
 {
  int ans=1;
  for(int i=0;i<b;i++)
  {
      ans*=(int)((long long )ans*a % n);
  }
 }

但是当b很大的时候很不理想

可以利用分治法(快速幂取模运算):
 模板代码:
int PowerMod(long long a, long long b, long long c)
{

    long long ans = 1;
    a = a % c;//缩小a
    while(b>0)
    {
        if(b%2==1)
            ans = (ans * a) % c;//如果指数是奇数需要先乘一个a
        b = b/2;
        a = (a * a) % c;//不断地将a平方并取模
    }
    return ans;
}//看似只是在不断的改变a,但是b不断除2 ,总会等于1,最后再根据a得出ans.//十分类似于二分查找






  • 0
    点赞
  • 0
    收藏
    觉得还不错? 一键收藏
  • 0
    评论

“相关推荐”对你有帮助么?

  • 非常没帮助
  • 没帮助
  • 一般
  • 有帮助
  • 非常有帮助
提交
评论
添加红包

请填写红包祝福语或标题

红包个数最小为10个

红包金额最低5元

当前余额3.43前往充值 >
需支付:10.00
成就一亿技术人!
领取后你会自动成为博主和红包主的粉丝 规则
hope_wisdom
发出的红包
实付
使用余额支付
点击重新获取
扫码支付
钱包余额 0

抵扣说明:

1.余额是钱包充值的虚拟货币,按照1:1的比例进行支付金额的抵扣。
2.余额无法直接购买下载,可以购买VIP、付费专栏及课程。

余额充值