栏目分类:
子分类:
返回
名师互学网用户登录
快速导航关闭
当前搜索
当前分类
子分类
实用工具
热门搜索
名师互学网 > IT > 软件开发 > 后端开发 > C/C++/C#

快速幂 ex ///快速乘 mul

C/C++/C# 更新时间: 发布时间: IT归档 最新发布 模块sitemap 名妆网 法律咨询 聚返吧 英语巴士网 伯小乐 网商动力

快速幂 ex ///快速乘 mul

快速幂(化幂为乘):快速计算 a * b % mod 的结果。(时间复杂度log n, 且可以在过程中取模)

快速乘(化乘为加) : log n

求a的b次方,取模mod(1<=a,b,mod<=1e18)

Input

多组输入,每组数据一行,3个正整数,分别为a,b,mod

Output

每组数据输出一行,为答案

Sample Input

2 10 10000000
5 100 1
0 2 37

Sample Output

1024
0
0

Hint

 这个题因为过程量较大, 会爆 longlong 所以要用快速乘 处理过程量。

 

快速乘模板:

对于 :a*53 = a*110101(二进制) 若 a = 5;

 

ll mul(ll a, ll b, ll mod)
{
    ll ans = 0;
    while(b)
    {
        if(b & 1)
        ans = (ans + a) % mod;
        b >>= 1;
        a = (a * 2) % mod;
    }
    return ans;
}
 快速幂

快速幂代码只有两点不同:

1 :    ll ans = 1; (因为是化幂为乘)

2 :    a =  ( a * a ) % mod;

ll ex(ll a, ll b, ll mod)
{
    ll ans = 1;
    while(b)
    {
        if(b & 1)
        ans = ans * a % mod; // 若数值过大。
    ans = mul(ans, a, mod);
        b >>= 1;
        a = a * a % mod;
    a = mul(a, a, mod);
    }
    return ans;
}

ac代码:

#include

#include

#include

#include

#include

using namespace std;

typedef  long long ll;

const int N = 110;

const int MOD = 1e9 + 7;

int t = 0;

int n = 0;

ll mul(ll a, ll b, ll mod)

{

    ll ans = 0;

    while(b)

    {

        if(b & 1)

        ans = (ans + a) % mod;

        b >>= 1;

        a = (a * 2) % mod;

    }

    return ans;

}

ll ex(ll a, ll b, ll mod)

{

    ll ans = 1;

    while(b)

    {

        if(b & 1)

        //ans = ans * a % mod;

    ans = mul(ans, a, mod);

        b >>= 1;

        //a = a * a % mod;

    a = mul(a, a, mod);

    }

    return ans;

}


 

int main()

{

    cin.tie(0);

    std::ios::sync_with_stdio(false);

    ll a, b, mod;

    while(~scanf("%lld%lld%lld", &a, &b, &mod))

    {

        cout << ex(a, b ,mod) << endl;

    }

}

转载请注明:文章转载自 www.mshxw.com
本文地址:https://www.mshxw.com/it/698183.html
我们一直用心在做
关于我们 文章归档 网站地图 联系我们

版权所有 (c)2021-2022 MSHXW.COM

ICP备案号:晋ICP备2021003244-6号