关闭

快速幂

123人阅读 评论(0) 收藏 举报
分类:

http://www.cnblogs.com/buptLizer/archive/2011/11/10/2244486.html

运算规则

模运算与基本四则运算有些相似,但是除法例外。其规则如下:

(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)


快速幂运算:

int pow(int a,int n)
{
    int rs=1;
    while(n)
    {
    if(n&1)
      rs=rs*a;
    a=a*a;
    n=n>>1;
  }
  return rs;
}
快速幂取余:

//求a^b%n,O(logb)
__int64 get_mi_mod(__int64 a,__int64 b,int n)
{
    if(0 == a)
        return 0;
    if(0 == b)
        return 1;
    __int64 rs=1;
    while(b)
    {
        if(b&1)
            rs=(rs*a)%n;
        a=(a*a)%n;
        b>>=1;
    }
    return rs;
}

0
0

查看评论
* 以上用户言论只代表其个人观点,不代表CSDN网站的观点或立场
    个人资料
    • 访问:170766次
    • 积分:4930
    • 等级:
    • 排名:第5753名
    • 原创:329篇
    • 转载:0篇
    • 译文:0篇
    • 评论:30条
    收藏网址